Technique for selecting the functions of the constitutive equations of creep and long-term strength with one scalar damage parameter

The strain-strength characteristics of aerostructures made of hardening materials under uniaxial tension in creep conditions are determined. The problem is reduced to a system of ordinary differential equations of the kinetic theory of creep with one scalar damage parameter. The approximate solutions of the problem are obtained with the help of the implicit Euler method and of the arc length method in combination with the explicit methods of the Runge–Kutta family for cylindrical St.45 steel samples and 3V titanium alloy plates.     

考虑微观结构特征长度演化的内变量黏塑性本构模型

在金属晶体材料高应变率大应变变形过程中,存在强烈的位错胞尺寸等微观结构特征长度细化现象,势必对材料加工硬化、宏观塑性流动应力产生重要影响。基于宏观塑性流动应力与位错胞尺寸成反比关系,提出了一种新型的BCJ本构模型。利用位错胞尺寸参数,修正了BCJ模型的流动法则、内变量演化方程,引入了考虑应变率和温度相关性的位错胞尺寸演化方程,建立了综合考虑微观结构特征长度演化、位错累积与湮灭的内变量黏塑性本构模型。应用本文模型,对OFHC铜应变率在10-4~103 s-1、温度在298~542 K、应变在0~1的实验应力-应变数据进行了预测。结果表明:在较宽应变率、温度和应变范围内,本文模型的预测数据与实验吻合很好;与BCJ模型相比,对不同加载条件下实验数据的预测精度均有较大程度的提高,最大平均相对误差从9.939%减小为5.525%。     

Viscoplastic dynamics of isotropic plates of variable thickness under explosive loading

A problem of viscoplastic dynamic bending of isotropic plates of variable thickness is formulated. A method for integrating the initial-boundary problem is developed. Numerical results are compared with a known analytical solution obtained within a rigid-plastic model; good agreement is demonstrated. The efficiency of the method developed is verified by numerical computations. It is shown that the final flexure of plates can be reduced severalfold by applying rational design. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 123–134, March–April, 2007.     

混凝土的一种标量损伤弹塑性本构模型

荷载作用下材料性能的劣化是混凝土结构的微观损伤机理,其宏观表现为结构刚度的折减和承载力的降低。论文推导了基于不可逆热力学过程的弹塑性标量损伤本构,给出时间离散的屈服准则。采用基于向后Euler法的应力更新算法——两步图形返回的最近点投影法,推导了满足迭代结果收敛假设的塑性参数及算法刚度张量,给出了空间梁单元本构积分算法的Jacobi矩阵。将模型用于混凝土简支梁的承载力试验模拟,与计算数据对比表明了模型和算法的合理性和有效性。     

Micromechanics-based constitutive modeling for unidirectional laminated composites

A micromechanics-based constitutive model is developed to predict the effective mechanical behavior of unidirectional laminated composites. A newly developed Eshelby’s tensor for an infinite circular cylindrical inclusion [Cheng, Z.Q., Batra, R.C., 1999. Exact Eshelby tensor for a dynamic circular cylindrical inclusion. J. Appl. Mech. 66, 563–565] is adopted to model the unidirectional fibers and is incorporated into the micromechanical framework. The progressive loss of strength resulting from the partial fiber debonding and the nucleation of microcracks is incorporated into the constitutive model. To validate the proposed model, the predicted effective stiffness of transversely isotropic composites under far field loading conditions is compared with analytical solutions. The constitutive model incorporating the damage models is then implemented into a finite element code to numerically characterize the elastic behavior of laminated composites. Finally, the present predictions on the stress–strain behavior of laminated composite plate containing an open hole is compared with experimental data to verify the predictive capability of the model.     

Rheological characterization and constitutive modeling of bread dough

Bread dough (a flour–water system) has been rheologically characterized using a parallel-plate, an extensional, and a capillary rheometer at room temperature. Based on the linear and nonlinear viscoelastic and viscoplastic data, two constitutive equations have been applied, namely a viscoplastic Herschel–Bulkley model and a viscoelastoplastic K–BKZ model with a yield stress. For cases where time effects are unimportant, the viscoplastic Herschel–Bulkley model can be used. For cases where transient effects are important, it is more appropriate to use the K-BKZ model with the addition of a yield stress. Finally, the wall slip behavior of dough was studied in capillary flow, and an appropriate slip law was formulated. These models characterize the rheological behavior of bread dough and constitute the basic ingredients for flow simulation of dough processing, such as extrusion, calendering, and rolling.     

Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.     

