Ansys复合材料结构分析总结

1、Ansys 复合材料结构分析总结说明：整理自Simwe说明：整理自Simwe论坛，复合材料版块，原创fea_stud，大家要感他呀目录1#复合材料结构分析总结（一）-一概述篇5#复合材料结构分析总结（二）-一建模篇10#复合材料结构分析总结（三）分析篇13#复合材料结构分析总结（四）-优化篇做了一年多的复合材料压力容器的分析工作，也积累了一些分析经验， 到了总结的时候了，回想起来，总最初采用I-deas，到MSC.Patran、Nastran，至ij最后选定Ansys为自己的分析工具，确 实有一些东西值得和大家分享，与从事复合材料结构分析的朋友门共同探讨。（一）概述篇复合材料是由一种以上具有不

2、同性质的材料构成，其主要优点是具有优异的材料性 能，在工程应用中典 型的一种复合材料为纤维增强复合材料，这种材料的特性表现为正 交各向异性，对于这种材料的模拟，很多 的程序都提供了一些处理方法，在I-Deas、Nastran、 Ansys 中都有相应的处理方法。笔者最初是用 I-Deas 下建立各项异性材料结合 三维实体结构单 元来模拟（由于研究对象是厚壁容器，不宜采用壳单元），分析结果还是非常好的，而且 l-Deas 强大的建模功能，但由于课题要求要进行压力容器的优化分析， 而且必须要自己写 优化程序，I-Deas 的二次开发功能开放性不是很强，所以改为MSC.Patran, Patran提

3、供了一种非常好的二次开发编程语言PCL （以后在MSC的版中专门给大家贴出这部分容），采用 Patran 结合 Nastran 的分析环境，建立了基于正交各项 异性和各项异性两 种分析模型，但最终发现，在得到的最后结果中，复合材料层之间的 应力结果始终不合理，而模型是没有问 题的（因为在I-Deas中，相同的模型结果是合理的），于是最后转向Ansys，刚开始接触Ansys，真有相 见恨晚的感觉，丰富的单元库，开放的二次开发环境（APDL语言），下面就重点写Ansys的容。在ANSY醉序中，可以通过各项异性单元（Solid 64 ）来模拟，另外还专门提供了一类层合单元（LayerElements

4、 ）来模拟层合结构（Shell 99, Shell 91, Shell 181, Solid 46 和 Solid 191 ）的复合材料。采用ANSYS序对复合材料结构进行处理的主要问题如下：（1 ） 选择单元类型 针对不同的结构和输出结果的要求，选用不同的单元类型。 Shell 99线性结构壳单元，用于较小或中等厚度复合材料板或壳结构，一般长度方向和厚度方向的比值大于10 ；Shell 91 非线性结构壳单元，这种单元支持材料的塑性和大应变行为；Shell 181 有限应变壳单元，这种单元支持几乎所有的包括大应变在的材料 的非线性行为；Solid 46三维实体结构单元，用于厚度较大的复合材料

5、层合壳或实体结大变形。（2）定义层属性配置主要是定义单层的层属性，对于纤维增强复合材料，在这里可以定义单层厚度、 纤维方向等。（3）定义失效准则 支持多种失效准则，不过我还是没有用他，而是自己写了通过应力结果采用二 次蔡胡准则程序来 判断的。（4）其他的一些建模技巧和后处理指导 在我的分析工作中，主要采用了三维实体结构单元。关于Solid 46单元（1）Solid 46是用于模拟复合材料厚壳或实体的8节点三维层合结构单元，单元节点有x,y和z方向三个结构自由度，单元允许最多250层不同的材料；（2）这种单元的定义包括：8 个节点、各层厚度、各层材料方向角和正交各 项异性材料属 性，其中每层可以

6、为面两个方向双线性的不等厚层；（3）在材料定义时，只需定义材料主方向和材料坐标系（单元坐标系）一 致的材料参数， 不一致的复合材料层通过定义材料方向角（该层材料主方向和材料坐标系所成的角度）由程序自动 转换；（4）通过选择不同的层直接在单元坐标下获取单元应力，包括三个方向的 应力和面剪切应力，而不需要通过应力应变的转换来获取；论坛问答：Q ANSYS 如何处理失效后的材料退化呢？A:ANSYS没有直接提供材料失效后的退化，但可以自己写程序让ANSYS八行。ANSYS可以用失效准则判断材料是否失效，之后刚度降低可以通过实验测得。再将实验数据输入到ANSYS中，对失效的单元重新进行分析。共同讨论！

