考虑细观变形特征的镁合金宏观本构模型及数值模拟

为了能有效描述镁合金宏观各向异性塑性行为,考虑了滑移、孪生、去孪生三种细观变形模式的特点,给出了相应的硬化函数;根据VonMises屈服准则,发展了一种镁合金宏观本构模型及其迭代算法。模型将变形模式的开启与晶粒取向相关联,同时针对镁合金孪生变形时引起的晶粒重新定向问题,描述了一种晶向偏转的方法。在此基础上编写了ABAQUS/UMAT材料用户子程序;利用开发的本构模型,开展了单轴拉伸、单轴压缩、单轴循环拉压加载条件下镁合金塑性行为的数值模拟,并对随机织构下的镁合金板材轧制过程进行了有限元仿真实验。模拟结果表明:单轴拉伸、单轴压缩和循环加载情形下的镁合金宏观硬化行为与实验结果基本吻合;轧制后镁合金板材表现出了应力-应变不均匀特性,多晶织构演化结果与实验结果基本一致。说明文中所提出的宏观本构模型、晶向偏转模型能够有效描述镁合金的宏观塑性行为和织构演化。     

用于板料成形反向模拟的特殊应力更新算法

提出了一种改进的反向模拟法,以最终构型为研究对象,采用Euler坐标系,基于虚功原理获得有限元列式. 改进的反向模拟法采用了一种基于塑性流动理论的本构方程,可以充分考虑应变历史对塑性变形的影响. 为了避免流动理论应力更新算法过程中关于未知量\Delta\lambda 的非线性方程的求解,引入等效应力思想,无需Newton-Raphson迭代直接计算未知量\Delta \lambda . 盒形件的拉深实例中,传统的基于塑性形变本构方程的反向模拟法和改进的基于塑性流动本构方程的反向模拟法计算结果,分别与基于增量有限元法的正向数值模拟求解器LS-DYNA计算结果进行对比. 通过获得的坯料轮廓、成形极限图、等效应变分布、计算效率等的比较,验证了所提出的基于塑性流动理论本构模型的应力更新算法的有效性.      

FeCrAl合金材料织构对其宏观力学本构关系的影响研究

FeCrAl合金具有优良的高温抗氧化性和耐辐照性能,是事故容错核燃料包壳的重要候选材料. 其在加工过程和热处理过程中易形成α纤维织构（<110>//RD）和γ纤维织构（<111>//ND）,会影响材料的宏观力学性能与深加工成形能力. 本研究针对具有不同织构的多晶FeCrAl合金,建立了代表性体元模型, 使用晶体塑性有限元方法,在ABAQUS/Explicit中模拟材料单轴加载下的宏观应力应变曲线,分析不同织构对FeCrAl合金宏观力学本构关系的影响. 计算结果表明,对于具有α织构、γ织构和晶粒无择优取向的材料,在轧向上的应力应变曲线差异较小. γ织构会引起材料强烈的各向异性,在轧面法向上的屈服强度远高于轧向和横向上的屈服强度,这是因为晶粒的<111>方向平行于加载方向,滑移系难以启动. 提高γ纤维织构的比例,将增大轧面法向上的屈服强度. 本研究可以为优化FeCrAl合金材料织构、加工条件和材料力学性能提供参考.     

FeCrAl合金材料织构对其宏观力学本构关系的影响研究

FeCrAl合金具有优良的高温抗氧化性和耐辐照性能,是事故容错核燃料包壳的重要候选材料. 其在加工过程和热处理过程中易形成α纤维织构（<110>//RD）和γ纤维织构（<111>//ND）,会影响材料的宏观力学性能与深加工成形能力. 本研究针对具有不同织构的多晶FeCrAl合金,建立了代表性体元模型, 使用晶体塑性有限元方法,在ABAQUS/Explicit中模拟材料单轴加载下的宏观应力应变曲线,分析不同织构对FeCrAl合金宏观力学本构关系的影响. 计算结果表明,对于具有α织构、γ织构和晶粒无择优取向的材料,在轧向上的应力应变曲线差异较小. γ织构会引起材料强烈的各向异性,在轧面法向上的屈服强度远高于轧向和横向上的屈服强度,这是因为晶粒的<111>方向平行于加载方向,滑移系难以启动. 提高γ纤维织构的比例,将增大轧面法向上的屈服强度. 本研究可以为优化FeCrAl合金材料织构、加工条件和材料力学性能提供参考.     

高强度钢板热成形本构理论与实验分析

热成形(热冲压)过程中硼钢的热、力、相变耦合关系是研究热成形理论的基础, 同时也是决定热成形工艺及数值模拟准确性的关键因素. 对热成形硼钢进行高温拉伸及淬火实验: 硼钢板材试样在奥氏体化(950℃)后保温一定时间, 然后在连续冷却的同时施加拉伸力, 记录此过程中力、位移、膨胀量及温度的变化. 通过对不同冷却速率及不同拉伸力情况下上述物理量的变化规律及微观组织性能的分析, 研究硼钢相变过程中的热、力、相变耦合关系. 建立了硼钢相变过程中的热、力、相变耦合模型. 通过引入混合定律对热成形过程中的多相材料热力学参数和力学性能进行等效分析; 对热成形应变组成及其形成机理进行了分析, 引入了相变体积应力及相变塑性应力等新概念. 硼钢高温流动应力采用修改的Norton-Hoff形式, 并通过实验确定了流动应力的材料常数. 在此基础上将热、力、相变耦合关系引入热成形本构方程中, 分别建立了高强度钢板热成形的全量形式及增量形式本构方程. 对U形零部件热成形过程进行了数值模拟, 并与实验结果进行比较, 结果证明建立的本构理论的有效性.      

