A nonclassical constitutive model for crystal plasticity and its application

I.IntroductionTheresearchoncrystalplasticitycanbedatedbackto1930'sill.AherHill[=],HillandRicely]builtupaperfectsystemofthegeometryandkineticsofCrystalplasticity,itsapplicatiollbecolllesmoreandmoreattractive,especiallyintheallalysisofpolycrystallinelllaterialssubjectedtoInultiaxiallynonproportionalcyclicloading.Ithasalreadybeenrealizedthattileillll,ol'talltforthepracticalapplicationis:tofindal.ealisticalldeasilyapplicableof'ystallineconstitutiverelationandaneffectivenumericalapproach.Tilecon…     

HILL屈服准则与晶体塑性模型对FCC单晶材料塑性各向异性描述能力的比较

采用宏观HILL模型和晶体塑性模型对面心立方单晶(FCC)材料的非均匀交形进行了数值模拟，意在比较两种不同尺度的模型对塑性各向异性的描述能力的差异。为了使两种模型具有可比性，对于FCC单晶材科，本文提出一种用晶体塑性模型来确定HILL模型中各向异性参数的标定方法。数值分析表明，两类模型对单晶体塑性各向异性的描述能力存在着差异。对FCC单晶材料，HILL模型对各向异性的预测能力没有晶体塑性模型细致，晶体塑性模型更能追踪塑性各向异性的变化。但两种模型对应力应变响应预测的趋势是一致的。对两种模型描述的差异，做了详细的分析。     

晶粒数量对多晶集合体初始各向异性的影响

Taylor类多晶晶体粘塑性模型被用于研究晶粒数量对随机分布多晶体拉伸塑性各向异性的影响。分别沿包含不同晶粒数量的多晶集合体的各方向进行单向拉伸数值模拟实验,得到多晶集合体各方向在一定等效应变下的等效应力,并用云图和等高线表示在多晶体的参考球面上。定义了描述多晶集合体各向异性程度的参考指标。讨论了三种确定晶体随机取向的方法。计算结果表明：晶粒数量有限的多晶集合体的应力应变响应仍有一定的各向异性,且随着晶粒数量增多,多晶集合体的各向异性程度降低;就所包含晶粒数相同的多晶集合体来说,在确定晶粒随机取向时,选取不同的方法对它的各向异性程度也有一定的影响。     

An endochronic constitutive model for rate-dependent and nonproportional cyclic plasticity

The total stress response of material is decomposed into a sum of an equilibrium stress response and a non-equilibrium overstress response. Correspondingly, the rate-independent intrinsic time and the rate-dependent intrinsic time are defined respectively. Additional hardening functions for describing the isotropic and anisotropic nonproportional effects are assumed to be related to the accumulation of plastic strain component along the normal of equilibrium stress trajectory, in which the effects of geometry of the loading path are included. An endochronic constitutive model for rate-dependent, nonproportional cyclic plasticity is formulated and applied to simulate the stress responses of stainless steel XCrNi18.9 for some typical loading programs at different loading rates. A comparison between predicted results and experimental ones by Haupt and Lion shows that the former are in agrreement with the latter.     

层合板冲击损伤的PD率效应本构模型分析

从近场动力学(简称PD)理论的PMB材料模型出发,结合Kelvin-Voigt粘弹性固体模型,建立PD率效应本构模型。采用LAMMPS软件模拟了环氧树脂板、纤维增强复合材料单向层板和多向铺层层合板受冲击的情况。通过分析各板的冲击损伤,探索纤维对板的增强作用。同时,分析了不同冲击速度下层合板上下表面的损伤程度,初步探讨了从低速碰撞到高速冲击过程中复合材料层合板的破坏机理及规律。     

A multi-dimensional composite model for plastic continua under polyaxial loading condition

By replacing a medium with reinforcing components oriented and distributed uniformly in a multi-dimensional space, a constitutive model is constructed. The components are extended/compressed compatibly with the strain and the resultant of load exerted on them to balance the stress. Their load-elongation relation can be determined from a conventional material test. Each component undergoes different elongation history depending on its own orientation during deformation, so that the model can simulate elasto-plastic behavior of materials under polyaxial loading conditions. The incremental constitutive matrix has been derived for application in numerical analysis and a yield criterion is also introduced. A few subsequent yield surfaces have been predicted and compared with experiments.     

Multiscale modelling of the induced plastic anisotropy in bcc metals

