塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

三维板料成形过程的显式有限元分析

本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问题。板材料服从Ｈｉｌ各向异性弹塑性准则，板料与模具之间的接触界面由主仆面接触搜寻法处理，接触力由罚参数法计算。文中给出了几个三维成形过程的计算实例。数值算例表明，本文方法具有较高的计算精度和计算效率，可在微机上分析中等复杂程度的成形过程     

考虑温度影响的混凝土弹塑性本构模型及其应用

基于弹塑性理论的框架,引入热力学中的自由能概念,建立了考虑温度影响的混凝土弹塑性本构模型。假设混凝土为弹塑性解耦材料,通过采用包含温度的自由能函数建立了弹性应力-应变关系的表达式;基于Drucker-Prager强度准则构造屈服函数,以塑性剪应变作为硬化参量,通过在硬化规律中引入温度变量来反映温度对混凝土塑性性质的影响,得到了塑性应力-应变关系的表达式。基于ABAQUS有限元软件平台,应用半隐式回映算法编制了上述弹塑性本构模型的UMAT子程序,完成了该模型的程序实现。通过对C20及C40混凝土单元在不同温度条件下单轴压缩与双轴压缩的数值模拟,并将模拟结果与已有文献试验结果进行了对比,二者符合较好,验证了该模型的正确性。应用该模型模拟混凝土材料的力学特性,开展了不同温度下集中荷载作用的混凝土简支梁及混凝土轴心受压柱力学响应的数值模拟。计算结果表明,随着温度的升高,相同荷载作用下混凝土构件的力学响应逐渐增大,反映了混凝土材料的力学性能随着温度的升高逐渐退化,这符合不同温度条件下混凝土构件力学响应的基本规律。     

弹塑性损伤本构关系的显式积分算法研究

首先引入弹塑性损伤本构关系,分别从材料软化与残余应变两个方面,描述伪脆性材料的非线性行为.针对结构动力分析中的强非线性问题,给出了弹塑性损伤本构关系的显式积分算法.算法中引入算子分解的思想,将弹塑性本构关系分成塑性与损伤两个模块.首先求解塑性模块,根据有效应空间塑性演化公式,采用前进欧拉算法,直接构造塑性演化的预测值,并且根据屈服函数的漂移构造了误差限公式,作为衡量显式算法精度的指标.将塑性模块求解的结果代入损伤模块,可以方便地求得损伤变量的演化,并最终得到更新后的应力.整个求解过程不需要迭代,可最大程度的算法稳定性.将论文建立的本构关系显式算法与结构分析显式算法结合,构造了结构显式分析方法,并模拟了两个经典算例,算例结果验证了论文方法的有效性.     

屈服准则及切线模量修正的弹塑性计算模型

为提高不同应力状态下金属弹塑性行为的模拟精度,采用含应力三轴度修正的von Mises屈服准则,材料弹塑性本构关系在等效应力$\!$-$\!$-$\!$等效应变曲线基础上,提出一个切线模量为主应力函数的理论修正项. 将这两个修正项结合,以子程序形式编程嵌入ABAQUS主程序,以此模拟几种不同形状试样的弹塑性行为,并将其他屈服准则在单一曲线假设下编程与之对比. 模拟结果与真实试验结果对比发现,对于屈服而言,含应力三轴度修正的vonMises屈服准则比其他屈服理论准确;对于弹塑性阶段计算而言,提出的切线模量为主应力函数这一假设比单一曲线假设更加接近真实试验.      

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

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

平面应变各向异性本构关系及在应力波传播模拟中的应用

为了研究平面应变条件下各向异性材料中应力波传播的特点,利用各向异性弹性Hooke定律、 Tsai-Hill屈服准则、经典塑性流动理论,引入修正的物态方程计及高压下的体积压缩非线性,建立了平面应 变条件下正交各向异性复合材料的弹塑性本构关系,并且分析了二维问题中材料变形引起的主轴旋转及客 观应力率修正问题。最后采用动态显式有限元方法自行编写程序模拟某种纤维增强复合材料碰撞过程中平 面应力波的传播,模拟结果显示,在平面应变条件下应力波在该材料的传播过程中表现出明显的二维效应、各 向异性特点及弹塑性特点。     

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

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