7、 Ansys确实没有直接提供材料失效后的退化的处理方法。我们在进行复合材料结构分析时，通常采 用单层模量退化的估算方法，这种估算方法就是将带有裂纹层的横向、 剪切模量与泊松系数全部用一组经过DF 因子退化的新值替代，为了考虑压缩强度的下降， 对单向复合材料的压缩强度也要 DF 因子退化（详细信息可以参考蔡为仑的 复合材料设计 一书），这样，我们就可以再结合Ansys的APDL来处理了。建模篇复合材料是一种各向异性材料，对于纤维增强复合材料又是一种正交各向异性材料，因 此，在进行复合材料结构建模的时候要特别注意的一个重要的问题，就是材料的方向性。 下面，就我个人的分析经验，对复合材料结构的建模作

8、一个总结。. 结构坐标系、单元坐标系、材料坐标系和结果坐标系 建立复合材料结构模型，存在一个结构坐标系，用于确定几何元素的位置，这个坐标可以是笛卡尔坐标系、 柱坐标系或者是球坐标系；单元坐标系是每个单元的局部坐标系，一般用来描述整个单元； 材料坐标系是确定材料属性方向的坐标系，一般没有专门建立的材料坐标系，而是参考其他坐标系，如整体结构坐标系，或单元坐标系，在Ansys程序中，材料坐标是由单元坐标唯一确定的， 要确定材料坐标，只要确定单元坐标就行了；结果坐标系是在进行结果输出时所使用的坐标系，也是一般参考其他坐标系。在Ansys程序中，关于坐标系有人做过专门的总结。见后。. 用于复合材料结构分

9、析的单元用于复合材料分析的单元主要有两类， 一类是层合单元，如 Shell 99, Shell 91, Shell181, Solid 46 和 Solid 191 ;另一类是各向异性单元，如Solid64 ;这些材料都有不同的 处理方法，层合单元，在一个单元可以包含多层信息，包括各层的材料、厚度和方向；各项 各向异性单元，在 一个单元，只能包含一种材料信息， 而且所得到的计算结果还要进行一些 处理，因此有一定的局限性。3. 单元坐标的一致性问题在进行复合材料结构建模的时候， 有些时候结构几何比较复杂， 很难用统一的坐标 来确定单元坐标系，即使对一些规则的几何（如圆桶） ，在用旋转方法生成几何

10、时，不同的 面法向也会带来单元坐标的不一致，这就使得材料输入的时候存在问题并使计算结果错误， 因此，在几何建模 时要特别注意这一问题， 笔者也没有得到一些复杂几何进行单元划分时保 持单元一致的合适方法。4. 一个实例5. 下面的命令流显示了不同的几何生成方法会产生不同的单元坐标方向 /PREP7 !*Create Material*MPTEMP,MPTEMP,1,0MPDATA,EX,1,2.068e8MPDATA,PRXY,1,0.29MPTEMP,MPTEMP,1,0*MPDATA,DENS,1,7.82e-6*Create Eleme nt Type*ET,1,SOLID95KEYOPT

11、,1,1,1KEYOPT,1,5,0 KEYOPT,1,6,0KEYOPT,1,11,0I*CSYS,1HS=80!*create two keypo ints along axial K,101,0,0,0,K,102,0,0,400,!*create keypo ints K,1,61,0,0,K,2,HS,0,0,K,5,100,0,0, K,11,61,0,178,K,12,HS,0,178, K,15,HS+10,0,178, K,111,61,0,178,K,112,HS,0,178,K,115,HS+10,0,178, K,21,61,0,2450, K,22,HS-4,0,24

12、50,K,25,HS+6,0,2450,I*!*create areas by keypo ints FLST,2,4,3FITEM,2,21FITEM,2,111FITEM,2,112FITEM,2,22A,P51XFLST,2,4,3FITEM,2,22FITEM,2,112FITEM,2,115FITEM,2,25A,P51XI*FLST,2,2,5,ORDE,2FITEM,2,1FITEM,2,-2FLST,8,2,3FITEM,8,101FITEM,8,102VROTAT,P51X, , , , , ,P51X, ,90,1,TYPE, 1MAT,1REAL,ESYS,0SECNUM

13、,MSHAPE,0,3DMSHKEY,1FLST,5,2,6,ORDE,2FITEM,5,1FITEM,5,-2CM,_Y,VOLUVSEL, , , ,P51XCM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2 运行上述命令流，查看一下单元坐标，再把命令流中下列部分FLST,2,4,3FITEM,2,21FITEM,2,111FITEM,2,112FITEM,2,22A,P51X改为：FLST,2,4,3FITEM,2,22 FITEM,2,21 FITEM,2,111 FITEM,2,112 A,

