Materially and geometrically non-linear woven composite micro-mechanical model with failure for finite element simulations

A computational micro-mechanical material model of woven fabric composite material is developed to simulate failure. The material model is based on repeated unit cell approach. The fiber reorientation is accounted for in the effective stiffness calculation. Material non-linearity due to the shear stresses in the impregnated yarns and the matrix material is included in the model. Micro-mechanical failure criteria determine the stiffness degradation for the constituent materials. The developed material model with failure is programmed as user-defined sub-routine in the LS-DYNA finite element code with explicit time integration. The code is used to simulate the failure behavior of woven composite structures. The results of finite element simulations are compared with available test results. The model shows good agreement with the experimental results and good computational efficiency required for finite element simulations of woven composite structures.     

复合材料层合板壳非线性力学的研究进展

复合材料层合板壳是由多种组分材料组合而成.与单一材料的板壳结构相比,它无明确的材料主方向,各层间材料间断和不连续,具有明显的几何非线性和材料非线性等新的特点.其失效模式也远比单一材料的情况复杂,具有如基体开裂、脱胶、分层、分层裂纹偏转、多分层以及分层传播等多种模式.各国学者基于不同的考虑,提出了多种方法研究复合材料层合板壳的失效.首先,在简要介绍了层合板壳线性力学基本理论的基础上,重点回顾了层合板壳结构非线性力学几种基本理论发展的过程,主要阐述了经典大挠度非线性理论、一阶剪切变形理论、高阶剪切变形理论、锯齿理论、广义分层理论的理论体系及基本公式,并对几种理论之间的联系和差异进行了总结;其次,介绍了当前层合结构非线性领域的研究进展,综述了典型复合材料板壳结构的失效机理及优化设计、复合材料板壳结构在复杂环境下的破坏机理、复合材料板壳结构的物理非线性、含脱层纤维增强复合材料板壳结构的破坏机理等各研究热点的最新研究成果;最后,对该领域未来的研究方向进行了展望.     

参数拟二次规划法及其在复合材料失效过程分析中的应用

本文介绍关于不连续线弹性静力学问题的多变量变分原理与相应的参数拟二次规划算法在复合材料层压板强度分析及纤维增强脆性基体材料单元细观失效过程分析中的应用。     

Modeling the effects of material non-linearity using moving window micromechanics

In the analysis of materials with random heterogeneous microstructure the assumption is often made that material behavior can be represented by homogenized or effective properties. While this assumption yields accurate results for the bulk behavior of composite materials, it ignores the effects of the random microstructure. The spatial variations in these microstructures can focus, initiate and propagate localized non-linear behavior, subsequent damage and failure. In previous work a computational method, moving window micromechanics (MW), was used to capture microstructural detail and characterize the variability of the local and global elastic response. Digital images of material microstructure described the microstructure and a local micromechanical analysis was used to generate spatially varying material property fields. The strengths of this approach are that the material property fields can be consistently developed from digital images of real microstructures, they are easy to import into finite element models (FE) using regular grids, and their statistical characterizations can provide the basis for simulations further characterizing stochastic response. In this work, the moving window micromechanics technique was used to generate material property fields characterizing the non-linear behavior of random materials under plastic yielding; specifically yield stress and hardening slope, post yield. The complete set of material property fields were input into FE models of uniaxial loading. Global stress strain curves from the FE–MW model were compared to a more traditional micromechanics model, the generalized method of cells. Local plastic strain and local stress fields were produced which correlate well to the microstructure. The FE–MW method qualitatively captures the inelastic behavior, based on a non-linear flow rule, of the sample continuous fiber composites in transverse uniaxial loading.     

Non-linear analysis of functionally graded fiber reinforced composite laminated plates,Part I: Theory and solutions

