基于余能原理的有限变形问题有限元列式

利用基面力概念,推导了一种基于余能原理的有限变形问题显式有限元列式,可应用于结构的大位移、大转动问题。以基面力为状态变量来表达单元的余能,将有限变形情况下的单元余能分解为变形余能部分和转动余能部分,利用Lagrange乘子法推导出余能原理有限元的控制方程,编制了相应的非线性有限元程序。通过算例分析,说明该列式和程序的可靠性和精确性。     

浅曲梁大挠度分析的余能有限元方法

以高玉臣提出的弹性大变形余能原理为基础,利用Lagrange乘子,放松平衡方程和力边界条件对余能泛函的约束,推导出广义的余能原理.根据极分解定理,将变形分为刚性转动和纯变形两部分,则余能也包含相应的两部分,一部分与刚性转动有关,而另一部分与纯变形有关.使用线弹性本构关系,建立了可用于几何非线性计算的有限元模型.应用更新的Lagrange列式法,给出了增量形式的有限元公式.数值计算结果表明,该方法可用于浅曲粱的几何大变形计算.     

Discrete formulation of mixed finite element methods for vapor deposition chemical reaction equations

The vapor deposition chemical reaction processes, which are of extremely extensive applications, can be classified as a mathematical model by the following governing nonlinear partial differential equations containing velocity vector, temperature field, pressure field, and gas mass field. The mixed finite element (MFE) method is employed to study the system of equations for the vapor deposition chemical reaction processes. The semidiscrete and fully discrete MFE formulations are derived. And the existence and convergence (error estimate) of the semidiscrete and fully discrete MFE solutions are demonstrated. By employing MFE method to treat the system of equations for the vapor deposition chemical reaction processes, the numerical solutions of the velocity vector, the temperature field, the pressure field, and the gas mass field can be found out simultaneously. Thus, these researches are not only of important theoretical means, but also of extremely extensive applied vistas.     

幂律全塑性罚函数随机有限元

用罚函数有限元方法解决平面应变下的体积下可压缩问题,应用摄动有限元理论发展幂律全塑性随机有限元,并以弹性模是E,泊松比u和节点坐标等的基本随机变量,推导有限元列式,给出单板受拉,梁受弯的算例。     

A new finite element formulation for both compressible and nearly incompressible fluid dynamics

A new finite element formulation designed for both compressible and nearly incompressible viscous flows is presented. The formulation combines conservative and non‐conservative dependent variables, namely, the mass–velocity (density * velocity), internal energy and pressure. The central feature of the method is the derivation of a discretized equation for pressure, where pressure contributions arising from the mass, momentum and energy balances are taken implicitly in the time discretization. The method is applied to the analysis of laminar flows governed by the Navier–Stokes equations in both compressible and nearly incompressible regimes. Numerical examples, covering a wide range of Mach number, demonstrate the robustness and versatility of the new method. Copyright © 2000 John Wiley & Sons, Ltd.     

The random variational principle in finite deformation of elasticity and finite element method

ＴＨＥＲＡＮＤＯＭＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＩＮＦＩＮＩＴＥＤＥＦＯＲＭＡＴＩＯＮＯＦＥＬＡＳＴＩＣＩＴＹＡＮＤＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＧａｏＨａｎｇ－ｓｈａｎ（高行山）（ＮｏｒｔｈｗｅｓｔｅｎＰｏｌｙｔｅｃｈｎｉｃａｌＵ...     

A variational principle and finite element formulation for multi-physics PLZT ceramics

This research presents an extended variational principle and a finite element formulation for multi-physics analysis of PLZT (lanthanum zirconate titanate) ceramics by including the photovoltaic and optothermic effects. The photo-induced electrical and thermal, and mechanical fields, as well as the heat source effect due to light illumination are all considered in the formulation. A generalized variation approach for advanced multi-physics PLZT problems is developed and then the relevant finite element formulation is established.     

