Plastic limit analysis of incompatible finite element method

This paper describes an incompatible finite element model satisfying the consistencycondition of energy to solve the numerical precision problem of finite element solution inperfectly plastic analysis.In this paper the reason and criterion of the application of themodel to plastic limit analysis are discussed,and an algorithm of computing plastic limitload is given.     

Rigid finite element and its applications in engineering

Rigid Finite Element Method (RFEM) was proposed to simulate the mechanical behavior of discontinuous structures such as rock and soil structures. The authors' work on the theory and applications of RFEM is summarized in this paper. Based on the theory of RFEM, the Elastic Body-Seams Model (EBSM) is proposed to take the deformation and damage of rock masses into account.     

PARAMETRIC VARIATIONAL PRINCIPLE BASED ELASTIC-PLASTIC ANALYSIS OF HETEROGENEOUS MATERIALS WITH VORONOI FINITE ELEMENT METHOD

The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.     

Variational principles of perforated thin plates and finite element method for buckling and post-buckling analysis

On the basis of the general theory of perforated thin plates under large deflections^{[1, 2]}, variational principles with deflectionw and stress functionF as variables are stated in detail. Based on these principles, finite element method is established for analysing the buckling and post-buckling of perforated thin plates. It is found that the property of element is very complicated, owing to the multiple connexity of the region. Project supported by National Natural Science Foundation of China.     

Finite element analysis of the duct flow of Bingham plastic fluids: an application of the variational inequality

The duct flow of Bingham plastic fluids is analysed with the variational inequality-based finite element method. The problem of tracking the yield surface is solvable through the regularization technique which can be easily incorporated into the existing finite element code. The existence theorem of this method was established through the theory of variational inequalities. A small positive constant is added to the second shear rate invariant, resulting in an apparent viscosity of finite magnitude in the unyielding plug zone. This makes the minimization of the non-differential variational integral possible. In order to achieve convergence at small regularization parameter, a zero-order continuation is employed. It is also shown that a fine tessellation of the flow domain is necessary for tracking the yield surfaces unambiguously. Two classes of duct flow, namely axial flows in eccentric annuli and in an L-shaped duct, were investigated. In both cases it was easy to show the presence of the mobile plugs around the duct centres from the axial velocity profiles; however, the stagnant plugs at the narrow side in eccentric annuli with large eccentricity and near the apex of right-angled corners in an L-shaped duct could only be identified from the calculation of the distributions of the second shear rate or shear stress invariant. © 1997 John Wiley & Sons, Ltd.     

Micro/macromechanical plastic limit analyses of composite materials and structures

The load-bearing capacities of ductile composite materials and structures are studied by means of a combined micro/macromechanics approach. Firstly, on the microscopic scale, the aim is to get the macroscopic strength domains by means of the homogenization theory of micromechanics. A representative volume element (RVE) is selected to reflect the microstructures of the composite materials. By introducing the homogenization theory into the kinematic limit theorem of plastic limit analysis, an optimization format to directly calculate the limit loads of the RVE is obtained. And the macroscopic yield criterion can be determined according to the relation between macroscopic and microscopic fields. Secondly, on the macroscopic scale, by introducing the Hill's yield criterion into the kinematic limit theorem, the limit loads of orthotropic structures such as unidirectional fiber-reinforced composite structures are worked out. The finite element modeling of the kinematic limit analysis is deduced into a nonlinear mathematical programming with equality-constraint conditions that can be solved by means of a direct iterative algorithm. Finally, some examples are illustrated to show the application of the present approach. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctoral Dissertation of China (No. 200025), the Fund of the Ministry of Education of China for Returned Oversea Scholars and the Basic Research Foundation of Tsinghua University.     

Plotting isoclinics for hybrid photoelasticity and finite element analysis

Displacement-based finite element method formulations are coupled with stress-based photoelasticity analysis. As the stress field is discontinuous at the interelement boundaries, the introduced smoothing procedure enables the generation of high-quality digital images acceptable for hybird experimental-numerical techniques. The proposed methods are applicable for the analysis of static and dynamic results of experimental photoelasticity.     

The random variational principle and finite element method

In this paper, we introduced the random materials, geometrical shapes, force and displacement boundary condition directly into the functional variational formulations and developed a unified random variational principle and finite element method with the small parameter perturbation method. Numerical examples showed that the methods have the advantages of the simple and convenient program implementation, and are effective for the random mechanics problems.     

A modified finite element method for determining equilibrium capillary surfaces of liquids with specified volumes

This paper describes a modified finite element method (MFEM) for determining the static equilibrium shape of the capillary surface of a liquid with a prescribed volume constrained by rigid boundaries with arbitrary shapes. It is assumed that the liquid is in static equilibrium under the influence of surface tension, adhesion, and gravity forces. This problem can be solved by employing the conventional FEM; however, a major difficulty arises due to the presence of the volume (integral) constraint and usually requires the use of the Lagrange multiplier method, the sequential unconstrained minimization technique, or the augmented Lagrange multiplier method. With the MFEM, the space variables defining the equilibrium surfaces (or curves) are expanded in terms of parametric interpolation functions, which are designed such that the boundary conditions and the integral constraint equation are automatically satisfied during each iteration of a direct numerical search process. Hence, there is no need to include Lagrange multipliers and/or penalty factors and the problem can be treated more simply as one involving unconstrained optimization. This investigation indicates that the MFEM is more efficient and reliable than the other methods. Results are presented for several case study problems involving liquid solder drops. Copyright © 2000 John Wiley & Sons, Ltd.     

齿轮精密锻造三维有限元分析

开发了齿轮精密锻造三维有限元分析软件系统 ,采用三维弹塑性大变形有限元分析方法 ,分析了精密锻造过程的力学响应 ,获得上模、精锻齿轮和下模应力与变形的详尽信息。