率相关本构方程积分新算法

提出一种积分率相关本构方程的隐式积分新算法,引入0～1范围内的缩放因子λ对本构方程进行间接求解,可以避免直接求解等效塑性应变或等效塑性应变率时,由于其数值过大或过小而造成的收敛困难或收敛失败,实现对率相关本构方程的快速准确求解.以B-P统一本构方程及双曲正弦本构方程为例,验证了算法的可行性.结果表明,新算法对于准静态变形条件下的无硬化本构方程也可以得出准确的解.     

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

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

海冰动力学中本构模型研究的若干进展

从20世纪70年代, 人们在不同尺度下建立了一系列的海冰动力学本构模型用于海冰数值模拟和预测. 将目前应用的海冰动力学本构模型分为弹塑性、黏塑性、各向异性和颗粒流体动力学中的黏弹塑性模型4类,并分别讨论了各类模型的特点和适用范围. 尤其对在极区及副极区大、中尺度下广泛应用的黏塑性及其改进的本构模型进行了深入地分析. 最后得出: 在大、中尺度下建立海冰动力学的黏弹塑性本构模型是提高海冰动力学计算精度的有效途径; 将小尺度下采用黏弹塑性本构模型的海冰颗粒流体动力学的计算时效进一步改进后, 可在大、中尺度下对海冰的重叠、堆积特性进行模拟; 进一步开展海冰动力学的尺度效应研究, 进行海冰本构模型的实验验证, 并建立不同尺度模型间的相互联系也应是海冰动力学本构模型研究的重要内容.      

任意边界形状的二维刚塑性成形过程的无网格分析方法

将无网格伽辽金方法引入到塑性成形过程模拟,结合刚塑性材料假设,提出了基于刚塑性理论的无网格伽辽金方法,推导了其刚度矩阵方程和求解列式,给出了模具形状任意的二维塑性成形问题摩擦力边界条件的施加方法以及无网格方法应用于任意边界形状的塑性成形问题时的坐标转换关系,建立了无网格方法模拟任意边界形状的塑性成形问题的步骤,并编写了相应的计算程序。应用建立的方法对典型塑性成形过程进行了无网格方法分析,通过与刚塑性有限元方法分析结果的比较,验证了本文所建立方法的可行性。     

考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

PREDICTION OF EARING IN TEXTURED ALUMINIUM SHEETS BASED ON CRYSTAL PLASTICITY

The phenomenon of earing is investigated in the present study based on the theory of crystal plasticity with the dynamic explicit finite element program developed. Firstly texture analysis is carried out of rolled aluminium alloy Al5052 by means of X-ray technique. Then from the texture coefficients an analytical expression for the orientation distribution function (ODF) is derived making use of the computer algebraic language Mathematica4.0, which makes it easier to discretize the ODF into a series of Eulerian angles representing the distribution of lattices and further the preferred orientation (texture) of crystals of the original sheets. For the polycrystal model, the material is described using crystal plasticity where each material point in grains with each grain modelled as an FCC crystal with 12 distinct slip systems. The modified Taylor theory of crystal plasticity is used and only the initial texture is taken into consideration during large plastic deformation. Numerical simulation of earing has been performed for an aluminium sheet with texture and one with crystals exhibiting random distribution to demonstrate the effect of texture of materials on their plastic anisotropy and formability. Project supported by the National Natural Science Foundation of China (No. 59875025).     

A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

Finite element analysis of localization in FCC polycrystalline sheets under plane stress tension

Localization phenomena in thin sheets subjected to plane stress tension are investigated. The sheet is modelled as a polycrystalline aggregate, and a finite element analysis based on rate-dependent crystal plasticity is developed to simulate large strain behaviour. Accordingly, each material point in the specimen is considered to be a polycrystalline aggregate consisting of a large number of FCC grains. The Taylor model of crystal plasticity theory is assumed. This analysis accounts for initial textures as well as texture evolution during large plastic deformations. The numerical analysis incorporates certain parallel computing features. Simulations have been carried out for an aluminum sheet alloy, and the effects of various parameters on the formation and prediction of localized deformation (in the form of necking and/or in-plane shear bands) are examined.     

Finite element modeling of tube hydroforming of polycrystalline aluminum alloy extrusions

