牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

一致切线刚度法在三维弹塑性有限元分析中的应用

本文提出了一致切线刚度法，并把它应用于三维弹塑性有限元分析问题。从而解决了增量迭代弹塑性有限元分析方法中长期存在的速度慢、精度低问题，一致切线刚度法满足加卸载互补准则，即没有应力漂移现象，具有一阶精度、二阶迭代收敛速度、计算量少和无条件稳定等优点，借助算例对一致切线刚度法和传统切线刚度法（包括路径相关和路径无关两种结构变量更新格式）从计算精度、迭代收敛速度和计算量等几方面进行了比较。     

屈服准则在有限元软件中实现的正确性验证

基于向后欧拉完全隐式积分方法,采用回映算法将Yld2000-2d屈服准则嵌入ABAQUS商业有限元软件.基于严格的力学推导,提出了能够有效验证屈服准则嵌入工作的正确性的方法.针对57540铝合金板,采用上述方法验证了Yld2000-2d屈服准则在ABAQUS有限元软件中的嵌入工作的正确性.     

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

一定初缺陷杆在轴向冲击下弹塑性动态屈曲有限元计算

本文用有限元方法分析了一定初缺陷杆受轴向冲击的弹塑性动态屈曲。由变形功相关的屈曲判据求出屈曲时间,计算了初缺陷及冲击载荷形状和大小对屈曲时间的影响。     

应力函数及其对偶理论在有限元中的应用

借助于Cosserat连续介质模型，探讨了应力函数和位移对避免有限元C$^{1}$ 连续性困难的互补性作用. 通过对应力函数对偶理论的深入分析，为将应力函数列式得到的 余能单元转化为具有一般位移自由度的势能单元提供了严格的理论基础，在此基础上， 给出应用应力函数构造有限元的一般方法.     

Hill（1979）屈服准则的展开及在金属薄板成形极限图中的应用

本文对于Hill(1979)屈服准则在平面应力条件下进行了展开,给出了五种特定条件下的等效应力和等效应变速率表达式,对这些特例的凸条件检验表明,其中的四种可做为合理的屈服函数,将这些屈服函数应用于金属薄板在拉伸范围内的成形极限图的M—K理论预测,获得了较为理想的结果。     

牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

框支剪力墙结构在地震力作用下弹塑性分析的有限元方法

建立了框支剪力墙结构弹塑性分析的有限元模型,采用有旋转自由度的高精度平面单元分析钢筋混凝土墙体在地震力作用下混凝土开裂、屈服和钢筋屈服的破损过程,为钢筋混凝土框支剪力墙的抗震设计提供较精确的分析方法。     

广义双剪粘弹塑性本构模型在渤海海冰动力学中的应用

在连续介质力学基础上建立了一个广义双剪粘弹塑性海冰动力学本构模型。该模型在海冰屈服前采用Kelvin-Vogit粘弹性模型,考虑中间主应力和静水压力对海冰屈服的影响选用广义双剪应力屈服准则作为海冰屈服判据,屈服后采用相关联的正则流动法则。采用该本构模型对渤海海冰动力过程进行了48小时数值模拟,讨论了辽东湾海冰的厚度、密集度、冰速和主应力的分布规律,其中海冰厚度分布与卫星遥感资料符合良好,从而有效地验证了该广义双剪粘弹塑性本构模型在海冰动力学中的可靠性。     

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.     

Finite element analysis of elastic-plastic plane stress problems based upon Tresca yield criterion

Summary A finite element technique for the elastic-plastic analysis of two dimensional structures subjected to conditions of plane stress and monotonically increasing loads is presented. The complete load deformation history as well as the propagation of the yield zones in the structure up to plastic collapse are studied. The material is assumed to be elastic-perfectly plastic and yielding is governed by Tresca yield condition. Plastic stress-strain relations for the sides and corners of Tresca yield condition are derived in terms of the components of the stresses and strains along a fixed reference coordinate system and the direction of the principal stress. The load is applied in small increments and the principal stress direction for each plastic element during a load increment is determined by an interpolation technique which leads to stresses that satisfy the yield condition. Numerical examples are given to illustrate the accuracy of the results obtained by the proposed method.

---

