各向异性两相材料尖劈奇性场的非协调元分析

提出了一个基于位移的、分析柱状各向异性两相材料尖劈端部邻域的奇性位移场和应力场问题的非协调元特征分析法. 该方法从柱状扇区的散度定理出发，将柱状扇区控制方程的弱式化为一个与虚功原理相同形式的方程，采用一种新的非协调元技术把所导出的``虚功原理'转化为标准一阶特征方程的求解问题. 非协调元法中，尖劈端部邻域的位移场假定没有采用奇异变换技术，有限元的单元形式是一维的. 将柱状各向异性两相材料尖劈视为``广义平面应变'问题，位移场与坐标z无关，只关注界面端的幂奇异性而不考虑对数奇异性. 运用该方法给出了柱状各向异性两相材料尖劈端部奇性应力指数、奇性位移角分布和应力角分布的算例. 所有的计算结果表明，该方法使用的单元少而且精度较高.     

The finite element method of singular perturbation problem

In this paper we construct new finite element subspace using polynomials of different degrees and the new finite element scheme is established. The convergence of the scheme and the stability of the reduced difference equation are proved.     

Superconvergence analysis of lower order anisotropic finite element

The convergence analysis of the lower order nonconforming element pro- posed by Park and Sheen is applied to the second-order elliptic problem under anisotropic meshes.The corresponding error estimation is obtained.Moreover,by using the interpo- lation postprocessing technique,a global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is derived.Numerical results are also given to verify the theoretical analysis.     

Discontinuous Galerkin finite element discretization of a strongly anisotropic diffusion operator

The discretization of a diffusion equation with a strong anisotropy by a discontinuous Galerkin finite element method is investigated. This diffusion term is implemented in the tracer equation of an ocean model, thanks to a symmetric tensor that is composed of diapycnal and isopycnal diffusions. The strong anisotropy comes from the difference of magnitude order between both diffusions. As the ocean model uses interior penalty terms to ensure numerical stability, a new penalty factor is required in order to correctly deal with the anisotropy of this diffusion. Two penalty factors from the literature are improved and established from the coercivity property. One of them takes into account the diffusion in the direction normal to the interface between the elements. After comparison, the latter is better because the spurious numerical diffusion is weaker than with the penalty factor proposed in the literature. It is computed with a transformed coordinate system in which the diffusivity tensor is diagonal, using its eigenvalue decomposition. Furthermore, this numerical scheme is validated with the method of manufactured solutions. It is finally applied to simulate the evolution of temperature and salinity due to turbulent processes in an idealized Arctic Ocean. Copyright © 2014 John Wiley & Sons, Ltd.     

A precise singular finite element for crack analysis

The multi-variable finite element algorithm based on the generalized Galerkin’smethod is more flexible to establish a finite element model in the continuum mechanies.Byusing this algorithm and numerical tests a new singular finite element for elasto-plasticfracture analysis has been formulated.The results of numerical tests show that the newelement possesses high accuracy and good performance.Some rules for formulating amulti-variable singular finite element are also discussed in this paper.     

A finite element implementation of the stress gradient theory

In this contribution, a finite element implementation of the stress gradient theory is proposed. The implementation relies on a reformulation of the governing set of partial differential equations in terms of one primary tensor-valued field variable of third order, the so-called generalised displacement field. Whereas the volumetric part of the generalised displacement field is closely related to the classic displacement field, the deviatoric part can be interpreted in terms of micro-displacements. The associated weak formulation moreover stipulates boundary conditions in terms of the normal projection of the generalised displacement field or of the (complete) stress tensor. A detailed study of representative boundary value problems of stress gradient elasticity shows the applicability of the proposed formulation. In particular, the finite element implementation is validated based on the analytical solutions for a cylindrical bar under tension and torsion derived by means of Bessel functions. In both tension and torsion cases, a smaller is softer size effect is evidenced in striking contrast to the corresponding strain gradient elasticity solutions.

     

Finite element implementation of a generalised non-local rate-dependent crystallographic formulation for finite strains

This work describes the finite element implementation of a generalised strain gradient and rate-dependent crystallographic formulation for finite strains and general anisothermal conditions based on a multiplicative decomposition of the deformation gradient. The implementation involved the development of both a novel finite element formulation to determine the spatial slip rate gradients at each material point, and an implicit numerical integration scheme at the constitutive level to update the stresses and solution dependent variables. The time-integration procedure uses a Newton–Raphson scheme with a single level of iteration to solve the incremental non-linear equations associated with the non-local constitutive formulation. Closed-form solutions for the relevant fourth-order Jacobian tensors are given. The proposed numerical scheme is formulated in a general form and hence should be applicable to most existing crystallographic models. The crystallographic formulation is then used to investigate the effect of the morphology and volume fraction of the reinforcing phase of a two-phase single crystal on its macroscopic behaviour.     