An energy stable monolithic Eulerian fluid‐structure finite element method

When written in an Eulerian frame, the conservation laws of continuum mechanics are similar for fluids and solids leading to a single set of variables for a monolithic formulation. Such formulations are well adapted to large displacement fluid‐structure configurations, but stability is a challenging problem because of moving geometries. In this article, the method is presented; time implicit discretizations are proposed with iterative algorithms well posed at each step, at least for small displacements; stability is discussed for an implicit in time finite element method in space by showing that energy decreases with time. The key numerical ingredient is the Characterics‐Galerkin method coupled with a powerful mesh generator. A numerical section discusses implementation issues and presents a few simple tests. It is also shown that contacts are easily handled by extending the method to variational inequalities. This paper deals only with incompressible neo‐Hookean Mooney‐Rivlin hyperelastic material in 2 dimensions in a Newtonian fluid, but the method is not limited to these; compressible and 3D cases will be presented later.     

A finite element formulation for laminar flow of a fluid with microstructure

A finite element formulation for the steady laminar flow of an incompressible fluid with microstructure has been developed. The particular fluids considered are commonly known as micropolar fluids, in which case suspended particulate microstructures are modelled by an ‘extended’ continuum formulation. The particle microspin is a new kinematic variable which is independent of the classical vorticity vector and thereby allows relative rotation between particles and the surrounding fluid. This formulation also gives rise to couple stresses in addition to classical force or traction stresses. The finite element formulation utilizes a variational approach and imposes conservation of mass through a penalty function. A general boundary condition for microspin has been incorporated whereby microspin at a solid boundary is constrained to be proportional to the fluid vorticity. The proportionality constant in this case can vary from zero to unity. Sample solutions are presented for fully developed flow through a straight tube and compared with an analytical solution. Results are also generated for flow through a constricted tube and compared with a Newtonian fluid solution.     

On the stability of bubble functions and a stabilized mixed finite element formulation for the Stokes problem

In this paper we investigate the relationship between stabilized and enriched finite element formulations for the Stokes problem. We also present a new stabilized mixed formulation for which the stability parameter is derived purely by the method of weighted residuals. This new formulation allows equal‐order interpolation for the velocity and pressure fields. Finally, we show by counterexample that a direct equivalence between subgrid‐based stabilized finite element methods and Galerkin methods enriched by bubble functions cannot be constructed for quadrilateral and hexahedral elements using standard bubble functions. Copyright © 2008 John Wiley & Sons, Ltd.     

Base force element method of complementary energy principle for large rotation problems

Using the concept of the base forces, a new finite element method （base force element method, BFEM） based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.     

Broadband energy exchanges between a dissipative elastic rod and a multi-degree-of-freedom dissipative essentially non-linear attachment

We analyze complex, multi-frequency, non-linear modal interactions in the damped dynamics of a viscously damped dispersive finite rod coupled to a multi-degree-of-freedom essentially non-linear attachment. We perform a parametric study to show that the attachment can be an effective broadband energy absorber and dissipater of shock energy from the rod. It is shown that strong targeted energy transfer from the rod to the attachment occurs when there is strong stiffness asymmetry in the attachment. For weak viscous dissipation, a clear understanding of dynamical transitions in the integrated rod-non-linear attachment system can be gained by wavelet transforming the time series and superimposing the resulting wavelet spectra in the frequency-energy plot (FEP) of the periodic orbits of the underlying Hamiltonian system. Two distinct NES configurations are analyzed in detail, and their damped responses are analyzed by the Hilbert-Huang transform (HHT). We show that the HHT is capable of analyzing even complex non-linear damped transitions, by providing the dominant frequency components (or equivalently, time scales) at which the non-linear phenomena take place, and clarifying the series of non-linear resonance captures between the rod and attachment dynamics that are responsible for the broadband energy exchanges in this system.     

On a new extended finite element method for dislocations: Core enrichment and nonlinear formulation

