基于显式有限元方法的两相介质弹塑性动力反应计算分析

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性...     

金属板条在圆柱形模中弯曲成形的大变形有限元分析

本文采用大变形弹塑性有限元法对金属板条在柱形模中的压弯成形过程进行了数值模拟,并与实验进行了比较。首先给出了纠正的拉格朗日有限元公式和基于弹塑性乘法分解的超弹性塑性本构关系。对接触摩擦问题的处理采用了罚函数法。通过对数值结果的分析得出了一些对弯曲工艺的设计有指导价值的结论。     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

非均质饱和多孔介质弹塑性动力分析的广义耦合扩展多尺度有限元法

针对非均质饱和多孔介质弹塑性动力问题分析提出了一种广义耦合扩展多尺度有限元方法。首先,提出了基于细尺度等效刚度阵的粗尺度单元数值基函数构造方法,并给出了构造数值基函数的一般公式,所构造的耦合数值基函数有效考虑了动力相关效应与固液之间的耦合效应。其次,针对弹塑性非线性问题迭代求解,给出了基于摄动方法的位移与孔隙压强降尺度计算修正方案。最后,针对材料的强非均质特征,利用多节点粗单元技术来提高多尺度有限元方法的计算精度。通过与基于精细网格的传统有限元分析结果对比,验证了本文所提出方法的有效性与高效性。     

塑性节点法的研究及其在薄壳结构中的应用

本文采用变分原理及离散塑性流动定律假设研究了弹塑性有限元的新方法—塑性节点法及其所表示的物理意义。结果表明:该方法的物理实质是所述的塑性流动定理的离散。在一般意义下,本文对其进行了修正,典型板壳例题和复杂壳体结构的弹塑性有限元分析显示,计算结果与文献资料和实验数据符合较好。证明用文中方法进行壳体结构的弹塑性分析是行之有效的和可靠的。     

Dispersion error reduction for acoustic problems using the smoothed finite element method (SFEM)

The smoothed finite element method (SFEM), which was recently introduced for solving the mechanics and acoustic problems, uses the gradient smoothing technique to operate over the cell‐based smoothing domains. On the basis of the previous work, this paper reports a detailed analysis on the numerical dispersion error in solving two‐dimensional acoustic problems governed by the Helmholtz equation using the SFEM, in comparison with the standard finite element method. Owing to the proper softening effects provided naturally by the cell‐based gradient smoothing operations, the SFEM model behaves much softer than the standard finite element method model. Therefore, the SFEM can significantly reduce the dispersion error in the numerical solution. Results of both theoretical and numerical experiments will support these important findings. It is shown clearly that the SFEM suits ideally well for solving acoustic problems, because of the crucial effectiveness in reducing the dispersion error. Copyright © 2015 John Wiley & Sons, Ltd.     

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.     

Mixed finite element and differential quadrature method for free and forced vibration and buckling analysis of rectangular plates

This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.     

Spectral strip-element method for beams treated with active constrained layer damping

For vibration analysis of beams fully treated with active constrained layer damping (ACLD), a new approach called spectral strip-element method (SSEM) based on the spectral finite element method (SFEM) is proposed. It can avoid difficulties in solving the characteristic equation with higher orders and unknown parameters for wave numbers when using the SFEM; simultaneously, advantages of a very few elements and high accuracy of the SFEM are kept. A numerical example shows that the proposed method is very effective and reliable, compared with the exact solutions resulted from the spectral transfer matrix method (STMM).     

Effective numerical methods for elasto-plastic contact problems with friction

Several effective numerical methods for solving the elasto-plastic contact problems with friction are presented. First, a direct substitution method is employed to impose the contact constraint conditions on condensed finite element equations, thus resulting in a reduction by half in the dimension of final governing equations. Second, an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation, which distinguishes two kinds of nonlinearities, and makes the solution unique. In addition, Positive-Negative Sequence Modification Method is used to condense the finite element equations of each substructure and an analytical integration is introduced to determine the elasto-plastic status after each time step or each iteration, hence the computational efficiency is enhanced to a great extent. Finally, several test and practical examples are presented showing the validity and versatility of these methods and algorithms. The Project Supported by National Natural Science Foundation of China.     

Coupling finite difference method with finite particle method for modeling viscous incompressible flows

A hybrid approach to couple finite difference method (FDM) with finite particle method (FPM) (ie, FDM-FPM) is developed to simulate viscous incompressible flows. FDM is a grid-based method that is convenient for implementing multiple or adaptive resolutions and is computationally efficient. FPM is an improved smoothed particle hydrodynamics (SPH), which is widely used in modeling fluid flows with free surfaces and complex boundaries. The proposed FDM-FPM leverages their advantages and is appealing in modeling viscous incompressible flows to balance accuracy and efficiency. In order to exchange the interface information between FDM and FPM for achieving consistency, stability, and convergence, a transition region is created in the particle region to maintain the stability of the interface between two methods. The mass flux algorithm is defined to control the particle creation and deletion. The mass is updated by N-S equations instead of the interpolation. In order to allow information exchange, an overlapping zone is defined near the interface. The information of overlapping zone is obtained by an FPM-type interpolation. Taylor-Green vortices and lid-driven shear cavity flows are simulated to test the accuracy and the conservation of the FDM-FPM hybrid approach. The standing waves and flows around NACA airfoils are further simulated to test the ability to deal with free surfaces and complex boundaries. The results show that FDM-FPM retains not only the high efficiency of FDM with multiple resolutions but also the ability of FPM in modeling free surfaces and complex boundaries.     

厚壁圆筒自紧问题的弹粘塑性有限元分析

在双剪应力强度理论的基础上，应用弹粘塑性模型求解厚壁圆筒自紧问题，算例表明，这一算法是可行有效的。     

