复合材料层合板多尺度交互渐进损伤分析

纤维增强复合材料强度的准确表征是复合材料力学性能研究的核心问题之一。该文以碳纤维增强树脂基复合材料层合板为研究对象,基于宏观-细观多尺度分析方法,根据复合材料的物理失效模式分别给出了基体和纤维的细观失效准则,同时考虑基体失效对复合材料层合板纤维轴向力学性能的影响。提出了新的刚度退化方式,可准确表征复合材料层合板的损伤演化过程,开展了复合材料层合板四点弯模型的多尺度交互渐进损伤分析和试验验证。结果表明:基于多尺度方法的复合材料层合板宏-细观交互渐进损伤分析结果与试验结果吻合较好,新的刚度退化方式可以准确模拟层合板的失效过程。     

含界面相三维全五向编织复合材料拉伸行为模拟及失效机制研究

基于三维全五向（Q5D）编织复合材料的细观结构模型,通过引入界面相单元,建立了含界面相Q5D编织复合材料单轴拉伸损伤失效分析模型。应用Python语言实现对ABAQUS的二次开发,将Linde等提出的失效准则和Von-Mises应力准则分别用于纱线和基体的渐进损伤判断,并确定材料的整体失效模式;对于界面相,采用Quads准则进行损伤判断。利用周期性位移边界条件,对含界面相Q5D编织复合材料的纵向拉伸应力-应变行为进行了渐进损伤数值模拟,详细讨论了在纵向拉伸载荷作用下材料的细观损伤起始、扩展和最终失效的演化过程,分析了材料的细观损伤失效机制,预测了材料的极限破坏强度,并研究了界面相性能对材料整体力学行为的影响规律。研究结果表明,数值模拟结果与实验值吻合较好,验证了渐进损伤模型的有效性,为该类材料的力学分析和优化设计奠定了基础。     

三维五向编织复合材料渐进损伤分析的数值方法

基于连续损伤理论, 推导了含损伤裂纹的正交各向异性材料的应力-应变关系。建立了考虑界面脱粘破坏的三维五向编织复合材料细观体胞模型, 并将有限元网格尺寸和单元裂纹尺寸引入损伤演化方程。采用Hashin准则和von-Mises准则分别判断纱线与基体的初始损伤, 结合Eshelby-Mori-Tanaka方法确定材料的刚度退化系数。利用ANSYS有限元软件对材料进行渐进损伤分析, 得到了材料在单向拉伸作用下的应力-应变曲线和极限强度。计算结果表明, 轴纱为材料的主要承力部分, 小编织角材料的破坏模式主要为纱线的拉伸断裂, 界面破坏情况较为严重。模拟计算结果与实验结果吻合较好。     

平纹编织SiC/SiC复合材料多尺度建模及强度预测

连续SiC纤维增强SiC基体复合材料（SiC/SiC）具有优异的高温力学性能、辐照稳定性及较低的氚渗透率,在核工程结构领域具有良好的应用前景,掌握其承载状态下的损伤演化和强度性能,对SiC/SiC复合材料的应用具有重要指导意义。本文基于平纹编织SiC/SiC复合材料的制备过程和组分材料分布的多尺度特性,考虑复合材料微观结构的局部近似周期性,建立了纤维丝尺度和纤维束尺度单胞模型。使用有限元分析软件对纤维丝尺度模型的弹性性能和强度性能进行预测,将这些性能参数代入纤维束尺度模型,引入Tsai-Wu失效准则,根据材料的不同失效模式并对失效单元进行方向性刚度折减,模拟了平纹编织SiC/SiC复合材料在单轴拉伸载荷下的渐进损伤过程。数值模拟曲线与试验曲线吻合较好,实现了对平纹编织SiC/SiC复合材料强度的有效预测。      

三维编织复合材料渐进损伤模拟及强度预测

采用考虑纤维束相互挤压的纤维束截面八边形单胞模型, 引入周期性边界条件, 对三维编织复合材料的渐进损伤过程进行数值模拟, 并预测了材料的拉伸强度。通过在应变能密度函数中引入损伤状态变量, 建立了含损伤材料的刚度矩阵, 运用基于不同失效模式下损伤状态变量的刚度渐进折减法表征材料积分点损伤, 通过数值结果与试验结果的对比, 分析了Hashin和Tsai-Wu两种准则作为判定纤维束起始损伤的适用性。分析表明: 基于引入不同失效模式的Tsai-Wu准则的数值模拟和试验结果吻合良好; Hashin准则不适合作为编织纤维束的损伤判据; 不同编织角材料的失效机制不同。      