This paper deals with the large amplitude vibration, non-linear bending and postbuckling of fiber reinforced composite laminated plates resting on an elastic foundation in hygrothermal environments. Two kinds of fiber reinforced laminated plates, namely, uniformly distributed and functionally graded reinforcements, are considered. The material properties of fiber reinforced laminated plates are estimated through a micromechanical model and are assumed to be temperature-dependent and moisture-dependent. The motion equations are based on a higher order shear deformation plate theory that includes plate-foundation interaction and the hygrothermal effect. A two-step perturbation technique is employed to determine the non-linear to linear frequency ratios of plate vibration, the load-deflection and load-bending moment curves of plate bending, and postbuckling equilibrium paths of laminated plates.     

材料非线性对复合材料层合板热自由边界效应的影响

本文用准三维有限元法研究了材料非线性对复合材料层合板热自山边界效应的影响,给出了修正型Hahn-Tsai非线性应力-应变关系的三维形式。由本文非线性分析方法得到的层间应力与以往由线性分析方法得到的层间应力做了比较,结果表明:材料非线性能显著降低层间剪应力的集中程度,但对层间正应力影响不太明显。     

Non-linear strain rate dependent micro-mechanical composite material model for finite element impact and crashworthiness simulation

The present study aims at implementation of a strain rate dependent, non-linear, micro-mechanics material model for laminated, unidirectional polymer matrix composites into the explicit finite element code LSDYNA. The objective is to develop an accurate and simple micro-mechanical, rate dependent material model, which is computationally efficient. Within the model a representative volume cell is assumed. The stress-strain relation including rate dependent effects for the micro-model is derived for both shell elements and solid elements. Micro-failure criterion is presented for each material constituent and failure mode. The implemented model can deal with problems such as impact, crashworthiness, and failure analysis under quasi-static loads. The developed material model has a wide range of applications such as jet engine jackets, armor plates, and structural crashworthiness simulation. The deformation response of two representative composite materials with varying fiber orientation is presented using the described technique. The predicted results compare favorably to experimental values.     

含椭圆孔非对称层合板的Green函数

运用作者建立的复合材料层合板问题的求解方法，得到了含椭圆孔和裂纹的无限大非对称层合板在广义集中载荷作用下的Ｇｒｅｅｎ函数．利用这些结果可研究含孔或裂纹的非对称层合板在集中载荷作用下的力学行为，还可结合边界元方法对复杂层合板结构建立有效的数值分析方法     

Geometrically non-linear analysis of laminated composite structures using a 4-node co-rotational shell element with enhanced strains

To demonstrate the solutions of linear and geometrically non-linear analysis of laminated composite plates and shells, the co-rotational non-linear formulation of the shell element is presented. The combinations of an enhanced assumed strain (EAS) in the membrane strains and assumed natural strains (ANS) in the shear strains improve the behavior of 4-node shell element. To secure computational efficiency in the incremental non-linear analysis, the present element uses the form of the resultant forces pre-integrated through the thickness. The transverse shear stiffness of the laminates is defined by an equilibrium approach instead of the shear correction factor. Numerical examples of this study show very good agreement with the references.     

A semi-analytical approach for the non-linear large deflection analysis of laminated rectangular plates under general out-of-plane loading

A semi-analytical approach for the geometrically non-linear analysis of rectangular laminated plates with general inplane and out-of-plane boundary conditions under a general distribution of out-of-plane loads is developed. The analysis is based on the elastic thin plate theory with geometrically non-linear von Kármán strains. The solution of the non-linear partial differential equations is reduced to an iterative sequential solution of non-linear ordinary differential equations using the multi-term extended Kantorovich method. The efficiency, accuracy, and convergence of the proposed method are examined through a comparison with other semi-analytical methods and with finite element analyses. The capabilities of the approach and its applicability to the non-linear large deflection analysis of plate structures are demonstrated through various numerical examples. Emphasis is placed on combinations of lamination, boundary, and loading conditions that cannot be analyzed using alternative semi-analytical methods.     

Fatigue Damage in Composite Materials

