有限元方法中的光滑积分伪弱形式

提出了一种光滑积分伪弱形式,将光滑积分拓展至被积函数非偏导项求解。结合光滑应变技术和伪弱形式,可实现有限元系统方程统一光滑积分求解,即对刚度矩阵和质量矩阵中的应变矩阵和形函数矩阵均可进行光滑积分处理,并转化为光滑子域的边界积分。光滑积分伪弱形式与光滑应变技术比较,增加了形函数矩阵不定积分处理过程,且没有降低有限元求解对形函数连续性的要求。不过,伪弱形式改变了单元积分的求解形式,连续质量矩阵求解也无需坐标映射和雅可比矩阵计算。以轴对称二维问题为研究对象,结果表明极度不规则三角形和四边形单元光滑积分伪弱形式在静态和动态有限元方程求解中也具有很好的精度。     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

基本解方法与边界节点法求解Helmholtz方程的比较研究

基本解方法和边界节点法是基于径向基函数的两种重要无网格边界离散数值技术。针对Helmholtz方程,本文比较研究这两种数值方法在不同计算区域问题上的计算精度、插值矩阵对称性、病态性及计算成本。数值试验结果表明,两种方法都可以有效求解边界数据准确的Helmholtz问题。在数值离散过程中,两种方法都可以通过调整配置点的位...     

三维Euler方程的分区和并行计算

三维全机绕流区域分解成多块子区域,多块区域之间采用迎风型通量守恒内边界耦合条件,分区计算总体区域,形成总体耦合流场的分区数值解。利用PVM并行环境,采用纯结点并行计算编程方式和“先进先出”的同步控制等待机制,对三维复杂流动跨音速流场相应分区实现了多区域并行计算。分析了影响并行效率的主要因素,将并行计算结果与串行计算结果和实验结果作了比较,讨论了多种区域分解数目的并行计算效率。在负载平衡程度较好时,可得到较高的并行效率。     

连续胶体材料化学-力学耦合有限元理论及其在ABAQUS中UEL程序实现

基于哈密顿原理,得到水凝胶的化学-力学耦合控制方程的等效积分形式和有限元形式。在整体坐标系下推导出用形函数表示的化学-力学耦合应变矩阵和单元刚度矩阵,并且得到在局部坐标系下的离散化形式。结合ABAQUS软件,编制了用户单元子程序UEL,通过数值算例验证了所开发单元的正确性,为在ABAQUS软件中实现各种耦合问题的有限元UEL编程提供了参考依据。     

两个新型的八结点块体单元—三维离散算子差分法

给出了弹性力学三维问题的离散算子差分法 ,讨论离散算子差分法在三维问题中的特点 ,意在为该方法的进一步发展提供依据 ,为应用弱形式进行数值求解的研究提供参考。本文从弹性力学平衡方程更为一般的弱形式出发 ,给出了含边界参数的弱形式方程。由该方程不仅可以得到有限元法 ,还可得到离散算子差分法。给出了两个八结点块体单元 ,虽然单元中位移函数是非协调的 ,不需特殊处理便可保证离散格式收敛 ,并对单元位移有十分好的反映能力。     

拟协调有限元弱形式的辛算法

基于弱形式的力学方程,阐述了弱形式广义方程是拟协调有限元的内在本质,用弱形式给出的微分方程和边界条件根本上是降低了函数光滑性,不过对工程问题而言,给出的有限元解比原始方程更接近真实解,其数值解就是广义协调方程的直接解,同时满足平衡和几何方程弱连续条件。进而就导出的对偶体系弱形式哈密尔顿方程,采用辛相似变换,利用平方约化法求解哈密尔顿矩阵特征值问题,使其哈密尔顿结构得到了保证。辛算法具有较强的有效性,可以解决常规有限元难以适应的领域,对计算力学发展有着重要的作用。     

基于点插值的配点型无网格法解Helmholtz问题

基于点插值法的思想,用三角函数作为基函数在局部支持域内构造具有Kroneckerδ函数性、单位分解性、高阶连续性、再生性和紧支性的形函数.用配点法离散微分方程,得到了具有稀疏带状性的系数矩阵,用GMERS方法求解代数方程组,分别研究了Helmholtz问题的边界层问题和波传播问题.通过数值算例可以发现,给出的数值结果非常接近于精确解,且随着节点的增加,其精确度越来越高,具有良好的收敛性.     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。     

通过Adomian分解法求解二维Helmholtz方程

提出基于Adomian分解法求解二维Helmholtz方程。通过Adomian分解法可以把Helmholtz微分方程和边界条件分别转换成递归代数公式和适用符号计算的简单代数公式。利用边界条件可以很容易得到方程的解析解表达式。Adomian分解法的主要特点在于计算简单快速,并且不需要进行线性化或离散化。最后给出数值实例以验证Adomian分解法求解二维Helmholtz方程的有效性。通过数值计算可以发现,基于Adomian分解法的计算结果非常接近精确解,并且该方法具有良好的收敛性。这表明Adomian分解法能够快速有效求解Helmholtz方程。