STAGS computational procedure for progressive failure analysis of laminated composite structures

Non-linear analyses are difficult simulations to perform and place increased demands not only on the computational systems but also on the analyst. The number of possible problems and/or difficulties increases significantly compared to those for linear elastic analyses. Combined material and geometric non-linearities challenge the analyst and the solution algorithms. Progressive failure and damage propagation for composite structures result in even more complexities due to the discrete, abrupt changes in local material stiffness. Analysts need to interrogate the computed solutions carefully based on their understanding of structural mechanics, material behavior, computational procedures and non-linear phenomena to distill correct physics from such simulations. This paper describes the computational strategy incorporated into the STAGS non-linear finite element analysis code with special emphasis on progressive failure analysis of laminated composite structures. Results for selected laminated composite structures are used to demonstrate the PFA capability.     

含孔平面编织混杂铺层层合板压缩破坏仿真

对ANSYS软件平台进行了二次开发，建立了平面编织混杂铺层层合板损伤累积有限元模 型，对含孔平面编织混杂铺层层合板的压缩破坏行为进行了仿真计算. 研究结果表明，所建 模型能够较好地仿真含孔平面编织层合板的损伤破坏过程，并能准确预测层合板的 损伤类型、损伤扩展以及破坏强度；同时提出并验证了适合该类型层合板的损伤判 据和衰减准则. 该模型所给结果形象直观，适合工程应用.     

水平荷载下型钢再生混凝土柱-钢梁组合框架受力性能非线性数值分析

在型钢再生混凝土柱-钢梁组合框架拟静力试验的基础上,采用Abaqus软件建立该组合框架有限元模型并进行非线性数值分析,获取组合框架的变形图、应力云图及荷载-位移骨架曲线,分析组合框架的受力破坏特征,验证了有限元模型的合理性,并对组合框架进行了参数分析。结果表明,数值模拟计算值和试验值对比误差较小,数值计算模型能够较好地模拟该组合框架的受力性能;组合框架符合强柱弱梁的破坏机制;另外,提高型钢强度或再生混凝土强度对组合框架承载力和刚度有利,但对其变形能力不利;框架承载力和刚度随着梁柱线刚度的增加而提高;增大轴压比对于组合框架的抗震性能和延性不利。研究结论可为该类绿色组合框架的工程应用提供参考。     

Woven fabric composite material model with material nonlinearity for nonlinear finite element simulation

The objective of the current investigation is to develop a simple, yet generalized, model which considers the two-dimensional extent of woven fabric, and to have an interface with nonlinear finite element codes. A micromechanical composite material model for woven fabric with nonlinear stress-strain relations is developed and implemented in ABAQUS for nonlinear finite element structural analysis. Within the model a representative volume cell is assumed. Using the iso-stress and iso-strain assumptions the constitutive equations are averaged along the thickness direction. The cell is then divided into many subcells and an averaging is performed again by assuming uniform stress distribution in each subcell to obtain the effective stress–strain relations of the subcell. The stresses and strains within the subcells are combined to yield the effective stresses and strains in the representative cell. Then this information is passed to the finite element code at each material point of the shell element. In this manner structural analysis of woven composites can be performed. Also, at each load increment global stresses and strains are communicated to the representative cell and subsequently distributed to each subcell. Once stresses and strains are associated to a subcell they can be distributed to each constituent of the subcell i.e. fill, warp, and resin. Consequently micro-failure criteria (MFC) can be defined for each constituent of a subcell and the proper stiffness degradation can be modeled if desired. This material model is suitable for implicit and could be modified for explicit finite element codes to deal with problems such as crashworthiness, impact, and failure analysis under static loads.     

Multi-level modeling of woven glass/epoxy composite for multilayer printed circuit board applications

The objective of this paper is to develop a hybrid homogenization method to predict the elastic properties of a common woven glass/epoxy composite substrate for multilayer circuit board applications. Comprehensive high resolution 3D finite element (FE) models of a quarter of the repeated unit cell (RUC) for the woven glass/epoxy composite were developed based on different micromechanical schemes. . Specifically, four different micromechanics schemes were investigated: self-consistent, Mori–Tanaka, three-phase approach and composite cylinder assemblage (CCA). The element based strain concentration matrices were determined and used to obtain the homogenized woven glass/epoxy composite properties via a specially developed MATLAB code. Attention was further devoted to the predictions of the homogenized elastic moduli of the multilayer printed circuit board (PCB). The results from our simulations, based on Mori–Tanaka and CCA, are in good agreement with existing experimental results, indicating that the newly proposed homogenization scheme can be used as a design tool to predict the overall properties of woven composite materials typically used in multilayer PCB applications.     

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.     

