张拉结构非线性分析两节点曲线单元有限元法

提出适合张拉结构几何非线性分析的两节点曲线单元有限元方法．假定索元的初始形状呈二次抛物线,根据单索的平衡条件、几何和物理关系建立了索元的位移函数;由拉格朗日应变的定义建立了可以考虑任意次高阶位移影响的索元轴向应变的精确表达式,并基于拉格明日描述方法和虚功原理得到了索元的非线性平衡方程与切线刚度矩阵．采用荷载增量法与Newton－Raphson法相结合的混合法进行了实例计算,结果表明：本文方法的精度明显优于两节点直线索单元,适合于大跨度索阿、索穹顶等张拉结构的几何非线性分析．     

焊接结构应力场计算中的一类五面体可变节点等参单元

本文应用三维有限元理论，构造了一类五面体６～１５节点等参单元的形函数，提出了在五面体单元内的Ｇａｕｓｓ积分处理方法，建立了相应的单元刚度矩阵。这类单元的节点可以在６～１５之间任意选择。将其与ＳＡＰ５中的六面体８～２１节点等参单元匹配，计算了对接板接头的理论应力集中系数。经与有关文献结果比较表明：这类单元适合于计算焊接结构的应力场，并且具有良好的精度。     

拱型结构在参、强激励下的非线性振动分析

利用数值解析法研究了拱型结构在参数激励与强迫激励联合作用下的非线性振动特性。得到了轴向力与固有频率的关系及轴向力对发生主共振,１／２亚谱共振的影响,由于１／２亚谐共振是高频激振低频响应,是最危险区域,应得到足够的重视,为工程设计提供了可靠的理论依据。     

平面4节点广义等参单元

在文献[1]建立的广义平面矩形单元的基础上,建立了平面4节点广义等参单元,并进行了深入研究,包括分片试验、单元的变形模态以及形函数的影响等.通过典型算例,对单元的性能进行了检验,数值结果表明,常规情况下建立的广义等参单元较传统单元具有更高的计算精度.     

非线性结构分析中三维过渡等参单元及其应力平滑方法

给出了三维四、六、八节点过渡等参元的模型，推导出了相应的应力平滑公式，编制了计算程序，并通过弹塑性问题和接触问题例题考核，结果良好。     

钢管混凝土柱-钢梁节点的力学性能分析

基于弹塑性有限元理论建立了钢管混凝土柱-钢梁节点荷载-位移全过程非线性有限元模型,在单元分析中采用改进的AUL表述推导得到梁柱单元刚度矩阵方程,同时考虑了材料的物理非线性和单元的几何非线性,并编制了非线性有限元程序NLFEACFST。采用该模型对相关研究者和作者进行的节点试验进行了分析,理论计算结果与试验结果比较表明,该模型具有很好的适用性和精度。在理论分析模型得到试验结果验证的基础上,对典型的中柱节点进行了荷载-位移全过程非线性特性分析,并对影响节点承载力和荷载-位移骨架曲线的因素进行了参数分析,为进一步从理论研究钢管混凝土框架结构的力学性能创造了条件。     

张力结构的非线性有限元分析

张力结构中的索结构分析通常采用非线性二节点直线或二次曲线单元，但是在大跨度结构，尤其是索弯顶结构的分析中，常用的单元已不能满足精度的要求，本文提出了考虑自重作用下有初始垂度的五结点非线性空间曲线元模型，放弃了一些特殊的假定，考虑了应变表达式中高阶量的影响，推导出了适合于弹性大位移几何非线性分析的Ｌａｇｒａｎｇｄ方程及其切线刚度矩阵，编制了相应的计算程序。     

等曲率井中有重钻柱屈曲的非线性有限元分析

建立了等曲率井中有重钻柱屈曲的平衡方程及对应的泛函表达式,用有限元法对等曲率井 中有重钻柱屈曲过程进行了分析,给出了钻柱正弦屈曲和螺旋屈曲临界载荷的定义. 力学模 型中考虑了重力、钻柱上端井斜角和井眼轨迹曲率半径对屈曲的影响. 分析结果表明：载荷 增大时,钻柱的下端先出现局部屈曲,随后屈曲向钻柱上部扩展,导致钻柱发生整体屈曲, 屈曲位移、井壁约束力线密度和钻柱弯矩都呈周期性变化;重力对等曲率井中钻柱的屈曲有 较强的抑制作用,其影响不可忽略;井眼轨迹曲率半径越小,钻柱上端井斜角越大则对钻柱 屈曲的抑制作用越强.     

复杂结构二十节点等参元建模法

提出用单一的二十节点等参元模拟平面、梁、板壳和板-梁结构的建模方法.通过调整板-梁接触面耦合节点的数目,分析接触面耦合节点的数目对接触面应力分布和内力的影响.结果表明,二十节点等参元完全具有模拟平面问题、梁和板壳的能力,能提供更多的相关信息;通过调整接触面耦合节点的数目和组成方式,二十节点等参元独立建模法可以模拟板-梁结构中的刚性联接,铰链联接和介于二者之间的弹簧联接的优越性能.这种建模方法不仅大为简化了复杂结构的建模过程,而且使有限元模型建立在更直观和真实的物理基础之上,避免引入传统建模方式所必须引入的联接单元,有利于有限单元法的推广和应用.     

焊接结构应力场计算中的一类五面体可变节点等参单元

本文应用三维有限元理论，构造了一类五面体６－１５节点等参单元的形函数，提出了在五面体单元内的Ｇａｕｓｓ积分处理方法，建立了相应的单元刚度矩阵，这类单元的节点可以在６－１５之间任意选择，将其与ＡＳＰ５中的六面体８－２１节点等参单元匹配，计算了对接板接头的理论应力集中系数，经有与有关文献结果比较表明：这类单元适合于计算焊接结构的应力场，并且具有良好的精度。     