冲压板材拉伸筋阻力的一种有效数值计算方法

将弹塑性有限变形的拟流动角点本构理论和厚向各向异性屈服函数引入弹塑性动力显式有限元列式,对板材通过拉伸筋的变形过程及拉伸筋阻力进行了数值模拟,并与有关实验结果进行了比较很好的一致性,表明了该计算方法的有效性。进而,数值研究了拉伸筋形状、界面摩擦状况以及材料的厚向各向异性对拉伸筋阻力的影响,为实际覆盖件冲压成形中拉伸筋的设置提供了重要的定量依据。     

基于对数应力率大应变本构关系的参数标定研究

在所有率型弹塑性本构模型中,只有对数应力率对应的本构模型能够满足自适应准则.基于对数应力率,采用实心圆轴扭转实验,对大应变弹塑性本构模型中的参数标定问题进行了讨论.推导出了考虑Swift效应时端部自由实心圆轴扭转变形的变形率、对数旋率、Kirchhoff应力及Kirchhoff应力的对数应力率.对于等向强化大应变弹塑性本构关系,给出了由实心圆轴扭转实验标定的、基于Kirchhhoff应力对数应力率的本构关系中塑性刚度函数的表达式.分析了扭转圆轴的Swift效应对塑性刚度函数的影响.结果表明,实心圆轴扭转的轴向伸长变形和径向变形对基于对数应力率大应变本构关系中的塑性刚度函数都有影响.当不考虑Swift效应时,所得塑性刚度函数表达式与不考虑Swift效应时基于Jaumann应力率的塑性刚度函数表达式相同.     

Model of plastic anisotropy evolution with texture-dependent yield surface

Model of evolution of plastic anisotropy due to crystallographic texture development, in metals subjected to large deformation processes, is presented. The model of single grain with the regularized Schmid law proposed by Gambin is used. Evolution of crystallographic texture during drawing, rolling and pure shear is calculated. Phenomenological texture-dependent yield surface for polycrystalline sheets is proposed. Evolution of this yield surface is compared with evolution of phenomenological higher order yield surfaces proposed by Hill and Barlat with Lian for drawing, rolling and pure shear processes. The change of the Hill yield surface and the Barlat–Lian yield surface is obtained by replacing material parameters present in these conditions by texture-dependent functions.     

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.     

Axisymmetric Thermoelastoplastic Stress State of Isotropic Solids of Revolution Under Impulsive Loading

Consideration is given to the solution of a dynamic problem for a solid of revolution with an arbitrary meridional section under impulsive thermomechanical loading inducing elastoplastic strains. The theory of small elastoplastic deformations is used. The constitutive equations are linearized by the variable-parameter method. The unloading process is described by a linear law. The solution technique involves the finite-element approximation in spatial coordinates and the finite-difference representation of time derivatives. Based on the principle of linear summation, recurrent relations are derived for successive evaluation of nodal displacements by an explicit scheme in time. Solution for cylinders and disks are presented to illustrate the influence of elastoplastic deformations on wave processes     

A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials

In this paper, a generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials is derived. The evolution equation for the active yield surface with reference to the memory yield surface is obtained by considering the continuous expansion of the active yield surface during the unloading/reloading process. The incremental constitutive relation based on the associated flow rule is then derived for a general yield function for pressure insensitive and sensitive materials. Detailed incremental constitutive relations for materials based on the Mises yield function, the Hill quadratic anisotropic yield function and the Drucker–Prager yield function are derived as the special cases. The closed-form solutions for one-dimensional stress–plastic strain curves are also derived and plotted for materials under cyclic loading conditions based on the three yield functions. In addition, the closed-form solutions for one-dimensional stress–plastic strain curves for materials based on the isotropic Cazacu–Barlat yield function under cyclic loading conditions are summarized and presented. For materials based on the Mises and the Hill anisotropic yield functions, the stress–plastic strain curves show closed hysteresis loops under uniaxial cyclic loading conditions and the Masing hypothesis is applicable. For materials based on the Drucker–Prager and Cazacu–Barlat yield functions, the stress–plastic strain curves do not close and show the ratcheting effect under uniaxial cyclic loading conditions. The ratcheting effect is due to different strain ranges for a given stress range for the unloading and reloading processes. With these closed-form solutions, the important effects of the yield surface geometry on the cyclic plastic behavior due to the pressure-sensitive yielding or the unsymmetric behavior in tension and compression can be shown unambiguously. The closed form solutions for the Drucker–Prager and Cazacu–Barlat yield functions with the associated flow rule also suggest that a more general anisotropic hardening theory needs to be developed to address the ratcheting effects for a given stress range.     