The phenomenon of fatigue is critical for designing structures including elements made of composite materials. The accurate prediction of the life and fatigue resistance of laminated composites is one of the subjects of inquiry in materials science. The ability of predicting the life of laminates is important for designing, operation, and safety analysis of a composite structure under specific conditions. To predict reliably the life of structures, it is necessary to know the mechanisms of cyclic deformation and damage. It is also necessary to develop a qualitative theory of fatigue failure that should be based on the concepts of solids mechanics. Developing such a theory requires to evaluate the microscopic parameters and the macroscopic variables of the material at the level of a laminate and the structure and to determine exactly the load modes acting on the system.     

Dynamic stability of uncertain laminated beams subjected to subtangential loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. Thus the understanding of the effect of uncertainties in laminated structures on their static and dynamic responses is highly important for a reliable design of such structures. In this research, we focus on the probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: exact Monte Carlo simulation, sensitivity-based Monte Carlo simulation, and probabilistic FEA. These methods were integrated into the existing finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems to study the stability analysis using the dynamic criterion.     

Non-linear free periodic vibrations of variable stiffness composite laminated plates

Steady-state free vibrations, with large amplitude displacements, of variable stiffness composite laminated plates (VSCL) are analysed. The intentions of this research are: (1)?to find out how the natural frequencies and (mode) shapes evolve with the displacement amplitude in this new type of laminated composite material; (2)?to describe modal interactions in VSCL due to energy interchanges under the coupling induced by non-linearity; (3)?to compare the VSCL with traditional, constant stiffness, laminated plates. The VSCL of interest here have curvilinear fibres and the numerical analysis carried out is based on a recently developed p-version finite element with hierarchic basis functions. The element follows first-order shear deformation theory and considers Von Kármán??s non-linear terms. The time domain equations of motion are first reduced using the linear modes of vibration and then transformed to the frequency domain via the harmonic balance method. These frequency domain equations are solved by an arc-length continuation method.     

基于新修正偶应力理论的Mindlin层合板稳定性分析

基于新的各向异性修正偶应力理论提出一个Mindlin复合材料层合板稳定性模型。该理论包含纤维和基体两个不同的材料长度尺度参数。不同于忽略横向剪切应力的修正偶应力Kirchhoff薄板理论,Mindlin层合板考虑横向剪切变形引入两个转角变量。进一步建立了只含一个材料细观参数的偶应力Mindlin层合板工程理论的稳定性模型。计算了正交铺设简支方板Mindlin层合板的临界载荷。计算结果表明该模型可以用于分析细观尺度层合板稳定性的尺寸效应。     

功能梯度材料构件三维分析的细观元模型

提出一种新颖的功能梯度构件分析的细观元法,给出了方法模型、基本算式及特点与功能。细观元法对构件的常规有限单元内部设置密集细观单元以反映材料特性梯度变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算。这种细观元法既能充分反映材料功能梯度及组分变化特性,而其计算单元与自由度又与常规有限元一样,是一种针对功能梯度构件分析的有效数值方法。算例表明了细观元法对不同情况下功能梯度构件分析的适应性与精度。     

Deformation of a laminated composite with a physically nonlinear reinforcement and microdamageable matrix

The structural theory of short-term microdamage is generalized to a laminated composite with a microdamageable matrix and physically nonlinear reinforcement. The basis for the generalization is the stochastic elasticity equations of a laminated composite with a porous matrix. Microvolumes in the matrix material meet the Huber-Mises failure criterion. The damaged-microvolume balance equation for the matrix is derived. This equation and the equations relating macrostresses and macrostrains of a laminated composite with porous matrix and physically nonlinear reinforcement constitute a closed-form system of equations. This system describes the coupled processes of physically nonlinear deformation and microdamage occurring in different composite components. Algorithms for computing the microdamage-macrostrain relationships and deformation diagrams are developed. Uniaxial tension curves are plotted for a laminated composite with linearly hardening reinforcement __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 47–56, November 2005.     

Short-term microdamage of laminated material with nonlinear matrix and microdamaged reinforcement

