一种新的求解对流占优问题的自适应网格细化方法

在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。     

数值流形单元法数学网格自适应

基于数值流形方法和有限覆盖技术,将有限元法的后验误差估计理论及h型网格自适应技术推广应用到数值流形单元法中,提出了数值流形单元法的后验误差估计方法和数学网格自适应技术,并编制了相应的程序。数值算例表明,经过网格自适应,可以在粗糙的初始网格基础上得到质量比较理想的网格,计算结果可达到用户要求的精度。     

改进的Z^2应力恢复过程与h型自适应有限元分析

建议了一种较为精确的边界应力求解方法，并用于改进Ｚｉｅｎｋｉｅｗｉｃｚ－Ｚｈｕ（Ｚ＾２）应力恢复过程。改进过程增加的计算量不大，但可有效地改善后验误估计精度。ｈ型自适应有限元分析结果表明，改进过程更有利于最优网格寻求工作。     

对流占优问题的无网格自适应布点方案

应用常规数值方法求解对流占优的对流扩散方程时会出现非物理的数值伪振荡现象.因此本文提出了一种基于无网格径向点插值法的自适应布点方案,并成功地解决了对流占优时的数值伪振荡问题.在自适应布点的实施过程中,该方案将无网格方法中的背景积分单元作为自适应控制的梯度计算单元,并将该控制单元场函数梯度的大小作为自适应的梯度控制指标,然后给定相应的梯度控制限,通过控制指标和梯度限的比较来指示高梯度区域进行自适应中心加点和梯度计算单元的分解.数值结果表明:这种基于无网格径向点插值法的自适应布点方案不仅能有效地消除对流占优时的数值伪振荡现象,而且它还具有计算精度高、数值稳定性好、算法实施简单、前后处理方便的优点.     

无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。     

求解对流占优问题的无网格自适应稳定化方案

应用标准的无网格方法求解对流占优问题时会出现非物理的数值伪振荡现象,采用MF-SUPG、MFGLS、MFSGS等稳定化方法可以有效地消除数值伪振荡.因此本文基于无网格径向点插值法提出了一种自适应布点方案,并分别与MFSUPG、MFGLS、MFSGS方法相结合.数值模拟表明:当扩散系数较小时,三种稳定化方法均可以有效地消除对流占优问题大部分区域的数值伪振荡,但稳定化后其解在边界处仍有振荡存在,而结合自适应方案后的三种稳定化方法均可以彻底地消除数值伪振荡,且具有计算精度高、稳定性好、算法实施简单、前后处理方便.     

积分背景网格对自适应EFG法拓扑优化的影响研究

采用应力能量范数作为误差指标,探讨了EFG法中积分背景网格对计算精度的影响,得到了合理划分背景网格的建议;建立了以节点密度为设计变量、以最小化柔度为优化日标的拓扑优化模型。采用以节点密度值为加点判据的自适应规则加点方案,开展了连续体结构的拓扑优化研究,该加点方案能有效地减少设计变量的个数,探讨了背景网格对拓扑优化结果的影响。算例结果表明,采用合适的背景网格不仅能进一步减少设计变量的个数,而且能够改善拓扑优化结果的光滑性,使计算效率和精度得到提高。     

粘聚裂纹扩展的强化有限元h型网格自适应模拟

非连续变形分析和非规则节点处理是基于单元细划的粘聚裂纹扩展网格自适应模拟的关键。首先,利用强化有限单元法中数学单元和物理单元分离的特点,通过引入过渡单元,将适用于非连续变形描述的数学模式覆盖法和方便处理非规则节点的物理模式重构法结合,提出了强化有限单元法的统一关联法则,并导出了相应的单元列式。其次,基于数学裂纹尖端影响域和裂尖单元尺寸,提出了基于强化有限单元法的粘聚裂纹扩展过程模拟的h型网格自适应策略。最后,通过两个算例验证了本文方法的合理性和有效性。     

含裂纹体的单位分解扩展无网格方法

不连续体的数值模拟尤其是动态裂纹的追踪问题一直是工程界研究的热点和难点问题。无网格方法仅仅需要结点信息,非常适合于求解这类问题。基于单位分解思想,在移动最小二乘近似函数(MLS)中根据裂纹面的不连续位移增加一个Heaviside函数,在裂尖则增加四个扩展函数描述渐进裂纹位移场;应用Galerkin方法推导了平衡方程的离散线性方程,并给出了求解裂纹问题应力强度因子的计算公式。与其他类型的扩展无网格相比,在裂尖处近似函数不需要使用可视准则,很容易生成r1/2奇异;另一个优势是影响域并没有因为裂纹的存在而改变,不会降低方程的稀疏性,求解效率较高。数值算例表明,该方法能方便有效地模拟不连续问题,具有十分广阔的应用空间。     

