一个改进的平面梁单元

根据有限单元法基本原理 ,提出了一个变截面平面梁单元 ,推导了其单元钢度矩阵。这一改进的梁单元用于分析梁高呈线性变化及二次抛物线变化的矩形截面梁 ,将得到准确解。文中给出了一个变截面悬臂梁算例 ,计算表明 ,这一改进的梁单元使变截面梁的分析大大简化     

一种精确积分的四结点板弯曲单元

本文从Ｍｉｎｄｌｉｎ／Ｒｅｉｓｓｎｅｒ理论出发,采用一种新的平行四边形母单元和相应的形函数推导四结点板弯曲单元刚度矩阵的精确积分解。弯曲应变和横向剪切应变分别采用不同的插值公式构成单元刚度矩阵。理论和算例分析表明本文方法克服了“闭锁”现象并能应用于很薄的板,单元刚度矩阵计算速度比采用数值积分计算的同类单元的快四倍。     

薄膜横向振动的有限元分析

构造了一种3节点三角形膜单元,以适用于平面薄膜横向振动的有限元分析.在给出单元形函数的基础上,根据最小势能原理建立了薄膜自由振动方程,并推导了单元刚度矩阵和单元质量矩阵.研究结果表明,单元刚度矩阵和单元质量矩阵形式简单,且自由度少;通过两个典型算例,证明3节点三角形膜单元的计算结果非常接近理论解,同时可以达到很高的精度...     

三维弹性接触问题的接触面单元法

基于虚功原理，推导了三维接触面单元的刚度矩阵，并引进预留单元的概念，避免了接触过程中由于接触面变化节点和单元需要重新编号的麻烦．文中算例验证了计算方法的有效性     

变截面高层框筒结构的矩阵传递法

以连续化数学模型和能量变分原理导出的变截面高层框简结构的微分方程组为基础；根据变截面框简结构的截面沿高度为阶形变化的特点，以每个相同截面的层作为一个计算单元，每个单元由微分方程组导出其单元矩阵、截面矩阵和传递矩阵。进而采用传递矩阵法进行分析。此法概念清楚，计算简单．并以算例说明其应用。     

斜拉桥分析的索膜单元理论及其应用

本文提出分析斜拉桥的索膜单元理论,建立起索膜单元概念与计算模型,推导出索膜单元刚度矩阵,说明了用索膜单元分析斜拉桥的方法细节,应用索膜单元对144+360+144米公路钢斜拉桥方案进行了分析计算,得出了主梁弯矩影响面与扭矩图,讨论了索梁锚固点与截面水平主惯性轴偏心对主梁弯矩的影响。     

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

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

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

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

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

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

反向Qm6非协调元

为通过强式分片试验,Qm6单元对Q6单元非协调部分的[G]矩阵进行了特殊的计算处理,但抗畸变性能下降,本文提出对有关处理反向进行,以恢复甚至提高抗畸变性能。分析了Qm6单元的原理,指出其实质是修改雅可比矩阵[J]的伴随矩阵[J*],在非协调部分[G]矩阵的计算时,把[J*]看成可变量,由Qm6的对应点向Q6方向进行反向搜索,查找有利的计算点。进行了典型和苛刻的算例测试,结果表明反向调整是有效的,调整系数取镜像值-1以及扩展到-2时,新单元的抗畸变性能优于原Q6和Qm6,其中取-2对消除剪切闭锁是最优点;除弱式分片试验外,总体性能和精度接近各类4节点四边形单元的最好水平。由于方法和原理简便,实现以及推广到三维问题都有显著优势。     

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

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

无单元法分析薄板自由振动问题

提出了分析四边简支和四边固定薄板自由振动问题的无单元法，推导了无单元法的插值函数，从变分原理出发导出薄板自由振动的质量矩阵和刚度矩阵，编制了相应的计算程序，通过计算实例与其它方法的结果进行比较．数值结果表明无单元法具有一系列优点。     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。     