A structural theory of short-term microdamage is proposed for a two-component laminated composite with microdamageable reinforcement and physically nonlinear matrix. The basis of the theory is the stochastic elasticity equations of a laminated composite with a porous reinforcement. Microvolumes in the reinforcement material meet the Huber-Mises failure criterion. The damaged-microvolume balance equation for the reinforcement is derived. This equation and the equations relating macrostresses and macrostrains of a laminated composite with porous reinforcement and physically nonlinear matrix constitute a closed-form system of equations. This system describes the coupled processes of physically nonlinear deformation and microdamage occurring in different composite components. Algorithms for computing the microdamage-macrostrain relationships and deformation diagrams are developed. Uniaxial tension curves are plotted for a laminated composite with linearly hardening matrix __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 3–12, December 2005.     

Large rotations in first-order shear deformation FE analysis of laminated shells

The paper deals with the geometrically non-linear analysis of laminated composite beams, plates and shells in the framework of the first-order transverse shear deformation (FOSD) theory. A central point of the present paper is the discussion of the relevance of five- and six-parameter variants, respectively, of the FOSD hypothesis for large rotation plate and shell problems. In particular, it is shown that the assumption of constant through-thickness distribution of the transverse normal displacements is acceptable only for small and moderate rotation problems. Implications inherent in this assumption that are incompatible with large rotations are discussed from the point of view of the transverse normal strain-displacement relations as well as in the light of an enhanced, accurate large rotation formulation based on the use of Euler angles. The latter one is implemented as an updating process within a Total Lagrangian formulation of the six-parameter FOSD large rotation plate and shell theory. Numerical solutions are obtained by using isoparametric eight-node Serendipity-type shell finite elements with reduced integration. The Riks-Wempner-Ramm arc-length control method is used to trace primary and secondary equilibrium paths in the pre- and post-buckling range of deformation. A number of sample problems of non-linear, large rotation response of composite laminated plate and shell structures are presented including symmetric and asymmetric snap-through and snap-back problems.     

考虑面内破坏过程的含分层损伤复合材料加筋层合板的动力响应

基于复合材料层合板一阶剪切理论，推导了复合材料层合板单元的刚度阵和质量阵列式；同时采用了Adams应交能法与Rayleigh阻尼模型相结合的方法，构造了相应的阻尼阵列式；为了防止在低阶模态中分层处出现的上、下子板不合理的嵌入现象，建立了含分层损伤复合材料加筋层合板动力分析的分层分析模型和虚拟界面联接模型。在上述模型和理论基础上，采用了Tsai提出的刚度退化准则和动力响应分析的精细积分法，对含分层损伤复合材料加筋层合板结构进行了动力响应和破坏分析。通过算例，分别讨论了外载频率、分层位置，以及破坏过程的刚度退化对含损伤复合材料加筋层合板动力响应特征的影响，得到了一些具有理论和工程价值的结论。     

Postbuckling of imperfect rectangular composite plates under inplane shear closed-form approximate solutions

In this paper, the postbuckling behavior of rectangular orthotropic laminated composite plates with initial imperfections under inplane shear load is investigated in a closed-form analytical manner. The plates under consideration are assumed to be infinitely long in the longitudinal direction. At the longitudinal edges, two different sets of boundary conditions are considered, specifically 1) simply supported edges and 2) fully clamped edges. Using Timoshenko-type shape functions for both the initial bifurcational buckling analysis and the subsequent Marguerre-type postbuckling studies, closed-form analytical solutions for the buckling loads and for the postbuckling state variables are derived. A comparison with geometrically non-linear finite element computations shows that the derived analysis approaches are suitable for postbuckling studies in load ranges not too far beyond bifurcational buckling as they are currently relevant for e.g., composite airframe structural analysis and design. Due to their strictly closed-form analytical nature, the presented analysis methods can be used conveniently in engineering practice for all application purposes where computational time is a crucial factor, especially for preliminary analysis and design or optimization procedures.