共查询到20条相似文献，搜索用时 0 毫秒
1.
涉及复杂材料弹塑性损伤问题数值计算研究时,不仅需要选择恰当预测损伤和破坏的本构模型,还需要有效和稳健的本构积分算法。首先,阐述了在热力学和连续介质力学框架下建立弹塑性损伤本构模型的基本步骤;其次,基于Lemaitre等向硬化弹塑性损伤耦合本构模型、相应的本构积分算法-完全隐式返回映射算法(Fully Return Mapping Algorithm)和一致切线模量,采用C++语言在Visual 6.0 环境下编制有限元本构求解程序,在塑性损伤修正步中求解返回映射方程时,选取一种简单的形式,只需迭代求解一个标量非线性方程,计算效率较高。最后,通过缺口圆棒数值算例初步验证了程序的正确性,并编制接口程序对计算结果进行可视化。研究结果表明积分算法的有效性及程序的正确性,Lemaitre等向硬化弹塑性损伤耦合本构模型能够较好地模拟韧性材料的破坏发展过程,可以求解类似的有限元边界值问题,为考虑损伤特性的韧性材料结构研究和设计奠定基础。 相似文献
2.
Feng Lu Zhang Guang Zhang Ke-shi J. Q. Xu Yu Hai-dong 《International journal for numerical methods in engineering》2004,61(3):406-432
Based on the concept of continuum damage mechanics, an anisotropic damage model for single crystals under the theory of crystal plasticity is presented. Damage and inelastic deformations are incorporated in the proposed model which is developed within the framework of thermodynamics with internal state variables. The dependence of the plastic anisotropy on the damage evolution has been considered. The anisotropic damage is characterized kinematically here through a second-order damage tensor which is physically based. The proposed model can successfully describe the interaction between the evolution of micro-structure of single crystals such as lattice orientation and the hardness development of each slip system and the process of material degradation. The Newton–Raphson iterative scheme is used to integrate the constitutive equations that work directly with the evolution equations for the elastic deformation gradient. The consistent algorithmic tangent stiffness for the present algorithm is formulated. The prescribed algorithm together with the consistent algorithmic tangent stiffness has been implemented into the ABAQUS finite element code by using user subroutine. Using the loading processes with homogeneous deformations and simulation of the classical tensile test of a notched bar illustrate the basic aspects of the model described. Numerical simulations show the validation and performance of the present model and algorithm. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
3.
G. Asensio C. Moreno 《International journal for numerical methods in engineering》2003,57(7):991-1014
The return mapping algorithm is one of the most efficient procedures to solve elasto‐plastic problems. However, a criticism that may be lodged against this method is the difficulty of the practical computation of the consistent tangent matrix when the return is non‐radial. Much research has been done to handle this matrix. In this paper, a unified approach is presented in such a way that a simple closed‐form expression gives the consistent tangent matrix for the classical constitutive relations (von Mises, Tresca, Mohr–Coulomb, Drucker–Prager). The basic ideas are in the properties of eikonal equations appearing in several fields as image treatment, short time computation in elastic waves and others. The same kinds of ideas can be extended to non‐classical models. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
4.
Ming‐Xiang Chen Qi Peng Junjie Huang 《International journal for numerical methods in engineering》2014,99(9):654-681
A simple and compact representation framework and the corresponding efficient numerical integration algorithm are developed for constitutive equations of isotropic elastoplasticity. Central to this work is the utilization of a set of mutually orthogonal unit tensor bases and the corresponding invariants. The set of bases can be regarded equivalently as a local cylindrical coordinate system in the three‐dimensional coaxial tensor subspace, namely, the principal space. The base tensors are given in the global coordinate system. Similar to the principal space approach, the proposed method reduces the problem dimension from six to three. In contrast to the conventional approach, the transformation procedure between the principal space and the general space is avoided and explicit computation of the principal axes is bypassed. With the proposed technique, the matrices, which need to be inverted during iteration, take a simple form for the great majority of constitutive equations in use. The tangent operator consistent with the proposed algorithm can be decomposed into the direct sum of two linear maps over the coaxial tensor subspace and the subspace orthogonal to it. Consequently, its closed form is derived in an extremely simple manner. Finally, numerical examples demonstrate the high quality performances of the proposed scheme. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
5.
6.
Artem S. Semenov Albrecht C. Liskowsky Herbert Balke 《International journal for numerical methods in engineering》2010,81(10):1298-1340
Return mapping algorithms for a rather general class of phenomenological rate‐independent models for ferroelectroelastic materials are presented. The fully coupled thermodynamically consistent three‐dimensional constitutive model with two internal variables (remanent polarization vector and remanent strain tensor) proposed by C. M. Landis in 2002 is used for the simulation of electromechanical hysteresis effects in polycrystalline ferroelectric ceramics. Based on the operator splitting methodology, the return mapping algorithm employs the closest point projection scheme to obtain an efficient and robust integration of the constitutive model. The consistent tangent operator is obtained in closed form by linearizing the return mapping algorithm, and is found to be non‐symmetric in the general case due to the dependence of the switching criterion on internal variables. Conditions that provide the symmetry of the consistent tangent matrix are analyzed. The compactness and generality of the received relations are achieved by means of using the thermodynamically based compact notations combining mechanical and electrical values. Both the cases scalar potential finite element (FE) formulation (primary variables: strain and electric field) and vector potential FE formulation (primary variables: strain and electric displacement) are considered. The accuracy and robustness of the algorithms are assessed through numerical examples. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
7.
8.
Karel Matou&#x; Antoinette M. Maniatty 《International journal for numerical methods in engineering》2004,60(14):2313-2333
Anisotropic, elasto‐viscoplastic behaviour in polycrystalline materials is modelled using a new, updated Lagrangian formulation based on a three‐field form of the Hu‐Washizu variational principle to create a stable finite element method in the context of nearly incompressible behaviour. The meso‐scale is characterized by a representative volume element, which contains grains governed by single crystal behaviour. A new, fully implicit, two‐level, backward Euler integration scheme together with an efficient finite element formulation, including consistent linearization, is presented. The proposed finite element model is capable of predicting non‐homogeneous meso‐fields, which, for example, may impact subsequent recrystallization. Finally, simple deformations involving an aluminium alloy are considered in order to demonstrate the algorithm. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
9.
J. Clausen L. Damkilde L. Andersen 《International journal for numerical methods in engineering》2006,66(6):1036-1059
A new return method for implicit integration of linear isotropic yield criteria is presented. The basic idea is to perform all the manipulations in the principal stress space and thereby achieve very simple formulae for calculating the plastic corrector stresses, based on the constant gradient of such criteria. The return formulae are in closed form and no iteration is required. The method accounts for three types of stress return: return to a single yield plane, to a discontinuity line at the intersection of two yield planes and to a discontinuity point at the intersection between three or more yield planes. The infinitesimal and the consistent elastoplastic constitutive matrix are calculated for each type of stress return, as are the conditions to ascertain which type of return is required. The method is exemplified with the Mohr–Coulomb yield criterion. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
10.
J. YANG G. Z. WANG F. Z. XUAN S. T. TU 《Fatigue & Fracture of Engineering Materials & Structures》2013,36(6):504-514
Constraint can be divided into two conditions of in‐plane and out‐of‐plane, and each of them has its own parameter to characterize. However, in most cases, there exists a compound change of both in‐plane and out‐of‐plane constraint in structures, a unified measure that can reflect both of them is needed. In this paper, the finite element method (FEM) was used to calculate the equivalent plastic strain (ɛp) distribution ahead of crack tips for specimens with different in‐plane and out‐of‐plane constraints, and the FEM simulations based on Gurson–Tvergaard–Needleman (GTN) damage model and a small number of tests were used to obtain fracture toughness for the specimens with different constraints. Unified measure and characterisation parameter of in‐plane and out‐of‐plane constraints based on crack‐tip equivalent plastic strain has been investigated. The results show that the area APEEQ surrounded by the ɛp isoline ahead of crack tips can characterize both in‐plane and out‐of‐plane constraints. Based on the area APEEQ, a unified constraint characterisation parameter Ap was defined. It was found that there exists a sole linear relation between the normalised fracture toughness JIC/Jref and regardless of the in‐plane constraint, out‐of‐plane constraint and the selection of the ɛp isolines. The unified JIC/Jref−reference line can be used to determine constraint‐dependent fracture toughness of materials. The FEM simulations with the GTN damage model (local approach) can be used in obtaining the unified JIC/Jref−reference line for materials with ductile fracture. 相似文献
11.
为了克服一般弹塑性损伤模型不能反映岩体结构、岩块强度、应力状态的影响以及非线性破坏特征等问题,该文基于广义的Hoek-Brown（HB）屈服准则,考虑损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立了岩体弹塑性损伤本构模型,给出了损伤变量定义及演化方程。针对该模型在数值求解过程中存在的奇异点问题,从主应力空间推导了弹塑性损伤模型的完全隐式返回映射求解算法,包括弹性预测、塑性修正和损伤修正三个步骤。通过ABAQUS软件的用户子程序接口Umat,实现了弹塑性损伤模型的数值求解过程。采用单轴、三轴压缩试验和隧道算例对模型算法进行验证和分析,结果表明,所建立的HB损伤本构模型能够很好地描述岩体材料的力学特性,在实际岩体工程的损伤模拟中效果令人满意,计算结果对工程有指导意义。 相似文献
12.
P. Ellsiepen S. Hartmann 《International journal for numerical methods in engineering》2001,51(6):679-707
For the numerical solution of materially non‐linear problems like in computational plasticity or viscoplasticity the finite element discretization in space is usually coupled with point‐wise defined evolution equations characterizing the material behaviour. The interpretation of such systems as differential–algebraic equations (DAE) allows modern‐day integration algorithms from Numerical Mathematics to be efficiently applied. Especially, the application of diagonally implicit Runge–Kutta methods (DIRK) together with a Multilevel‐Newton method preserves the algorithmic structure of current finite element implementations which are based on the principle of virtual displacements and on backward Euler schemes for the local time integration. Moreover, the notion of the consistent tangent operator becomes more obvious in this context. The quadratical order of convergence of the Multilevel‐Newton algorithm is usually validated by numerical studies. However, an analytical proof of this second order convergence has already been given by authors in the field of non‐linear electrical networks. We show that this proof can be applied in the current context based on the DAE interpretation mentioned above. We finally compare the proposed procedure to several well‐known stress algorithms and show that the inclusion of a step‐size control based on local error estimations merely requires a small extra time‐investment. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
13.
M. Kobayashi M. Mukai H. Takahashi N. Ohno T. Kawakami T. Ishikawa 《International journal for numerical methods in engineering》2003,58(10):1523-1543
This paper describes the implicit integration and consistent tangent modulus of an inelastic constitutive model with transient and steady strain rates, both of which are time‐ and temperature‐dependent; the transient rate is influenced by the evolution of back stress decomposed into parts, while the steady rate depends only on applied stress and temperature. Such a non‐unified model is useful for high‐temperature structural analysis and is practical owing to the ease in determining material constants. The implicit integration is shown to result in two scalar‐valued coupled equations, and the consistent tangent modulus is derived in a quite versatile form by introducing a set of fourth‐rank constitutive parameters into the discretized evolution rule of back stress. The constitutive model is, then, implemented in a finite element program and applied to a lead‐free solder joint analysis. It is demonstrated that the implicit integration is very accurate if the multilinear kinematic hardening model of Ohno and Wang is employed, and that the consistent tangent modulus certainly affords quadratic convergence to the Newton–Raphson iteration in solving nodal force equilibrium equations. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
14.
L. S. Miers J. C. F. Telles 《International journal for numerical methods in engineering》2008,76(7):1090-1107
A meshless procedure, based on boundary integral equations, is proposed to analyze elastoplastic problems. To cope with non‐linear problems, the usual boundary element method introduces domain discretization cells, often considered a ‘drawback’ of the method. Here, to get rid of the standard element and cell, i.e. boundary and domain discretization, the orthogonal moving least squares (also known as improved moving least squares) method is used. The algorithm adopted to solve these particular inelastic non‐linear problems is a well‐established, criterion‐independent implicit procedure, previously developed by the authors. Comparative results are presented at the end to illustrate the effectiveness of the proposed techniques. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
15.
Takuya Anjiki Koichi Hashiguchi 《International journal for numerical methods in engineering》2022,123(1):291-303
The extended overstress model is formulated based on the subloading surface model with the smooth elastic-inelastic transition, which is called the subloading-overstress model. Therein, the rigorous translation rules of the elastic-core and the similarity-center, the limitation in the expansion of the subloading surface, and so forth are incorporated. The model possesses the basic structure capable of describing the monotonic/cyclic loading behaviors at the general rate of deformations from the quasi-static to the impact loading. The experiments were conducted using the spheroidal graphite cast iron under the various loading conditions. It was verified that the experimental results can be simulated accurately by the subloading-overstress model. Further, the complete implicit stress integration algorithm based on the return-mapping projection is formulated for the present subloading-overstress model and implemented into Abaqus through UMAT. Then, the deformation analyses of the R-notched cylinder were performed by the present algorithm. Consequently, the performability of the present algorithm is verified by the analyses of the boundary-value problem under the cyclic loadings. 相似文献
16.
Xiaojia Wang Qingchun Meng Weiping Hu 《Fatigue & Fracture of Engineering Materials & Structures》2019,42(6):1373-1386
Fatigue damage of butt‐welded joints is investigated by a damage mechanics method. First, the weld‐induced residual stresses are determined by using a sequentially coupled thermo‐mechanical finite element analysis. The plastic damage of material is then calculated with the use of Lemaitre's plastic damage model. Second, during the subsequent fatigue damage analysis, the residual stresses are superimposed on the fatigue loading, and the weld‐induced plastic damage is considered as the initial damage via an elasto‐plastic fatigue damage model. Finally, the fatigue damage evolution, the relaxation of residual stress, and the fatigue lives of the joints are evaluated using a numerical implementation. The predicted results agree well with the experimental data. 相似文献
17.
R. Fincato S. Tsutsumi 《International journal for numerical methods in engineering》2018,113(11):1729-1754
18.
M. Kadkhodayan L. C. Zhang 《International journal for numerical methods in engineering》1995,38(14):2413-2431
This paper proposed a new efficient method, the consistent DXDR method, for analysing general elastic-plastic problems. Two important factors in computational plasticity, the convergence and stability, were addressed. The method was constructed through a natural combination of two vital components: the DXDR algorithm which had been proved to be a powerful equation solver developed by the authors and a stable consistent algorithm for the integration of the constitutive equations of plasticity. Numerical examples demonstrated in detail the efficiency, accuracy, and stability attainable in solving various engineering problems. 相似文献
19.
Zhi‐Qiao Wang Guan‐Suo Dui 《International journal for numerical methods in engineering》2008,75(12):1435-1456
This paper presents alternative forms of hyperelastic–plastic constitutive equations and their integration algorithms for isotropic‐hardening materials at large strain, which are established in two‐point tensor field, namely between the first Piola–Kirchhoff stress tensor and deformation gradient. The eigenvalue problems for symmetric and non‐symmetric tensors are applied to kinematics of multiplicative plasticity, which imply the transformation relationships of eigenvectors in current, intermediate and initial configurations. Based on the principle of plastic maximum dissipation, the two‐point hyperelastic stress–strain relationships and the evolution equations are achieved, in which it is considered that the plastic spin vanishes for isotropic plasticity. On the computational side, the exponential algorithm is used to integrate the plastic evolution equation. The return‐mapping procedure in principal axes, with respect to logarithmic elastic strain, possesses the same structure as infinitesimal deformation theory. Then, the theory of derivatives of non‐symmetric tensor functions is applied to derive the two‐point closed‐form consistent tangent modulus, which is useful for Newton's iterative solution of boundary value problem. Finally, the numerical simulation illustrates the application of the proposed formulations. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
20.
T.‐S. Cao P. Montmitonnet P.‐O. Bouchard 《International journal for numerical methods in engineering》2013,96(9):561-583
In an effort to implement Gurson‐type models into a mixed velocity–pressure finite element formulation with the MINI‐element P1 + ∕ P1, the algorithm proposed by Aravas (IJNME, 1987) to integrate the pressure dependent plasticity as well as the formulations of consistent tangent moduli have been analyzed. This work firstly reviews and clarifies the mathematical basis of the formulations used by Aravas (IJNME, 1987) and demonstrates the equality of the tangent moduli formulations proposed by Govindarajan and Aravas (CNME, 1995) and Zhang (CMAME, 1995), which are widely used in the literature. A unified formulation to calculate the tangent moduli is proven, and its accuracy is also investigated by the finite difference method. The implementation of the Gurson–Tvergaard–Needleman model is then detailed for the mixed velocity–pressure finite element formulation, which employs the MINI‐element P1 + ∕ P1. Due to the particularity of this element, one needs to calculate two tangent moduli instead of one. The formulas for calculating the ‘linear tangent modulus’ and the ‘bubble tangent modulus’ are then detailed. Finally, comparison tests are carried out with ABAQUS (Dassault System, Simulia Corp., Providence, RI, USA) in order to validate the present implementation for both homogeneous and heterogeneous deformations. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献