Numerical implementation of pseudo-arclength algorithm in nonlinear finite element analysis

A numerical procedure to calculate the pre-buckling and postbuckling response of general structures is presented. This procedure is based on the pseudo-arclength algorithm suggested by E. Riks et al., which has some numerical difficulties during implementation of large applied analysis programs. To overcome these difficulties, a scheme based on rank-1 modification of the matrix is proposed. Some example show this procedure behaves well in passing through the limit point and is rather efficient. recommended by Prof. Chen Yaosong and Prof. Wu Jike     

基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现

依据有限元理论,结合内聚力模型法则,推导出二维二次粘结界面单元在大位移情况下的数值格式,得到用形函数表示的单元位移模式、载荷向量和刚度矩阵,并进行了离散化。基于ABAQUS软件的自定义扩展模块,编制了相应的用户单元子程序UEL,通过数值算例验证了该程序的准确性和有效性。这一成果能为在ABAQUS软件中开展相关数值研究,以及开发其他类型的内聚力界面有限单元提供思路和参考。     

Improved Carroll's hyperelastic model considering compressibility and its finite element implementation

Acta Mechanica Sinica - In engineering component design, material models are increasingly used in finite element simulations for an expeditious and less costly analysis of the design prototypes. As...     

轴对称弹性应变梯度理论公式推导及有限元实现

用张量运算推导了弹性应变梯度轴对称问题的基本公式.建立了应变梯度轴对称不协调元的弱连续条件,进一步建立了满足弱连续条件的应变梯度轴对称18-DOF三角形单元(BCIZ+ART9),其中BCIZ满足线性应变C0连续,用于计算应变ε;ART9满足常曲率C1弱连续,用于计算应变梯度η0数值结果表明该单元通过C0-1分片检验并能体现材料的尺度效应.     

A novel singular finite element of mixed-mode crack problems with arbitrary crack tractions

A novel singular finite element is presented to study cracked plates with arbitrary traction acting on crack surfaces. Firstly, the analytical solution around crack tips is determined using the symplectic dual approach. Subsequently, the solution is used to develop a novel singular finite element, which depicts accurately the characteristic of singular stresses field near crack tips. And the novel element can be applied to solve cracked plates, and both Mode I and Mode II stress intensity factors can be determined directly and accurately. Lastly, two numerical examples are given to illustrate the present method.     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

各向异性位势问题边界元法中几乎奇异积分的解析算法

导出了一种解析积分算法,精确计算了二维各向异性位势问题边界元法中近边界点的几乎奇异积分。对线性单元,几乎奇异积分可用解析公式直接计算。对二次单元,可将其细分为几个线性单元,采用该解析公式间接近似计算。当内点离积分单元较远时,仍然保持常规高斯数值积分模式;而当内点离其较近时,高斯积分结果失效,采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元比线性元计算结果更精确。     

智能结构有限元动力模型的建立及主动振动控制和抑制

采用一种新的压电板单元，建立了含有分布压电传感元件和执行元件结构（智能结构）的有限元动力模型。利用两种反馈控制律，研究了智能结构振动控制与抑制的问题，并提出了智能结构主动振动控制和抑制的一种方法。最后，提供了数值示例，说明本文提出方法的应用。     

A singular finite element for Stokes flow: The stick–slip problem

Abrupt changes in boundary conditions in viscous flow problems give rise to stress singularities. Ordinary finite element methods account effectively for the global solution but perform poorly near the singularity. In this paper we develop singular finite elements, similar in principle to the crack tip elements used in fracture mechanics, to improve the solution accuracy in the vicinity of the singular point and to speed up the rate of convergence. These special elements surround the singular point, and the corresponding field shape functions embody the form of the singularity. Because the pressure is singular, there is no pressure node at the singular point. The method performs well when applied to the stick–slip problem and gives more accurate results than those from refined ordinary finite element meshes.     

On the implementation of anisotropic yield functions into finite strain problems of sheet metal forming

In this article the implementation of anisotropic yield functions into finite element investigations of orthotropic sheets with planar anisotropy is discussed within a plane-stress context. Special attention is focused on the proper treatment of the orientation of the anisotropic axes during deformation into the finite-strain range. As an example problem the hydrostatic bulging of a membrane is considered in conjunction with a recently proposed anisotropic yield function. It is shown that the aspects of the plane-stress assumption, which do not come into consideration in isotropic analyses, can play an important role on the accuracy of the solution when the rotation of the orthotropic axes enters the computation directly due to the presence of material anisotropy. When the material anisotropy is considered and when the deformation of the workpiece is not limited to the plane of the undeformed sheet (such as cup drawing, hydrostatic bulging, etc.), the numerical experiments indicate that the only correct formulation is the one based on numerically imposing the requirement that for the plane-stress application, the in-plane material axes have to remain in the plane of the sheet during the deformation.