Homogenization of inelastic solid materials at finite strains based on incremental minimization principles. Application to the texture analysis of polycrystals

The paper presents new continuous and discrete variational formulations for the homogenization analysis of inelastic solid materials undergoing finite strains. The point of departure is a general internal variable formulation that determines the inelastic response of the constituents of a typical micro-structure as a generalized standard medium in terms of an energy storage and a dissipation function. Consistent with this type of finite inelasticity we develop a new incremental variational formulation of the local constitutive response, where a quasi-hyperelastic micro-stress potential is obtained from a local minimization problem with respect to the internal variables. It is shown that this local minimization problem determines the internal state of the material for finite increments of time. We specify the local variational formulation for a distinct setting of multi-surface inelasticity and develop a numerical solution technique based on a time discretization of the internal variables. The existence of the quasi-hyperelastic stress potential allows the extension of homogenization approaches of finite elasticity to the incremental setting of finite inelasticity. Focussing on macro-deformation-driven micro-structures, we develop a new incremental variational formulation of the global homogenization problem for generalized standard materials at finite strains, where a quasi-hyperelastic macro-stress potential is obtained from a global minimization problem with respect to the fine-scale displacement fluctuation field. It is shown that this global minimization problem determines the state of the micro-structure for finite increments of time. We consider three different settings of the global variational problem for prescribed displacements, non-trivial periodic displacements and prescribed stresses on the boundary of the micro-structure and develop numerical solution methods based on a spatial discretization of the fine-scale displacement fluctuation field. Representative applications of the proposed minimization principles are demonstrated for a constitutive model of crystal plasticity and the homogenization problem of texture analysis in polycrystalline aggregates.     

Dynamic analysis of beam-soil interaction systems with material and geometrical nonlinearities

In this paper a Hybrid Domain Boundary Element Method is developed for the geometrically nonlinear dynamic analysis of inelastic Euler-Bernoulli beams of arbitrary doubly symmetric simply or multiply connected constant cross-section, resting on viscous inelastic Winkler foundation. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse dynamic loading and bending moments in both directions as well as to axial loading, while its edges are subjected to the most general boundary conditions. A displacement based formulation is developed and inelastic redistribution is modelled through a distributed plasticity (fibre) approach. A uniaxial hysteretic law is considered for the evolution of the plastic part of the normal stress following the phenomenological hysteresis model, while hysteretic force-displacement model is also employed in order to describe the inelastic behaviour of the Winkler springs. Numerical integration over the beam cross sections is performed in order to resolve the hysteric parts of the stress resultants. Application of the boundary element technique yields a system of nonlinear Differential-Algebraic Equations, which are written in state-space form and solved by an incremental–iterative solution strategy. Numerical examples are worked out confirming the accuracy and the computational efficiency of the proposed beam formulation, as well as the significant influence of material and geometrical nonlinearities in the response of beam-soil interaction systems.     

基于介观力学信息的颗粒材料损伤——愈合与塑性宏观表征

本文在二阶计算均匀化框架下提出了颗粒材料损伤--愈合与塑性的多尺度表征方法. 颗粒材料结构在宏观尺度模型化为梯度Cosserat连续体,在其有限元网格的每个积分点处定义具有离散颗粒介观结构的表征元. 建立了表征元离散颗粒系统的非线性增量本构关系. 表征元周边介质作用于表征元边界颗粒的增量力与增量力偶矩以表征元边界颗粒的增量线位移与增量转动角位移、当前变形状态下表征元离散介观结构弹性刚度、以及凝聚到表征元边界颗粒的增量耗散摩擦力表示. 基于平均场理论与Hill定理,导出了基于介观力学信息的梯度Cosserat连续体增量非线性本构关系. 在等温热动力学框架下定义了表征颗粒材料各向异性损伤--愈合和塑性的损伤、愈合张量因子与综合损伤、愈合效应的净损伤张量因子和塑性应变. 此外,定义了损伤和塑性耗散能密度与愈合能密度,以定量比较材料损伤、愈合、塑性对材料失效的效应. 应变局部化数值例题结果显示了所建议的颗粒材料损伤--愈合--塑性表征方法的有效性.      

Dynamic analysis of laminated cross-ply composite non-circular thick cylindrical shells using higher-order theory

Here, the dynamic analysis of laminated cross-ply composite non-circular thick cylindrical shells subjected to thermal/mechanical load is carried out based on higher-order theory. The formulation accounts for the variation of the in-plane and transverse displacements through the thickness, abrupt discontinuity in slope of the in-plane displacements at the interfaces, and includes in-plane, rotary inertia terms, and also the inertia contributions due to the coupling between the different order displacement terms. The strain–displacement relations are accurately accounted for in the formulation. The shell responses are obtained employing finite element approach in conjunction with direct time integration technique. A detailed parametric study is carried out to bring out the effects of length and thickness ratios, eccentricity parameters and number of layers on the thermal/mechanical response characteristics of non-circular shells.     