用可退化有限单元进行平面连续体拓扑优化

提出一种基于准则法的平面连续体结构拓扑优化方法。该方法根据单元在结点处的厚度或应力值决定删除单元或将矩形单元退化为三角单元 ,使迭代过程更加平稳。在一定的计算量下 ,提高了计算精度 ,减少了迭代次数 ,使优化结构的边界更加光顺。该方法可同时应用于等厚度和变厚度平面连续体的拓扑优化问题。本文推导了变厚度矩形有限单元的刚度矩阵 ,并给出了相应的拓扑优化迭代方法     

降落伞膜索非线性有限元模拟

用膜单元和索单元模拟降落伞织物绳索系统,基于完全拉格朗日格式的非线性有限元方法编程计算降落伞的结构动力学特性。采用增量与迭代混合方法改善非线性计算的收敛特性并结合HHT隐式时间推进方法减小整体迭代计算量。使用修正应力应变张量导数的方法模拟膜单元单向应力状态并针对膜单元和索单元分别进行了非线性有限元计算验证。最后针对C-9型降落伞建立三维有限元模型,根据设定流速对伞衣施加均匀压强载荷,将模拟展开的结果与使用相同模型、不同方法商业软件的文献进行对比,显示了隐式非线性有限元方法模拟降落伞膜索系统大变形动力学的能力。     

不同材料交界面上接触应力的有限元分析

工程中有许多不同材料组成的结构。对这类复合结构，正确计算材料交界面上的接触应力是非常重要的。通常的位移有限元法得到的应力是不连续且欠准确的。本文提出了一种新的用于不同材料交界面处的接触单元。应用虚位移原理和结合考虑材料交界面上的约束条件，文章建立了接触单元的刚度矩阵和等效荷载向量。该单元具有算法简单，程序容易实现的优点。计算实例证实了该单元的实用性和可靠性。     

无单元法计算钢筋混凝土筏板

对无单元法的基本理论作了介绍,并从变分原理出发推导Ｗｉｎｋｌｅｒ地基上钢筋混凝土筏板的无单元法刚度矩阵,编制了相应的计算程序,通过计算实例验证了本文方法的有效性。     

考虑约束扭转的薄壁梁单元刚度矩阵

推导了薄壁空间梁单元刚度矩阵 ,考虑了双向弯曲及截面约束扭转对杆件轴向变形的影响 ;计算了截面的翘曲变形 ,以及二次剪应力对翘曲变形的影响 ,可适用于任意截面 (包括开口、闭口和混合剖面 )的薄壁杆件。计算结果表明 ,考虑约束扭转的薄壁梁单元刚度矩阵有相当好的精确度 ,可以用于薄壁杆件的静动力分析。     

二维弹性接触问题的接触面单元法

本文基于虚功原理，推导了二维接触面单元的刚度矩阵，并引进预留单元的概念，避免了接触过程中由于接触面变化，节点和单元需要重新编号的麻烦。采用位移和应力联合控制的增量法控制加载过程，文中给出考题验证计算方法的有效性，并给出本文方法在汽检轮机中根与轮缘接触问题中的应用实例。     

加权残值法在钢筋混凝土拱桥非线性有限元分析中的应用

本文用圆弧梁离散拱肋;用圆柱拖带坐标、三次位移插值函数及平截面假定来描述单元位形;用加权残值配点法来消除曲梁单元的剪力与膜力闭锁。按基于连续介质力学的Ｕ．Ｌ．列式建立单元增量平衡方程,以考虑几何非线性。假定钢筋为理想弹塑性材料。按三参数各向同性强化塑性模型,建立混凝土的弹塑性本构矩阵。将拱单元分段分块,根据钢筋及砼的本构特性,建立拱单元及梁段单元的弹塑性刚度矩阵,以考虑材料非线性。用编制的程序对两座模型拱桥进行计算,计算结果与模型测试结果接近