发明名称 一种基于细观力学退化模型的复合材料结构失效分析方法
摘要 本发明涉及一种基于细观力学退化模型的复合材料结构失效分析方法,以复合材料三维基本材料性能为输入参数,通过细观力学方法计算得到复合材料损伤后的材料性能,考虑了常见的七种失效模式,包括纤维拉伸,纤维压缩,基体拉伸,基体压缩,纤维-基体剪切,分层拉伸和分层压缩,最终得到复合材料渐进损伤分析中的材料退化模型,并用来进行复合材料结构的失效分析,不仅能够相当准确地预测结构失效强度,同时还能够预测结构的失效模式及位置。本发明基于细观力学退化模型的复合材料结构失效分析方法中,材料退化模型是通过理论分析得出,将复合材料的宏观失效与细观机理联系起来,与现有的基于经验或试验的材料退化模型相比,大大减少了时间和试验成本。
申请公布号 CN103698199B 申请公布日期 2015.11.18
申请号 CN201310447115.0 申请日期 2013.09.26
申请人 北京航空航天大学 发明人 张建宇;周龙伟;山美娟;陈玉丽;赵丽滨
分类号 G01N3/00(2006.01)I;G06F19/00(2011.01)I 主分类号 G01N3/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明;顾炜
主权项 一种基于细观力学退化模型的复合材料结构失效分析方法,其特征在于包括以下步骤:步骤A,根据复合材料结构几何参数,建立三维有限元模型,施加边界条件及初始载荷;步骤B,基于三维有限元模型,进行应力分析,获得复合材料结构各单层的应力分布结果;步骤C,将复合材料结构各单层内的应力结果带入失效准则,进行材料失效判断;步骤D,如果材料没有发生失效,则保持材料性能不变,增大施加在有限元模型上的载荷,转入步骤B;步骤E,如果材料发生失效,进行结构失效判断;步骤F,如果结构没有失效,根据基于细观力学的材料退化模型进行材料刚度退化,增大施加在有限元模型上的载荷,转入步骤B;(F1)首先按照给定的失效准则判断失效部位材料的失效模式;(F2)若发生纤维拉伸失效,其退化系数d<sub>ft</sub>和损伤模量<img file="FDA0000789515480000011.GIF" wi="80" he="87" />分别为:d<sub>ft</sub>=c<sub>m</sub>E<sub>m</sub>/E<sub>11</sub><img file="FDA0000789515480000012.GIF" wi="269" he="98" />其中E<sub>11</sub>为复合材料纵向模量,E<sub>m</sub>为基体模量,c<sub>m</sub>为基体体积含量;(F3)若发生纤维压缩失效,其退化系数d<sub>fc</sub>和损伤模量<img file="FDA0000789515480000013.GIF" wi="82" he="91" />分别为:<img file="FDA0000789515480000014.GIF" wi="1956" he="201" /><img file="FDA0000789515480000015.GIF" wi="281" he="96" />其中S<sub>11</sub>,S<sub>22</sub>,S<sub>12</sub>,S<sub>66</sub>是单层的正轴柔度矩阵元素,E<sub>11</sub>为复合材料纵向模量,角度θ随着压缩载荷的增加在0到π/2之间变化,角度满足:<img file="FDA0000789515480000016.GIF" wi="535" he="155" />(F4)若发生基体拉伸失效,其退化系数d<sub>mt</sub>和损伤模量<img file="FDA0000789515480000017.GIF" wi="300" he="87" />及损伤泊松比<img file="FDA0000789515480000018.GIF" wi="100" he="84" /><img file="FDA0000789515480000021.GIF" wi="77" he="84" />分别为:d<sub>mt</sub>=0<img file="FDA0000789515480000022.GIF" wi="1376" he="93" />其中E<sub>22</sub>为复合材料横向模量,G<sub>12</sub>和G<sub>23</sub>分别为复合材料剪切模量,ν<sub>12</sub>及ν<sub>23</sub>分别为复合材料泊松比;(F5)若发生基体压缩失效,其退化系数d<sub>mc</sub>和损伤模量<img file="FDA0000789515480000023.GIF" wi="296" he="94" />及损伤泊松比<img file="FDA0000789515480000024.GIF" wi="99" he="84" /><img file="FDA0000789515480000025.GIF" wi="79" he="82" />分别为:d<sub>mc</sub>≈0<img file="FDA0000789515480000026.GIF" wi="1404" he="86" />其中E<sub>22</sub>为复合材料横向模量,G<sub>12</sub>和G<sub>23</sub>分别为复合材料剪切模量,ν<sub>12</sub>及ν<sub>23</sub>分别为复合材料泊松比;(F6)若发生分层拉伸失效,其退化系数d<sub>dt</sub>和损伤模量<img file="FDA0000789515480000027.GIF" wi="296" he="86" />及损伤泊松比<img file="FDA0000789515480000028.GIF" wi="99" he="84" /><img file="FDA0000789515480000029.GIF" wi="77" he="79" />分别为:d<sub>dt</sub>=0<img file="FDA00007895154800000210.GIF" wi="1316" he="93" />其中E<sub>33</sub>为复合材料厚度方向模量,G<sub>13</sub>和G<sub>23</sub>分别为复合材料剪切模量,ν<sub>13</sub>及ν<sub>23</sub>分别为复合材料泊松比;(F7)若发生分层压缩失效,其退化系数d<sub>dc</sub>和损伤模量<img file="FDA00007895154800000211.GIF" wi="294" he="90" />及损伤泊松比<img file="FDA00007895154800000212.GIF" wi="98" he="80" /><img file="FDA00007895154800000213.GIF" wi="74" he="83" />分别为:d<sub>dc</sub>≈0<img file="FDA00007895154800000214.GIF" wi="1335" he="94" />其中E<sub>33</sub>为复合材料厚度方向模量,G<sub>13</sub>和G<sub>23</sub>分别为复合材料剪切模量,ν<sub>13</sub>及ν<sub>23</sub>分别为复合材料泊松比;(F8)若发生纤维‑基体剪切失效,其退化系数d<sub>fm1</sub>、d<sub>fm2</sub>和损伤模量<img file="FDA00007895154800000215.GIF" wi="328" he="90" />及损伤泊松比<img file="FDA00007895154800000216.GIF" wi="300" he="86" />分别为:<img file="FDA00007895154800000217.GIF" wi="693" he="165" /><img file="FDA0000789515480000031.GIF" wi="1808" he="99" />其中,G<sub>12</sub>,G<sub>13</sub>和G<sub>23</sub>为复合材料剪切模量,ν<sub>12</sub>,ν<sub>13</sub>及ν<sub>23</sub>分别为复合材料泊松比,<img file="FDA0000789515480000032.GIF" wi="197" he="78" />和<img file="FDA0000789515480000033.GIF" wi="84" he="78" />为表达式分别为:<img file="FDA0000789515480000034.GIF" wi="808" he="209" /><img file="FDA0000789515480000035.GIF" wi="689" he="210" />其中,c<sub>f</sub>为纤维体积含量,G<sub>m</sub>为基体剪切模量;步骤G,如果结构发生失效,计算结束,则此时所施加的载荷即为复合材料结构的失效强度。
地址 100191 北京市海淀区学院路37号