等参元逆变换算法在渗流—位移耦合场分析中的应用

在渗流 -位移耦合场分析中应用等参元逆变换算法解决了两个场之间的数据传递问题 ,所提出的逆变换算法是基于解析性质的新算法 ,具有简便、高效和高精度的特点 ,便于在实际工程中应用。     

基于加权残量法的多变量等参元简明列式及其应用

根据作者提出的建立多变量有限元模型的系统化方法,应用加权残量法的基本原理,给出了建立多变量等参元模型的简明列式。然后应用该列式,构造了简明有效的平面四结点多变量等参元MQ4和MQ5前者不含内部自由度,后者含有内部自由度。文中结合若干典型算例对所构造的单元性能进行了考核,并与若干典型协调与非协调等参元的数值结果进行了比较。数值结果表明,无论是在规则网格还是在不规则网格情况下,本文构造的多变量等参元,不仅能通过分片检验的要求,而且表现出了良好的性能。     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

张拉结构非线性分析的五节点等参单元

本文针对张拉结构的特点，提出了一种五节点等参数单元有限元模型，采用四次多项式作为位移插值函数及单元初始形状函数，并假定索是理想柔性的且满足虎克定律，基于修正的Ｌａｇｒａｎｇｉａｎ坐标描述法，建立了非线性有限元基本方程和切线刚度矩阵，利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法进行了实例计算。结果表明：本文方法精度极高，可供大跨度索网，索穹顶，拉线塔等张拉结构分析，设计时采用。     

求解大柔度梁/绳索的ANCF单元及隐式迭代格式

为模拟大柔度梁/绳索结构的变形和大范围运动,基于绝对节点坐标方法ANCF(Absolute nodal coordi-nate formulation)和HHT(Hilber-Hughes-Taylor)积分方法,建立了ANCF单元的隐式动力学迭代格式.得到了简洁的节点等效力向量,且进一步导出了切线刚度矩阵的全部公式,...     

连续体结构拓扑优化的节点独立连续映射法

提出基于节点拓扑变量的独立连续映射连续体结构拓扑优化方法,将常规ICM中基于单元拓扑变量描述改进为基于节点拓扑变量描述,通过单元形函数插值保证了拓扑变量场的连续性.引入拓扑变量和设计变量倒数关系,根据复合函数链式求导法则,推导了各类结构响应量的敏度分析公式.以位移约束问题为例,遵循ICM法建模方式建立了连续体结构拓扑优化模型.基于数字图像处理过滤技术,提出两种改进的过滤方式以消除数值不稳定性现象.最后通过二维拓扑优化数值算例证明方法和模型的可行性和有效性.     

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

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

索-梁组合结构的建模及面内动力特性分析

利用哈密顿变分原理以及结构动静态构型的影响,建立了索-梁组合结构的约化运动学控制方程.考虑到边界条件和耦合连接条件,我们研究了体系的面内特征值问题.根据求解得到的面内特征值方程,并通过分段函数的引入,结构的模态函数可以被直接确定.随后,我们研究了参数垂跨比f,刚度比K和质量比m对面内固有频率的影响.研究发现从结构的频率谱图中可以看出频率跳跃现象是存在的,另外,频率穿越现象也是十分明显.随后,考虑到局部模态和整体模态,结合之前确定的特征值方程及分段振型函数,我们研究了索-梁组合结构可能的模态形状.最后,我们讨论了索-梁组合结构可能发生的内共振形式,比如面内1:1内共振形式以及1:2内共振形式.研究表明梁的静态构型不仅直接影响到耦合力连接条件,还将影响索-梁组合结构频率的确定.     

复杂板壳结构声振耦合特性的理论与实验研究

0引言在国家民用工业及国防工业减振降噪重大应用需求的牵引下,本文通过理论建模、实验验证和数值计算研究了汽车、高速列车、舰艇/潜艇及航空航天飞行器中常用典型复杂板壳结构的声振耦合特性[1-13],建立了相对完善可靠的结构声振耦合特性理论表征体系,分析了关键结构参数对结构声振耦合特性的影响,揭示了弯曲波在结构中的传播规律及结构的声辐射/传声特性,提山了轻质、高强度、声辐射小及隔声性能优良的复杂板壳结构的创新优化设计概念,初步建立了综合结构质量、力学刚度和声振耦合特性的优化设计理论和判据,为典型板壳/板腔结构在民用工业及国防工业中的应用奠定了理论基础、实验依据并提供了技术支撑.     