基于响应面法的复合材料连接结构可靠性分析

复合材料的力学性能参数比金属材料呈现出更大的分散性。本文提出采用加权线性响应面法对复合材料单钉连接结构进行可靠性分析。采用三维有限元模型进行复合材料单钉连接的确定性渐进失效分析。单向铺层的力学性能参数和强度参数均考虑为输入随机变量,用于可靠性分析。详细研究了不同失效准则和不同刚度退化方法的影响。结果表明,Olmedo-Santiuste失效准则与Camanho-Matthews刚度退化方法的组合能够给出与试验均值最接近的强度分析结果。另外,还仔细研究了响应面参数f的影响。结果表明,取f=1的线性加权响应面能够为复合材料单钉连接可靠性分析提供较好的统计特性逼近值。     

复合材料弹塑性多尺度分析模型与算法

对材料非线性多尺度分析的计算模型与算法进行研究.在构建周期分布单胞分析算法的基础上,发展针对复合材料结构材料非线性多尺度分析的一般有限元方法.方法的特点是将所建立的单胞分析过程作为有限元分析的子程序嵌入到总程序系统当中,完成对应的高斯点应力计算,因而使所发展的方法具有实现方便的特点.给出数值计算结果,验证了方法与所发展的多尺度有限元分析程序的正确与有效性.     

含材料非线性的复合材料单钉接头累积损伤分析

发展了静拉伸复合材料接头层合板三维逐渐损伤模型,考虑了单层复合材料在材料1-2面及3-1面上具有明显非线性剪切应力-应变关系的叠层非线性效应,结合有限元技术即应力分析、失效判定准则及损伤过程中材料性能退化等,对接头层合板损伤扩展进行了模拟,结果表明考虑材料非线性的影响与实验结果吻合更好.     

复合材料柱管撞击分析中的有限元方法

建立了复合材料柱管结构在撞击状态下的有限元分析模型,考虑了撞击大变形引起的几何非线性、撞击接触的非线性、纤维树脂复合材料的各向异性及材料的物理非线性．给出了利用该有限元方法分析玻璃环氧柱管撞击吸能特性,并与试验进行了比较,得到了较好的吻合结果,表明该有限元方法是有效的．此外,还讨论了纤维的铺设角度对复合材料管的吸能性能的影响．     

黏弹夹芯板的有限元建模及实验研究

从有限元分析和数值模拟及实验验证的角度研究了黏弹夹芯板的频率依赖振动特性。夹芯板中间层为黏弹性材料，其刚度和阻尼的频率依赖性行为直接影响系统的模态频率和阻尼，并导致振动模式求解的复杂化。采用三阶七参数Biot模型描述黏弹性材料频率相关的黏弹性行为。开发了三层四节点28自由度的夹芯板单元，基于经典板理论和哈密顿原理建立了黏弹夹芯板的有限元动力学方程。通过引入辅助耗散坐标，将Biot模型和黏弹夹芯板的有限元动力学模型结合起来，并将其转化为常规二阶线性系统形式，极大简化了求解非线性振动特性的过程。对一边固定、另三边自由的黏弹夹芯板进行了前三阶固有频率和损耗因子的预测，并与实验结果对比。数值模拟结果和实验结果吻合良好，说明所提有限元方法是正确有效的。     

钢－混凝土组合梁非线性有限元分析

提出了部分剪力连接钢-混凝土组合梁非线性有限元分析的简化模型，该模型同时考虑了钢、混凝土材料的非线性本构关系和剪力连接件的滑移非线性对组合梁刚度和强度的影响，推导了考虑滑移的剪力连接件单元刚度矩阵。利用该模型计算了连续组合梁的极限荷载、挠度、应力以及滑移，与已有的理论计算结果和实验结果吻合，证明本方法分析钢-混凝土组合梁非线性的可靠性和实用性。     

Prediction of stress–strain and fracture behaviour of an 8-Harness satin weave ceramic matrix composite

