二维有限元线法自适应分析的若干新进展

有限元线法(FEMOL)是一种优良的半解析、半离散方法,将其比拟为广义一维问题,遂可将一维有限元中十分成功的单元能量投影(EEP)超收敛算法以及基于该法的自适应求解方法推广到二维有限元线法分析中,至今已在二维Poisson方程和弹性力学平面问题中取得了令人满意的进展.该文旨在报道这些进展和成果.该文简要介绍了线法的EE...     

二维变分不等式问题的自适应有限元分析

有限元后处理超收敛计算的EEP(单元能量投影)法以及基于该法的自适应有限元分析已在一维变分不等式问题的求解中取得显著成功。以此为基础,该文对二维变分不等式问题成功地实现了自适应有限元分析。该文提出二维区域二分法和二维C 检验技术,有效地提升了松弛迭代的收敛速率,进而应用EEP 超收敛公式计算超收敛解答,用其检验误差并指导网格细分。该文给出的典型数值算例表明该文算法高效、稳定、精确,解答可逐点以最大模度量满足精度要求,堪称为数值精确解。     

自适应有限元线法在二维无穷域问题中的应用

无穷域问题广泛存在于实际工程中,半解析、半离散的数值计算方法—有限元线法（Finite ElementMethod of Lines,简称FEMOL）对其具有较好的适应性。在已有的映射型FEMOL无穷单元理论的基础上,基于单元能量投影（Element Energy Projection,简称EEP）法的自适应FEMOL被应用于二维无穷域问题的求解。用户只需输入稀疏的初始网格和误差限,算法即自动生成优化的FEMOL网格,该网格上常规单元和无穷单元的FEMOL解均按最大模度量满足给定误差限。文中首先介绍二维FEMOL的原理策略、无穷单元的构建,然后概述基于EEP法的自适应FEMOL算法,并讨论其对无穷域问题的适用性,之后对圆柱绕流的Poisson方程问题、带孔无穷大板单向拉伸的弹性力学平面问题、受圆形均布荷载半空间体的三维轴对称问题进行了自适应分析,最终不仅给出了满足误差限的函数（位移）解,也给出了具有优良性态的导数（应力）解,从而为无穷域问题的求解提供了一种高效可靠的新途径。     

膜结构极小曲面找形的一种自适应有限元分析

找形分析是膜结构设计中的关键环节,但在数学上,膜结构的极小曲面找形分析是一个高度非线性问题,一般无法求得其解析解,因此数值方法成为重要工具。近年来,基于单元能量投影法（EEP法）的一维非线性有限元的自适应分析已经取得成功,基于EEP法的二维线性有限元自适应分析也被证实是有效、可靠的。在此基础上,该文提出一种基于EEP法的二维非线性有限元自适应方法,并成功将之应用于膜结构的找形分析。其主要思想是,通过将非线性问题用Newton法线性化,引入现有的二维线性问题的自适应求解技术,进而实现二维有限元自适应分析技术从线性到非线性的跨越,将非线性有限元的自适应分析求解从一维问题拓展到二维问题。该方法兼顾求解的精度和效率,对网格自适应地进行调整,最终得到优化的网格,其解答可按最大模度量逐点满足用户设定的误差限。该文综述介绍了这一进展,并给出数值算例用以表明该方法的可行性和可靠性。     

二维自由振动问题的自适应有限元分析初探

自由振动反映结构动力特性,是抗震分析和结构设计的重要基础。近年来,基于单元能量投影(EEP)法的自适应有限元分析已在一系列线弹性及非线性问题中取得成功,而有限元线法(FEMOL)自适应分析在二维自由振动问题中的应用也被证实是有效的。在此基础上,该文进一步提出二维自由振动问题的自适应有限元分析方法。通过将特征值问题线性化,合理引入二维线性问题的EEP超收敛计算和自适应求解技术,该法可得到满足精度要求的自振频率和按最大模度量满足用户给定误差限的振型。该文以弹性薄膜为例,介绍了这一进展,并给出数值算例以表明该方法的有效性和可靠性。     

一维EEP自适应技术新进展:从线性到非线性

有限元后处理中超收敛计算的EEP(单元能量投影)法以及基于该法的自适应分析方法对线性ODE(常微分方程)问题的求解已经获得了全面成功,也推动了非线性ODE问题自适应求解的研究。经过研究,已经实现了一维有限元自适应分析技术从线性到非线性的跨越,该文意在对这方面的进展作一简要综述与报道。该文提出一种基于EEP法的一维非线性有限元自适应求解方法,其基本思想是通过线性化,将现有的线性问题自适应求解方法直接引入非线性问题求解,而无需单独建立非线性问题的超收敛计算公式和自适应算法,从而构成一个统一的、通用的非线性问题自适应求解算法。该文给出的数值算例表明所提出的算法高效、稳定、通用、可靠,解答可逐点按最大模度量满足用户给定的误差限,可作为先进高效的非线性ODE求解器的核心理论和算法。     

二维有限元线法超收敛解答计算的EEP法

