三维五向及全五向编织复合材料的三单胞结构模型_张超.pdf 9页

第４５卷第２期南京航空航天大学学报Ｖｏｌ．４５Ｎｏ．２２０１３年４月ＪｏｕｒｎａｌｏｆＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓ＆ＡｓｔｒｏｎａｕｔｉｃｓＡｐｒ．２０１３三维五向及全五向编织复合材料的三单胞结构模型张超许希武严雪（南京航空航天大学机械结构力学及控制国家重点实验室，南京，２１００１６）摘要：针对三维五向及全五向编织复合材料，基于四步法编织工艺携纱器的运动规律，研究了各单胞区域纱线的运动轨迹和空间走向，给出了材料内胞、面胞和角胞的纱线拓扑结构。通过给定一系列的位置节点和节点间线性插值函数的方法来描述纱线空间轨迹。考虑纱线间的相互挤压状况，在假设编织纱和轴纱横截面形状分别为扁平椭圆形和正方形的基础上，系统建立了三维五向及全五向编织复合材料的三单胞参数化结构模型。给出了两种编织材料的纤维体积含量计算公式并进行了对比分析，明确了三维全五向编织复合材料的轴向性能优势。分析了三维五向编织材料的表面多胞结构及三维全五向编织材料的内胞截面图，并与实际试件结构进行对比，结果表明本文所建模型合理准确。关键词：三维五向；三维全五向；编织复合材料；三单胞；结构模型；参数化中图分类号：ＴＢ３３２文献标志码：Ａ文章编号：１００５－２６１５（２０１３）０２－０１７０－０９ＴｈｒｅｅＵｎｉｔ－ＣｅｌｌＳｔｒｕｃｔｕｒｅＭｏｄｅｌｓｏｆ３－ＤＦｉｖｅ－ＤｉｒｅｃｔｉｏｎａｌａｎｄＦｕｌｌＦｉｖｅ－ＤｉｒｅｃｔｉｏｎａｌＢｒａｉｄｅｄＣｏｍｐｏｓｉｔｅｓＺｈａｎｇＣｈａｏ，ＸｕＸｉｗｕ，ＹａｎＸｕｅ（ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙｏｆＭｅｃｈａｎｉｃｓａｎｄＣｏｎｔｒｏｌｏｆＭｅｃｈａｎｉｃａｌＳｔｒｕｃｔｕｒｅｓ，ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓ＆Ａｓｔｒｏｎａｕｔｉｃｓ，Ｎａｎｊｉｎｇ，２１００１６，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｂａｓｅｄｏｎｔｈｅｆｏｕｒ－ｓｔｅｐｂｒａｉｄｉｎｇｐｒｏｃｅｓｓａｎｄｔｈｅｍｏｔｉｏｎｌａｗｏｆｙａｒｎｃａｒｒｉｅｒｓ，ｔｈｅｓｐａｔｉａｌｍｏｖｅ－ｍｅｎｔｔｒａｃｅｓｏｆｂｒａｉｄｉｎｇｙａｒｎｓａｒｅｓｔｕｄｉｅｄｆｏｒ３－Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌａｎｄｆｕｌｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓ－ｉｔｅｓ．Ａｎｄｔｈｅｔｏｐｏｌｏｇｙｇｅｏｍｅｔｒｉｃａｌｍｏｄｅｌｓｆｏｒｔｈｅｉｎｔｅｒｉｏｒ，ｓｕｒｆａｃｅａｎｄｃｏｒｎｅｒｒｅｇｉｏｎｓｏｆｔｈｅｐｒｅｆｏｒｍａｒｅｅｓｔａｂｌｉｓｈｅｄ．Ｔｈｅｓｐａｔｉａｌｙａｒｎｐａｔｈｉｓｆｉｔｔｅｄｂｙａｓｅｒｉｅｓｏｆｐｏｓｉｔｉｏｎｎｏｄｅｓａｎｄｌｉｎｅａｒｉｎｔｅｒｐｏｌａｔｉｏｎｆｕｎｃｔｉｏｎｍｅｔｈｏｄ．Ｃｏｎｓｉｄｅｒｉｎｇｔｈｅｓｑｕｅｅｚｉｎｇｓｉｔｕａｔｉｏｎｂｅｔｗｅｅｎｙａｒｎｓａｎｄｗｉｔｈｔｈｅａｓｓｕｍｐｔｉｏｎｔｈａｔｔｈｅｃｒｏｓｓ－ｓｅｃ－ｔｉｏｎｏｆｂｒａｉｄｉｎｇｙａｒｎｉｓｅｌｌｉｐｔｉｃａｌ，ｔｈｅａｘｉａｌｙａｒｎａｎｄｔｈｅａｘｉａｌｓｔａｔｉｃｙａｒｎａｒｅｓｑｕａｒｅ，ａｎｄｔｈｅｐａｒａｍｅｔｅｒ－ｉｚｅｄｔｈｒｅｅｕｎｉｔ－ｃｅｌｌｓｔｒｕｃｔｕｒｅｍｏｄｅｌｓｏｆ３－Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌａｎｄｆｕｌｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓａｒｅｆｏｕｎｄｅｄ．Ｔｈｅｆｏｒｍｕｌａｓｏｆｆｉｂｅｒｖｏｌｕｍｅｆｒａｃｔｉｏｎａｒｅｇｉｖｅｎａｎｄｔｈｅａｘｉａｌａｄｖａｎｔａｇｅｓｏｆ３－Ｄｆｕｌｌｆｉｖｅ－ｄｉｒｅｃ－ｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓａｒｅｓｈｏｗｎｔｈｒｏｕｇｈａｃｏｍｐａｒａｔｉｖｅａｎａｌｙｓｉｓ．Ｔｈｅｓｕｒｆａｃｅｍｕｌｔｉ－ｃｅｌｌｓｔｒｕｃｔｕｒｅｍｏｄｅｌｏｆ３－Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓａｎｄｔｈｅｃｒｏｓｓ－ｓｅｃｔｉｏｎｏｆｉｎｔｅｒｉｏｒｕｎｉｔ－ｃｅｌｌｍｏｄｅｌｏｆｆｕｌｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓａｒｅａｎａｌｙｚｅｄａｎｄｃｏｍｐａｒｅｄｗｉｔｈｔｈｅａｃｔｕａｌｓｐｅｃｉｍｅｎｓｔｒｕｃｔｕｒｅ，ｗｈｉｃｈｖｅｒｉｆｉｅｄｔｈｅｒａｔｉｏｎａｌｉｔｙａｎｄｅｆｆｅｃｔｉｖｅｎｅｓｓｏｆｔｈｅｓｔｒｕｃｔｕｒｅｍｏｄｅｌｓ．Ｋｅｙｗｏｒｄｓ：３－Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎ；３－Ｄｆｕｌｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎ；ｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓ；ｔｈｒｅｅｕｎｉｔ－ｃｅｌｌ；ｓｔｒｕｃｔｕｒｅｍｏｄｅｌ；ｐａｒａｍｅｔｅｒｉｚｅｄ三维编织复合材料纱束空间多向延伸并相互材料层间强度低、易分层的致命弱点；且具有抗冲交织形成整体网状结构，从根本上克服了层合复合击、抗损伤、耐磨损、缺口欠敏感、近净截面设计制基金项目：国家自然科学基金（１１２７２１４６）资助项目。收稿日期：２０１２－０４－２３；修订日期：２０１２－０８－０７通信作者：许希武，男，教授，博士生导师，Ｅ－ｍａｉｌ：ｘｗｘｕ＠ｎｕａａ．ｅｄｕ．ｃｎ。 第２期张超，等：三维五向及全五向编织复合材料的三单胞结构模型１７１造等诸多层合复合材料不具备的优点，在以航空航复合材料的三单胞结构模型，计算分析了两种编织天为代表的高新技术领域具有十分广阔的应用前材料的纤维体积含量，并以内部编织角和纱线横截景。面尺寸作为输入参数，实现了三单胞结构模型的参三维四向编织复合材料显著改善了材料厚度数化，为基于三单胞模型的材料力学性能有限元分方向的力学性能，但同时其面内力学性能也有明显析奠定了基础。下降。为了强化编织成形方向材料的面内力学性１四步法三维五向及全五向编织能，研究发展了三维五向编织复合材料。三维五向工艺编织复合材料是以三维四向编织工艺为基础，织造过程中在各行编织纱中间沿轴向加入了不参与编四步法三维五向编织工艺以四步法三维四向织的第五向纱（称为轴向纱）而形成的一种新的整编织工艺为基础，在每行相邻编织纱携纱器中间加体编织结构材料。由于轴向纱的引入，提高了材料入轴向纱携纱器；而三维全五向编织工艺则是在三的轴向力学性能，为编织结构材料作为主承力构件维五向编织工艺的基础上，在相邻行之间继续加入提供了基础。但现有的三维五向编织结构中存在一定数量的轴向不动纱携纱器。每一机器循环中，较多无轴纱占据的空位，限制了其轴向力学性能的编织纱携纱器运动四步，且每步运动距离相等，为进一步提升。文献［１］提出了三维全五向编织结构相邻一个编织纱携纱器位置。轴向纱携纱器只沿的概念，通过在三维五向编织预成形件的所有空位行向运动，不沿列向运动。轴向不动纱携纱器在整处增加轴向不动纱，以消除轴向空位，提高材料的个编织过程中均保持不动。纤维体积含量。制成的三维全五向编织复合材料本节仅给出三维全五向编织工艺中携纱器的能够有效地弥补三维四向、五向编织材料轴向力学四步运动规律，对于三维五向编织工艺，只需去除性能的不足。轴向不动纱携纱器，其他携纱器运动规律完全相目前，关于三维五向编织复合材料的研究已经同。如图１所示，○代表编织纱携纱器，×代表轴取得了一定的研究成果，包括几何模型、内胞实体向纱携纱器，代表轴向不动纱携纱器。第一步，［２－５］［６－９］相邻行的编织纱和轴向纱携纱器交错移动一个位模型的建立，弹性性能的预测，强度及损伤［１０－１２］置；第二步，相邻列的编织纱携纱器交错移动一个性能的分析等方面，而对于三维全五向编织复［１］位置；第三步与第一步运动相反；第四步与第二步合材料的研究还处于起步阶段。刘振国分析了现有三维四向、五向编织复合材料力学性能的缺运动相反。一个编织循环后，机器底盘携纱器的排陷，从材料细观结构角度提出了改进途径，并提出列阵形恢复到初始位置，紧接着的“打紧”工序使纱了三维全五向编织复合材料的概念。随后，Ｚｈａｎｇ线间紧密接触，编织结构趋于稳定，预制件沿成形［１３］方向延伸了一个编织花节，长度为ｈ。不断地重复等考虑了纱束间的挤压状况，建立了三维全五［１４］上述的四步运动再加上“打紧”工序，最终编织成具向编织材料的实体内胞模型。最近，皮秀标等分析了三维全五向编织材料的纱线运动规律及空有预订截面形状的预制件。间交织形态，建立了材料的多单胞几何模型。矩形编织预制件中编织纱总数为Ｎ１，轴向纱编织材料的细观结构决定了其宏观力学性能，总数为Ｎ２，轴向不动纱总数为Ｎ３，分别为建立合理准确的细观实体结构模型是对材料进行Ｎ１＝ｍｎ＋ｍ＋ｎ（１）有限元力学分析的基础。至今，关于三维五向及全Ｎ２＝ｍｎ（２）五向编织复合材料的细观实体建模还仅仅局限于Ｎ３＝（ｍ－１）（ｎ－１）（３）内胞结构，而关于材料面胞、角胞的三维建模还未式中：ｍ和ｎ分别为主体纱的行数和列数。见报道，限制了对材料损伤和强度等具有“局部失对于三维五向及全五向编织预制件，纱线总数效特性”的力学问题的深入研究。本文从四步法三分别为Ｎ和Ｎｆ，那么维五向及全五向编织工艺出发，分析了携纱器的运Ｎ＝Ｎ１＋Ｎ２（４）动规律，研究了各单胞区域纱线的面内及空间运动Ｎｆ＝Ｎ１＋Ｎ２＋Ｎ３（５）轨迹，给出了这两种编织材料内胞、面胞和角胞的２纱线运动规律及三单胞几何模型纱线拓扑结构。在假设编织纱线和轴纱线横截面分别为扁平椭圆形和正方形的基础上，以空间样条２．１纱线的面内运动规律曲线描述纱线轨迹，建立了三维五向及全五向编织编织过程中，携纱器牵引着编织纱线和轴向纱 １７２南京航空航天大学学报第４５卷图２（ｂ）为所有纱线的Ｘ－Ｙ面内运动轨迹，◇表示轴向纱，表示轴向不动纱，实线表示编织纱，轴向纱和轴向不动纱均匀地分布于编织纱中间。对于三维五向编织结构，表示的轴向不动纱并不存在，形成无轴纱占据的空位。图２纱线面内运动规律２．２单胞选取本文采用三单胞划分思想，认为三维五向及全图１三维四步法编织工艺过程五向编织复合材料由周期性的内胞、面胞和角胞构成，三单胞的取法如图２（ｂ）所示。三类单胞的取线在机器底盘上做周期性的交错运动，轴向不动纱［４］向均平行于试件横截面边界，在后续基于三单胞携纱器保持不动。编织纱携纱器总体上沿“Ｚ”字模型进行有限元力学分析时，这种取法有利于边界形连续运动，穿越主体纱区域，抵达边界后停动一条件的施加和单胞应力应变场的分析。需要指出步，再返回主体纱区域。经过若干个循环后，携纱的是，只有ｍ和ｎ同时为偶数时，这种单胞取法才器将遍历所有边界，回到起始位置。能完整描述织物的整体结构。如图２（ａ）所示，“打紧”工序后在纱线张力作２．３纱线的空间运动规律及三单胞几何模型用下，内部编织纱线处于拉伸状态而保持直线走三维五向及全五向编织工艺的一个机器循环向，纱线在水平面内的投影为编织纱携纱器运动过包含两个运动步骤，即位于Ｘ－Ｙ平面内携纱器的程中相邻位置点中点的连线；对于表面区域编织纱交错四步运动和Ｚ方向“打紧”运动。对于所有携线，在外表面停动一步（Ｊ－Ｊ′）后再返回内部区域，纱器，每运动一步，纱线在空间均沿Ｚ方向延伸纱线面内运动轨迹为直折线段；对于角部区域纱ｈ／４高度。线，在外角部区域需停动两次（Ｎ－Ｎ′，Ｏ－Ｏ′），纱线２．３．１内部纱线空间运动规律及内胞纱线拓扑结面内运动轨迹为直双折线段。对于轴向纱线运动，构以携纱器Ｑ为例，携纱器始终只在ＱＲ内运动，ＱＲ如图３（ａ）所示，内部编织纱携纱器Ａ由点Ａ中点Ｓ即为轴向纱在水平面内的投影；轴向不动运动到点Ｅ经历四步，经“打紧”工序后，纱线空间纱在编织过程中并不运动，其在水平面内投影为携轨迹为直线，位于与试件横截面成±４５°的两组正纱器所在的位置点。交平面内，且与Ｚ方向夹角为γ，定义为内部纱线 第２期张超，等：三维五向及全五向编织复合材料的三单胞结构模型１７３编织角。如图３（ｂ）所示，轴向纱携纱器Ｑ的四步运动始终只在ＱＲ范围内，并不参与编织，经“打紧”工序后，轴向纱空间轨迹为平行于Ｚ轴的直线，并通过轴向纱携纱器位置点Ｑ和Ｒ的中点。轴向不动纱携纱器在整个编织过程中并不运动，轴向不动纱在空间保持平行于Ｚ方向的直线状态延伸。轴向纱和轴向不动纱空间轨迹在面胞、角胞区域与内胞区域完全相同。三维全五向编织材料的内胞纱线拓扑结构如图３（ｃ）所示，对于三维五向编织材料则没有轴向不动纱所代表的拓扑纱线。图４表面编织纱线空间运动规律和面胞几何模型图３内部纱线空间运动规律和内胞几何模型２．３．２表面编织纱线空间运动规律及面胞纱线拓扑结构图５角部编织纱线空间运动规律和角胞几何模型如图４（ａ）所示，表面编织纱携纱器Ｆ经历五步运动由内部区域进入表面后再返回内部区域，携与其他区域纱线自然连接过渡。直线段ＷＸ与Ｚ纱器位置点分别为Ｆ，Ｉ，Ｊ，Ｊ′，Ｉ′和Ｋ，其中Ｊ－Ｊ′停方向夹角为β，定义为角部纱线编织角。三维全五动一步，这导致表面纱线的空间拓扑结构与内部纱向编织复合材料的角胞纱线拓扑结构如图５（ｂ）所线不同。经“打紧”工序后，假设表面纱线空间轨迹示，对于三维五向材料则没有轴向不动纱所代表的［１１］为两条直线段ＴＵ和ＵＶ组成，并与其他区域拓扑纱线。纱线自然连接过渡。直线段ＴＵ与Ｚ方向夹角为３三单胞结构模型θ，定义为表面纱线编织角。三维全五向编织材料的面胞纱线拓扑结构如图４（ｂ）所示，对于三维五３．１基本假设向编织复合材料则没有轴向不动纱所代表的拓扑三维编织复合材料实际细观结构非常复杂，本纱线。文中三维五向及全五向编织复合材料三单胞结构２．３．３角部编织纱线空间运动规律及角胞纱线拓模型的建立基于如下基本假设：扑结构（１）内部编织纱线空间沿直线延伸，考虑到其如图５（ａ）所示，角部编织纱携纱器经历８步同时受到编织纱、轴向纱及轴向不动纱（对于全五由点Ｌ运动到点Ｐ，纱线由内部区域进入角部区向编织结构）的挤压，假设其截面形状为扁平椭圆域后再返回内部区域，其中Ｎ－Ｎ′，Ｏ－Ｏ′各停动一形，长半轴为ａ，短半轴为ｂ。步，这导致角部纱线的空间拓扑结构与内部、表面（２）表面及角部编织纱线空间由直线段组成，纱线均不相同。在“打紧”工序后，假设角部纱线空纱束截面形状也是扁平椭圆形，考虑到外部编织纱［１１］间轨迹为３条直线段ＷＸ，ＸＹ和ＹＺ组成，并束与内部编织纱束挤压状态的差异，引入截面形状 １７４南京航空航天大学学报第４５卷缩小比例系数φｆ和φｃ。角胞厚度Ｔｃ＝２槡（２＋ｒ）ｂ＋槡φｃｂ（１５）（３）轴向纱和轴向不动纱空间为平行于Ｚ轴角部编织角的直线走向，截面形状均为正方形，两组平行对边ｔａｎβ＝（２＋ｒ）ｂ／２（（３／８）ｈ）＝ｔａｎγ／３（１６）与试件横截面边界成±４５°，边长分别ｒｂ，ｒｆｂ，ｒ，ｒｆ纱线编织角是三维编织复合材料重要的结构为截面尺寸因子。参数。对比分析三单胞模型各区域的纱线编织角，（４）表面和角部区域编织纱线向内部延伸过程由式（１３，１６）看出，表面编织角和角部编织角两者中的横截面渐进变化。相等，小于内部编织角。这说明外部编织纱线空间（５）编织工艺相对稳定，一定长度范围内保证构型与内部纱线不同，外部编织纱线更靠近于材料编织结构均匀一致。轴向，这导致外部纱线的挤压比内部区域严重。内部编织纱、轴向纱及轴向不动纱的截面形状３．３纤维体积含量如图６所示，基于上述基本假设，认为“打紧”工序以下先计算三维五向编织材料三单胞结构的后纱线紧密接触，则有纤维体积含量：［４］２（６）ａ＝ｂｃｏｓγ槡ｒ＋４ｒ＋３对于内部单胞，有ｒ＝ｒ（２＋ｒ）２２ｆ（７）内胞体积Ｕｉ＝８ｂｈ（１７）编织纱体积Ｙｉｂ＝４πａｂｈ／ｃｏｓγ（１８）轴纱体积Ｙ２２ｉａ＝４ｒｂｈ（１９）纤维体积含量Ｖｉｆ＝（Ｙｉｂεｉ＋Ｙｉａεａ）／Ｕｉ（２０）对于表面单胞，有面胞体积２Ｕｆ＝槡２（２＋ｒ）（２槡２＋槡２ｒ＋槡φｆ）ｂｈ（２１）编织纱体积３πａｂｈ１＋φｆφｆ（２２）图６纱线横截面形状Ｙｓｂ＝＋２（２ｃｏｓγｃｏｓθ）２２轴纱体积Ｙｓａ＝２ｒｂｈ（２３）３．２三单胞结构参数纤维体积含量三维全五向编织结构是在三维五向编织结构３πａｂｈεｉ＋φｆεｆφｆεｆ的基础上在结构空位处沿Ｚ方向引入轴纱，基于＋＋Ｙｓａεａ２（２ｃｏｓγｃｏｓθ）Ｖｆｆ＝（２４）本文假设，轴纱的引入并不改变五向编织材料单胞Ｕｆ的结构参数，所以这两种编织材料的结构参数表达对于角部单胞，有式相同。根据结构建模的基本假设，结合图３～５角胞体积所示的各区域单胞的纱线拓扑几何特征，三单胞实２２Ｕｃ＝（２槡２＋槡２ｒ＋槡φｆ）ｂｈ（２５）体模型的结构参数如下：编织纱体积对于内胞，结构参数为πａｂｈ１＋φｃ３φｃＹｃ＝＋（２６）内胞宽度Ｗｉ＝２槡２（２＋ｒ）ｂ（８）２（２ｃｏｓγｃｏｓβ）２２内胞厚度Ｔｉ＝２槡２（２＋ｒ）ｂ（９）轴纱体积Ｙｃａ＝ｒｂｈ（２７）内胞高度ｈ＝４（２＋ｒ）ｂ／ｔａｎγ（１０）纤维体积含量对于面胞，结构参数为πａｂｈεｉ＋φｃεｃ３φｃεｃ＋＋Ｙｃａεａ２（２ｃｏｓγｃｏｓβ）面胞宽度Ｗｆ＝２槡２（２＋ｒ）ｂ（１１）Ｖｃｆ＝（２８）Ｕｃ面胞厚度Ｔｆ＝２槡（２＋ｒ）ｂ＋槡φｆｂ（１２）在各区域单胞纤维体积含量的表达式中：εｉ，εｆ和εｃ表面编织角分别为内部、表面及角部编织纱线的纱线填充系ｔａｎθ＝（２＋ｒ）ｂ／２（（３／８）ｈ）＝ｔａｎγ／３（１３）数；εａ为轴纱的纱线填充系数。对于角胞，结构参数为在三维全五向编织材料内胞结构中，含４根轴角胞宽度Ｗｃ＝２槡（２＋ｒ）ｂ＋槡φｃｂ（１４）向纱，也含有４根完整的轴向不动纱；在面胞结构 第２期张超，等：三维五向及全五向编织复合材料的三单胞结构模型１７５中，含有两根轴向纱，并含有１根完整的轴向不动设区间［ａ，ｂ］被ｎ＋１个位置节点ａ＝ｔ０＜ｔ１＜纱；在角胞结构中，含有１根轴向纱，并含有１／４根…＜ｔｎ－１＜ｔｎ＝ｂ分割，向量函数Ｓ在ｎ个分割区间完整的轴向不动纱。假设轴向纱和轴向不动纱的［ｔｉ，ｔｉ＋１］（ｉ＝０，１，…，ｎ－１）上均是ｔ的函数，其表纱线填充系数相同，则对于三维全五向编织材料各达式为单胞结构纤维体积含量的计算，只需在三维五向编烄Ｓ０（ｔ）ｔ０≤ｔ≤ｔ１织材料计算公式（１７～２８）的基础上，将式（１９）的系Ｓ１（ｔ）ｔ１≤ｔ≤ｔ２Ｓ（ｔ）＝烅（３０）数由４改为８，式（２３）的系数由２改为３，式（２７）的系数由１改为１．２５，其他计算公式完全相同。烆Ｓｎ－１（ｔ）ｔｎ－１≤ｔ≤ｔｎ三维五向及全五向编织复合材料总纤维体积在所有的分割区间上，插值函数均采用线性形式含量为各区域单胞纤维体积含量与其占整个结构Ｓｉ（ｔ）＝ａｉ＋ｂｉｔ（３１）百分比之积的和，即每个分割区间上Ｓｉ（ｔ）的表达式都有２个未知Ｖ＝ＶｉＶｉｆ＋ＶｆＶｆｆ＋ＶｃＶｃｆ（２９）系数，ｎ个区间上共２ｎ个未知系数需要确定．向量式中，Ｖｉ，Ｖｆ，Ｖｃ分别为内部、表面和角部各区域占函数Ｓ必须通过所有的ｎ＋１个节点，而节点的位整个复合材料结构的体积分数。置是给定的，这便提供了２ｎ个方程，函数Ｓ就唯一在携纱器排列为６×６编织方式且纤维束截面确定。尺寸确定的条件下，本文给出了两种典型编织角的３．５三单胞实体结构模型五向及全五向材料三单胞和整体结构的纤维体积目前，对三维编织复合材料的力学性能进行有含量，如表１所示。对于五向编织结构，内胞纤维限元分析预测时，研究者往往认为材料尺寸足够体积含量最少，而对于全五向编织结构，内胞纤维大，内胞区域占材料的主体从而忽略面胞、角胞的体积含量最大。相对于三维五向编织材料，由于轴影响［８－１２］。对于表现复合材料混合效应的刚度预向不动纱的引入，全五向编织材料的纤维体积含量测，这种假设还不致产生较大的误差，但材料强度、有较大提升，达到６０％以上。当携纱器排列增加、损伤等问题具有“局部失效特性”，即材料失效起始材料截面尺寸增大时，内胞区域占整个结构比例将于局部（尤其是外表面），这时若忽略材料的面胞和大幅上升，三维全五向编织材料的轴向力学性能优角胞，仅仅基于周期性边界条件，研究内胞响应来势将更为显著。代表整个复合材料结构，其前提假设并不合理，所表１典型编织角下三单胞模型及整体结构的纤维体积含量得结果也不精确。故需要建立合理准确的三单胞％实体模型，以深入研究材料的力学性能。五向编织全五向编织建立实体模型时，通过给定纱线的截面形状，纤维体积含量结合纱线的空间轨迹来确定各个纱线的实体构型，γ＝１５°γ＝３５°γ＝１５°γ＝３５°在所有纱线的空间组合体中取单胞空间尺寸的立内胞４４．４７４４．４７６２．４９６２．４９方体，立方体内部所围的纱线结构即是单胞的实体面胞５０．９２４８．６４５８．９５５６．６７模型。三维五向及全五向编织复合材料的三单胞角胞５７．２０５３．１４６１．２３５７．１６实体模型如图７所示，考虑到编织纱束在表面、角整体结构４９．０８４７．５１６０．６９５９．１２部区域与内部区域挤压状态的差异，两种编织材料三单胞结构模型建立时，表面和角部截面形状缩小３．４纱线空间轨迹描述系数φｆ和φｃ均设为０．９。在靠近内胞一侧的位置根据本文基本假设，经“打紧”工序后，编织纱节点上纱线截面取内胞截面的尺寸，而在面胞、角胞线轨迹在内部区域为空间直线，在表面及角部区域内部的位置节点上，纱线截面取乘以了截面缩小系为空间分段直线，轴向纱及轴向不动纱轨迹为平行数后的值，这样，在面胞、角胞编织纱线的整个空间于Ｚ轴的直线。实体建模时，将纱线空间轨迹看轨迹上，纱线截面尺寸在各个位置节点上是不同的。作是通过纱线横截面中心点的空间向量Ｓ函数，其实体成形时，软件会自动插值各位置节点处的截面值是沿纱线中心线方向距离ｔ的函数，通过给定一形状，以保证编织纱线由外部区域向内部区域延伸系列的位置节点和节点间插值函数来确定纱线过程中横截面渐进变化。轴向纱和轴向不动纱的截轨迹。面形状和大小在外部区域与内部区域保持不变。 １７６南京航空航天大学学报第４５卷图７三单胞的实体结构模型４参数化结构建模三维编织复合材料的细观结构复杂，建立其三维实体模型时结构参数众多。若针对每一特定结构参数都建立相应的实体模型，准备工作量过大。本文通过建立三单胞模型各参数间的关联约束方程，实现了三单胞结构模型的参数化。４．１单胞实体模型随输入参数的变化经分析三单胞模型各结构参数间关系，三维五向及全五向参数化单胞模型的输入参数设定为内部编织角γ、编织纱线椭圆截面短半轴长ｂ、轴纱截面尺寸因子ｒ及编织纱截面形状缩小系数φｆ和，三单胞模型的其他结构参数均可由输入参数结φｃ合关联约束方程确定。当任一输入参数发生改变时，参数化的三单胞实体模型均会即时作出相应更新，从而实现快速建模。图８三维全五向参数化三单胞模型本节仅给出其他输入参数确定时，三维全五向编织结构内胞、面胞及角胞模型随内部编织角γ的ｔａｎα＝Ｗｆ／ｈ＝（槡２／２）ｔａｎγ（３２）变化情况，如图８所示。图１０给出了本文所建的三维五向编织复合材料多面胞结构模型与实际预制件结构表面照片的４．２表面、内部单胞结构模型分析对于三维编织复合材料试件，测量内部编织角对比情况。可以看出，多面胞模型准确有效地反映比较困难，实际应用中均是测量材料编织角α（简了材料的表面纱线结构特征。称编织角），定义为试件表面相邻同向编织纱线排图１１为三维全五向编织复合材料试件横截面列形成的表面纱线纹路与编织成形方向Ｚ的夹图，轴向纱及轴向不动纱截面形状近似为正方形，角，如图９所示。编织纱近似为椭圆形。不存在没有被轴纱占据的 第２期张超，等：三维五向及全五向编织复合材料的三单胞结构模型１７７图１２内胞模型的截面图图９编织预成形件表面纱线拓扑模型中均未考虑纱线的空间弯曲、扭转等细节特征，建立更加准确的三维五向编织材料的单胞模型还需做更加细致和深入的工作。５结论（１）基于三维五向及全五向编织工艺携纱器运动规律，分析了纱线的面内及空间运动轨迹，建立了两种材料内胞、面胞和角胞的纱线拓扑结构。图１０多面胞模型与预制件表面照片的对比（２）通过给定位置节点和节点间线性插值函数的方法来拟合纱线空间轨迹，在假设编织纱线和轴空位，纱线结构紧密稳定且未出现明显弯曲变形。纱线横截面分别为扁平椭圆形和正方形的基础上，从中取出一个内胞，与图１２（ａ）中本文所建三维全系统建立了两种材料的三单胞参数化结构模型。五向编织结构内胞模型截面图进行对比，两者非常（３）相对于三维五向编织复合材料，全五向材接近。因为缺少角胞纱线结构的实验观测结果，这料的纤维体积含量得到较大提升，大幅增强了材料里角胞模型未与实验进行对比。从以上已有对比的轴向力学性能分析可以看出，本文模型在纱束截面形状、空间构（４）所建参数化模型在纱束截面形状、空间构型及交织状况等方面均与实际情况非常吻合，说明型及交织状况等方面均与实际情况非常吻合，说明了所建三单胞结构模型的合理准确性。了三单胞结构模型的合理准确性。文献［１］指出，三维五向编织材料结构空位的（５）三维五向编织结构由于轴向空位的存在引存在引起轴向纱和编织纱横截面形状的改变，纱束起纱线在截面和中心线方向产生弯曲变形。全五内纤维丝产生弯曲，纱束的中心线也会产生弯曲。向编织复合材料不存在空位结构，使编织结构更加三维五向编织结构中的轴向纱截面多为三角形和紧密，纱线间的挤压更均匀。［１５］扇形相互过渡变化。编织纱线和轴向纱线的弯曲、扭转会对材料的整体力学性能产生一定的影参考文献：响，但在本文所建的三维五向内胞结构模型（图１２［１］刘振国．三维全五向编织预型件的概念［Ｊ］．材料工（ｂ））中，轴向纱截面为正方形，且轴向纱和编织纱程，２００８（Ｓ１）：３０５－３１２．空间保持直线走向，即模型中没有引入纱线的截面ＬｉｕＺｈｅｎｇｕｏ．Ｃｏｎｃｅｐｔｏｆｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌａｌｌｆｉｖｅ－变形和中心线弯曲。且目前已有的实体建模研究ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｐｅｒｆｏｒｍｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＭａｔｅｒｉａｌｓＥｎｇｉｎｅｅｒｉｎｇ，２００８（Ｓ１）：３０５－３１２．［２］ＷｕＤＬ．Ｔｈｒｅｅ－ｃｅｌｌｍｏｄｅｌａｎｄ５Ｄｂｒａｉｄｅｄｓｔｒｕｃｔｕｒａｌｃｏｍｐｏｓｉｔｅｓ［Ｊ］．ＣｏｍｐｏｓｉｔｅｓＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，１９９６，５６（３）：２２５－２３３．［３］徐正亚．三维五向编织复合材料细观结构及工艺分析［Ｄ］．天津：天津工业大学，２００７．［４］徐焜，许希武．三维五向矩形编织复合材料的细观结［１４］构模型［Ｊ］．南京航空航天大学学报，２００８，４０（２）：图１１三维全五向编织复合材料纱线截面形状 １７８南京航空航天大学学报第４５卷１６８－１７３．ｓｏｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｏｆｔｈｒｅｅ－ａｘｉａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅＸｕＫｕｎ，ＸｕＸｉｗｕ．Ｍｉｃｒｏｓｔｒｕｃｔｕｒｅｍｏｄｅｌｏｆ３－Ｄ５－ｍａｔｅｒｉａｌｓ［Ｊ］．ＣｏｍｐｏｓｉｔｅｓＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，ｄｉｒｅｃｔｉｏｎａｌｒｅｃｔａｎｇｕｌａｒｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．Ｊｏｕｒ－２００６，６６（１５）：２９５４－２９６４．ｎａｌｏｆＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏ－［１１］徐焜．三维多向编织复合材料力学性能及渐进损伤研ｎａｕｔｉｃｓ，２００８，４０（２）：１６８－１７３．究［Ｄ］．南京：南京航空航天大学，２００８．［５］李金超，张一帆，孙菲，等．三维五向编织复合材料的［１２］张超，许希武，郭树祥．含界面脱粘三维五向编织复合力学性能分析Ⅰ：细观结构模型［Ｊ］．复合材料学报，材料单向拉伸损伤失效机理研究［Ｊ］．航空材料学报，２００９，２６（１）：１５０－１５５．２０１１，３１（６）：７３－８０．ＬｉＪｉｎｃｈａｏ，ＺｈａｎｇＹｉｆａｎ，ＳｕｎＦｅｉ，ｅｔａｌ．ＭｅｃｈａｎｉｃｓＺｈａｎｇＣｈａｏ，ＸｕＸｉｗｕ，ＧｕｏＳｈｕｘｉａｎｇ．Ｄａｍａｇｅａｎｄａｎａｌｙｓｉｓｏｆｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｆａｉｌｕｒｅｍｅｃｈａｎｉｓｍａｎａｌｙｓｉｓｏｆ３Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉ－ｃｏｍｐｏｓｉｔｅｓⅠ：Ｍｉｃｒｏｓｔｒｕｃｔｕｒａｌｍｏｄｅｌ［Ｊ］．ＡｃｔａＭａ－ｄｅｄｃｏｍｐｏｓｉｔｅｓｗｉｔｈｉｎｔｅｒｆａｃｅｄｅｂｏｎｄｉｎｇｕｎｄｅｒｕｎｉｄｉ－ｔｅｒｉａｅＣｏｍｐｏｓｉｔａｅＳｉｎｉｃａ，２００９，２６（１）：１５０－１５５．ｒｅｃｔｉｏｎａｌｔｅｎｓｉｏｎ［Ｊ］．ＪｏｕｒｎａｌｏｆＡｅｒｏｎａｕｔｉｃａｌＭａｔｅｒｉ－［６］ＬｉＤＳ，ＬｕＺＸ，ＣｈｅｎＬ，ｅｔａｌ．Ｍｉｃｒｏｓｔｒｕｃｔｕｒｅａｎｄａｌｓ，２０１１，３１（６）：７３－８０．ｍｅｃｈａｎｉｃａｌｐｒｏｐｅｒｔｉｅｓｏｆ３Ｄａｎｄ５－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄ［１３］ＺｈａｎｇＦ，ＬｉｕＺＧ，ＷｕＺ，ｅｔａｌ．Ａｎｅｗｓｃｈｅｍｅａｎｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌｏｆＳｏｌｉｄｓａｎｄｍｉｃｒｏｓｔｒｕｃｔｕｒａｌｍｏｄｅｌｆｏｒ３Ｄｆｕｌｌ５－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉ－Ｓｔｒｕｃｔｕｒｅｓ，２００９，４６（１７／１８）：３４２２－３４３２．ｄｅｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＡｅｒｏｎａｕｔｉｃｓ，［７］卢子兴，刘子仙．三维五向编织复合材料的弹性性能２０１０，２３（１）：６１－６７．［Ｊ］．北京航空航天大学学报，２００６，３２（４）：４５５－４６０．［１４］皮秀标，钱坤，曹海建，等．三维全五向编织复合材料ＬｕＺｉｘｉｎｇ，ＬｉｕＺｉｘｉａｎ．Ｅｌａｓｔｉｃｐｒｏｐｅｒｔｉｅｓｆｏｒ３ｄｉ－的细观结构分析［Ｊ］．宇航材料工艺，２０１１，４１（６）：３９－ｍｅｎｓｉｏｎ５ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．Ｊｏｕｒｎａｌ４３．ｏｆＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕ－ＰｉＸｉｕｂｉａｏ，ＱｉａｎＫｕｎ，ＣａｏＨａｉｊｉａｎ，ｅｔａｌ．Ｍｉｃｒｏ－ｔｉｃｓ，２００６，３２（４）：４５５－４６０．ｓｔｒｕｃｔｕｒｅａｎａｌｙｓｉｓｏｆ３Ｄｆｕｌｌ５－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄ［８］ＬｉＤＳ，ＦａｎｇＤＮ，ＪｉａｎｇＮ，ｅｔａｌ．Ｆｉｎｉｔｅｅｌｅｍｅｎｔｃｏｍｐｏｓｉｔｅｓ［Ｊ］．ＡｅｒｏｓｐａｃｅＭａｔｅｒｉａｌｓａｎｄＴｅｃｈｎｏｌｏ－ｍｏｄｅｌｉｎｇｏｆｍｅｃｈａｎｉｃａｌｐｒｏｐｅｒｔｉｅｓｏｆ３Ｄｆｉｖｅ－ｄｉｒｅｃ－ｇｙ，２０１１，４１（６）：３９－４３．ｔｉｏｎａｌｒｅｃｔａｎｇｕｌａｒｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．Ｃｏｍｐｏｓ－［１５］陈利，徐正亚．三维五向编织复合材料中纱线截面形ｉｔｅｓ：ＰａｒｔＢ，２０１１，４２（６）：１３７３－１３８５．状实验分析［Ｊ］．复合材料学报，２００７，２４（４）：１２８－［９］ＸｕＫ，ＸｕＸＷ．Ｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｏｆｍｅｃｈａｎｉｃａｌ１３２．ｐｒｏｐｅｒｔｉｅｓｏｆ３Ｄｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓＣｈｅｎＬｉ，ＸｕＺｈｅｎｇｙａ．Ｅｘｐｅｒｉｍｅｎｔａｌａｎａｌｙｓｉｓｏｎｔｈｅ［Ｊ］．ＭａｔｅｒｉａｌｓＳｃｉｅｎｃｅａｎｄＥｎｇｉｎｅｅｒｉｎｇＡ，２００８，４８７ｙａｒｎ′ｓｓｅｃｔｉｏｎｏｆｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｆｉｖｅ－ｄｉｒｅｃｔｉｏｎａｌ（１／２）：４９９－５０９ｂｒａｉｄｅｄｃｏｍｐｏｓｉｔｅｓ［Ｊ］．ＡｃｔａＭａｔｅｒｉａｅＣｏｍｐｏｓｉｔａｅ［１０］ＭｉｒａｖｅｔｅＡ，ＢｉｅｌｓａＪＭ，ＣｈｉｍｉｎｅｌｌｉＡ，ｅｔａｌ．３Ｄｍｅ－Ｓｉｎｉｃａ，２００７，２４（４）：１２８－１３２．