14、P51X 再看一下单元坐标。ANSYS 坐标系总结工作平面（Working Plane ）工作平面是创建几何模型的参考（X,Y）平面，在前处理器中用来建模（几何和网格）总体坐标系在每开始进行一个新的ANSYS八析时，已经有三个坐标系预先定义了。它们位于模型的总体原点。三种类型为：CS,0:总体笛卡尔坐标系CS,1:总体柱坐标系 CS,2:总体球坐标系 数据库中节点坐标总是以总体笛卡尔坐标系，无论节点是在什么坐标系中创建的。 局部坐标系 局部坐标 系是用 户定 义的 坐标系。局部 坐标 系可 以通过 菜单路 径 WorkplaneLocal CSCreate LC 来 创建。激活的坐标系是分析中

15、特定时间的参考系。缺省为总体笛卡尔坐标系。当创建了一个新的坐标系时，新坐标系变为激活坐标系。 这表明后面的激活坐标系的命令。 菜单中激活坐标系的 路径 WorkplaneChange active CS to 。节点坐标系 每一个节点都有一个附着的坐标系。节点坐标系缺省总是笛卡尔坐标系并与总体笛卡尔坐标 系平行。节点力和 节点边界条件（约束）指的是节点坐标系的方向。时间历程后处理器/POST26中的结果数据是在节点坐标系下表达的。而通用后处理器/POST1中的结果是按结果坐标系进行表达的。例如：模型中任意位置的一个圆，要施加径向约束。 首先需要在圆的中心创建一个柱坐标系并分配一个坐标系（例如C

16、S,11）。这个局部坐标系现在成为激活的坐标系。然后选择圆上的 所有节点。通过使用 Prep7Move/ModifyRotate Nodal CS to active CS,选择节点的节点坐标系的朝向将沿着激活坐标系的方向。未选择节点保持不变。 节点坐标系的显示通过菜单路径PltctrlsSymbolsNodal CS这些节点坐标系的X方向现在沿径向。约束这些选择节点的X方向，就是施加的径向约束。 注意：节点坐标系总是笛卡尔坐标系。可以将节点坐标系旋转到一个局部柱坐标下。这种情 况下，节点坐标系的X方向指向径向，Y方向是周向（theta ）。可是当施加theta方向非零位移时，ANSYS 总是

17、定义它为一个笛卡尔丫位移而不是一个转动（丫位移不是theta位移）。单元坐标系 单元坐标系确定材料属性的方向（例如，复合材料的铺层方向）。对后处理也是很有用的，诸如提取梁和壳单元的膜力。单元坐标系的朝向在单元类型的描述中可以找到。结果坐标系/Postl 通用后处理器中（位移，应力，支座反力）在结果坐标系中报告，缺省平行于总体笛 卡尔坐标系。这 意味着缺省情况位移，应力和支座反力按照总体笛卡尔在坐标系表达。无论节点和单元坐标系如何设定。 要恢复径向和环向应力， 结果坐标系必须旋转到适当的坐标系下。这可以通过菜单路径Post1Options for output实现。/POST26时间历程后处理器

18、中的结果总是以节点坐标系表达。显示坐标系显示坐标系对列表圆柱和球节点坐标非常有用（例如，径向，周向坐标）。建议不要激活这个坐标系进行显示。屏幕上的坐标系是笛卡尔坐标系。 显示坐标系为柱坐标系，圆弧将显示为 直线。这可能引起混乱。因此在以非笛卡尔坐标系列表节点坐标之后将显示坐标系恢复到总 体笛卡尔坐标系。 分析篇下面就我对碳纤维增强复合材料压力容器分析过程中所做的工作，从复合材料材料参数转 化、复合材料强度准则、结构刚强度分析几方面写些我的心得，与大家共同探讨。复合材料材料参数的转化 单向纤维增强复合材料（也称单向板）是指纤维按照同一方向平行排列的复合材料，是构成层合板和壳的基本元素，可认为是一

19、种正交各向异性材料，也是一种横观各向同性材料（存 在一个各向同性面）， 在进行有限元计算时，必须知道复合材料的弹性特性参数，并由弹性特性参数来计算正交各向异性材料的9个参数（在ANSY皐序中定义材料时所需3个弹性模量、3个泊松系数和3个剪切模量），单向复合材料特性的计算有许多种方法，主要的方法有Halpin-Tai的弹 性力学方法，这种方法根据弹性理论将复杂的纤维与树脂间的关系用一 组方程来表示，通过求解方程组，解得弹性参数，我们使用的9个弹性参数的计算是通过单 向复合材料的刚度 矩阵转化得到，下面是用APDL 语言编写的材料转化程序。MAT_PAR_COMP!*this macro is u