A recently developed finite element method for the modeling of dislocations is improved by adding enrichments in the neighborhood of the dislocation core. In this method, the dislocation is modeled by a line or surface of discontinuity in two or three dimensions. The method is applicable to nonlinear and anisotropic materials, large deformations, and complicated geometries. Two separate enrichments are considered: a discontinuous jump enrichment and a singular enrichment based on the closed-form, infinite-domain solutions for the dislocation core. Several examples are presented for dislocations constrained in layered materials in 2D and 3D to illustrate the applicability of the method to interface problems.     

Materially and geometrically non-linear woven composite micro-mechanical model with failure for finite element simulations

A computational micro-mechanical material model of woven fabric composite material is developed to simulate failure. The material model is based on repeated unit cell approach. The fiber reorientation is accounted for in the effective stiffness calculation. Material non-linearity due to the shear stresses in the impregnated yarns and the matrix material is included in the model. Micro-mechanical failure criteria determine the stiffness degradation for the constituent materials. The developed material model with failure is programmed as user-defined sub-routine in the LS-DYNA finite element code with explicit time integration. The code is used to simulate the failure behavior of woven composite structures. The results of finite element simulations are compared with available test results. The model shows good agreement with the experimental results and good computational efficiency required for finite element simulations of woven composite structures.     

三晶体中三晶交线附近弹性应变能集中的有限单元分析

提出应变能集中这一概念,对三晶体中三晶交线附近的弹性应变能的集中特性进行了计算,发现当加载方向平行于三晶交线时,三晶交线附近有弹性应变能集中,但是集中程度不是很高,当加载方向垂直于三晶交线时,根据三晶体中组元晶体的相对位置的不同和取向的不同,三晶交线附近或者没有弹性应变能集中,或者有弹性应变能集中,而且在一些情况下与平行于三晶交线方向加载的情况相比,弹性应变能的集中程度要高很多。     

Geometrically non-linear analysis of arbitrary elastic supported plates and shells using an element-based Lagrangian shell element

In the present paper, the ELF (element-based Lagrangian formulation) 9-node ANS (assumed natural strain) shell element was combined with the spring element for geometrically non-linear analysis of plates and shells sustained by arbitrary elastic edge supports that are subjected to variation in loading.This particular spring element serves as tool for modeling an arbitrary elastic edge support with 6 DOF (degrees of freedom). The elastic edge support was modeled by combining different spring models. The ANS method was used to overcome shear and membrane locking problems inherent in some thin plate and shell problems. In the formulation of the ELF characteristic arrays, the expression of element strains was adopted in the framework of the element natural coordinates. The non-linear analysis results of idealized edge supports were validated against the reference solutions available in the literature. As a result of the numerical test, the combination of the ELF 9-node shell element and spring element shows an exceptional performance for non-linear analysis of plates and shells under elastic edge supports.     

五零能模式材料设计的非线性有限元计算方法

五零能模式材料是一种新型人工超材料,特征为其弹性模量矩阵的6个特征值中5个为零,可用等效体积模量来描述,表现出类似流体的性质,可被应用于声学隐声斗篷的设计中.然而,根据Norris A N[1]提出的理论,设计五零能模式材料时,与应用变换声学方法设计一般声学人工超材料不同,要求其满足一非线性偏微分方程约束.论文利用非线性有限元的完全拉格朗日方法,推导了这一偏微分方程的弱形式,并给出了相应的非线性有限元计算列式,以及迭代求解的具体算法.最后,给出了五零能摸式材料设计的二维和三维坐标变换数值算例.     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

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.     

Rigid finite element and limit analysis

According to the lower bound theorem of limit analysis the Rigid Finite Element Method (RFEM) is applied to structural limit analysis and the linear programmings for limit analysis are deduced in this paper. Moreover, the Thermo-Parameter Method (TPM) and Parametric Variational principles (PVP) are used to reduce the computational effort while maintaining the accuracy of solutions. A better solution is also obtained in this paper. The project supported by National Natural Science Foundation of China