弹性平面问题通用程序

本文给出了以有限元法解弹性平面问题的通用程序。使用该程序计算弹性平面问题不受边界形状、边界位移、边界力以及板厚变化的限制。计算时首先划好网格，根据网格填写初始数据表并穿孔，之后利用本程序在 ７０９机上按一般上机操作便可以算题。如果利用其它ＡＬＧＯＬ ６０语言机算题也可以考虑使用。     

边界元法与有限元法耦合的对称—迭代求解

本文针对二维弹性静力问题的边界元法与有限元法的耦合方法,首先论述了建立在配点离散方法之上的边界元刚度矩阵的不对称性,继而提出了边界元法与有限元法耦合的对称—迭代求解方法,并讨论了迭代的收敛性。计算结果表明,该方法有很高的精度。     

一种新的用于二维弹性静力学的快速多极边界元法

快速多极边界元法(fastmultipole BEM)是近几年发展起来的边界元新型算法。本文提出了一种新型的适合二维弹性静力学问题的快速多极边界元格式，并用于含有多个夹杂的二维复合材料的应力分析。数值结果表明这种方法非常适合解决大规模问题。     

直接边界元法奇异积分的研究

在用直接边界元法解决弹性问题时,当场点与加载点重合时,边界积分方程将出现奇异.为了减少计算误差,有必要求出奇异积分的解析式.应用高等数学基本理论推导出二维弹性问题直接边界元法奇异积分的解析式,采用线性单元离散边界,求出了奇异的对角线子矩阵元素的解析表达式.根据理论推导的结果,编制了相应的计算程序.可用于分析弹性二维问题的位移场和应力场.     

平面弹塑性问题边界元法与有限元法耦合的对称—迭代求解

在一些弹塑性问题中,塑性区域比较小,大部分区域是弹性变形区域,若采用边界元与有限元耦合的方法,则能充分发挥两种数值计算的优点.在以往的边界元法与有限元法耦合的论文中,边界元子域的刚度矩阵都取为对称阵,尽管易于耦合,但会产生明显的误差.若细分边界元网格,又会大大增加计算时间.本文讨论用作者提出的对称-迭代求解方法处理弹塑性问题中边界元法与有限元法的耦合.     

直齿圆柱齿轮三维接触应力边界元分析

提供了一种有摩擦的三维弹性接触问题边界元分析方法及其前处理技术，并将此法用于一对处于啮合状态的直齿轮，得到了满意的结果。该程序使用方便，只需输入齿轮的基本参数，即可得到齿轮不同啮合位置的网格、应力分布。     

混合有限元法在弹性接触问题中的应用

讨论了利用Hellinger-Reissner二类变量广义变分原理解决定边界接触问题的混合有限元法,并且给出这种计算方法的数值实例。数值结果表明该方法在定边界弹性接触问题中的应用行之有效。     

直接边界元法及其在弹性力学问题中的应用

文章通过弹性力学问题的基本解将域内微分方程变换成边界上的积分方程,然后在边界上离散;由已知边界位移和边界应力直接求出未知边界位移和边界应力,并得出据以计算整个问题域的位移场和应力场。最后运用此方法求解一个弹性力学问题并与有限元法的计算结果进行了比较。     

一种计算轴对称磁场的边界元—有限元混合法的研究

由静磁场的基本方程和有限元的原理出发，研究了一种新的计算带非饱和磁介质轴对称磁场的边界元－有限元混合法，编制了计算电磁聚焦成象系统聚焦磁场的程序，这种方法解决了单纯使用有限元法的边办封闭问题，为成像系统及其它磁透镜的工程设计提供了一种有效方法。     

动节点对自适应接触边界元法

Kamiya提出的弹性接触问题边界元分离解法和申光宪、陈一鸣、刘德义等引入最优罚因子的边界元分离解法，关于计算中接触状态的判别比较模糊，所得接触区与真实接触区相差很大。对此，本文进行了严密的分析，提出一种新的动节点对自适应分离接触迭代判别方法及实现方案并详细描述了数据结构，同时，还编制了自适应动节点对分离接触迭代方法的边界元程序(2-DECPAOBEM)，证实了本文新方案、算法的有效性，从而获得更精确的真实接触区及边界元解。     

自然边界元与有限元的直接耦合及其并行化实现

首先介绍了有限元和自然边界元的直接耦合法,然后以一个含3个裂纹等截面杆的扭转问题为例编写了并行求解程序,并同串行程序运行时间作了比较,从实践上说明了并行计算在有限元和自然边界元的直接耦合法中的可行性与有效性。     