Finite element modeling of tube hydroforming requires information about the anisotropy of the extruded aluminum tube. Unlike sheet metals, the complex geometry of extruded tubes makes it difficult, except in extrusion direction, to directly measure material properties. Therefore, polycrystalline models provide a good alternative for calculating the anisotropy of the tube in all directions and under various loading conditions. Using a rate-independent single crystal yield surface and rigid plasticity, a Taylor-type polycrystalline model was developed and implemented into ABAQUS/Explicit finite element (FE) code using VUMAT. The constitutive model was then used to calculate the crystallographic texture evolution during the hydroforming of an extruded aluminum tube. Initial crystallographic texture measured using orientation imaging microscopy (OIM) and uniaxial tensile test data obtained along the extrusion direction were input to this FEA model. In order to efficiently and practically simulate the tube hydroforming process using the polycrystalline model, sensitivity to the number of grain orientation, total simulation time, and number of finite elements were studied. Predicted results agreed very well with experimentally measured strain obtained from tube hydroforming process.     

镁合金晶体塑性本构模型与非均匀孪生变形分析

微结构演化对镁合金材料力学性能有着显著的影响,为了揭示镁合金宏观塑性各向异性特性与非均匀孪生变形的关系,开展了不同路径下的单轴加载试验以及采用含滑移、孪生机制的晶体塑性本构模型对试验条件下的镁合金变形行为进行数值模拟研究。文中本构模型描述了滑移与孪生变形机制以及晶格转动的机制,同时研究采用三维微结构代表性有限元模型,其包含晶粒尺寸、晶向和晶界倾角等微结构参数。研究结果表明,轧制镁合金具有强烈的宏观塑性各向异性行为,并对这种镁合金塑性各向异性行为的模拟结果以及多晶织构的模拟演化结果与试验测量进行对比,结果都基本吻合。对孪生非均匀变形模拟分析表明,镁合金宏观塑性各向异性行为与滑移、孪生变形机制的不同启动组合紧密相关,同时多晶体内应力的非均匀分布受到孪生变形的严重影响。而不同晶粒尺寸的晶粒所发生的孪生变形有比较大的差异,造成孪晶变体在晶粒内的分布极不均匀。本研究可为通过微结构的合理配置来设计和控制材料的力学性能提供理论依据.     

NUMERICAL SIMULATION OF NON-PLANAR PLASTIC ANISOTROPY OF COLD ROLLING POLYCRYSTALLINE MATERIALS

The plastic anisotropy of sheet metal is usually caused by preferred orientation of grains, developed by mechanical deformation and thermal treatment. In the present study, a Taylor-like polycrystal model suggested by Asaro and Needleman is applied to investigate the evolution of the anisotropic behavior of a face centered cubic (FCC) polycrystalline metal, which is considered having {111} (110) slip systems, by stretching it along an arbitrary direction after it has undergonea plane-strata compression that rationally simulates the cold rolling process of FCC polycrystalline pure aluminium. By using the Taylor-like polycrystal model, pole figures are obtained to describe the texture development of polycrystalline aggregate after plane-strain compression, and then the plastic anisotropy of polycrystalline aggregate is evaluated by stretching the polycrystalline aggregate in different direction in term of yield stress. According to the results, the contours of longitudinal flow stress in three-dimensional orientation space are given and analyzed. Experiment results similar to the prediction of planar anisotropy can be found inthe literature written by Takahashi et al. that in directly show the correctness of the prediction of non-planar plastic anisotropy by this analysis.     

On the modeling and simulation of induced anisotropy in polycrystalline metals with application to springback

The presence of initial, and the development of induced, anisotropic elastic and inelastic material behavior in polycrystalline metals, can be traced back to the influence of texture and dislocation substructural development on this behavior. As it turns out, via homogenization or other means, one can formulate effective models for such structure and its effect on the macroscopic material behavior with the help of the concept of evolving structure tensors. From the constitutive point of view, these quantities determine the material symmetry properties. Most importantly, all dependent constitutive fields (e.g., stress) are by definition isotropic functions of the independent constitutive variables, which include these evolving structure tensors. The evolution of these tensors during loading results in an evolution of the anisotropy of the material. From an algorithmic point of view, the current approach leads to constitutive models which are quite amenable to numerical implementation. To demonstrate the applicability of the resulting constitutive formulation, we apply it to the case of metal plasticity with combined hardening involving both deformation- and permanently induced anisotropy. Comparison of simulation results based on this model for the bending tension of aluminum-alloy sheet-metal strips with corresponding experimental ones show good agreement.     

On the r-value of textured sheet metals

The r-value of a sheet metal is a measure of plastic anisotropy frequently used for prediction of performance in deep-drawing. It has also figured prominently in the literature for validation of theories where the predicted angular dependence of r is compared with the measured dependence. As plastic anisotropy in sheet metals is caused mainly by the preferred orientations of grains within the polycrystalline metal, it is natural to ask how r would depend on the orientation distribution function (ODF) w which defines the crystallographic texture of the polycrystal. In this paper a general formula relating r to w is derived for textured sheet metals whose plastic flow behavior is governed by a plastic potential f(σ, w), the anisotropic part of which depends linearly on the texture coefficients; here σ denotes the deviator of the Cauchy stress. Specific forms of this formula for orthorhombic sheets of cubic and of hexagonal metals are explicitly given.