Implicit stress-update algorithm for isotropic Cam-clay model based on the subloading surface concept at finite strains

A stress-update algorithm for the subloading surface Cam-clay model is developed. The model is re-formulated to be compatible with the multiplicative hyperelasto-plasticity for finite strains. The algorithm aims at synthesizing the phenomenological capability of the subloading surface model and the numerical advantage of the return mapping. A closed-form representation of the algorithmic tangent is also derived. The proposed algorithm is implemented into a finite-element code and is shown to be capable of producing accurate results even for large incremental steps. Numerical examples of boundary-value problems demonstrate the robustness of the algorithm.     

Finite elasto-viscoplastic modeling of polymers including asymmetric effects

Asymmetric effects between compression and tension are a pronounced behavior for glassy polymers such as polycarbonate. For its simulation an elasto-viscoplastic framework is formulated within a geometrically nonlinear theory. Here a new approach within the concept of stress mode dependent weighting functions is used, where each material parameter is additively decomposed into a sum of weighted stress mode-related quantities. The characterization of the stress modes is obtained in the octahedral plane of the deviatoric stress space in terms of the mode angle, such that stress mode dependent scalar weighting functions can be constructed. The constitutive equations are formulated for large strains in terms of logarithmic Hencky strains and its work conjugated Hill stresses. The resulting evolution equations are updated using a semi-implicit Euler scheme, and the algorithmic tangent operator is derived for the finite element equilibrium iteration. The numerical implementation is also used to identify the material parameters thus resulting into a good agreement with experimental data. Furthermore, the model is used to simulate the cold drawing processes for a dumbbell-shaped specimen in tension and a perforated strip in compression and tension.     

非饱和多孔介质中混合元法的化学-热-渗流-力学耦合的本构模拟

在文献[1]中建立的多孔介质中化学-热-渗流-力学(CTHM)本构模型基础上,针对文献[2]建立的非饱和多孔介质中热-渗流-力学耦合分析的混合有限元方法,发展了非饱和多孔介质中混合元的化学-热-渗流-力学(CTHM)耦合本构模拟算法。采用非关联流动多重屈服准则模拟非饱和多孔介质的材料非线性行为。推导了u-pw-pa-T形式的包含了耦合率本构方程积分的向后欧拉映射算法和一致性弹塑性切线模量矩阵(单元刚度矩阵)的混合元一致性算法。本文给出了临界状态线(CSL)和状态边界面(SBS)两个屈服准则的一致性算法。数值结果显示了本文所发展的混合元耦合本构模拟算法在模拟由热、化学、力学荷载共同引起的多孔介质中化学-热-渗流-力学(CTHM)耦合行为的能力和有效性。     

On the plastic bifurcation and post-bifurcation of axially compressed beams

The paper presents two new results in the domain of the elastoplastic buckling and post-buckling of beams under axial compression. (i) First, the tangent modulus critical load, the buckling mode and the initial slope of the bifurcated branch are given for a Timoshenko beam (with the transverse shear effects). The result is derived from the 3D J₂ flow plastic bifurcation theory with the von Mises yield criterion and a linear isotropic hardening. (ii) Second, use is made of a specific method in order to provide the asymptotic expansion of the post-critical branch for a Euler-Bernoulli beam, exhibiting one new non-linear fractional term. All the analytical results are validated by finite element computations.     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

THIRD ORDER SHEAR DEFORMATION MODEL FOR LAMINATED SHELLS WITH FINITE ROTATIONS： FORMULATION AND CONSISTENT LINEARIZATION

The paper presents an approach for the formulation of general laminated shells based on a third order shear deformation theory. These shells undergo finite (unlimited in size) rotations and large overall motions but with small strains. A singularity-free parametrization of the rotation field is adopted. The constitutive equations, derived with respect to laminate curvilinear coordinates, are applicable to shell elements with an arbitrary number of orthotropic layers and where the material principal axes can vary from layer to layer. A careful consideration of the consistent linearization procedure pertinent to the proposed parametrization of finite rotations leads to symmetric tangent stiffness matrices. The matrix formulation adopted here makes it possible to implement the present formulation within the framework of the finite element method as a straightforward task.     

Associative coupled thermoplasticity at finite strain with temperature-dependent material parameters

The paper deals with associative coupled thermoplasticity at finite strains. J2 plasticity model is based on hyperelastic formulation with multiplicative decomposition of deformation gradient into elastic and plastic parts. As a novel aspect, temperature dependence of all material parameters is introduced. For such model, variational procedure is carried out and discretization by the mixed finite element method is performed. Consistent linearization is carried out and algorithmic elasto-plastic tangent moduli are given. Verification of algorithm is provided by two examples.     

A note on nonassociated plastic flow rules

The analytical properties of the constitutive equations in plasticity with a nonassociated flow rule are investigated. Under the assumption of small deformations the directional stiffness (and compliance) rule is considered and the relevant spectral properties of the tangent stiffness tensor are assessed. It is shown that the directional stiffness may be larger than the elastic. It may also be negative in the case of a formally perfectly plastic material and so the nonassociative flow rule represents (spurious) softening in terms of an associated flow rule. The issue of uniqueness at finite strains is briefly addressed, whereby use is made of the tangent stiffness tensor relating the velocity gradient to the first Piola-Kirchhoff stress rate. The relevant spectral properties, which generalise those from the small deformation case, are found explicit. A sufficient condition for uniqueness is given in terms of a critical (upper bound) value of the hardening modulus.     

Simulation of rate dependent plasticity for polymers with asymmetric effects

Polycarbonate is an amorphous polymer which exhibits a pronounced strength-differential effect between compression and tension. Also strain rate and temperature influence the mechanical response of the polycarbonate. The concept of stress mode dependent weighting functions is used in the proposed model to simulate the asymmetric effects for different loading speeds. In this concept, an additive decomposition of the flow rule is assumed into a sum of weighted stress mode related quantities. The characterization of the stress modes is obtained in the octahedral plane of the deviatoric stress space in terms of the mode angle, such that stress mode dependent scalar weighting functions can be constructed. The resulting evolution equations are updated using a backward Euler scheme and the algorithmic tangent operator is derived for the finite element equilibrium iteration. The numerical implementation of the resulting set of constitutive equations is used in a finite element program for parameter identification. The proposed model is verified by showing a good agreement with the experimental data. After that the model is used to simulate the laser transmission welding process.     

Effect of strain hardening in elastic–plastic transition behavior in a hemisphere in contact with a rigid flat

An axisymmetrical hemispherical asperity in contact with a rigid flat is modeled for an elastic–plastic material on the lines of the Kogut–Etsion Model (KE Model) and the Jackson–Green Model (JG Model). The present work extends the previous KE and JG works, accounting for the effect of realistic material behavior in terms of the varying yield strengths and the isotropic strain hardening behavior. The predicted results show that the transition behavior of the materials from the elastic–plastic to the fully plastic case is influenced by the yield strength and the tangent modulus (E_t) and such transition do not take place at specific values of interference ratios as suggested by the KE model. New empirical relations are proposed to determine the contact load and the contact area based on the analysis. Numerical results from the finite element modeling are also validated with an experimental ball on flat configuration approach.