Übersicht Es wird eine finite element-Methode vorgeschlagen zur Berechnung zweidimensionaler elastischplastischer Baukonstruktionen, die einem ebenen Spannungszustand mit monoton wachsender Last unterworfen sind. Das Belastungs-Verformungs-Verhalten und die Ausbreitung der Fließbereiche bis zum plastischen Zusammenbruch werden untersucht. Das Material soll der Fließbedingung von Tresca genügen. Spannungs-Verformungs-Berechnungen für die Ränder und Ecken des Tresca-Bereiches werden in Komponenten der Spannungen und Dehnungen längs der Achsriclitungen eines festen Bezugssystems und der Hauptspannungsrichtungen ausgedrückt. Die Belastung wird in kleinen Stufen aufgebracht, und die Hauptspannungsrichtungen werden für jedes plastische Teilchen während des Lastaufbringens durch Interpolation bestimmt. Die so erhaltenen Spannungen genügen der Fließbedingung. Durch numerische Beispiele wird die Genauigkeit des vorgeschlagenen Verfahrens demonstriert.

---

---

On leave from Northwestern University, Evanston, Illinois     

Dynamic finite element with diagonalized consistent mass matrix and elastic-plastic impact calculation

There are some common difficulties encountered in elastic-plastic impact codes such as EPIC^[1,2], NONSAP^[3] etc. Most of these codes use the simple linear functions usually taken from static problems to represent the displacement components. In such finite element formulation, the strain and stress components are constants in every element. In the equations of motion, the stress components in general appear in the form of their space derivatives. Thus, if we use such form functions to represent the displacement components, the effect of internal stresses to the equations of motion vanishes identically. The usual practice to overcome such difficulties is to establish as self-equilibrium system of internal forces acting on various nodal points by means of transforming equations of motion into variational form of energy relation through the application of virtual displacement principle. The nodal acceleration is then calculated from the total force acting on this node from all the neighbouring elements. The transformation of virtual displacement principle into the variational energy form is performed on the bases of continuity conditions of stress and displacement throughout the integrated space. That is to say, on the interface boundary of finite element, the assumed displacement and stress functions should be conformed. However, it is easily seen that, for linear form function of finite element calculation, the displacement continues everywhere, but not the stress components. Thus, the convergence of such kind of finite element computation is open to question. This kind of treatment has never been justified even in approximation sense. Furthermore, the calculation of nodal points needs a rule to calculate the mass matrix. There are two ways to establish mass matrix, namely lumped mass method and consistent mass method [4]. The consistent mass matrix can be obtained naturally through finite element formulation, which is consistent to the assumed form functions. However, the resulting consistent mass matrix is not in diagonalized form, which is inconvenient for numerical computation. For most codes, the lumped mass matrix is used, and in this case, the element mass is distributed in certain assumed proportions to all the nodal points of this element. The lumped mass matrix is diagonalized with diagonal terms composed of the nodal mass. However, the lumped mass assumption has never been justified. All these difficulties are originated from the simple linear form functions usually used in static problems.In this paper, we introduce a new quadratic form function for elastic-plastic impact problems. This quadratic form function possesses diagonalized consistent mass matrix, and non-vanishing effect of internal stress to the equations of motion. Thus with this kind of dynamic finite element, all above-said difficulties can be eliminated.     

Diagonalized consistent mass matrix and the dynamical finite element analysis of elastic-plastic impact in axisymmetrical problems

In this paper, the diagonalized consistent mass matrix is found for the triangular ring element in axisymmetrical problems. The results of this work eliminate the feeling of uncertainty and arbitrariness of lumped mass method on the one hand and the difficulty of computation due to non-diagonalized character of consistent mass method on the other. This paper gives also the foundations of the finite element analysis of elastic-plastic axisymmtrical impact problems.     

On the thermodynamics of elastic-plastic materials with temperature-dependent moduli and yield stresses

Summary The subject of this article is the thermodynamics of perfect elastic-plastic materials undergoing unidimensional, but not necessarily isothermal, deformations. The first and second laws of thermodynamics are employed in a form in which only the following quantities appear: the temperature , the elastic strain _e, the plastic strain _p, the elastic modulus (gq), the yield strain (gq), the heat capacity (_e, _p,), the latent elastic heat _e(_e, _p, ), and the latent plastic heat _p(_e, _p, ). Relations among the response functions , , , _e, and _p are derived, and it is shown that a set of these relations gives a necessary and sufficient condition for compliance with the laws of thermodynamics. Some observations are made about the existence and uniqueness of energy and entropy as functions of state.Dedicated to Clifford Truesdell on the occasion of his 60th birthdayThis research was supported by the U.S. National Science Foundation.     

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontinuous deformation analyses (3D DDA) has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block deformation. In this paper, 3D DDA is coupled with tetrahedron finite elements to tackle these two problems. Tetrahedron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topology shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Validation is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demonstrates the robustness and versatility of the coupled method.