改进的Z~2应力恢复过程与h型自适应有限元分析

建议了一种较为精确的边界应力求解方法,并用于改进Ｚｉｅｎｋｉｅｗｉｃｚ－Ｚｈｕ（Ｚ２）应力恢复过程。改进过程增加的计算量不大,但可有效地改善后验误差估计精度。ｈ型自适应有限元分析结果表明,改进过程更有利于最优网格寻求工作     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

MESHLESS METHOD FOR 2D MIXED-MODE CRACK PROPAGATION BASED ON VORONOI CELL

A meshless method integrated with linear elastic fracture mechanics (LEFM) is presented for 2D mixed-mode crack propagation analysis. The domain is divided automatically into sub-domains based on Voronoi cells, which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path. For each increment of the crack extension, the meshless method is applied to carry out a stress analysis of the cracked structure. The J-integral, which can be decomposed into mode I and mode II for mixed-mode crack, is used for the evaluation of the stress intensity factors (SIFs). The crack-propagation direction, predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method. The meshless results agree well with the experimental ones, which validates the accuracy and efficiency of the method.     

ADAPTIVE MESHLESS METHOD BASED ON LOCAL FIT TECHNOLOGY

An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement.New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.     

任意梯度分布功能梯度涂层平面裂纹分析

提出可以分析任意梯度功能梯度材料的分层模型,并采用该模型研究功能梯度涂层平面裂纹问题.采用Fourier变换和传递矩阵法将该混合边值问题化为奇异积分方程组,通过数值求解获得应力强度因子.考察了分层模型的有效性,以及材料梯度变化形式、结构几何尺寸和材料梯度参数对裂纹应力强度因子的影响,发现结构几何尺寸、材料梯度变化形式、...     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题,推导了相应的离散方程. 与 其它基于网格的方法相比,有限点法采用移动最小二乘法构造形函数,只需要节点信息,不 需要划分网格,用配点法离散控制方程,可以直接施加边界条件,不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

功能梯度材料平面问题的辛弹性力学解法

将辛弹性力学解法推广用于功能梯度材料平面问题的分析,考虑沿长度方向弹性模量为指数函数变化而泊松比为常数的矩形域平面弹性问题,给出了具体的求解步骤. 提出了移位Hamilton矩阵的新概念,建立起相应的辛共轭正交关系;导出了对应特殊本征值的本征解,发现材料的非均匀特性使特殊本征解的形式发生明显的变化.      

A NEW IHIERARCHICAL BOUNDARY ELEMENT METHOD AND ITS ADAPTIVE PROCESSES FOR PLATE BENDING PROBLEM

In this paper a new hierarchical boundary element method is introducedfor solving the problem of plate bending.Exact solutions of the governing equationsare used inside the domain together with independent deflections and rotations on theboundary.A generalized variational principle is employed to achieve the boundaryelement formulation.Further,the adaptive processes of the method are also discussed.By virtue of the error estimate technique proposed by Zienkiewicz et al.areasonable error indicator and adaptive scheme are suggested.Numerical examplesillustrate the high accuracy of the new elements and show the excellent efficiency ofthe adaptive computation described in this paper.     

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.     

MESHLESS METHOD OF DUAL RECIPROCITY HYBRID RADIAL BOUNDARY NODE METHOD FOR ELASTICITY

Combining the radial point interpolation method （RPIM）, the dual reciprocity method （DRM） and the hybrid boundary node method （HBNM）, a dual reciprocity hybrid radial boundary node method （DHRBNM） is proposed for linear elasticity. Compared to DHBNM, RPIM is exploited to replace the moving least square （MLS） in DHRBNM, and it gets rid of the deficiency of MLS approximation, in which shape functions lack the delta function property, the boundary condition can not be applied easily and directly and it＇s computational expense is high. Besides, different approximate functions are discussed in DRM to get the interpolation property, in which the accuracy and efficiency for different basis functions are compared. Then RPIM is also applied in DRM to replace the conical function interpolation, which can greatly improve the accuracy of the present method. To demonstrate the effectiveness of the present method, DHBNM is applied for comparison, and some numerical examples of 2-D elasticity problems show that the present method is much more effective than DHBNM.     

ON THE PLANE PIEZOELECTRIC PROBLEM OF A LOADED CRACK TERMINATING AT A MATERIAL INTERFACE

The plane problem of a crack terminating at the interface of a bimaterial piezoelectric,and loaded on its faces, is treated. The emphasis is placed on how to transform this problem into anon-homogeneous Hilbert problem. To make the derivation tractable, the concept of the axial conju-gate is introduced and related to the complex conjugate. The angle between the crack line and the inter-face may be arbitrary. Numerical results are given to illustrate the stress singularity at crack tip.