编织型层合结构复合材料力学性能研究
摘要
复合材料技术是一项具有战略意义的国防关键技术，在一定程度上，先进复合材料的研究水平和应用程度是一个国家科技发展水平的代表，特别是在飞机制造业，各种先进的飞机无不与先进的材料技术紧密。随着复合材料在飞机上应用比例的加大，在复合材料制造领域，先进的数字化设计制造技术和计算机辅助工程(CAE)技术等得到了广泛应用。
本研究利用本课题组特色技术，以国家/行业标准为基准，以大型商用有限元软件ABAQUS软件作为建模计算平台，采用Python语言脚本参数化设计语言，对某飞机蒙皮CKFR-PMCs编织型层合结构复合材料的性能进行了研究。完成了铺层结构代表性体积单元(RVE)和相关力学性能试验的试样数值模型的建立，模拟了包括拉伸/压缩试验、剪切试验、弯曲性能试验、层间剪切强度试验、纵横剪切试验在内的复合材料基本力学试验。应用Python语言编写脚本程序，实现了编织型层合结构复合材料有限元模型中纤维体积分数含量的控制以及各向异性材料局部坐标与整体坐标下输出变量的转换。对层合结构复合材料各结构层的碳纤维/芳纶纤维材料、基体材料的性能进行了反推预测，对40%纤维含量复合材料性能的进行了验证性计算。
本文对编织型层合结构复合材料的力学性能预测给出了较为完整的思路与方法。对于不同性能的纤维、基体材料，不同纤维体积分数含量，不同纤维铺设铺层角度，不同加载方向的纤维增强复合材料和编织复合材料的力学性能预测具有一定的通用性。对后面开展编织型层合结构复合材料的分层多尺度研究和跨尺度研究提供了数值上的准备与方法上的参考。
Composite technology is a strategic key technology for national defense. To someextent, the research and application of the advanced composite represent thedevelopment level of a nation's science and technology. Especially in aircraft industry,various kinds of advanced aircraft are closely related to advanced materials technology.As the application of composite increases in aircraft industry, advanced digital designmanufacturing technology and computer aided engineering (CAE) technology havebeen widely applied in composite manufacturing field.
This research use our special skill to study the efficiency of the braided laminatedcomposite of a CKFR-PMCs, with the national/industry standard as a benchmark, thelarge commercial finite element software ABAQUS as modeling computing platforms,the Python language script parameter to design language. The study has established thelayer structure representative volume unit (RVE) and related sample numerical modelsof mechanical properties test, including the basic mechanical tests of composite such asextrusion/compression test, shear test, flexural performance test, interplay shearstrength test,and longitudinal transverse shear test.The study applies Python languageto program and achieves to control the portion of the fiber volume fraction in braidedlaminated composite finite element models and the conversion of the output variablesof the anisotropic material from the local coordinate system to global coordinatesystem. What’s more, the study also makes an inverse prediction of the performance ofthe carbon/Kevlar fiber and matrix materials in the layers of braided laminatedcomposite and calculates the performance of the composite containing40%fiber.
This paper gives a more complete idea and method to predict the mechanicalproperty of the braided laminated composite. The results also can be applied into theprediction of different properties of the fiber, matrix materials, different fiber volumefraction content, and different fiber laying layer angles, different loading directions offiber reinforced materials and the mechanical performance of the braided laminatedcomposite. It will also offer some numerical preparation and method to carry outfurther multi-scale and cross-scale research on the braided laminated composite.
本研究利用本课题组特色技术，以国家/行业标准为基准，以大型商用有限元软件ABAQUS软件作为建模计算平台，采用Python语言脚本参数化设计语言，对某飞机蒙皮CKFR-PMCs编织型层合结构复合材料的性能进行了研究。完成了铺层结构代表性体积单元(RVE)和相关力学性能试验的试样数值模型的建立，模拟了包括拉伸/压缩试验、剪切试验、弯曲性能试验、层间剪切强度试验、纵横剪切试验在内的复合材料基本力学试验。应用Python语言编写脚本程序，实现了编织型层合结构复合材料有限元模型中纤维体积分数含量的控制以及各向异性材料局部坐标与整体坐标下输出变量的转换。对层合结构复合材料各结构层的碳纤维/芳纶纤维材料、基体材料的性能进行了反推预测，对40%纤维含量复合材料性能的进行了验证性计算。
本文对编织型层合结构复合材料的力学性能预测给出了较为完整的思路与方法。对于不同性能的纤维、基体材料，不同纤维体积分数含量，不同纤维铺设铺层角度，不同加载方向的纤维增强复合材料和编织复合材料的力学性能预测具有一定的通用性。对后面开展编织型层合结构复合材料的分层多尺度研究和跨尺度研究提供了数值上的准备与方法上的参考。
Composite technology is a strategic key technology for national defense. To someextent, the research and application of the advanced composite represent thedevelopment level of a nation's science and technology. Especially in aircraft industry,various kinds of advanced aircraft are closely related to advanced materials technology.As the application of composite increases in aircraft industry, advanced digital designmanufacturing technology and computer aided engineering (CAE) technology havebeen widely applied in composite manufacturing field.
This research use our special skill to study the efficiency of the braided laminatedcomposite of a CKFR-PMCs, with the national/industry standard as a benchmark, thelarge commercial finite element software ABAQUS as modeling computing platforms,the Python language script parameter to design language. The study has established thelayer structure representative volume unit (RVE) and related sample numerical modelsof mechanical properties test, including the basic mechanical tests of composite such asextrusion/compression test, shear test, flexural performance test, interplay shearstrength test,and longitudinal transverse shear test.The study applies Python languageto program and achieves to control the portion of the fiber volume fraction in braidedlaminated composite finite element models and the conversion of the output variablesof the anisotropic material from the local coordinate system to global coordinatesystem. What’s more, the study also makes an inverse prediction of the performance ofthe carbon/Kevlar fiber and matrix materials in the layers of braided laminatedcomposite and calculates the performance of the composite containing40%fiber.
This paper gives a more complete idea and method to predict the mechanicalproperty of the braided laminated composite. The results also can be applied into theprediction of different properties of the fiber, matrix materials, different fiber volumefraction content, and different fiber laying layer angles, different loading directions offiber reinforced materials and the mechanical performance of the braided laminatedcomposite. It will also offer some numerical preparation and method to carry outfurther multi-scale and cross-scale research on the braided laminated composite.
引文
[1]梁军，陈晓峰，庞宝君.多向编织复合材料的力学性能研究[J].力学进展，1999，29(2)
[2]马文锁，赵允岭，冯伟.三维编织复合材料理论研究进展[J].材料科学与工程学报，2006，24(4)：631-636
[3]詹世革，孟庆国，方岱宁.航空航天纺织结构复合材料力学性能研究进展与展望[J].中国基础科学(综述评述)，2008，2：5-10
[4]杜善义.先进复合材料与航空航天[J].复合材料学报，2007，24(1)：1-11
[5]师昌绪.材料大词典.北京:化学工业出版社，1994：2-8
[6] Robert I M. Advanced composite structures research in Australia[J]. CompositeStructures，2002，57：3-10
[7]益小苏.航空先进复合材料的应用技术及其数字化与CAE[R].中国科学院技术科学论坛第二十三次学术报告会议(CAE自主创新发展战略).上海：中国科学院，2006
[8]益小苏.先进复合材料技术研究与发展[M].北京：国防工业出版社，2006
[9]益小苏，李宏运，连建民等.先进复合材料的CAE虚拟技术[J].航空制造技术，2008，14：34-39
[10]钟群鹏.失效分析基础知识[M].理化检验(物理分册)，2005，2(1)：35-37
[11] Mouritz A P. Review of applications for advanced three-dimensional fiber textilecomposites[J]. Composites: Part A，1999，30：1445-1461
[12] Mouritz A P，Bannister M K，Falzon P J，et al. Review of applications foradvanced three-dimensional fibre textile composites[J]. Compos Part A，1999，30(12)：1445
[13] Wang Y Q，Wang A S D. Spatial distribution of yarns and mechanical properties in3D braided tubular composites[J]. Appl Compos Mater，1997，4(2)：121
[14] Jan Z，Michal S. Homogenization of balanced plain weave composites withimperfect microstructure：Part I-Theoretical formulation[J]. Int J Solids Struct，2004，41(22/23)：6549
[15] Chiu C H，Tsai K H，Huang W J. Crush-failure modes of2D triaxially braidedhybrid composite tubes[J]. Compos Sci Technol，1999，59：1713
[16] Quek S C，Waas A，Shahwan K W，et al. Compressive response and failure ofbraided textile composites：Part2Computations[J]. Int J Non Linear Mech，2004，39(4)：649
[17]庞宝君，杜善义，韩杰才.三维四向编织复合材料细观组织及分析模型[J].复合材料学报，1999，16(3)：135
[18]马文锁，冯伟.二维编织复合材料几何结构的平面群分析[J].北京科技大学学报.2007，29(2)：226-230
[19]于平，周平等.二维平纹编织复合材料压缩力学行为研究[J].计算力学学报，2008，25(2)：224-228
[20]陈陆平，潘敬哲，钱令希.复合材料纤维/基体界面失效问题的参变量有限元数值模拟[J].复合材料学报，1993，10(1)：381-393
[21]张双寅.复合材料层间剪切强度[J].工程力学，1985，3(4)：26-27
[22]范天佑.断裂理论基础[M].北京：科学技术出版社，2003，76-78
[23]王铎，杜善义，王殿富等.断裂力学[M].黑龙江：哈尔滨工业大学出版社，1989，110-111．
[24]赵艳华，徐世烺. I-II复合型裂纹脆性断裂的最小J2准则[J].工程力学，2002，19(4)：94-98
[25]张磊，谢宗蕻，苏霓.离位增韧复合材料层压板准静态压痕试验研究[J].航天器环境工程，2009，12(26)：128-130
[26]钟明辉.编织复合材料[J].纤维复合材料，1992，2：42-46
[27] ISH IKAWA T，CHOU T W. Stiffness and strength behavior of woven fabriccomposite[J]. Mat Sci，1982，17(11)：32l1-l3220
[28]. ISH IKAWA T，CHOU T W. One two-dimensional micromechanical analysis ofwoven fabric composite[J]. A IAA，1983，21：1714-1721
[29]. ISH IKAWA T， MATSUSH IMA M， HAYASHI Y， et al. Experimentalconfirmation of theory of elastic modulus of fabric composites[J]. CompositeMaterials，1985，19(9)：443-458
[30] HUANG Z M. The mechanical properties of composites reinforced with wovenand braided fabrics[J]. Composites Science and Technology，2000，60：479-498
[31]燕瑛，楼畅，成传贤等.机织复合材料力学性能的细观分析与实验研究[J].复合材料学报，2001，18(2)：109-113
[32]卢子兴，胡奇.三维编织复合材料压缩力学性能的实验研究[J].复合材料学报，2003，20(6)：67-72
[33] GUDEN M，YILDIRIM U，HALL W I. Effect of strainrate on the compressionbehavior of a woven glass fiber/SC-15composite[J]. Polymer Testing，2004，23：719-725
[34] AKIL O，YILDIRIM U，GUDEN M，et al. Effent of strain rate on the compressionbehavior of a woven fabric S2-glass fiber reinforced vinylester composite[J]. PolymerTesting，2003，22：883-887
[35]刘书山.复合材料应力分析的均匀化方法[J].力学学报，1997，29(3)，306-313
[36]刘书山.复合材料弹性性能尺度效应[J].大连理工大学学报，2004，44(2)：200-205
[37]宋迎东，雷友锋.通用单胞法模型预测纤维增强钛基复合材料弹性性能[C].第十届航空发动机结构强度与振动会议征文集.安徽歙县，2000，10：137-141
[38]雷友锋，宋迎东，高德平.纤维增强金属基复合材料弹塑性行为分析[J].机械科学与技术，2003，22(6)：982-986
[39]雷友锋，宋迎东，高德平.纤维增强金书基复合材料细观几何结构对宏观弹性性能的影响[J].宇航材料工艺，2003，33(1)：43-48
[40]宋迎东，雷友锋，孙志刚等.一种新的纤维增强复合材料细观力学模型[J].南京航空航天大学学报，2003，35(4)：435-440
[41]曹礼群，崔俊芝.整周期复合材料弹性结构的有限元计算[J].计算数学，1998，20(3)：279-290
[42]崔俊芝.整周期复合材料弹性结构的双尺度渐近分析方法[J].应用数学学报，1999，21(1)：38-46
[43] Blassiau S，Thionnet A，Bunsell A R. Micromechanics of load transfer in aunidirectional carbon fiber-reinforced epoxy composite due to fibre failures：Part3.Multiscale reconstruction of composite behaviour[J]. Composite Structures，2007，1-12
[44] Naimark O B，Davydova M，Plekhov O A. Nonlinear and structural aspects oftransitions from damage to fracture in composites and structures[J]. Computers andStructures，2000，76：67-75
[45] Tang X D，Whitcomb J D，Kelkar A D. Progressive failure analysis of2×2braidedcomposites exhibiting multiscale heterogeneity[J]. Composites Science andTechnology，2006，66：2580-2590
[46]刘昊，陈宪战，傅德莲. CAE技术现状及其在石油物探机械装备设计中的应用[J].物探装备，2010，20(2)：92-95
[47] Glimm J， Sharp H D. Multiscale science： a challenge for the twenty-firstcentury[J]. SIAM News，1997，30(8)：1~7
[48]曹礼群.材料物性的多尺度关联与数值模拟[J].材料科学研究，2002，24(6)：23-30
[49]沈观林，胡更开.复合材料力学[M].北京：清华大学出版社，2006，3-4
[50]郑晓霞，郑锡涛，缑林虎.多尺度方法在复合材料力学分析中的研究进展[J].力学进展，2010，40(1)：41-47
[51]曹金凤，王旭春，孔亮. Python语言在ABAQUS中的应用[M].北京：机械工业出版社.2011：2-8
[52]石亦平，周玉蓉. ABAQUS有限元分析实例详解[M].北京：机械工业出版社.2006：2-10
[53] Wesley J Chun. Python核心编程(第二版)[M].北京：人民邮电出版社.2008：3-9
[2]马文锁，赵允岭，冯伟.三维编织复合材料理论研究进展[J].材料科学与工程学报，2006，24(4)：631-636
[3]詹世革，孟庆国，方岱宁.航空航天纺织结构复合材料力学性能研究进展与展望[J].中国基础科学(综述评述)，2008，2：5-10
[4]杜善义.先进复合材料与航空航天[J].复合材料学报，2007，24(1)：1-11
[5]师昌绪.材料大词典.北京:化学工业出版社，1994：2-8
[6] Robert I M. Advanced composite structures research in Australia[J]. CompositeStructures，2002，57：3-10
[7]益小苏.航空先进复合材料的应用技术及其数字化与CAE[R].中国科学院技术科学论坛第二十三次学术报告会议(CAE自主创新发展战略).上海：中国科学院，2006
[8]益小苏.先进复合材料技术研究与发展[M].北京：国防工业出版社，2006
[9]益小苏，李宏运，连建民等.先进复合材料的CAE虚拟技术[J].航空制造技术，2008，14：34-39
[10]钟群鹏.失效分析基础知识[M].理化检验(物理分册)，2005，2(1)：35-37
[11] Mouritz A P. Review of applications for advanced three-dimensional fiber textilecomposites[J]. Composites: Part A，1999，30：1445-1461
[12] Mouritz A P，Bannister M K，Falzon P J，et al. Review of applications foradvanced three-dimensional fibre textile composites[J]. Compos Part A，1999，30(12)：1445
[13] Wang Y Q，Wang A S D. Spatial distribution of yarns and mechanical properties in3D braided tubular composites[J]. Appl Compos Mater，1997，4(2)：121
[14] Jan Z，Michal S. Homogenization of balanced plain weave composites withimperfect microstructure：Part I-Theoretical formulation[J]. Int J Solids Struct，2004，41(22/23)：6549
[15] Chiu C H，Tsai K H，Huang W J. Crush-failure modes of2D triaxially braidedhybrid composite tubes[J]. Compos Sci Technol，1999，59：1713
[16] Quek S C，Waas A，Shahwan K W，et al. Compressive response and failure ofbraided textile composites：Part2Computations[J]. Int J Non Linear Mech，2004，39(4)：649
[17]庞宝君，杜善义，韩杰才.三维四向编织复合材料细观组织及分析模型[J].复合材料学报，1999，16(3)：135
[18]马文锁，冯伟.二维编织复合材料几何结构的平面群分析[J].北京科技大学学报.2007，29(2)：226-230
[19]于平，周平等.二维平纹编织复合材料压缩力学行为研究[J].计算力学学报，2008，25(2)：224-228
[20]陈陆平，潘敬哲，钱令希.复合材料纤维/基体界面失效问题的参变量有限元数值模拟[J].复合材料学报，1993，10(1)：381-393
[21]张双寅.复合材料层间剪切强度[J].工程力学，1985，3(4)：26-27
[22]范天佑.断裂理论基础[M].北京：科学技术出版社，2003，76-78
[23]王铎，杜善义，王殿富等.断裂力学[M].黑龙江：哈尔滨工业大学出版社，1989，110-111．
[24]赵艳华，徐世烺. I-II复合型裂纹脆性断裂的最小J2准则[J].工程力学，2002，19(4)：94-98
[25]张磊，谢宗蕻，苏霓.离位增韧复合材料层压板准静态压痕试验研究[J].航天器环境工程，2009，12(26)：128-130
[26]钟明辉.编织复合材料[J].纤维复合材料，1992，2：42-46
[27] ISH IKAWA T，CHOU T W. Stiffness and strength behavior of woven fabriccomposite[J]. Mat Sci，1982，17(11)：32l1-l3220
[28]. ISH IKAWA T，CHOU T W. One two-dimensional micromechanical analysis ofwoven fabric composite[J]. A IAA，1983，21：1714-1721
[29]. ISH IKAWA T， MATSUSH IMA M， HAYASHI Y， et al. Experimentalconfirmation of theory of elastic modulus of fabric composites[J]. CompositeMaterials，1985，19(9)：443-458
[30] HUANG Z M. The mechanical properties of composites reinforced with wovenand braided fabrics[J]. Composites Science and Technology，2000，60：479-498
[31]燕瑛，楼畅，成传贤等.机织复合材料力学性能的细观分析与实验研究[J].复合材料学报，2001，18(2)：109-113
[32]卢子兴，胡奇.三维编织复合材料压缩力学性能的实验研究[J].复合材料学报，2003，20(6)：67-72
[33] GUDEN M，YILDIRIM U，HALL W I. Effect of strainrate on the compressionbehavior of a woven glass fiber/SC-15composite[J]. Polymer Testing，2004，23：719-725
[34] AKIL O，YILDIRIM U，GUDEN M，et al. Effent of strain rate on the compressionbehavior of a woven fabric S2-glass fiber reinforced vinylester composite[J]. PolymerTesting，2003，22：883-887
[35]刘书山.复合材料应力分析的均匀化方法[J].力学学报，1997，29(3)，306-313
[36]刘书山.复合材料弹性性能尺度效应[J].大连理工大学学报，2004，44(2)：200-205
[37]宋迎东，雷友锋.通用单胞法模型预测纤维增强钛基复合材料弹性性能[C].第十届航空发动机结构强度与振动会议征文集.安徽歙县，2000，10：137-141
[38]雷友锋，宋迎东，高德平.纤维增强金属基复合材料弹塑性行为分析[J].机械科学与技术，2003，22(6)：982-986
[39]雷友锋，宋迎东，高德平.纤维增强金书基复合材料细观几何结构对宏观弹性性能的影响[J].宇航材料工艺，2003，33(1)：43-48
[40]宋迎东，雷友锋，孙志刚等.一种新的纤维增强复合材料细观力学模型[J].南京航空航天大学学报，2003，35(4)：435-440
[41]曹礼群，崔俊芝.整周期复合材料弹性结构的有限元计算[J].计算数学，1998，20(3)：279-290
[42]崔俊芝.整周期复合材料弹性结构的双尺度渐近分析方法[J].应用数学学报，1999，21(1)：38-46
[43] Blassiau S，Thionnet A，Bunsell A R. Micromechanics of load transfer in aunidirectional carbon fiber-reinforced epoxy composite due to fibre failures：Part3.Multiscale reconstruction of composite behaviour[J]. Composite Structures，2007，1-12
[44] Naimark O B，Davydova M，Plekhov O A. Nonlinear and structural aspects oftransitions from damage to fracture in composites and structures[J]. Computers andStructures，2000，76：67-75
[45] Tang X D，Whitcomb J D，Kelkar A D. Progressive failure analysis of2×2braidedcomposites exhibiting multiscale heterogeneity[J]. Composites Science andTechnology，2006，66：2580-2590
[46]刘昊，陈宪战，傅德莲. CAE技术现状及其在石油物探机械装备设计中的应用[J].物探装备，2010，20(2)：92-95
[47] Glimm J， Sharp H D. Multiscale science： a challenge for the twenty-firstcentury[J]. SIAM News，1997，30(8)：1~7
[48]曹礼群.材料物性的多尺度关联与数值模拟[J].材料科学研究，2002，24(6)：23-30
[49]沈观林，胡更开.复合材料力学[M].北京：清华大学出版社，2006，3-4
[50]郑晓霞，郑锡涛，缑林虎.多尺度方法在复合材料力学分析中的研究进展[J].力学进展，2010，40(1)：41-47
[51]曹金凤，王旭春，孔亮. Python语言在ABAQUS中的应用[M].北京：机械工业出版社.2011：2-8
[52]石亦平，周玉蓉. ABAQUS有限元分析实例详解[M].北京：机械工业出版社.2006：2-10
[53] Wesley J Chun. Python核心编程(第二版)[M].北京：人民邮电出版社.2008：3-9