有限元线法(FEMOL)是一种优良的半解析、半离散方法,但其解答存在解析方向和离散方向的精度不相称的弱点。本文提出将二维有限元线法比拟为广义一维问题的概念,遂可将新近提出的一维有限元超收敛计算的单元能量投影(EEP)法推广到二维有限元线法分析中。经有限元线法后处理中EEP超收敛计算而获得的解答,继承和保留了一维有限元中的出色表现,不但使任意一点的位移和应力的解答在两个方向具有相当的精度,而且都具有超收敛性质。文中以二维Poisson方程问题为例,具体给出了有限元线法EEP超收敛的公式,并给出了数值算例,用以表明本法的可行性和有效性。     

一维变分不等式问题的自适应有限元分析新探

结构工程中的弹性薄膜接触和杆件弹塑性扭转等问题是典型的变分不等式问题,对其高效精确求解,特别是满足给定精度要求下的自适应求解,是挑战性课题。该文作者新近成功实现了一维变分不等式问题的自适应有限元分析,该文对此进展作一报道。对于变分不等式的有限元求解,该文提出区域二分法和C检验技术,极大提升了松弛迭代的收敛速度,一般4次~5次线性解即可得到收敛的有限元解答,进而采用作者提出的EEP(单元能量投影)超收敛公式计算超收敛解答,用其检验误差并指导网格细分,逐步得到堪称为数值精确解的解答,亦即得到按照最大模度量逐点满足精度要求的解答。该文给出的数值算例表明所提出的算法具有高效、可靠、精确的优良特性。     

有限元线法EEP超收敛计算简约格式的再简约

有限元线法(FEMOL)是一种优良的半解析、半离散方法。为提高离散方向的精度,一维有限元中超收敛计算的单元能量投影(EEP)法,已成功地应用于有限元线法,导出一套简约格式的计算公式,该简约格式对位移及其导数都能给出比有限元线法解高出一阶的超收敛精度。该文对这套简约格式做进一步分析,发现其中的导数计算公式可以进一步简化,遂得到更加简约的导数计算公式。数值试验验证了该文方法的有效性。     

二维自适应技术新进展:从有限元线法到有限元法

二维有限元线法(FEMOL)的自适应分析已经取得成功,而且表现出色.然而,为了进一步推广应用领域,提高效率和效能,将其先进的自适应技术在最常用的有限元法(FEM)当中实现,便成为必然追求.经过近年的研究,已经基本实现了二维自适应分析技术从FEMOL到FEM的跨越,该文意在对这方面的进展作一简要综述与报道.从FEMOL出...     

有限元线法求解非线性模型问题——Ⅲ.薄膜的固有振动

本文是有限元线法(FEMOL)求解非线性模型问题的系列工作之三,对薄膜的固有振动这一特征值模型问题作了分析求解。文中,首先用FEMOL对特征值问题的泛函进行半离散,得到相应的常微分方程(ODE)特性值问题;然后,利用若干ODE变换技巧将问题转换成标准的非线性ODE问题;最后,采用一个新近研究出的有效算法,对各阶特征对进行了方便有效、精确可靠的求解。文中出示的算例展示了该法的功效。     

中厚圆柱壳自由振动的动力刚度法分析

该文阐述了将动力刚度法应用于中厚圆柱壳的自由振动分析。从考虑横向剪切变形和转动惯量的中厚壳理论出发,将圆柱壳的振动分解为一系列确定环向波数下的一维振动问题。用常微分方程求解器COLSYS求解该一维问题的动力刚度,通过Wittrick-Williams算法及导护型牛顿法求得该环向波数下结构的频率和振型。由于求解动力刚度时使用COLSYS对控制方程进行了精确求解,所以该文方法是精确方法。数值算例验证了中厚圆柱壳壳段固端频率计数J0计算方法的可靠性。综合表明:应用动力刚度法对中厚圆柱壳自由振动进行分析是可靠、精确的。     

中厚椭球壳自由振动动力刚度法分析

介绍精确动力刚度法分析中厚椭球壳自由振动具体实施方法,据环向波数不同将中厚椭球壳自由振动分解为一系列确定环向波数的一维振动;利用控制方程Hamilton形式建立动力刚度关系,用常微分方程求解器COLSYS求解控制方程获得单元动力刚度,用Wittrick-Williams算法求得该环向波数下椭球壳自振频率。数值算例给出中厚圆球壳及椭球壳不同边界条件的自振频率,验证动力刚度法高效、可靠、精确。     

壳体的有限元线法分析(Ⅱ)——数值算例

本文介绍了用前文所构造的任意曲面壳体的四边形有限元线法[1]单元所作的几个数值算例。算例表明该单元具有精度高、网格适应性好、厚薄通用的特点，采用p收敛技术可顺利克服闭锁，同时获得高精度的位移和内力，是求解壳体结构的一种有竞争力的半解析方法。     

时域自适应精细算法求解自激振动问题

应用时域自适应精细算法求解自激振动问题,将各物理量在离散时段内展开,可进行自适应递推计算,对非线性计算不做附加假设,避免了非线性迭代。数值分析表明,时域自适应精细算法可有效地用于自激振动问题的求解。