An RVE-based multiscale modeling method for constitutive relations

This article reports an efficient method to characterize constitutive responses based on multiscale modeling for fluid flow in heterogeneous media based on the concept of representative volume element (RVE). Between different scales, it is considered as the basic principles for down-scaling information the conservation of velocity and of the strain rate tensor. Within this context, we formulate (i) the problem to be solved at the micro-scale, (ii) the up-scaling procedure which involves homogenization rules, and (iii) the generalized principle of multiscale virtual power. The complete theory for constitutive modeling is revisited and shown that when employing multiscale analysis among the suitable variational arguments we are able to obtain, in a straightforward manner, new constitutive behavior between kinematic motions and actions. Some examples of application of fluid flow in heterogeneous media with obstacles are presented to show the consequences of the proposed approach.     

Characteristic solutions of scalar newton equations in one space dimension

Newton equations are dynamical systems on the space of fields. The solutions of a given equation which are curves of characteristic fields for its force are planar and have constant angular momentum. Separable solutions are characteristic with angular momentum equal to zero. A Newton equation is separable if and only if its characteristic equation is homogeneous. Separable equations correspond to invariants of homogeneous ordinary differential equations, and those associated with a given homogenous equation correspond to its generalized dilation symmetries. A Newton equation is compatible with the characteristic condition if and only if its characteristic equation is linear. Such equations correspond to invariants of linear ordinary differential equations. Those associated with a given linear equation correspond to the central force problems on its solution space. Regardless of compatibility, any Newton equation with a plane of characteristic fields has non-separable characteristic solutions.     

Stress-strain state of flexible cylindrical panels with variable geometric parameters

International Applied Mechanics -     

Multi-cycle deformation of silicone elastomer: observations and constitutive modeling with finite strains

Observations are reported on a medical grade of silicone elastomer in uniaxial tensile tests up to breakage of specimens, short-term relaxation tests, and cyclic tests with monotonically increasing maximum elongation ratios. Experimental data in cyclic tests demonstrate the fading memory phenomenon: stress–strain diagrams for two specimens with different deformation histories along the first n?1 cycles and coinciding loading programs for the other cycles become identical starting from the nth cycle. A constitutive model is developed in cyclic viscoplasticity of elastomers with finite strains, and its adjustable parameters are found by fitting the experimental data. Ability of the stress–strain relations to predict the mechanical response in cyclic tests with various deformation programs is confirmed by numerical simulation.     

Multiscale constitutive modeling and numerical simulation of fabric material

A multiscale model for a fabric material is introduced. The model is based on the assumption that on the macroscale the fabric behaves as a continuum membrane, while on the microscale the properties of the microstructure are accounted for by a constitutive law derived by modeling a pair of overlapping crimped yarns as extensible elasticae. A two-scale finite element method is devised to solve selected boundary-value problems.     

Elliptic yield cap constitutive modeling for high porosity sandstone

Field and laboratory investigators have observed thin, tabular zones of localized compressional deformation without shear in high porosity sandstone. These ‘compaction bands’ display greatly reduced porosity, and may affect the withdrawal of fluids from reservoirs. Studies addressing band formation as a type of strain localization predict the onset of the bands in a range of constitutive parameters roughly consistent with experiments, but are highly dependent on the constitutive relation used. In particular, the hardening modulus in shear and the slope of the yield surface in a plot of shear stress versus mean compressive stress are critical to localization predictions. Previous yield cap constitutive models employed a single deformation criterion, linking hydrostatic and shear response. In this work, we propose an elliptic yield cap model employing separate inelastic deformation parameters along each axis of the ellipse. The two deformation parameters allow the proposed surface to change in aspect ratio as it deforms, and allow a negative hardening modulus in shear without a negative hydrostatic modulus. Some cases with simplified modeling are shown for illustrative purposes, followed by a comparison with existing models. The proposed model displays similar strain behavior to the other models, but predicts localization under less restrictive conditions.     

Multi-dimensional constitutive modeling of SMA during unstable pseudoelastic behavior