This paper presents a new framework to predict the qualitative and quantitative variation in local plastic anisotropy due to crystallographic texture in body-centered cubic polycrystals. A multiscale model was developed to examine the contribution of mesoscopic and local microscopic behaviour to the macroscopic constitutive response of bcc metals during deformation. The model integrated a dislocation-based hardening scheme and a Taylor-based crystal plasticity formulation into the subroutine of an explicit dynamic FEM code (LS-DYNA). Numerical analyses using this model were able to predict not only correct grain rotation during deformation, but variations in plastic anisotropy due to initial crystallographic orientation. Optimal results were obtained when {1 1 0}〈1 1 1〉, {1 1 2}〈1 1 1〉, and {1 2 3}〈1 1 1〉 slip systems were considered to be potentially active. The predicted material heterogeneity can be utilised for research involving any texture-dependent work hardening behaviour, such as surface roughening.     

Investigation on elasto-plastic constitutive model coupled with damage for localization phenomena

On the basis of existing plasticity-based damage model for plasticity coupled with damage for localization analysis, constitutive parameter identification was carried out through a series of numerical tests at local level. And then improvements were made on the expressions of the evolution laws of damage. Strain localization phenomena were simulated with a typical double-notched specimen under tensions. Numerical results indicate the validity of the proposed theory.     

Criterion for judging satisfaction of consistency conditions in rate-dependent constitutive relations

Based on irreversible thermodynamics, the criterion for judging the satisfaction of consistency conditions in rate-dependent constitutive relationship is deduced by introducing four basic hypotheses. Formulas for solving internal variables are given. It makes the rate-dependent model applicable no matter whether the consistency conditions can be satisfied or not. This research was funded by the National Natural Science Foundation of China (No. 59739180) and the Foundation of National Key Laboratory of Coastal and Offshore Engineering of Dalian University of Technology (No. 9701).     

Study on the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy based on a new crystal plasticity constitutive model

In this paper, a crystal plasticity based constitutive model (Yu et al., 2013) is extended to describe the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy by considering the internal heat production. Two sources of internal heat productions are included in the proposed model, i.e., the mechanical dissipations of inelastic deformation and the transformation latent heat in the NiTi shape memory alloy. With an assumption of uniform temperature field in the alloy specimen, a simplified evolution law of temperature field is obtained by the first law of thermodynamics and the heat boundary conditions. An explicit scale-transition rule is adopted to extend the proposed single crystal model to the polycrystalline version. The capability of the extended polycrystalline model to describe the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy is verified by comparing the predictions with the corresponding experimental ones. The comparison demonstrates that the proposed constitutive model considering the internal heat production predicts the rate-dependent cyclic deformation of super-elastic NiTi shape memory alloy fairly well.     

A dislocation density-based single crystal constitutive equation

Single crystal constitutive equations based on dislocation density (SCCE-D) were developed from Orowan’s strengthening equation and simple geometric relationships of the operating slip systems. The flow resistance on a slip plane was computed using the Burger’s vector, line direction, and density of the dislocations on all other slip planes, with no adjustable parameters. That is, the latent/self-hardening matrix was determined by the crystallography of the slip systems alone. The multiplication of dislocations on each slip system incorporated standard 3-parameter dislocation density evolution equations applied to each slip system independently; this is the only phenomenological aspect of the SCCE-D model. In contrast, the most widely used single crystal constitutive equations for texture analysis (SCCE-T) feature 4 or more adjustable parameters that are usually back-fit from a polycrystal flow curve. In order to compare the accuracy of the two approaches to reproduce single crystal behavior, tensile tests of single crystals oriented for single slip were simulated using crystal plasticity finite element modeling. Best-fit parameters (3 for SCCE-D, 4 for SCCE-T) were determined using either multiple or single slip stress–strain curves for copper and iron from the literature. Both approaches reproduced the data used for fitting accurately. Tensile tests of copper and iron single crystals oriented to favor the remaining combinations of slip systems were then simulated using each model (i.e. multiple slip cases for equations fit to single slip, and vice versa). In spite of fewer fit parameters, the SCCE-D predicted the flow stresses with a standard deviation of 14 MPa, less than one half that for the SCCE-T conventional equations: 31 MPa. Polycrystalline texture simulations were conducted to compare predictions of the two models. The predicted polycrystal flow curves differed considerably, but the differences in texture evolution were insensitive to the type of constitutive equations. The SCCE-D method provides an improved representation of single-crystal plastic response with fewer adjustable parameters, better accuracy, and better predictivity than the constitutive equations most widely used for texture analysis (SCCE-T).     

混凝土WW三参数率相关动态本构模型