二维编织陶瓷基复合材料偏轴拉伸力学性能预测

基于二维编织C/SiC复合材料的细观结构,建立了碳纤维丝/热解碳界面/SiC基体和纤维束/表层SiC基体两种尺度下的细观单胞模型,通过有限元方法计算碳纤维丝/热解碳界面/SiC基体模型的等效弹性常数和强度,然后代入纤维束/表层SiC基体模型中计算,并引入Tsai-Wu失效准则,考虑不同失效模式的损伤,建立了二维编织C/SiC复合材料的渐进损伤模型,模拟了其偏轴拉伸应力-应变行为。针对该模型,阐述了二维编织C/SiC复合材料单胞模型在复杂应力状态下其纤维束的损伤过程。数值模拟结果与实验数据吻合较好,验证了模型的有效性,为该种材料的力学性能分析提供了一种有效方法。     

含脱粘损伤的复合材料加筋板压缩破坏渐进损伤数值分析

本文提出了基于连续损伤力学的复合材料层合加筋壁板渐进损伤分析模型。该模型采用界面单元模拟筋条和壁板之间的连接界面,连接界面和复合材料层板分别采用Quads准则和Hashin准则作为失效判据,基于内嵌连续损伤状态变量的材料刚度退化方案,采用非线性有限元方法,研究了复合材料加筋壁板在压缩载荷下的破坏过程,分析了结构相应失效模式的细观失效机理。数值分析结果与实验数据吻合较好,证明该方法的合理性和有效性,并详细探讨了界面单元关键参数和层板铺设角度对加筋壁板结构力学响应的影响规律,得到了一些富有价值的结论。     

碳纤维编织复合材料层合板抗侵彻性能研究

通过弹道冲击实验开展了碳纤维编织复合材料层合板的抗侵彻性能研究，进行了动态响应分析和损伤模式分析。建立了基于Hashin失效和Yeh分层失效准则的渐进损伤模型，运用ABAQUS有限元软件模拟了碳纤维编织复合材料层合板的侵彻失效过程，采用Lambert-Jonas公式拟合了柱状弹侵彻层合板弹道极限曲线，对比分析了碳纤维编织复合材料层合板侵彻实验与数值模拟的弹道极限速度及损伤形貌。结果表明，层合板侵彻损伤模式主要为分层、纤维断裂和基体开裂失效，弹道极限速度数值模拟结果与实验结果吻合较好。     

简谐激励下复合材料加筋板的基体微裂纹损伤演化

研究了在简谐激励作用下复合材料加筋板基体微裂纹损伤的演化行为及其对加筋板动力特性的影响。基于平均微裂纹密度和断裂力学方法, 建立了复合材料加筋板基体微裂纹演化的刚度退化准则。由于该准则考虑了载荷作用周期数的影响, 从而能够更合理地分析周期性动载荷作用下基体微裂纹损伤演化规律。采用Mindlin一阶剪切理论和复合材料模态阻尼模型, 建立了复合材料加筋板动力分析的有限元方法, 研究了在简谐激励作用下, 含分层损伤复合材料加筋板振动过程中诱发的基体微裂纹损伤的演化、 刚度退化, 频率折减和动力响应。      

基于微观力学的复合材料气瓶爆破强度研究

为了研究复合材料储氢气瓶在内压和温度载荷作用下的复杂失效行为,提出一种基于微观力学的气瓶渐进失效分析方法.基于微观力学失效准则,分别判断纤维和基体的损伤与失效,并采用材料刚度退化的方法模拟复合材料失效后的力学行为.此外,建立了复合材料微观力学分析模型,并提取了相应宏-微观分析参数.整个渐进失效分析过程由ABAQUS用户子程序(UMAT)实现,分析得到了复合材料气瓶在热力耦合载荷下的复杂失效模式与最终爆破强度.有限元模拟与试验结果吻合良好,说明该方法可以为复合材料储氢气瓶的设计和优化提供理论指导.     

A new finite point generation scheme using metric specifications