A computationally economic finite-element-based approach has been developed to predict the stress–strain and fracture behaviour of an 8-Harness satin woven ceramic matrix composite with strain-induced damage. The finite element analysis utilises a solid element to model the behaviour of the homogenised orthotropic uni-directional tow and its matrix. The underpinning models of the tow and matrix, (Tang et al., 2009) capture the physics of the interactions between fibres and matrix; and, in this way, permit modelling that bridges the length scales of the fibres and full-scale components. The non-linear multi-axial stress–strain behaviour of the composite has been discretised by multi-linear elastic curves; and the latter has been used as input to a user defined subroutine, UMAT, in the commercial finite element package, ABAQUS. A partial unit cell model has been constructed of the 8-Harness satin weave composite of carbon fibres embedded in an amorphous carbon matrix, HITCO C/C. Predictions of the global stress–strain curve, which include the effects of fibre waviness, have been made for two failure modes: the first by deformation localisation, and the second by dynamic tow failure on fibre fracture, triggered by instantaneous pull-out deactivation. Comparisons have been made between the predictions and experimental data that exhibit two classes of fracture behaviour: brittle and quasi-ductile. The predicted results, both with and without tow waviness, compare well with the experimental data; however, the predictions for waviness are slightly better. The two extremes of experimental behaviour have been found to correspond with the two tow fracture criteria modelled.     

Virtual testing applied to transverse multiple cracking of tows in woven ceramic composites

The present paper proposes a method of virtual testing with a view to investigating the local response of tows within textile ceramic matrix composite (CMC) under various loading conditions. The method was developed on 2D woven SiC/SiC composites. It capitalizes on knowledge on mechanical damage phenomenology and data established in previous works. It is applied to isolated transverse tows subjected to uniaxial loading by parallel longitudinal tows. The transverse tows contain heterogeneities like matrix voids, fibres and interphases. Mesh for finite element analysis is constructed from micrographs of composite cross section. Cracks were introduced into the mesh for simulation of multiple cracking. Transverse tow tensile behavior and data on distributions of flaw populations were derived from finite element computations of stress-state. Results were compared to experimental observations.     

复合材料波纹梁盒段耐撞性的数值模拟

运用MSC/DYTRAN有限元软件,结合参数等效的方法对复合材料波纹梁盒段的碰撞过程进行了数值模拟,获得了波纹梁盒段碰撞时的破坏过程、载荷-时间曲线、能量吸收能力等数据,将计算结果与M.A.McCarthy的试验结果进行对比,可以发现两者比较吻合。从而说明采用参数等效的方法可以得到复合材料结构耐坠撞性设计中需要的数据。     

Finite element analysis of V-ribbed belts using neural network based hyperelastic material model

A three-dimensional finite element model was built to study V-ribbed belt pulley contact mechanics. The model consists of a pulley and a segment of V-ribbed belt in contact with the pulley. A material model for the belt, including the rubber compound and the reinforcing cord is developed. Rubber is modeled as hyperelastic material. The hyperelastic strain energy function is approximated by neural network trained by rubber test data. Reinforcing cord is modeled as elastic rebar. The material model developed is implemented in the commercial finite element code ABAQUS to simulate the V-ribbed belt-pulley system. A study is then conducted to investigate the effect of belt pulley system parameters on the contact mechanics. The effects of temperature and aging on belt materials are also investigated. The information gained from the analysis can be applied to optimize V-ribbed belt and pulley design.     

用于正交各向异性织物仿真的弹簧质点模型

针对正交各向异性织物,本文提出了一种含附加节点的三角形弹簧质点模型．通过冗余自由度的引入和缩聚,能够实现有限元模型和弹簧质点模型三角形单元刚度矩阵的精确相等;据此,针对一般正交各向异性材料常数和单元形状,推导出弹簧质点模型中弹簧刚度系数及附加节点位置的解析表达式,有效解决了传统弹簧质点模型在正交各向异性织物仿真中存在的计算精度不高、弹簧参数选取随意等缺陷;并进一步完善了基于弹簧质点模型的柔性织物空间变形算法．数值模拟表明本文提出的弹簧质点模型和变形算法在进行织物仿真计算时,具有精度好、通用性强等优点．     

复合材料层合板低速冲击响应和损伤分析

通过在有限元软件Abaqus/Explicit中编写用户材料子程序VUMAT,建立了一种基于能量演化的复合材料低速冲击渐进损伤模型．该模型以三维Hashin准则来预测层内损伤的起始,以一种简化的损伤变量来模拟层内损伤的演化,将具有双线性牵引-分离本构的零厚度界面单元插入层间来模拟分层损伤．采用该模型对14.7 J冲击能量下的纤维增强复合材料低速冲击过程进行了仿真分析,计算得到的冲击力-时间曲线、能量-时间曲线和损伤分布与试验结果吻合较好．根据数值模拟结果,分析了纤维、基体和分层损伤的扩展规律,为低速冲击下的复合材料结构设计提供了理论依据．     

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

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