机动法作影响线与反力-位移互等定理

根据反力-位移互等定理直接导出了影响线的机动作法,并给出了一种概念更为直接明了的计算和绘制超静定力影响线的方法.     

基于带状单元法的电磁弹性多层板复频散特性分析

超声波在弹性介质中传播时理论上可以分为传播、非传播及衰减模式,且分别与实、虚及复波数相对应,可用复频散特性来进行描述.这一性质在电磁弹性多层板结构中同样存在.该文利用带状单元法对其进行了分析和计算.首先,利用能量守恒定律并结合介质的本构方程、平衡方程及相应边界条件构建了结构的带状单元模型,然后利用该模型获得了电磁弹性多层板的频散方程,最后以一个由压电材料BaTiO3和压磁材料CoFe2O4构成的三层板为数值分析对象,给出了波在其中传播时的复频散特性,并利用计算结果证明了当前相关文献中材料常数引用的不合理性.     

基于力等价结构的复杂超静定结构求解新方法研究

超静定结构内力计算是结构力学课程学习的重点和难点,当超静定次数较多时,采用传统方法计算时难度极大且不易理解。本文总结提出了力等价结构思想,力等价结构与原结构在内力、变形上完全等价,基于该思想提出了复杂超静定结构内力计算的新方法,并将其应用到位移法、混合法求解多次超静定的复杂结构算例中。通过与传统方法对比,验证了所提力等价结构思想的可行性。所提方法可大大降低复杂超静定结构基本未知量的数目,简化计算过程,提高求解效率。

     

基于单元分解的结构断裂扩展有限元法

研究并提出了一种基于单元分解的扩展有限元方法,并将其用于求解结构断裂问题。首先,将不含加强节点的四边形单元剖分为四个三角形子域,通过加权平均获得单元中心处局部应变值;其次,基于三角形子域局部应变进一步构造单元刚度矩阵稳定项。最后,将该单元分解法应用到扩展有限元法的分析中。与传统扩展有限元方法相比,该方法可有效减少积分点数量、避免复杂等参变换且能保证裂纹尖端应力强度因子的求解精度。

     

???????????????????????

采用谱单元方法推导带阻尼梁的传递函数矩阵. 采用一维连续梁的解析解作为动力形状函 数. 与有限元法相比,谱单元方法的自由度和单元数目小且计算精度高. 传递函数表示为梁 的几何和物理参数的超越隐函数,不同于用模态参数表示的传统传递函数. 提出了采用遗传 算法的结构物理参数识别方法. 以变截面悬臂梁为数值算例,显示该方法的有效性.     

组合结构罚有限元分析与罚因子的影响

该文讨论了不同刚度组合结构有限元罚单元罚因子的选取，提出了罚因子选取的一般原则，在实例中，提出了罚单元的估计式及罚因子的取值范围，这对一般情况的数值计算也具有一定的参考价值。     

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

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

基于扩展有限元的应力强度因子的位移外推法

针对平面裂纹问题,阐述了扩展有限元法的单元位移模式、推导了扩展有限元法的控制方程、介绍了特殊单元的数值积分技术.基于最小二乘法,建立了应力强度因子位移外推法的计算公式.利用MATLAB编写计算程序,对平面裂纹问题用扩展有限元法进行了计算.基于扩展有限元法的计算结果,分别利用位移外推法和相互作用积分法,对平面裂纹的应力强度因子进行了计算.计算结果表明,位移外推法比相互作用积分法能更方便和准确地计算平面裂纹的应力强度因子.     

基于等效单元概念的梁结构损伤定位方法

研究目的是为梁结构提供一种基于等效单元概念的新损伤定位方法.区分刚度损伤和质量损伤两种情况, 分别对连续梁和简支梁给出了具体的损伤定位方法. 利用的具体指标是支座反力和跨中位移. 文中的新方法都可以得到以小区间表示的损伤位置结果, 有效性得到了数值模拟算例和模型实验的验证.     

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

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

基于自然单元法的极限上限分析

自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.