内时理论在板料成形失稳中的应用

应用Valanis提出的内时本构方程，研究了板料成形的拉伸失稳问题，推导出单向和双向拉伸应力状态下的内时本构方程，据此分析了分散性失稳和集中性失稳。该文推导出应用于拉伸失稳分析时内时理论的近似表达式，它对应于经典塑性理论解，同时给出了内时理论的完整迭代数值解。结果表明内时理论具有很好的适用性。     

An admissible form of an inelastic constitutive equation—II. Examples of elastoplastic constitutive equation

Some examples of elastoplastic constitutive equation are presented using the general theory reported in the preceding paper (Part I). Some examinations of them are given to show that the theory is self-consistent and useful especially for anisotropic materials or materials with anisotropy resulting from plastic deformation. Mises' and Yoshimura's yield functions and a kind of quadratic function are adopted as the yield function. Formulae of r-value after arbitrary pre-straining are given which are of paramount importance in the field of press-forming of sheet metals. Several examples of stress-strain curves for various loading paths are also given.     

SECTIONAL FINITE ELEMENT ANALYSIS OF COUPLED DEFORMATION BETWEEN ELASTOPLASTIC SHEET METAL AND VISCO-ELASTOPLASTIC BODY

The present paper is devoted to developing a new numerical simulation method for the analysis of viscous pressure forming (VPF), which is a sheet flexible-die forming (FDF) process. The pressure-carrying medium used in VPF is one kind of semisolid, flowable and viscous material and its deformation behavior can be described by the visco-elastoplastic constitutive model. A sectional finite element model for the coupled deformation analysis between the visco- elastoplastic pressure-carrying medium and the elastoplastic sheet metal is proposed. The resolution of the Updated Lagrangian (UL) formulation is based on a static explicit approach. The frictional contact between sheet metal and visco-elastoplastic pressure-carrying medium is treated by the penalty function method. Coupled deformation between sheet metal and visco-elastoplastic pressure-carrying medium with large slip is analyzed to validate the developed algorithm. Finally, the viscous pressure bulging (VPB) process of DC06 sheet metal is simulated. Good agreement between numerical simulation results and experimental measurements shows the validity of the developed algorithm.     

Forming of aluminum alloys at elevated temperatures – Part 2: Numerical modeling and experimental verification

The temperature-dependent Barlat YLD96 anisotropic yield function developed previously [Forming of aluminum alloys at elevated temperatures – Part 1: Material characterization. Int. J. Plasticity, 2005a] was applied to the forming simulation of AA3003-H111 aluminum alloy sheets. The cutting-plane algorithm for the integration of a general class of elasto-plastic constitutive models was used to implement this yield function into the commercial FEM code LS-Dyna as a user material subroutine (UMAT). The temperature-dependent material model was used to simulate the coupled thermo-mechanical finite element analysis of the stamping of an aluminum sheet using a hemispherical punch under the pure stretch boundary condition. In order to evaluate the accuracy of the UMAT’s ability to predict both forming behavior and failure locations, simulation results were compared with experimental data performed at several elevated temperatures. Forming limit diagrams (FLDs) were developed for the AA3003-H111 at several elevated temperatures using the M-K model in order to predict the location of the failure in the numerical simulations. The favorable comparison found between the numerical and experimental data shows that a promising future exists for the development of more accurate temperature-dependent yield functions to apply to thermo-hydroforming process.