A new approach to generate finite point meshes on 2D flat surface and any bi‐variate parametric surfaces is suggested. It can be used to generate boundary‐conforming anisotropic point meshes with node spacing compatible with the metric specifications defined in a background point mesh. In contrast to many automatic mesh generation schemes, the advancing front concept is abandoned in the present method. A few simple basic operations including boundary offsetting, node insertion and node deletion are used instead. The point mesh generation schemeis initialized by a boundary offsetting procedure. The point mesh quality is then improved by node insertion and deletion such that optimally spaced nodes will fill up the entire problem domain. In addition to the point mesh generation scheme, a new way to define the connectivity of a point mesh is also suggested. Furthermore, based on the connectivity information, a new scheme to perform smoothing for a point mesh is proposed toimprove the node spacing quality of the mesh. Timing shows thatdue to the simple node insertion and deletion operations, the generation speed of the new scheme is nearly 10 times faster than a similar advancing front mesh generator. Copyright © 2000 John Wiley & Sons, Ltd.     

Automatic generation of anisotropic quadrilateral meshes on three‐dimensional surfaces using metric specifications

A new algorithm for constructing full quadrilateral anisotropic meshes on 3D surfaces is proposed in this paper. The proposed method is based on the advancing front and the systemic merging techniques. Full quadrilateral meshes are constructed by systemically converting triangular elements in the background meshes into quadrilateral elements.By using the metric specifications to describe the element characteristics, the proposed algorithm is applicable to convert both isotropic and anisotropic triangular meshes into full quadrilateral meshes. Special techniques for generating anisotropic quadrilaterals such as new selection criteria of base segment for merging, new approaches for the modifications of the background mesh and construction of quadrilateral elements, are investigated and proposed in this study. Since the final quadrilateral mesh is constructed from a background triangular mesh and the merging procedure is carried out in the parametric space, the mesh generator is robust and no expensive geometrical computation that is commonly associated with direct quadrilateral mesh generation schemes is needed. Copyright © 2002 John Wiley & Sons, Ltd.     

Moving Particles Through a Finite Element Mesh

We present a new numerical technique for modeling the flow around multiple objects moving in a fluid. The method tracks the dynamic interaction between each particle and the fluid. The movements of the fluid and the object are directly coupled. A background mesh is designed to fit the geometry of the overall domain. The mesh is designed independently of the presence of the particles except in terms of how fine it must be to track particles of a given size. Each particle is represented by a geometric figure that describes its boundary. This figure overlies the mesh. Nodes are added to the mesh where the particle boundaries intersect the background mesh, increasing the number of nodes contained in each element whose boundary is intersected. These additional nodes are then used to describe and track the particle in the numerical scheme. Appropriate element shape functions are defined to approximate the solution on the elements with extra nodes. The particles are moved through the mesh by moving only the overlying nodes defining the particles. The regular finite element grid remains unchanged. In this method, the mesh does not distort as the particles move. Instead, only the placement of particle-defining nodes changes as the particles move. Element shape functions are updated as the nodes move through the elements. This method is especially suited for models of moderate numbers of moderate-size particles, where the details of the fluid-particle coupling are important. Both the complications of creating finite element meshes around appreciable numbers of particles, and extensive remeshing upon movement of the particles are simplified in this method.     

An alternative pseudoperipheral node finder for resequencing schemes

Many resequencing algorithms for reducing the bandwidth, profile and wavefront of sparse symmetric matrices have been published. In finite element applications, the sparsity of a matrix is related to the nodal ordering of the finite element mesh. Some of the most successful algorithms, which are based on graph theory, require a pair of starting pseudoperipheral nodes. These nodes, located at nearly maximal distance apart, are determined using heuristic schemes. This paper presents an alternative pseadoperipheral node finder, which is based on the algorithm developed by Gibbs, Poole and Stockmeyer. This modified scheme is suitable for nodal reordering of finite meshes and provides more consistency in the effective selection of the starting nodes in problems where the selection becomes arbitrary due to the number of candidates for these starting nodes. This case arises, in particular, for square meshes. The modified scheme was implemented in Gibbs-Poole-Stockmeyer, Gibbs-King and Sloan algorithms. Test problems of these modified algorithms include: (1) Everstine's 30 benchmark problems; (2) sets of square, rectangular and annular (cylindrical) finite element meshes with quadrilateral and triangular elements; and (3) additional examples originating from mesh refinement schemes. The results demonstrate that the modifications to the original algorithms contribute to the improvement of the reliability of all the resequencing algorithms tested herein for the nodal reordering of finite element meshes.     