线性等参单元逆变换的几何二分法及其修正

在证明线性等参单元的一个几何特性的基础上,提出了线性等参单元逆变换的几何二分法。几何二分法简单易行,避免求解高维高次非线性方程组。在几何二分法基础上,提出一种修正措施以消除畸形单元对所求结果的影响,几何二分法及其修正措施克服了以往方法存在的缺点,两个实例验证了由几何二分法及其修正措施进行等参逆变换的正确性。     

拉索穹顶结构非线性分析的混合有限元增量法

拉索穹顶结构是由受压桅杆和拉索组成的新型柔性大跨度空间组合结构,几何上表现为极强的非线性特性,计算困难,本文应用有限元法,结合拉索穹顶结构特征,假定拉索和桅杆的受力满足虎克宣定律,建立了可以直接考虑拉索垂度影响的两节点索单元模型,并与两节点直杆单元相结合,基于修正的拉格朗日描述方法和虚功原理建立了拉索穹顶结构非线性分析的混合有限元增量方程。采用荷载增量法与Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法相结合的求解     

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

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

含两类附加函数的扩展等参有限元法

基于扩展有限元的基本思想,提出一类指数型间断函数来模拟。由于裂纹或节理等非连续结构所导致的位移不连续现象,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减,同时,在非连续结构末端引入能反映其奇异场特性的三角基函数。本文用弱解形式推导了扩展有限元格式,并论证了两类附加函数在单元公共边上能够保持位移连续性这一要求。最后,编制了二维4节点和三维8节点的扩展等参有限元程序,并分别给出了算例,结果表明在模拟裂纹追踪时,扩展有限元法可行且有效。     

平面Cosserat模型有限元分析的4和8节点单元与分片试验研究

经典连续体理论不包括物质内部尺度,当考虑应变软化问题时,有限元结果对网格具有很强的依赖性。与经典连续介质力学理论不同,Cosserat连续体模型在传统平动自由度的基础上添加了一独立的旋转自由度,在本构模型中引入了内尺度参数。本文研究了基于Cosserat理论的平面4和8节点等参元以及8（4）节点线、角位移混合插值等参单元,给出Cosserat单元分片试验的实施过程。最后将单元运用到小孔应力集中问题的分析当中,通过计算结果与理论解的比较,表明了4和8节点以及8（4）节点等参元的适用性,为问题的非线性分析打下基础。     

一种新的覆冰导线舞动非线性有限元分析方法

覆冰导线舞动严重威胁输电线路的安全稳定运行,论文针对覆冰导线的大幅舞动问题,基于Hamilton原理建立了耦合三个平动自由度和一个扭转自由度以及偏心覆冰作用下导线的非线性舞动方程,提出了一种新的覆冰导线舞动分析非线性有限元模型,将相邻跨导线和绝缘子等效为线性弹簧单元,同时考虑覆冰导线非线性气动力和大幅舞动的几何非线性,采用Newmark-β时间积分结合修正的Newton-Raphson非线性迭代求解舞动有限元方程.论文数值解与D形覆冰导线舞动实测值无论在振幅还是频率方面均吻合得非常好,证明了论文新提出的导线舞动有限元模型的准确性.论文研究表明:舞动是一种以上下运动为主的低频振动,通常发生在一阶上下自振频率附近,其振幅和频率由输电线路的物理参数和风载荷唯一确定,与初始运动状态无关.通过与现有方法对比,论文在覆冰导线的舞动模型和气动力处理方面有较大的创新和改进,计算精度更高,效率较传统有限元法有显著提升,能为输电线路舞动研究和防舞抑舞控制提供一种新的方法与思路.     

A simple finite element method for Stokes flows with surface tension using unfitted meshes

This paper presents a simple finite element method for Stokes flows with surface tension. The method uses an unfitted mesh that is independent of the interface. Due to the surface force, the pressure has a jump across the interface. Based on the properties of the level set function that implicitly represents the interface, the jump of the pressure is removed, and a new problem without discontinuities is formulated. Then, classical stable finite element methods are applied to solve the new problem. Some techniques are used to show that the method is equivalent to an easy‐to‐implement method that can be regarded as a traditional method with a modified pressure space. However, the matrix of the resulting linear system of equations is the same as that of the traditional method. Optimal error estimates are derived for the proposed method. Finally, some numerical tests are presented to confirm the theoretical results. Copyright © 2015 John Wiley & Sons, Ltd.     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

预应力混凝土平面杆系结构的有限元方法

建立了基于有限元方法的考虑材料和几何非线性的任意截面预应力混凝土平面杆系结构的数值分析模型,可用于模拟预应力混凝土大跨度梁、单向偏压细长柱等的非线性全过程结构响应。引入修正的Rodriguez截面模型确定截面切线刚度,其中混凝土的贡献通过截面边界顶点定义的梯形单元来实现;在此基础上利用传统的平面非线性杆单元导出了标准有限元公式。通过两个算例验证了该模型的可靠性和适用性。     

附加自由阻尼结构系统的有限元分析和拓扑优化设计方法研究

给出了一种适用于车身等复杂附加阻尼结构系统的有限元建模和动态特性分析方法,包括附加自由阻尼薄壁结构和附加阻尼材料粘弹性特性的有限元建模方法以及附加阻尼结构动态特性的有限元分析方法。在此基础上,又以强迫振动响应最小为优化设计目标,给出了一个附加阻尼结构的拓扑优化设计方法,包括优化设计问题的列式和适用的求解算法。通过对一个简单薄壁构件和一个车身地板上的附加阻尼结构的拓扑优化设计,验证了提出的有限元分析和拓扑优化设计方法的有效性。