20、sed to calculate material parameters of compositeI*E1=1.81E8E2=1.03E7V21=0.28V12=E2*V21/E1V23=0.5V32=0.5G12=7.17E6RM=C0S（ARG1 ）RN=SIN（ARG1）RM2=RM*RMRM4=RM2*RM2RN2=RN*RNRN4=RN2*RN2RMN=RM*RNRMN2=RMN*RMN!* caculate stiffness matrice of uni direct ional composite materialW=（1.0+V23）*（1.0-V23-2.0*V21*V12

21、）VV=1.0/VVQ11=（1.0-V23*V32）*W*E1Q22=（1.0-V21*V12）*W*E2Q33=Q22Q12=V21*(1.0+V23)*VV*E2Q13=Q12Q23=(V23+V21*V12)*W*E2Q44=(1.0-V23-2.0*V21*V12)*W*E2*0.5Q55=G12Q66=Q55!* calculate equivale nt stiffness of composite material*HQ1 仁 Q11*RM4+2.0*(Q12+2.0*Q66)*RMN2+Q22*RN4HQ12=(Q11+Q22-4.0*Q66)*RMN2+Q12*(RM4+R

22、N4)HQ13=Q13*RM2+Q23*RN2HQ23=Q13*RN2+Q23*RM2HQ16=-RMN*RN2*Q22+RM2*RMN*Q11-RMN*(RM2-RN2)*(Q12+2.0*Q66)HQ22=Q11*RN4+2.0*(Q12+2.0*Q66)*RMN2+Q22*RM4HQ33=RN2*Q13+RM2*Q23HQ33=Q33HQ26=-RMN*RM2*Q22+RMN*RN2*Q11+RMN*(RM2-RN2)*(Q12+2.0*Q66)HQ36=(Q13-Q23)*RMNHQ44=Q44*RM2+Q55*RN2HQ45=(Q55-Q44)*RMNHQ55=Q55*RM2+Q44

