薄板弯曲单元的精化矩阵

本文提出一个改进薄板弯曲单元计算精度的一般性方法：用精化矩阵对原有的单元刚度阵进行修正。实践表明，这种方法可以有效地改变单元的刚柔特性，从而达到提高单元精度的目的。文中证明了，单刚中加进精化矩阵不会破坏原来单元的收敛性质。     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形分原理和全局拉拉格朗日方法，推导出几何非线性精化三角形薄板单元，对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象，数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形变分原理和全局拉格朗日方法，推导了几何非线性精化三角形薄板单元。对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象。数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

8—21节点块体精化不协调元

本文基于文（２）提出的精化直接刚度法的变分原理，构造了８－２１节点元的精化不协调位移模式，并装入ＳＡＰ５程序中，该单元能通过分片试验数据结果表明：采用精化直接刚度法，可以有效地改装单元特性，使新单元在具原单元的简单、高效、可靠性好的同时，在精度，适应性等方面有大幅度的提高。     

精化杂交压电固体单元的研究

基于力、电耦合问题的三类交量广义交分原理，提出了广义杂交压电单元列式。为了进一步改进单元的性能和保证单元能够通过分片检验，通过引入非协调模式、放松电学方程约束条件和单元间的弱连续性条件，建立了新的、修正的广义交分原理，在此基础上成功地引入了应力、应交的正交化插值模式，从而建立了精化杂交压电单元法，它继承了常规精化杂交单元的全部优点。文中所推导的八节点精化杂交压电固体单元列式完全避免了矩阵求逆运算，较广义杂交压电单元和杂交应力压电单元均显著提高了计算效率。数值算例表明，与同类型其他单元相比，该单元明显具有更好的对歪斜网格的适应性。     

一种提高薄板稳定分析精度的方法

在薄板稳定分析中，九参数三角形薄板单元因其形状简单，使用方便，在实际工程中得到了广泛应用。本文基于参数调正的几何刚度矩阵，对九参数三角形薄板单元的一致刚度矩阵进行了修正，为提高计算精度提供了一种有效方法。     

复合材料层合板精化高阶理论及其精化三角形板单元

提出一种新的精化高阶理论，该理论满足层间位移、应力连续条件，由此建立了三角形精化板单元。该单元满足单元间Ｃ１类弱连续条件，其收敛性得到保证，且具有简单、高效率的优点。     

大型结构分析程序系统单元库更新

阐述了大型结构分析程序系统单元库更新的概念及其意义，报告了装入ＡＰ－１程序系统单元库的精化四节点及八节点平面等参元的情况，给出了更新前、后单元通过考题库测试的结果及与其它程序系统结果的比较。数值结果表明通过更新单元来提高程序系统的分析精度与计算效率从而增强其竞争力是可行的。     

大型结构分析程序系统单元库更新

阐述了大型结构分析系统单元库更新的概念及含意义，报告了装入ＡＰ－１程序系统单元库的精化四节点及八节点平面等参元的情况，给出了更新前、后单元通过考题库测试的结果及与其它程序系统结果的比较。数值结果表明通过更新单元来提高程序系统的分析精度与计算效率从而增强其竞争力是可行的。     

埃尔米特梁单元的分块对角与高阶质量矩阵

埃尔米特梁单元常用的集中质量矩阵,是由挠度自由度对应的一致质量矩阵元素通过行求和或节点积分构造。然而,数值结果表明该集中质量矩阵在求解包含自由端的梁振动问题时,会出现频率精度掉阶现象。本文首先从保障质量矩阵最优收敛性的数值积分精度出发,分别针对三次和五次梁单元,发展了质量矩阵的梯度增强节点积分方案。利用梯度增强节点积分方案,可以得到具有分块对角形式的单元质量矩阵,而其组装的整体质量矩阵除边界节点外仍然呈现对角形式。对于两种单元,其分块对角质量矩阵分别具有4阶最优精度和6阶次优精度。再者,将标准一致质量矩阵和具有同阶精度的梯度增强节点积分质量矩阵进行优化组合,建立了具有超收敛特性的高阶质量矩阵。最后,通过数值算例系统验证了三次和五次单元的分块对角与高阶质量矩阵的频率计算精度。     

基于重叠划分的自由网格四边形单元计算方法

提出了一种基于重叠划分的自由网格四边形单元计算方法。这一方法将四边形单元引入到自由网格计算方法中，不仅提高了计算的精度，同时还保留了自由网格计算方法的特点。方法首先对分析域内自动生成的每一个节点建立一套临时三角形单元，利用这些临时三角形单元组合生成四边形单元，以节点为单位进行计算。由于各矩阵的计算与组集均以节点为中心进行处理，因而特别适合于并行计算环境。在详细介绍自由网格四边形单元计算方法的基础上，利用数值算例证实了这一方法改善计算精度方面的有效性。     

一种提高四边形四节点平面壳单元计算精度的新方法