沥青路面弹塑性动力响应分析

采用弹塑性理论,建立了沥青路面弹塑性动力响应分析的三维有限元模型,利用有限元法分析了沥青路面的弹性和弹塑性动力响应、以及弹塑性状态下层间接触对沥青路面力学性能的影响．结果表明：在相同条件下,沥青路面为弹塑性状态时得到的弯沉和最大主应变均比弹性状态时大;卸载后,弹塑性状态时存在残余变形,说明沥青路面的弹塑性动力学响应分析得到的结果和路面实际情况较符合;沥青路面在弹塑性状态下,层间完全光滑时其弯沉是完全连续时的6 倍,上面层最大竖向应变是层间连续时的3.7 倍,下面层处最大竖向应变是层间连续时的2.3 倍;卸载后,层间完全光滑时,面层A 点与B 点均存在残余应变;随着层间摩擦系数的增大,路面弯沉值减少,说明在弹塑性状态下,层间接触状态对沥青路面的动力响应有较大影响．     

Geometrical nonlinear elasto-plastic analysis of tensegrity systems via the co-rotational method

An efficient finite element formulation is presented for geometrical nonlinear elasto-plastic analyses of tensegrity systems based on the co-rotational method. Large displacement of a space rod element is decomposed into a rigid body motion in the global coordinate system and a pure small deformation in the local coordinate system. A new form of tangent stiffness matrix, including elastic and elasto-plastic stages is derived based on the proposed approach. An incremental-iterative solution strategy in conjunction with the Newton-Raphson method is employed to obtain the geometrical nonlinear elasto-plastic behavior of tensegrities. Several numerical examples are given to illustrate the validity and efficiency of the proposed algorithm for geometrical nonlinear elasto-plastic analyses of tensegrity structures.     

A method of solving the fuzzy finite element equations in monosource fuzzy numbers

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｒｅｈａｖｅｂｅｅｎｍａｎｙｒｅｓｅａｒｃｈｐａｐｅｒｓａｂｏｕｔｔｈｅｆｕｚｚｙｓｔｏｃｈａｓｔｉｃｓｔｒｕｃｔｕｒｅ[1～ 3],ａｎｄｆｕｚｚｉｎｅｓｓａｎｄｒａｎｄｏｍｎｅｓｓａｒｅｔｗｏｉｍｐｏｒｔａｎｔｆａｃｔｏｒｓｉｎｅｎｇｉｎｅｅｒｉｎｇ .Ｂａｓｅｄｏｎｔｈｅｏｐｅｒａｔｉｏｎｒｕｌｅｓｏｆｆｕｚｚｙｎｕｍｂｅｒｓａｎｄｉｎｔｅｒｖａｌｎｕｍｂｅｒｓ,ｔｈｅｆｕｎｃｔｉｏｎｕｎｄｅｒｆｕｚｚｙａｎｄｒａｎｄｏｍｆａｃｔｏｒｓｃａｎｂｅｔｒａｎｓｐｏｓｅｄｉ…     

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.     

Solving material nonlinear structural problems by a finite analytic method in local coordinates

The purpose of this paper is to develop a finite analytic (FA) numerical solution for the elasto-plastic problem of the total theory. Schemes for the FA method in local coordinates for solving non-linear governing equations in the form of Navier equations are derived, which can be utilized to solve the problem in a domain of arbitrary geometry. Numerical illustration shows that the schemes are effective and practical.     

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method (FEM) based on a proper orthogonal decomposition (POD) theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save memory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element (FE) scheme is shown to be unconditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations (FDEs).     

A finite element method for the one‐dimensional extended Boussinesq equations

A new finite element method for Nwogu's (O. Nwogu, ASCE J. Waterw., Port, Coast., Ocean Eng., 119 , 618–638 (1993)) one‐dimensional extended Boussinesq equations is presented using a linear element spatial discretisation method coupled with a sophisticated adaptive time integration package. The accuracy of the scheme is compared to that of an existing finite difference method (G. Wei and J.T. Kirby, ASCE J. Waterw., Port, Coast., Ocean Eng., 121 , 251–261 (1995)) by considering the truncation error at a node. Numerical tests with solitary and regular waves propagating in variable depth environments are compared with theoretical and experimental data. The accuracy of the results confirms the analytical prediction and shows that the new approach competes well with existing finite difference methods. The finite element formulation is shown to enable the method to be extended to irregular meshes in one dimension and has the potential to allow for extension to the important practical case of unstructured triangular meshes in two dimensions. This latter case is discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

FEM-FDM coupled liquefaction analysis of a porous soil using an elasto-plastic model

The phenomenon of liquefaction is one of the most important subjects in Earthquake Engineering and Coastal Engineering. In the present study, the governing equations of such coupling problems as soil skeleton and pore water are obtained through application of the two-phase mixture theory. Using au-p (displacement of the solid phase-pore water pressure) formulation, a simple and practical numerical method for the liquefaction analysis is formulated. The finite difference method (FDM) is used for the spatial discretization of the continuity equation to define the pore water pressure at the center of the element, while the finite element method (FEM) is used for the spatial discretization of the equilibrium equation. FEM-FDM coupled analysis succeeds in reducing the degrees of freedom in the descretized equations. The accuracy of the proposed numerical method is addressed through a comparison of the numerical results and the analytical solutions for the transient response of saturated porous solids. An elasto-plastic constitutive model based on the non-linear kinematic hardening rule is formulated to describe the stress-strain behavior of granular materials under cyclic loading. Finally, the applicability of the proposed numerical method is examined. The following two numerical examples are analyzed in this study: (1) the behavior of seabed deposits under wave action, and (2) a numerical simulation of shaking table test of coal fly ash deposit.