Analytical prediction of damage due to large mass impact on thin ply composites

This article presents analytical models for predicting large mass impact response and damage in thin-ply composite laminates. Existing models for large mass impact (quasi-static) response are presented and extended to account for damage phenomena observed in thin-ply composites. The most important addition is a set of criteria for initiation and growth of bending induced compressive fibre failure, which has been observed to be extensive in thin ply laminates, while it is rarely observed in conventional laminates. The model predictions are compared to results from previous tests on CFRP laminates with a plain weave made from thin spread tow bands. The experiments seem to confirm the model predictions, but also highlight the need to include the effects of widespread bending induced fibre failure into the structural model.     

Characterising and modelling variability of tow orientation in engineering fabrics and textile composites

Variability of tow orientation is unavoidable for biaxial engineering fabrics and their composites. Since the mechanical behaviour of these materials is strongly dependent on the fibre direction, variability should be considered and modelled as exactly as possible for more realistic estimation of their forming and infusion behaviour and their final composite mechanical properties. In this study, a numerical code, ‘VariFab’, has been written to model realistic full-field variability of the tow directions across flat sheets of biaxial engineering fabrics and woven textile composites. The algorithm is based on pin-jointed net kinematics and can produce a mesh of arbitrary perimeter shape, suitable for subsequent computational analysis such as finite element forming simulations. While the shear angle in each element is varied, the side-length of all unit cells within the mesh is constant. This simplification ensures that spurious tensile stresses are not generated during deformation of the mesh during forming simulations. Variability is controlled using six parameters that can take on arbitrary values within certain ranges, allowing flexibility in mesh generation. The distribution of tow angles within a pre-consolidated glass–polypropylene composite and self-reinforced polypropylene and glass fabrics has been characterised over various length scales. Reproduction of the same statistical variability of tow orientation as in these experiments is successfully achieved by combining the VariFab code with a simple genetic algorithm.     

Optimization of a regularized distortion measure to generate curved high‐order unstructured tetrahedral meshes

We present a robust method for generating high‐order nodal tetrahedral curved meshes. The approach consists of modifying an initial linear mesh by first, introducing high‐order nodes, second, displacing the boundary nodes to ensure that they are on the computer‐aided design surface, and third, smoothing and untangling the mesh obtained after the displacement of the boundary nodes to produce a valid curved high‐order mesh. The smoothing algorithm is based on the optimization of a regularized measure of the mesh distortion relative to the original linear mesh. This means that whenever possible, the resulting mesh preserves the geometrical features of the initial linear mesh such as shape, stretching, and size. We present several examples to illustrate the performance of the proposed algorithm. Furthermore, the examples show that the implementation of the optimization problem is robust and capable of handling situations in which the mesh before optimization contains a large number of invalid elements. We consider cases with polynomial approximations up to degree ten, large deformations of the curved boundaries, concave boundaries, and highly stretched boundary layer elements. The meshes obtained are suitable for high‐order finite element analyses. Copyright © 2015 John Wiley & Sons, Ltd.     

COMBINATION OF ADAPTIVITY AND MESH SMOOTHING FOR THE FINITE ELEMENT ANALYSIS OF SHELLS WITH INTERSECTIONS

A compatible hierarchical adaptive scheme is proposed which allows to control both density and geometrical properties of meshes with four-node linear finite shell elements. The algorithm produces a sequence of meshes with two aims, nearly equal distribution of the local error in each element and a mesh with regular elements, thus internal element angles near 90° and length ratios of adjacent element sides near unity. This goal is achieved in an efficient manner imposing a combination of a local smoothing algorithm with the adaptive mesh generation. New created nodes are positioned on the real shell surface and shell boundaries which may be given e.g. by CAD data. Also the shell directors are determined from the normals on the real geometry. Shell intersections are detected automatically as common curves of two adjacent shell parts. As a shell continuum cannot be assumed for these intersections and thus simple standard adaptive schemes fail, shell intersections have to be treated in a way similar to shell boundaries. For some numerical examples the developed algorithms are demonstrated and the resulting meshes are shown. © 1997 by John Wiley & Sons, Ltd.