有限元整体刚度矩阵的一种快速解算方法

根据有限元整体刚度矩阵的特点，提出了快速解算方程的一种方法，通过对程序的运行，可以实现运用常规内存快速计算４０００个节点的有限元计算，并可实现快速解算大型稀疏对称线性方程。     

PDETOOL工具箱在磁共振技术中的应用

有限元法是数值求解偏微分方程的一种典型方法,它在许多科学计算领域中的应用已愈来愈广泛。实际计算中,当偏微分方程定解问题的定义域是复杂几何形状或有复杂定解条件时,直接编写有限元原始程序具有相当难度,例如计算磁共振电阻抗成像技术（MREIT）中的正问题。基于有限元的MATLAB／PDETOOL工具箱库函数,对MREIT中的椭圆方程定解问题进行数值求解,还给出了求解规则区域的定解问题时,不同节点所需的运算时间。数值结果表明方法不仅能有效地求解此类问题而且节点的增加对于时间的影响并不大。     

节理岩体损伤模型的有限元分析

将岩体中固有的初始节理损伤考虑成宏观损伤场，通过引入损伤张量模拟岩体的多裂隙性。损伤张量由节理组的实测数据确定，并引入有效应力来计算节理岩体的力学效应。编制了损伤模型的平面非线性有限元程序，并采用损伤模型和常规有限元模型的有限元程序对某矿山边坡进行对比计算分析。结果表明，损伤模型的有限元解析能反映节理岩体的各向异性。     

缓坡方程的有限元解

波浪在缓变海底传播中,用缓坡方程求解波浪场可得到较理想的结果,对缓坡方程的有限元解法进行了论述,推导出了有限元近似方程,开发了计算程序;计算结果与其他方法的结果吻合良好。     

An eigenfunction expansion method for the elastodynamic response of an elastic solid with mixed boundary surfaces

The method of eigenfunction expansion is one of the most elegant methods for solving elastodynamic problems. The solution obtained from it is more concise than that obtained from the integral transform technique. Traditional eigenfunction expansion method is used for the elastodynamic problems with displacement and traction boundary conditions. In this paper, the method is generalized to study the elastodynamic response of an elastic solid with mixed boundary surfaces, and the exact analytical solution is derived. The dynamic response of a finite-length solid aluminum cylinder with two mixed end boundaries is numerically evaluated. The result computed from the analytical solution agrees very well with that obtained from finite element method （FEM）.     

平面相变热传导问题等效热容法的有限元解

用有限元法计算带相变的平面热传导问题,应用单元相变发生的热焓判据,使程序能够自启动,不必给出单元初始状态输入数据,相界面位置计算按线性插值处理,使程序有较高的计算效率,算例表明,等效热容法有限元程序热焓判据合理,相界面位置计算准确,并能成功地处理相变焓,计算出温度分布。     

基于接触摩擦的多片式钢板弹簧悬架预应力的计算

采用基于接触摩擦的有限元方法,在Ansys中使用面-面接触单元建立钢板弹簧有限元模型.对一组多片式钢板弹簧悬架进行预应力的计算,并与用传统预应力计算方法得到的结果和行业推荐值进行比较分析,证明新的方法更能反映簧片实际的受力情况.由此对簧片的自由曲率半径和自由弧高进行修正,既能满足设计要求,又能提高整副钢板弹簧的使用寿命.     

渗流随机有限元程序设计及煤层顶板涌水预测

基于随机有限元理论和渗流力学原理,通过推导渗流随机有限元列式,得到随机渗流单元水导矩阵,将可靠度计算引入采动岩体随机渗流之中.根据渗流随机有限元列式编写出可用于岩层渗流计算的随机有限元程序,并对俄霍布拉克煤矿的实际渗流问题进行了数值计算.结果表明,程序计算较好地反映了实际岩层的渗流特性;保护关键层和防水层对预防矿井的顶板突水具有重要作用.计算结果可为俄霍布拉克的开采提供合理化建议,也为随机有限元在渗流方面的深入研究提供依据.     

有限元法计算工业管道电场分布的自动剖分方法

给出了一种在计算工业管道电场分布时的计算机自动剖分方法，该方法是在标准区域内进行单元剖分，再坐标变换成实际区域。给出了用计算机Ｃ语言实现自动剖分的结果，并进行了讨论。