Fiber push-out study of a copper matrix composite with an engineered interface: Experiments and cohesive element simulation

The fiber push-out test is a basic method to probe the mechanical properties of the fiber/matrix interface of fiber-reinforced metal matrix composites. In order to estimate the interfacial properties, parameters should be calibrated using the measured load–displacement data and theoretical models. In the case of a soft matrix composite, the possible plastic yield of the matrix has to be considered for the calibration. Since the conventional shear lag models are based on elastic behavior, a detailed assessment of the plastic effect is needed for accurate calibration. In this paper, experimental and simulation studies are presented regarding the effect of matrix plasticity on the push-out behavior of a copper matrix composite with strong interface bonding. Microscopic images exhibited significant local plastic deformation near the fibers leading to salient nonlinear response in the global load–displacement curve. For comparison, uncoated interface with no chemical bonding was also examined where the nonlinearity was not observed. A progressive FEM simulation was conducted for a complete push-out process using the cohesive zone model and inverse fitting. Excellent coincidence was achieved with the measured push-out curve. The predicted results confirmed the experimental observations.     

Extended finite element method in plasticity forming of powder compaction with contact friction

In this paper, a new computational technique is presented based on the eXtended Finite Element Method (X-FEM) in pressure-sensitive plasticity of powder compaction considering frictional contact. In X-FEM, the need for mesh adaption to discontinuity interface is neglected and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of the elements. The technique is applied by employing additional functions, which are added to approximate the displacement field of the elements located on the interface. The double-surface cap plasticity model is employed within the X-FEM framework in numerical simulation of powder material. The plasticity model includes a failure surface and an elliptical cap, which closes the open space between the failure surface and hydrostatic axis. The moving cap expands in the stress space according to a specified hardening rule. The frictional behavior of contact between two bodies is modeled by using the X-FEM technique and applying the Heaviside enrichment function. The application of X-FEM technique in simulation of pressure-sensitive material is presented in an incremental manner and the role of sub-elements in simulation of contact treatment is demonstrated. Finally, several numerical examples are analyzed with special reference to plasticity forming of powder compaction.     

Approximate MPA-based method for performing incremental dynamic analysis

In performance based earthquake engineering, it is important to accurately predict the seismic demand and capacities of structures. One recent estimation method is incremental dynamic analysis (IDA), which requires a series of nonlinear response history analyses (RHA) of the structure under various ground motions, each scaled to multiple levels of intensity, selected to cover entire range of structural response from elasticity, to yield and finally global dynamic instability. The implementation of IDA requires intensive computation and detailed knowledge of the nonlinear RHA of structures. In response to the complexity of IDA, an approximate method based on modal pushover analysis (MPA-based IDA) was developed. In MPA-based IDA, seismic demands are computed using the nonlinear RHA of the equivalent SDF systems instead of using nonlinear RHA of MDF systems. The objective of this study is to develop a simpler MPA-based IDA procedure that can avoid nonlinear RHA of equivalent SDF systems. For this purpose, MPA-based IDA employs the empirical equation of the inelastic displacement ratio (C _R ), defined as the ratio of peak displacement of the inelastic SDF system to that of the corresponding elastic SDF system given the strength ratio R, and that of the collapse strength ratio (R _c), which is the ratio of collapse intensity to yield strength. The proposed procedure is verified by comparing the seismic demands and capacities of 6-, 9-, and 20-story steel moment frames as determined by the proposed method and exact IDA.     

Direct Extraction of Cohesive Fracture Properties from Digital Image Correlation: A Hybrid Inverse Technique

The accuracy of an adopted cohesive zone model (CZM) can affect the simulated fracture response significantly. The CZM has been usually obtained using global experimental response, e.g., load versus either crack opening displacement or load-line displacement. Apparently, deduction of a local material property from a global response does not provide full confidence of the adopted model. The difficulties are: (1) fundamentally, stress cannot be measured directly and the cohesive stress distribution is non-uniform; (2) accurate measurement of the full crack profile (crack opening displacement at every point) is experimentally difficult to obtain. An attractive feature of digital image correlation (DIC) is that it allows relatively accurate measurement of the whole displacement field on a flat surface. It has been utilized to measure the mode I traction-separation relation. A hybrid inverse method based on combined use of DIC and finite element method is used in this study to compute the cohesive properties of a ductile adhesive, Devcon Plastic Welder II, and a quasi-brittle plastic, G-10/FR4 Garolite. Fracture tests were conducted on single edge-notched beam specimens (SENB) under four-point bending. A full-field DIC algorithm was employed to compute the smooth and continuous displacement field, which is then used as input to a finite element model for inverse analysis through an optimization procedure. The unknown CZM is constructed using a flexible B-spline without any “a priori” assumption on the shape. The inversely computed CZMs for both materials yield consistent results. Finally, the computed CZMs are verified through fracture simulation, which shows good experimental agreement.     