根据一致粘塑性模型理论,结合所做的混凝土单轴动力特性试验,对常用的混凝土本构模型Willam—Warnke三参数模型进行改进,引入应变率的影响,推导了Willam—Warnke材料的动态率相关本构模型;并与试验的应力应变曲线进行比较,结果表明本文模型能够较好地反应混凝土的动力特性;最后,通过对一受冲击荷载的简支梁进行计算,并与线弹性及率无关塑性结果进行比较,结果表明考虑了应变率的影响后梁的动力特性发生了较大的变化。     

水泥砂浆的一个热粘弹性率型损伤本构模型

利用SHPB实验系统及自行研制的混凝土类材料快速高温加热设备,对水泥砂浆试件进行了不同 温度(20~600℃)和3种冲击速度下的实验,得到了不同温度和冲击速度下水泥砂浆试件的应力应变关系曲 线。基于ZWT粘弹性本构模型,并且考虑高温下水泥砂浆损伤演化规律都服从Weibull分布,提出了一个水 泥砂浆的热粘弹性率型损伤本构模型。通过数据拟合,获得了本构模型的相关参数,结果表明:理论预测和实 验结果吻合良好。     

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.     

A hierarchical model for rate-dependent polycrystals

A hierarchical model of a polycrystalline aggregate of rigid viscoplastic grains is formulated, and a robust and efficient computational algorithm for its solution is proposed. The polycrystalline aggregate is modeled as a binary tree. The leaves of the binary tree represent grains, and higher tree nodes represent increasingly larger sub-aggregates of grains. The root of the tree represents the entire polycrystalline aggregate. Velocity and traction continuity are enforced across the interface between the children of each non-leaf node in the binary tree. The hierarchical model explicitly models intergranular interactions but is nevertheless comparable in computational effort to the mean field models of polycrystal plasticity. Simulations of tensile, compressive, torsional, and plane strain deformation of copper lead to predictions in good agreement with experiments, and highlight the interconnection between grain deformations and intergranular constraints. It is inferred from the results that a hybrid mean field/hierarchical model represents a computationally efficient methodology to simulate polycrystal deformation while accounting for intergranular interactions.     

单晶有限变形的热力耦合计算模型

以弹性变形梯度作为基本变量,结合热力学理论构造了单晶有限变形的热、力耦合计算模型。该模型考虑了温度、变温速率以及塑性耗散等条件对单晶有限变形的影响,相对于传统的以弹性变形梯度为基本变量的晶体塑性模型,算法能够体现温度效应的影响。采用隐式的积分方法对建立的控制方程进行计算以保证求解过程的稳定。以1100Al单晶为例计算了不同升温、降温速率,以及不同应变率影响下的材料应力-应变的响应。结果表明,模型能较好地反映变温过程中,单晶各向异性性质的演化以及应力、应变之间关系的变化。     

A constitutive equation of elastic-visco-plasticity for polymers

The rheological behavior and a constitutive relationship of elastic-visco-plasticity for polymers are investigated in this paper. Several sets of experiments have been carried out to determine the material constants and to test the validity of constitutive formulation. It is shown that the theoretical profiles are in good agreement with experimental results. The rheological characteristic and the strain-rate effect of model are analyzed by computer simulation.The project is supported by the National Natural Science Foundation of China     

A novel constitutive model for semiconductors: The case of silicon

The photovoltaic industry relies heavily on solar-grade silicon multicrystals. Understanding their mechanical behavior requires the development of adequate constitutive models accounting for the effects of both high dislocation densities and complex loading situations in a wide range of temperature, strain rate, and impurity contents. The traditional model of Alexander and Haasen poses several limitations. We introduce in this work a novel constitutive model for covalent single crystals and its implementation into a rate-dependent crystal plasticity framework. It is entirely physically based on the dislocation generation, storage and annihilation processes taking place during plastic flow. The total dislocation density is segmented according to the dislocation mobility potential and their character. A dislocation multiplication law for the yield region more accurate than the one of Alexander and Haasen is proposed. The influence of additional dislocation sources created on forest trees, usually disregarded in models for semiconductors, is assessed. The dislocation velocity law combines three potentially rate-limiting mechanisms: the standard double kink mechanism, jog dragging and the influence of localized obstacles. The model is valid at finite strains, in multiple slip conditions and captures accurately the high temperature- and strain rate sensitivity of semiconductors. The experimental stress overshoot is qualitatively reproduced only when jog dragging is accounted for. Localized obstacles are shown not to have any significant effect on dislocation motion in silicon. The broader case of extrinsic semiconductors is discussed and the influence of dissolved oxygen on the upper yield stress of silicon monocrystals is successfully reproduced.