A multi-dimensional continuum-level constitutive model of shape memory alloys exhibiting pseudoelastic behavior has been developed in this paper. The current model is an extension of the one-dimensional model previously developed by the authors, which consists of a constitutive relation and a transformation evolution rule (kinetic relation). The constitutive relation is constructed based on the gradient of a transformation potential function (effective stress), and the kinetics of transformation are expressed based on a set of transformation surfaces in stress-temperature space. The constitutive relation allows for the reorientation of the transformation strain tensor according to the current state of stress. The model is used to simulate the localized deformation and transformation front propagation in NiTi strip under quasi-static extension. Special attention has been paid to the multi-axiality of the stress state at the transformation front during the forward and reverse transformations. It is shown that the ability of the model to consider the reorientation of the transformation strain tensor results in full recovery of transformation strain upon unloading, and a uniform distribution of stress in the transformed areas in spite of the localization of transformation. It also produces a stress gradient that dictates the direction of the transformation front propagation during both the forward and reverse transformations.     

Efficient and robust constitutive integrators for single-crystal plasticity modeling

Small-scale deformation phenomena such as subgrain formation, development of texture, and grain boundary sliding require simulations with a high degree of spatial resolution. When we consider finite-element simulation of metal deformation, this equates to thousands or hundreds of thousands of finite elements. Simulations of the dynamic deformations of metal samples require elastic–plastic constitutive updates of the material behavior to be performed over a small time step between updates, as dictated by the Courant condition. Further, numerical integration of physically-based equations is inherently sensitive to the step in time taken; they return different predictions as the time step is reduced, eventually approaching a stationary solution. Depending on the deformation conditions, this converged time step becomes short (10⁻⁹ s or less). If an implicit constitutive update is applied to this class of simulation, the benefit of the implicit update (i.e., the ability to evaluate over a relatively large time step) is negated, and the integration is prohibitively slow. The present work recasts an implicit update algorithm into an explicit form, for which each update step is five to six times faster, and the compute time required for a plastic update approaches that needed for a fully-elastic update. For dynamic loading conditions, the explicit model is found to perform an entire simulation up to 50 times faster than the implicit model. The performance of the explicit model is enhanced by adding a subcycling algorithm to the explicit model, by which the maximum time step between constitutive updates is increased an order of magnitude. These model improvements do not significantly change the predictions of the model from the implicit form, and provide overall computation times significantly faster than the implicit form over finite-element meshes. These modifications are also applied to polycrystals via Taylor averaging, where we also see improved model performance.     

On the constitutive modeling of coupled thermomechanical phase-change problems

This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient and robust numerical model, able for the numerical simulation of industrial solidification processes. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. All the material properties have been considered to be temperature dependent. The constitutive behavior has been defined by a thermoviscous/elastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. The constitutive model proposed accounts for a continuous transition between the initial liquid state, the intermediate mushy state and the final solid state taking place in a solidification process. In particular, a pure viscous deviatoric model has been used at the initial fluid-like state. A thermomecanical contact model, including a frictional hardening and temperature dependent coupled potential, is derived within a fully consistent thermodinamical theory. The numerical model has been implemented into the computational finite element code COMET developed by the authors. Numerical simulations of solidification processes show the good performance of the computational model developed.     

爆炸力学数值模拟中本构建模问题的讨论

就爆炸力学数值模拟中的本构建模问题的有关内容作了讨论,包括:波传播研究与材料的动态特性研究的相互依赖、容变律与畸变律的解耦与耦合、率型本构关系与失效准则、应变率效应与温度效应的等效性、加载本构关系与卸载本构关系以及流变过程与损伤演化过程的耦合等,并展望了今后的有关发展动向。     

Multi-cycle deformation of semicrystalline polymers: Observations and constitutive modeling

Experimental data are reported on isotactic polypropylene in multi-cycle uniaxial tensile tests where a specimen is stretched up to some maximum strain and retracted down to the zero minimum stress, while maximum strains increase with number of cycles. Fading memory of deformation history is observed: when two samples are subjected to loading programs that differ along the first n − 1 cycles only, their stress–strain diagrams coincide starting from the nth cycle. Constitutive equations are developed in cyclic viscoelasticity and viscoplasticity of semicrystalline polymers, and adjustable parameters in the stress–strain relations are found by fitting the experimental data. Results of numerical simulation demonstrate that the model predicts the fading memory effect quantitatively. To confirm that the observed phenomenon is typical of semicrystalline polymers, experimental data are presented in tensile cyclic tests with large strains on low density polyethylene and compressive cyclic tests on poly(oxymethylene).