《研究复合材料三维编织预成型件(11页).doc》由会员分享，可在线阅读，更多相关《研究复合材料三维编织预成型件(11页).doc（10页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-研究复合材料三维编织预成型件-第 10 页研究复合材料三维编织预成型件内/间层的剪切变形先进纺织复合材料的教育部重点实验室，复合材料研究所，天津工业大学，天津300387，中华人民共和国纺织和制衣部门，生物与农业工程系，美国加州大学戴维斯分校【摘要】这项研究提出了具有不同的面料密度的三维角联锁预成型件的面内剪切和层间剪切行为。对三维织物预成型件进行画框剪切试验，分析了剪切应力与剪切角度的非线性曲线和变形机理。设计了一个新的测试方法来确定的层间剪切性能表征。经过层间剪切试验后的样本，通过调查拉出的纱线和中间结构发现变形和破坏机制。结果表明织物密度对三维联锁预成型件面内剪切和层间剪切性能有重要的
2、影响，并且织物密度的增加，剪切行为减小。织物密度越低，可变形性越好。层间剪切破坏模式是从织物上引出的粘结剂纱线。希望该研究可以为建立理论模型提供试验基础。连续纤维增强复合树脂基材料引起了很多重视，这都是由于它们所具有的优势，例如高性能，加工周期短，维修和焊接的可能性14。虽然层压复合材料具有优异的面内力学性能，但是层间复合材料的应用范围因厚度受到限制，这是由于差的层间性质。三维纺织结构复合材料具有厚度优势，好的破坏误差和有利的影响，抗疲劳优点59。作为三维纺织结构加强的复合材料的一种，三维角联锁织物已被广泛地应用于工程领域，归于它在传统织机简单和有效地加工1012。另外，三维角联锁织物最吸引的
3、优点是具有近终成形能力的制造复合材料13。三维角联锁预成型件有卓越的机械性能和好的可成形性（图1）。随着预成型技术的发展，可以生产出形状复杂和不同尺寸的结构件。图1三维角联锁织物的半球成形在复合材料生产的结构集成制造中，三维角联锁是根据最终复合材料产品形状预成型，该形状可以是复杂的 1416。对于三维织物，平面内的行为和层间的行为是最重要的变形，并且剪切 行为材料变形的主要模式1719。研究三维角联锁织物层内和层间的剪切行为是有价值的，因为在生产中它们被广泛应用，尤其是成形工艺。二维织物的面内剪切行为已有比较好的研究。Zhu等人20,21通过实验测试仔细研究了二维织物面内表征特性，并且发现了发
4、现纱线的减少是起皱的一个关键。Hivet等 22,23使用相框测试方法研究了二维织物剪切性能，并指出在实验过程中，剪切结果对纱线的张力敏感。拉伸力随剪切角增大而增大。Lomov等24,25通过相框试验提出对在三种不同预张力状态的非平衡2/2斜纹玻璃/ PP织物的剪切测试，并且研究在纱线方向的拉力载荷对织物的剪切抗力的影响及可重复性的方法。基于二维Lin等人 26建立了织物的几何形状来模拟的面内剪切的有限元模型，仿真结果与实验相同。Cao等。 27 比较了相框剪切试验结果，这些结果来于用于制定标准的测试设置获得准确的和适当的材料特性的七个不同实验室的。Chen等28 开发了有限元模型来预测层压板
5、面内和层间的剪切性能。然而，三维角联锁织物的面内/层间的剪切行为很少有报道。Charmetant等 29建立了半球模型来仿真三维织物成形。在这篇文章，报道了一份仔细研究了关于不同织物密度的三维角联锁织物面内剪切和层间剪切行为。记录了剪切应力和剪切角度曲线和面内剪切测试的起皱位置，并且它们相互比较，分析了内层剪切试验的应力-位移曲线。另外，面内剪切非线性曲线的三个阶段被表征。呈现层间剪切破坏形态并且被比较从而在剪切测试过程中获得三维机织物的结构效应。它可以为研究成形性和理论分析提供基础。图2是三维角联锁织物的图片和半球横截面。织物样本是用玻璃纤维细丝制成的。它们的规格列在表1。三维角联锁织物结构
6、的草图示意图如图3 来图解说明具体的预成型结构。对于这种架构， 由三组不同的纱线系统组成，即经纱， 纬纱和粘合纱，它提出了一个层到层的角联锁 结构，其中经纱和纬纱都几乎直，粘合纱显示出不同的起伏，连接卷曲的纬纱层的上层和下层，使它们接合，形成一个稳定的织物结构10。非卷曲 经纱和纬纱敷设在0/90序列里，彼此无编织。粘合纱的线密度小于 该经线和纬线的线密度，它只是在预成型体起到了部分的连接作用。这种结构特征保证了沿经线和纬线方向的高刚度和高强度。图2三维角联锁变形照片；（a）是表面（b）是横截面表1三维角联锁变形的结构参数图3三维角联锁变形的示意图相框测试是一种有效的方式表征 织物的内部层剪切
7、性能。三维织物的剪切试验 根据相框的大小，样本的特征如图4所描述。为了防止大的变形时由夹具施加在预成型件的压力，剪切变形的中心区域是100mmX100mm，且四角落部份被切掉。剪切试验在岛津1kNE万能试验机器上以10 mm / min的十字头速度进行的。三个样本每个结构都被测试。在实验过程中，相框试验被改善是为了确保纯剪切载荷，例如，以最小化边缘约束和促进织物的夹持，等等30-34。图4用于剪切试验的织物样本的示意图问题：（1） 图像帧的臂由滚动轴承连接而不是轴连接。框架之间的摩擦是剪切的结果的一个重要因素，由于在剪切的初始阶段，纱线之间小的剪切力。之前摩擦 - 位移曲线和改进后的摩擦 -
8、位移曲线如图5所示。可以看出，框架的摩擦保持稳定且小，改善后的值在0.040.08 N范围内，摩擦不会影响试验结果。图5改善前后摩擦力的比较（2） 纱线的定位得到了提高。剪切结果和样本的定位高度相关，尤其是对带有较低卷曲的织物。如图6所示，定位不当，在剪切力测试时，可引起大的纱线张力纱线或者纱线松弛。因此，样品中的纱线应与框架平行。四个框架子啊中央位置开槽（如图7），槽的宽度等于 ,夹持部件中织物的宽度。当织物放在槽中，纱线不能移动。没有一个纱线出现偏差，纯剪切载荷可以得到。图6剪切变形后处于两个不同位置的纱线状态图7相框法的示意图（3） 一个新型的板是用来防止在测试过程中织物滑动。不同类型的
9、板用于相框试验中，最佳的折中是在粘合橡胶（如图8）的铝板。用超级胶水将两个板粘合到织物的每一面，然后将其拧入相框。由于它良好的弹性，橡胶层使纱线无滑动。这种技术可以导致在应变区域内好均匀性，对织物行为有轻微的影响。这些板也可以用于其它的测试。图8板的示意图层间剪切行为是重要的变形方式，尤其是对于三维织物37。设计一个新的测试装置来表征三维角联锁织物的跨层剪切性能，根据美国材料与试验协会：C273标准。在拉伸试验机的加载装置的示意图如图9所示。这是一台通用的岛津1kNE测试机器。图10给出层间剪切测试的样本图。测试的样本尺寸是分别沿经向和纬纱方向是 50 mm X40 mm，并且十字头速度为1m
10、m/min。首先用超级胶水将该样本被胶合在尺寸大小为80 mm X 60 mm的铝板上，两板彼此之间必须对准，用于固定在设备上。然后，胶合后的样本插入到测试装置上的槽内。最后，凹槽被拧在设备上。图9跨层试验装置的示意图图10跨层样本图（a）织物被胶于两板之间和（b）铝板该试验装置由两部分组成，分别是顶部部分和底部部分。这两部分分别连接样本的左边和右边（图9）。这个测试装置装有两个传感器（位移传感器和力传感器）。力和位移数据可以通过一台电脑获得。下面的方程是用来计算剪切强度： （1）其中P是剪切载荷，l和b是样本长度和宽度。在冲压制造中，三维预型件和二维织物的面内剪切变形仍然是主要的变形模式。剪
11、切变形是受到局部起皱的限制，直到纱线达到所谓的锁定角3840。基于夹具的变形构造，框架的剪切角由下面的公式来计算： （2） （3）其中Lf是框架的长度，d是十字头的位移 ，和分别是目前的角度和框架的剪切角度。相框测试采用三维预型件，章节2.2呈现，已被评为几个二维增强体（看参考文献27）。在考虑下，这个测试方法对三维角度相互联锁的预型件验证过41。图11显示出了不同阶段的三维织物样品的剪变形。在开始，经纱和纬纱是正交的（图11a ），没有剪切。然后，经纱和纬纱围绕织点旋转，在锁定角之前，剪切变形阻力主要是由于经纱和纬纱之间摩擦（图11b）。在剪切过程中，纱线之间的间隙逐渐消失和相邻的纱线相互接
12、触，但纱线的宽度几乎不降低（图11 c）所示。在锁定角后，大的剪切变形过程中，在横向压缩下，纱线的宽度开始减小，为纱线在起皱之前提供更多的剪切空间。剪切载荷迅速增加（图11 d）。进一步剪切导致纱线局部起皱，如图所12。起皱出现在样品沿载荷方向的中心区域，因为压缩来自左边和右边的相框接缝。图11三维变形过程的示意图图12剪切变形中起皱的产生，如图15所示。锁定剪切角与织物密度密切相关。当剪切角度达到一定值时，P3纤维体积 比其它两种高，这是首先起皱原因 。另外，在织物纱的纱线已被压缩在编织过程中特别是对P3，P3的纱线 宽度为1mm，而P1纱线宽度是1.4毫米和P2的纱线宽度是 1.7毫米。因
13、此，大的织物密度的剪切应力在压缩阶段迅速增加，织物是不容易变形的。纱线宽度和剪切角之间的关系可以被描述为如下 38,42： （4）其中Wo是初始纱线宽度，W是剪切过程中的纱线宽度。图13三维角联锁织物的面内剪切应力-剪切角曲线图14三维角联锁织物的纤维体积分数-剪切角曲线图15织物样本的剪切锁定角在锁定剪切角之前，纱线宽度变化不明显和粘合纱线宽度减小，由于相邻的经纱（或纬纱）的压缩。在这个阶段中，横向 压缩小，纱线之间仍有空隙（如图16）。当剪切角到达锁定剪切角， 经纱（或纬纱）的宽度迅速减小，这意味着纱线消失间隙，纱线上更大的侧横向压力（图16 b）。进一步剪切导致的平面外局部起皱。图16剪
14、切过程中的三维变形的微观结构；（a）是r=30，（b）是r=42因此，可以将维角联锁织物大的剪切变形划分为三个典型阶段，如图17. 在第一阶段中，纱线的宽度可以几乎无变化，摩擦是到剪切变形的主要阻力。在第二阶段，纱线宽度迅速减少和接结纱被压缩到扁平。在第三阶段，局部起皱 发生在织物上，同时纱线的宽度保持恒定， 因为织物达到最大纤维体积分数。图17三维角联锁织物剪切过程的三个阶段在三维角联锁增强体情况下，之前详细描述的机制受存在的穿过厚纱线的影响。接结纱阻碍面内纱线的旋转，增加局部压力和纱线的弯曲（见图16）。43里展示了用相框法获得的剪切行为对所使用的夹具敏感。因此，通过相框的剪切行为的比较应
15、当考虑用相同的实验框架。对于这篇研究，三维预型件用相同的测试装置测量的的剪切响应和一些在44试验研究的二维织物比较。为了观察织物构造的影响，比较如图18所示。表2说明了用于比较的二维织物的参数和相关的参考。三维预型件呈现一个小的初始非线性区域，它的张力取决于接合纱线，而结合纱线影响纱线的移动和摩擦。同样的原因决定了接下来三维织物的剪切阻力更快速增长。图18二维织物和三维变形的相框测试法表2二维织物的结构参数44跨层剪切是三维织物典型的机械性能。跨层剪切强度主要来源于从上层到低层交织的接合纱线45。图19显示的三维角度联锁织物微观结构的失效模式，分为两部分。跨层失效模式表明结合纱线从织物里拉出，
16、整体结构被破坏。三纱系统不再彼此交织（图19 b）。另外，在实验过程，一些经纱和纬纱中从织物的中间层脱落。当测试机器工作，样本一侧沿着横梁移动。当设备左边从边缘到样本中间，经纱和纬纱被剥离当接合纱线从织物里抽出。样本从织物边缘开始受损。图20显示了样本不同的剪切变形阶段。结构变松在变形后，部分接合纱线从样本推出，经纱和纬纱由接结纱连接的可自由移动(图20 b)。当应力达到的剪切强度时，粘合纱线可以从三维织物自由提取。对于经纱和纬纱的脱落没有约束。最后，粘合纱线完全从织物离开，结构完全受损(图20 c)。图19三维织物跨层剪切的失效形态；（a）是俯视图（b）是前视图图20三维角联锁织物的跨层剪切
17、；（a）是变形前，（b）是变形后，（c）是失效模式图21描述了3种样品剪切应力-位移曲线，在其中所有的曲线非线性增加。应力达到最大后，仅位移的增加可导致在织物扯分成两部分。三条曲线的趋势是相似的。剪切强度和刚度都随着织物密度增加而增加。此外，当织物密度的变化，剪切的应力成倍的增加。接合纱线从样本边缘拉出，经纱先从边缘抽出，然后纬纱也脱离。沿着经纱方向的剪切应力大于纬纱方向的剪切应力。这是因为接合纱线平行于经纱，与纬纱交织。这使得接结纱更难以从织物中抽出。但是，接结纱在纬纱方向只能起次要的作用。另外，沿着纬纱方向的最大值得曲线形状比沿着经纱方向的最大值得曲线形状更光滑。图21三维角联锁织物的跨层
18、剪切的剪切应力-位移曲线（a）是经纱方向（b）是纬纱方向图22给出了3种织物样品的剪切强度。可以得到沿经线的剪切强度比较沿纬纱的剪切强度大。当织物密度的增加，剪切强度在经纱方向上增加，和纬纱方向增加的方式相似。因为对于大密度织物的密度，在单位区域有更多交织的次数，所以接合纱线需要更多的拉力才能拉出。虽然在两个方向的最大应力不同，但最大应力产生的位移区域相似。对于经纱，范围是710.5毫米，对于纬纱对应的是68毫米。这可以是归于产于织物结构的不同的变形机理。图22三维角联锁织物的跨层剪切强度（a）是经纱方向（b）是纬纱方向三维织物的面内剪切，跨层剪切行为也受接合纱线的限制。接结纱减少层的滑动，并
19、在跨层剪切变形中扮演一个关键的角色。已被证实，三维织物具有优异的层间剪切性能46,47。三维预型件的跨层剪切机制和二维织物不同。二维织物的剪切行为主要来自织物层之间的摩擦，然而，三维织物的剪切性能依靠和纬纱交错的接合纱线48,49。可以得到织物结构对三维预型体的剪切行为有重要的影响。 三维预型体比二维织物有个呢过大的剪切载荷，接合纱线限制经纱和纬纱的滑动。在三维角联锁织物的冲压操作中，内/间层抗剪性能是两个重要因素。现有的研究重点在实验和面内较大的剪切的表征和层间剪切变形。应用修改后的相框试验法，获得更准确和可靠的结果。织物的密度对面内的剪切性能有重要的影响。基于实验，三维织物的大剪切行为机制
20、被进一步研究。设计一个新的测试装置来表征三维织物的跨层剪切特性。分析了失效模式和跨层剪切的变形机制。结果证明了对于跨层剪切变形，该测试方法是合适的。为了进一步研究大的剪切变形，基于本文描述的剪切测试，可以建立一个微观的理论模型。同时，可以观察变形后的三维角联锁织物的微观结构。另外，研究其他与内/间层抗剪行为相关的架构因素，为大的变形机制提供一个更普遍的结论。参考文献1 Chen XG, Taylor LW, Tsai LJ. An overview on fabrication of three-dimensional woven textile preforms for composites
21、. Text Res J 2011;81:93244.2 Bilisik K. Multiaxis 3D woven preform and properties of multiaxis 3D woven and 3D orthogonal woven carbon/epoxy composites. J Reinf Plast Compos 2010;29:117386. 3 Boisse P. Mesomacro approach for composites forming simulation. J Mater Sci 2006;41:65918.4 Wang B, Wu LZ, J
22、in X, Du SY, Sun YG, Ma L. Experimental investigation of 3D sandwich structure with core reinforced by composite columns. Mater Des 2010;31:158165165. 5 Bilisik K. Multiaxis three-dimensional weaving for composites: a review. Text Res J 2012;82:72543. 6 Ansar M, Wang XW, Zhou CW. Modeling strategies
23、 of 3D woven composites: a review. Compos Struct 2011;93:194763. 7 Liu YW, Lv LH, Sun BZ, Hu H, Gu BH. Dynamic response of 3D biaxial spacer weft-knitted composite under transverse impact. J Reinf Plast Compos 2006;25:162941.8 Hufenbach W, Boehm R, Thieme M, Winkler A, Maeder E, Rausch J, et al. Pol
24、ypropylene/glass fibre 3D-textile reinforced composites for automotive applications. Mater Des 2011;32:146876. 9 Li DS, Fang DN, Zhang GB, Hu H. Effect of temperature on bending properties and failure mechanism of three-dimensional braided composite. Mater Des 2012;41:16770. 10 Jin LM, Niu ZL, Jin B
25、C, Sun BZ, Gu BH. Comparisons of static bending and fatigue damage between 3D angle-interlock and 3D orthogonal woven composites. J Reinf Plast Compos 2012;31:93545. 11 Lin H, Zeng X, Sherburn M, Long AC, Clifford MJ. Automated geometric modelling of textile structures. Text Res J 2012;82:1689702. 1
26、2 Ivanov DS, Lomov SV, Bogdanovich AE, Karahan M, Verpoest I. A comparative study of tensile properties of non-crimp 3D orthogonal weave and multi-layer plain weave e-glass composites. Part 2: Comprehensive experimental results. Compos. A 2009;40:114457. 13 Li CD, Chen ZF, Zhu JX, Liu Y, Jiang Y, Gu
27、an TR, et al. Mechanical properties and microstructure of 3D orthogonal quartz fiber reinforced silica composites fabricated by silicasol-infiltration-sintering. Mater Des 2012;36:28995. 14 El Abed B, Msahli S. Numerical simulation of woven fabric wrinkling. J Text I 2011;102:7786. 15 Lomov SV, Bogd
28、anovich AE, Ivanov DS, Mungalov D, Karahan M, Verpoest I. A comparative study of tensile properties of non-crimp 3D orthogonal weave and multi-layer plain weave e-glass composites. Part 1: Materials, methods and principal results. Composites Part A 2009;40:113443. 16 Lomov SV, Perie G. Modeling thre
29、e-dimensional fabrics and three-dimensional reinforced composites: challenges and solutions. Text Res J 2011;81:2841. 17 Zhu B, Yu TX, Tao XM. Large shear deformation of e-glass/polypropylene woven fabric composites at elevated temperatures. J Reinf Plast Compos 2009;28:261530. 18 Zhu B, Yu TX, Tao
30、XM. An experimental study of in-plane large shear deformation of woven fabric composite. Compos Sci Technol 2007;67:25261. 19 Vanclooster K, Lomov SV, Verpoest I. Investigation of interply shear in composite forming. Int J Mater Form 2008;1:95760. 20 Zhu B, Yu TX, Teng J, Tao XM. Theoretical modelin
31、g of large shear deformation and wrinkling of plain woven composite. J Compos Mater 2009;43:12538.21 Zhu B, Yu TX, Tao XM. Large deformation and slippage mechanism of plain woven composite in bias extension. Composites Part A2007;38:18218. 22 Hivet G, Duong AV. Acontribution to the analysis of the i
32、ntrinsic shear behavior of fabrics. J Compos Mater 2010;45:695716. 23 Launay J, Hiveta G, Duong AV, Boisse P. Experimental analysis of in plane shear behaviour of woven composite reinforcements. Influenceof tensions. In: 10th ESAFORM conference on material forming, Pts Aand B; 2007. p. 10338. 24 Van
33、clooster K, Lomov SV, Verpoest I. Simulation of multi-layered composites forming. Int JMater Form 2010;3:6958. 25 Willems A, Lomov SV, Verpoest I, Vandepitte D. Picture frame shear tests on woven, textile composite reinforcements with controlled pretension. In: 10th ESAFORM conference on material fo
34、rming, Pts Aand B; 2007. p. 9991004. 26 Lin H, Clifford MJ, Long AC, Sherburn M. Finite element modelling of fabric shear. Model Simul Mater Sci 2009;17:116. 27 Cao J, Akkerman R, Boisse P, Chen J, Cheng HS, de Graaf EF, et al. Characterization of mechanical behavior of woven fabrics: experimental m
35、ethods and benchmark results. Composites Part A2008;39:103753. 28 Chen QQ, Boisse P, Park CH, Saouab A, BrardJ. Intra/inter-ply shear behaviors of continuous fiberreinforced thermoplastic composites in thermoforming processes. Compos Struct 2011;93:1692703. 29 Charmetant A, Orliac JG, Vidal-SallE, B
36、oisse P. Hyperelastic model for large deformation analyses of 3D interlock composite preforms. Compos Sci Technol 2012;72:135260. 30 Zhu B, Yu TX, Zhang H, Tao XM. Experimental investigation of formability of woven textile composite preform in stamping operation. Int JMater Form 2008;1:96972. 31 Boi
37、sse P, Zouari B, Daniel JL. Importance of in-plane shear rigidity in finite element analyses of woven fabric composite preforming. Composites Part A 2006;37:220112. 32 Lomov SV, Verpoest I. Model of shear of woven fabric and parametric description of shear resistance of glass woven reinforcements. C
38、ompos Sci Technol 2006;66:91933. 33 Harrison P, Abdiwi F, Guo Z, Potluri P, Yu WR. Characterising the sheartension coupling and wrinkling behaviour of woven engineering fabrics. Composites Part A 2012;43:90314. 34 Gasser A, Boisse P, Hanklar S. Mechanical behaviour of dry fabric reinforcements. 3D s
39、imulations versus biaxial tests. Compos Mater Sci 2000;32:43948. 35 Lomov SV, Barburski M, Stoilova T, Verpoest I, Akkerman R, Loendersloot R, et al. Carbon composites based on multiaxial multiply stitched preforms. Part 3: Biaxial tension, picture frame and compression tests of the preforms. Compos
40、 A Appl Sci Manuf 2005;36:1188206. 36 Sun HY, Pan N. Shear deformation analysis for woven fabrics. Compos Struct 2005;67:31722. 37 De Luycker E, Morestin F, Boisse P, Marsal D. Simulation of 3D interlock composite preforming. Compos Struct 2009;88:61523. 38 Liu L, Chen JL, Gorczyca JL, Sherwood JA.
41、Modeling of friction and shear in thermostamping of composites Part II. JCompos Mater 2004;38:193147. 39 Hivet G, Launay J, Gasser A, Daniel JL, Boisse P. Mechanical behavior of woven composite reinforcements while forming. JThermoplast Compos 2002;15:54555. 40 Willems A, Lomov SV, Verpoest I, Vande
42、pitte D. Drape-ability characterization of textile composite reinforcements using digital image correlation. Opt Laser Eng 2009;47:34351. 41 Carvelli V, Pazmino J, Lomov SV, Verpoest I. Deformability of anon-crimp 3D orthogonal weave e-glass composite reinforcement. Compos Sci Technol 2012;73:918. 4
43、2 Gorczyca JL. Modeling of friction and shear in thermostamping of composites Part I. JCompos Mater 2004;38:191129. 43 Willems A, Lomov SV, Verpoest I, Vandepitte D. Optical strain fieldsin shear and tensile testing of textile reinforcements. Compos Sci Technol 2008;68:80719. 44 Zhang YF, Chen L, Su
44、n F, Sun Y, Li YL, Tang BM, et al. Experimental research on shear deformation of tackif ied woven fabrics. Acta Mater Compos Sin 2009;26:2934. 45 Mehrabadi FA. Analysis of pure mode III and mixed mode (IIIplus II) interlaminar crack growth in polymeric woven fabrics. Mater Des 2013;44:42937. 46 Pand
45、ya KS, Veerraju C, Naik NK. Hybrid composites made of carbon and glass woven fabrics under quasi-static loading. Mater Des 2011;32:40949. 47 Desplentere F, Lomov SV, Woerdeman DL, Verpoest I, Wevers M, Bogdanovich A. Micro-CT characterization of variability in 3D textile architecture. Compos Sci Tec
46、hnol 2005;65:192030. 48 Bel S, Hamila N, Boisse P, Dumont F. Finite element model for ncf composite reinforcement preforming: importance of inter-ply sliding. Composites Part A 2012;43:226977. 49 Ullah H, Harland AR, Silberschmidt VV. Experimental and numerical analysis of damage in woven GFRP composites under large-deflectionbending. Appl Compos Mater 2012;19:76983. Y. Zhang et al. /Materials and Design 49 (2013)151159 159