平面壳单元是由平面应力单元和平板弯曲单元叠加组合而成,具有简单的理论表达,但是它在计算曲面壳体结构时误差较大。为了进一步提高平面壳单元的计算精度,本文提出了一种计算平面壳单元刚度矩阵的新方法。通过该方法在高斯积分点建立多个单元局部坐标系,并保证每个局部坐标系都位于单元在高斯点处的切平面上,从而可以有效适应曲面壳体形状,达到进一步提高平面壳单元计算精度的目的。为了在这种新坐标系下计算单元刚度矩阵,给出了求解形函数对局部坐标的导数、局部到自然坐标系积分转换的雅可比、以及局部到整体坐标系的转换矩阵的新型计算方法。通过将这些新坐标系以及新计算方法运用到平面壳单元DKQ24中,可以有效提高平面壳单元尤其是在计算曲面壳体时的精度。计算结果表明,本文方法和平面壳单元相结合,不仅具有平面壳单元简单的理论表达式,还能得到满意的精度。另外,本文方法还可以应用到其他类型的平面壳单元,为提高其他类型平面壳单元的计算精度提供了一种新的途径和思路,具有广阔的应用前景。     

无单元法在薄板稳定问题中的应用

用无单元法研究了薄板的弹性稳定问题，从滑动最小二乘法和变分原理出发导出了薄板的无单元法几何刚度矩阵，编制了相应的计算程序，并给出了算例，结果表明，方法合理可行，且精度高于有限元。     

单变量有限元的新思考：精化直接刚度法

对单变量有限元和多变量有限元进行了综合分析和比较。其中包括有限元法的单元的可靠性问题：收敛性，坐标不变性，伪零能模式；单元的性能；计算精度，对畸变网格的适应性及解除闭锁、解的晃动现象等问题，进一步对单变量有限元提出新的思考，有精化直接刚度法发展变量有限元，使之能兼备单变量有限元与多变量有限元之长，克服两者之短。     

直接基于单元平衡的变截面梁反应分析方法

固定形状的单元位移插值函数不能合理地近似变截面梁内部的位移变化,从而影响了传统梁单元用于计算变截面梁的精度.采用直接基于单元平衡的思想给出了计算变截面梁反应的有限元方法,解决了单元位移插值函数局限性所带来的问题.导出了变截面梁单元的单元刚度矩阵、单元等效节点荷载和单元一致质量矩阵.在此基础上,利用编制的程序进行了算例验证与分析.算例验证了本文理论的正确性,表明本文方法具有很高的计算精度.     

非线性杆单元荷载矩阵积分算子直接集成算法

在文献[1]的基础上，将积分算子直接集成单元刚度矩阵的方法扩展应用于非线性杆单元的荷载矩阵集成．该方法将分布荷载产生的内力函数分解为由未知参数表示的线形项，和与分布荷载特性相关的已知函数项的组合，计算积分算子，由积分算子的线形组合得到分布荷戟产生的杆件内力函数，从而得到荷载矩阵．该方法没有采用任何假定，具有广泛的普适性，能用于求解任意刚度分布杆件上作用任意分布荷载或集中荷载的单元荷载矩阵．对于刚度非连续、荷载分布非连续的情况，该方法避免了高斯数值积分的误差，具有更高的计算精度。     

高层框筒和筒中筒结构的整体稳定分析

本文将有限元法用于高层筒 体结构的整体稳定分析。导出了截面含双对称轴的矩形筒的筒单元刚度矩阵和单元几何刚度矩阵。可用于框筒和筒中筒结构的整体稳定计算。计算结果验证了本文方法的精度。     

计算框架结构之动态特性的改进的有限单元方法

本文根据模态综合法之理论与方法,在经典的平面杆单元之形函数中,增加了两个固定界面主模态函数,从而推导出了适合于框架结构动力分析的改进的有限单元模型、单元特性矩阵以及由此单元特性矩阵组集得到结构的总体特性矩阵的方法。通过算例表明:这样的有限单元方法,用于求解框架结构的固有特性时在相同的精度之情况下,结构自由度可以减少近一倍。从而有效地节约所用的计算机的内存与机时。     

钢筋混凝土剪力墙单元非线性分析模型及其应用

本文简要评述了几种常见的钢筋混凝土剪力墙非线性宏观单元模型后，着重对多垂直杆剪力墙非线性单元模型的几个重要问题如剪力变形的考虑方法，单元刚度矩阵的形式，垂直拉压杆及剪切弹簧的恢复力模型等进行了探讨与改进，最后给出了一个算例，并与试验结果比较，表明非线性宏观墙单元模型具有较好的计算精度。     

移动荷载作用下地基动力分析的有限元方法

通过对地基动力问题的基本方程进行变换,把基本方程变换到随荷载移动的运动坐标系中,通过加权残数法推导了相应的单元刚度矩阵,从而建立了移动问题的有限元格式,并发现移动荷载问题的单元刚度矩阵是对相应静力问题单元刚度矩阵的修正,在静力单元刚度矩阵的主对角元素上增加与移动速度有关的项,即可得到移动问题有限元的单元刚度矩阵,这样就将动力学问题转化为“拟静力”问题处理。文中用移动问题有限元方法计算了地基的动力响应,并与解析解进行了对比,以说明本方法具有较好的精度。