23、*RN2HQ66=(Q11+Q22-2*Q12)*RMN2+Q66*(RM2-RN2)*(RM2-RN2)QQ1 仁 HQ11QQ12=HQ12QQ22=HQ22QQ13=HQ13QQ23=HQ23QQ33=HQ33QQ44=(HQ44*HQ55-HQ45*HQ45)/HQ55QQ55=(HQ44*HQ55-HQ45*HQ45)/HQ44QQ66=HQ66Q(1)=QQ11Q(2)=QQ12Q(3)=QQ13Q(4)=QQ22Q(5)=QQ23Q(6)=QQ33Q(7)=QQ66Q(8)=QQ44Q(9)=QQ55!*QQQ=Q(1)*(Q (4) *Q (6)-Q (5) *Q(5)-Q(2

24、)*(Q (2) *Q (6)-Q (3)*Q (5) )+Q(3)*(Q (2) *Q (5)-Q(3)*Q(4)S 仁(Q (4)*Q( 6)-Q(5)*Q(5)/QQQS2=-(Q (2)*Q (6)-Q(3)*Q(5)/QQQS3=(Q (2)*Q (5)-Q(3)*Q(4)/QQQS4=(Q(1)*Q(6)-Q(3)*Q(3)/QQQS5=-(Q(1)*Q (5)-Q(2)*Q(3)/QQQS6=(Q(1)*Q (4)-Q(2)*Q(2)/QQQS7=1/Q S8=1/Q(8)S9=1/Q(9)EEX=1/S1EEY=1/S4EEZ=1/S6VXY=-S2*EEXVXZ=-S3*EE

25、XVYZ=-S5*EEYGXY=1/S7GYZ=1/S8GXZ=1/S9/EOF复合材料强度准则复合材料结构的受力及应力应变情况非常复杂，并要考虑各种应力应变的耦合和相互影响， 复合材料强度破坏准则基于结构的宏观破坏，一般来说复合材料的二次蔡-吴强度破坏准则较为精确。有兴趣的朋友可以参考科学出版的蔡为仑先生的复合材料设计这一本书。复合材料结构刚强度分析一般说来，复合材料结构总是受到空间力的作用，其应力分布是三维的，因此，复合材料结构的刚强度分析一般不宜采用复合材料的板壳理论(这种理论仅考虑板壳面的应力和横向剪切应力，而忽略法向应力)，同时，对于简单的结构(如板、壳)，可以得到弹性力学的一般解，

26、而对于大多数结 构来说，则必须用数值的方法计算，三维有限元分析是最常用的方法。采用 ANSY 皐序对复合材料进行刚强度分析的步骤如下：建立结构的几何模型由于复合材料分析单元一般都是六面体单元，因此，在建立几何时要特别考虑到网格 划分的方便。建立材料模型根据复合材料材料参数建立单向复合材料材料模型，我所采用的是碳纤维增强复合材 料，有两种建立方 法。a. 若选择单元为各向异性单元， 则根据单向复合材料的刚度矩阵或柔度矩阵建立各向异 性材料模型；b. 若选择层合单元，则可以建立相关的材料模型， 如单向复合材料则可以建立正交各向 异性材料模型选择单元类型并设置相关属性根据结构特征和计算要求，选择不同

27、的单元类型并设置单元属性(各种单元的选择依 据请参考概述篇或ANSY 醐助文件)网格划分在建立的几何实体上进行网格划分，对于复合材料，选择六面体三维实体单元，定义 单元属性，分别指定不同的材料属性，并保证材料坐标一致，运用有限元网格生成器进行网 格划分。（5）定义边界条件根据实际情况定义边界条件。（6）分析设定并提交计算设定分析类型及相关一些参数（7）结果后处理复合材料结构的分析结果在进行后处理时，非常重要的一点是选择合适的并与计算时所用的坐标一致的 结果坐标系，如对于回转体结构选择计算时的柱坐标。另外，对于用各向异性单元（Solid64 ）来模拟的计算结果在结果处理时必须保证应力应变关系的一

28、致，主要是在不同种复合材 料层间或者同一种复合材料不同铺层方向的层之间界面的应力应变情况， ANSYSH 处理中所得到的结果不完全 是正确的，应该根据法向应力联系，面应变连续的准则 来进行处理。与传统材料相比，复合材料具有可设计性 大的灵活性，复合材料的力学性能和机械性能， 环境条件要求，通过组分材料的选择匹配、 的达到预期目的，以满足工程设备的使用性能，复合材料结构分析总结（四）一一优化篇与传统材料相比，复合材料具有可设计性 大的灵活性，复合材料的力学性能和机械性能， 环境条件要求，通过组分材料的选择匹配、 的达到预期目的，以满足工程设备的使用性能，复合材料结构的多层次性为复合材料及其结 构

29、设计带来了极都可按照结构的使用要求和铺层设计及界面控制等材料设计手段，最大限度因此，在工程实践中对复合材料结构进行优将复合材料结构优化设计有很重要的现实意义， 下面以我所研究的复合材料压力容器为例 化以及在 ANSYSF 的实现过程给大家作一个介绍。将复合材料结构优1. 问题描述本文所涉及的复合材料压力容器是带有金属胆外缠碳纤维增强复合材料的复合容器，优化问题是：以金属胆壁厚、复合材料各缠绕层厚度和缠绕角为设计变量，在满足压力容器强度（金属胆层和复合材料层均满足强度要求）和重量要求的条件下，使压力容器的刚度最大。2. 优化模型根据纤维增强复合材料特性， 压力容器环向缠绕复合材料有利于提高容器刚

30、度，轴向平铺复合材料有利于提高容器刚度， 因此，模型采用 3 种缠绕角的方案，即靠近金属胆为环向 （90 度）缠绕，中间为缠绕，外部为轴向平铺（0 度），以各层的厚度（金属层和三层复合材料）和中间缠绕层的角度为优化参数，在压力容器强度约束的条件下，以压力容器一阶固有频率为优化目标。其数学模型如下：Maximize: f X ,其中 X (捲必,怡,人)丁 (hi, h?, hh, )TSubject to:z “ sh h2hs H,( h, h2, h3 0),0，90$(X) 1.2,S2(X) 1.5,c(X) Co其中，f为复合材料压力容器的一阶固有频率，S1和S2分别为金属胆的安全系

31、数和各复合材料层的强度比，通过有限元程序求得，为中间层复合材料缠绕角，hj、h2和hs分别为金属胆厚度、90度缠绕层厚度和度缠绕层厚度，H为hi、h2和hs的极限值，当总厚度确定后，0度缠绕层厚度由、h2、hs及总厚度确定，c为复合容器重量，C0为全压力容器重量上限。3. 优化算法基于 ANSYS 勺优化，可以直接使用 ANSYS 提供的优化模块， 根据上述优化模型， 建立优 化计算文 件，选择合适的优化算法，进行计算。同时，也可以通过APDL语言(甚至可以通过外部编程环境，如VC+, FORTRAN八)来自己编制优化算法，本文就是通过自己编制优化算法来实现的，采用的优化算法是复形调优法。算法描述如下：复形调优法是求解约束条件下n 维极值问题的重要方法， 通过构造复合形，计算各顶点的目标函数值，并