A generalized approach to formulate the consistent tangent stiffness in plasticity with application to the GLD porous material model

It has been shown that the use of the consistent tangent moduli is crucial for preserving the quadratic convergence rate of the global Newton iterations in the solution of the incremental problem. In this paper, we present a general method to formulate the consistent tangent stiffness for plasticity. The robustness and efficiency of the proposed approach are examined by applying it to the isotropic material with J₂ flow plasticity and comparing the performance and the analysis results with the original implementation in the commercial finite element program ABAQUS. The proposed approach is then applied to an anisotropic porous plasticity model, the Gologanu–Leblond–Devaux model. Performance comparison between the consistent tangent stiffness and the conventional continuum tangent stiffness demonstrates significant improvement in convergence characteristics of the overall Newton iterations caused by using the consistent tangent matrix.     

Micromechanical model for cohesive materials based upon pseudo-granular structure

A micromechanical model for cohesive materials is derived by considering their underlying microstructure conceptualized as a collection of grains interacting through pseudo-bonds. The pseudo-bond or the inter-granular force–displacement relations are formulated taking inspiration from the atomistic-level particle interactions. These force–displacement relationships are then used to derive the incremental stiffnesses at the grain-scale, and consequently, obtain the sample-scale stress–strain relationship of a representative volume of the material. The derived relationship is utilized to study the stress–strain and failure behavior including the volume change and “brittle” to “ductile” transition behavior of cohesive materials under multi-axial loading condition. The model calculations are compared with available measured data for model validation. Model predictions exhibit both quantitative and qualitative consistency with the observed behavior of cohesive material.     

平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体离散元法的增量平衡方程；而后推导了多刚体离散元的弹塑性弹簧系数矩阵，建立了多刚体离散元内力屈服面塑性铰法的增量求解格式，成功地进行了平面钢框架的弹塑性、大位移极限承载力分析。计算结果与其他数值方法或实验结果吻合良好，显示了多刚体离散元方法进行结构极限承载力分析这一复杂问题的优越性     

Moderately large flexural vibrations of composite plates with thick layers

In this paper moderately large amplitude vibrations of a polygonally shaped composite plate with thick layers are analyzed. Three homogeneous and isotropic layers with a common Poisson’s ratio are perfectly bonded and their arbitrary thickness and material properties are symmetrically disposed about the middle plane. Mindlin–Reissner kinematic assumptions are implemented layerwise, and as such model both the global and local response. Geometric nonlinear effects arising from longitudinally constrained supports are taken into account by Berger’s approximation of nonlinear strain–displacement relations. Overall cross-sectional rotations are defined and subsequently a correspondence of this complex problem to the simpler case of a homogenized shear-deformable nonlinear plate with effective stiffness and hard hinged boundary conditions is found. The nonlinear steady-state response of composite plates subjected to a time-harmonic lateral excitation is investigated and the phenomena of nonlinear resonance are studied and evaluated.     

Analysis of material instabilities in inelastic solids by incremental energy minimization and relaxation methods: evolving deformation microstructures in finite plasticity

We propose an approach to the definition and analysis of material instabilities in rate-independent standard dissipative solids at finite strains based on finite-step-sized incremental energy minimization principles. The point of departure is a recently developed constitutive minimization principle for standard dissipative materials that optimizes a generalized incremental work function with respect to the internal variables. In an incremental setting at finite time steps this variational problem defines a quasi-hyperelastic stress potential. The existence of this potential allows to be recast a typical incremental boundary-value problem of quasi-static inelasticity into a principle of minimum incremental energy for standard dissipative solids. Mathematical existence theorems for sufficiently regular minimizers then induce a definition of the material stability of the inelastic material response in terms of the sequentially weakly lower semicontinuity of the incremental variational functional. As a consequence, the incremental material stability of standard dissipative solids may be defined in terms of the quasi-convexity or the rank-one convexity of the incremental stress potential. This global definition includes the classical local Hadamard condition but is more general. Furthermore, the variational setting opens up the possibility to analyze the post-critical development of deformation microstructures in non-stable inelastic materials based on energy relaxation methods. We outline minimization principles of quasi- and rank-one convexifications of incremental non-convex stress potentials for standard dissipative solids. The general concepts are applied to the analysis of evolving deformation microstructures in single-slip plasticity. For this canonical model problem, we outline details of the constitutive variational formulation and develop numerical and semi-analytical solution methods for a first-level rank-one convexification. A set of representative numerical investigations analyze the development of deformation microstructures in the form of rank-one laminates in single slip plasticity for homogeneous macro-deformation modes as well as inhomogeneous macroscopic boundary-value problems. The well-posedness of the relaxed variational formulation is indicated by an independence of typical finite element solutions on the mesh-size.     

Development of the large increment method for elastic perfectly plastic analysis of plane frame structures under monotonic loading

The displacement-based finite element method dominates current practice for material nonlinear analysis of structures. However, there are several characteristics that may limit the effectiveness of this approach. In particular, for elastoplastic analysis, the displacement method relies upon a step-by-step incremental approach stemming from flow theory and also requires significant mesh refinement to resolve behavior in plastic zones. This leads to computational inefficiencies that, in turn, encourage the reconsideration of force-based approaches for elastoplastic problems.One of these force algorithms that has been recently developed is the large increment method. The main advantage of the flexibility-based large increment method (LIM) over the displacement method is that it separates the global equilibrium and compatibility equations from the local constitutive relations. Consequently, LIM can reach the solution in one large increment or in a few large steps, thus, avoiding the development of cumulative errors. This paper discusses the extension of the large increment methodology for the nonlinear analysis of plane frame structures controlled by an elastic, perfectly plastic material model. The discussion focuses on the power of LIM to handle these nonlinear problems, especially when plastic hinges form in the frame and ultimately as the structure approaches the collapse stage. Illustrative planar frame examples are presented and the results are compared with those obtained from a standard displacement method.     

A nonlinear theory for doubly curved anisotropic sandwich shells with transversely compressible core

In the present paper, an advanced geometrically nonlinear shell theory of doubly curved structural sandwich panels with transversely compressible core is presented. The model is based on the adoption of the Kirchhoff theory for the face sheets and a second/third order power series expansion for the core displacements. The theory accounts for dynamic effects as well as for initial geometric imperfections. In the v. Kármán sense, large displacement theory is employed with respect to the transverse direction while the displacement gradients with respect to the tangential directions are assumed to be small. The equations of motion are derived by means of Hamilton’s principle and hold valid for all types of elastic and elastic–plastic material models. The theory is illustrated by an analysis of the elastic buckling and postbuckling behavior of flat and curved sandwich panels using an extended Galerkin scheme. Owing to the assumed transverse flexibility of the core, both the global and the local (face wrinkling) instability modes can be addressed.     

A time-domain high-order spectral finite element for the simulation of symmetric and anti-symmetric guided waves in laminated composite strips

A time domain spectral finite element is developed for improving the efficiency of numerical simulations of guided waves in laminated composite strips. The finite element relies on a new generalized laminate mechanics model formulated to represent symmetric and anti-symmetric Lamb waves. The laminate mechanics incorporate third-order polynomial terms for the approximation of axial and transverse displacement fields through the thickness and consider the displacements of the upper and lower surfaces as degrees of freedom. The laminate theory formulation is easily expanded to a high-order layerwise model. Based on the resultant governing equations of the laminate section, a new finite element with 8 nodal degrees of freedom is formulated; its nodes are collocated with Gauss–Lobatto–Legendre integration points in order to improve computational efficiency. Stiffness and mass matrices are assembled and the transient response is predicted using the explicit central differences time integration scheme. The transient response of Aluminum, Carbon Fiber Reinforced Polymer laminated and sandwich strips is investigated. Numerical results are validated against a semi-analytical solution. The accuracy and computational efficiency of the introduced element regarding the prediction of symmetric and anti-symmetric wave propagation is also quantified.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

Nonlinear analysis of beams of variable cross section,including shear deformation effect

In this paper the analog equation method (AEM), a BEM-based method, is employed for the nonlinear analysis of a Timoshenko beam with simply or multiply connected variable cross section undergoing large deflections under general boundary conditions. The beam is subjected in an arbitrarily concentrated or distributed variable axial loading, while the shear loading is applied at the shear center of the cross section, avoiding in this way the induction of a twisting moment. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, the axial displacement and to two stress functions and solved using the AEM. Application of the boundary element technique yields a system of nonlinear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. Numerical examples with great practical interest are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. The influence of the shear deformation effect is remarkable.     

Nonlinear analysis of shear deformable beam-columns partially supported on tensionless three-parameter foundation

In this paper, a boundary element method is developed for the nonlinear analysis of shear deformable beam-columns of arbitrary doubly symmetric simply or multiply connected constant cross-section, partially supported on tensionless three-parameter foundation, undergoing moderate large deflections under general boundary conditions. The beam-column is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM-based method. Application of the boundary element technique yields a system of nonlinear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member as well as the shear forces along the span induced by the applied axial loading. Numerical examples are worked out to illustrate the efficiency, wherever possible, the accuracy and the range of applications of the developed method.