CN103455713B - 一种设计平面编织复合材料力学模量的方法 - Google Patents
一种设计平面编织复合材料力学模量的方法 Download PDFInfo
- Publication number
- CN103455713B CN103455713B CN201310363526.1A CN201310363526A CN103455713B CN 103455713 B CN103455713 B CN 103455713B CN 201310363526 A CN201310363526 A CN 201310363526A CN 103455713 B CN103455713 B CN 103455713B
- Authority
- CN
- China
- Prior art keywords
- cell space
- space unit
- woven composite
- warp
- yarn
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Woven Fabrics (AREA)
Abstract
一种设计平面编织复合材料力学模量的方法,它有四大步骤:步骤一、根据纤维束的平面编织方式即周期性和重复性,选择最小的重复性单元作为代表性体积元,由此确定其胞体单元;步骤二、根据外载荷施加方式以及步骤一中的胞体单元,对胞体单元内织布的纤维束进行受力分析,从而,建立平面编织复合材料胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能U*,并利用最小势能原理,求解胞体单元纤维织布的内力;步骤三、根据卡式定理或单位载荷法求解胞体单元纤维织布的变形,再根据应力-应变本构方程,得到平面编织复合材料织布的力学模量;步骤四、根据混合定理,得到平面编织复合材料的弹性模量。
Description
技术领域
本发明提供一种设计平面编织复合材料力学模量的方法,属于复合材料设计领域。
背景技术
平面编织复合材料作为一种新型轻质高效纺织结构复合材料,具有工艺自动化程度水平高、经济性好、结构整体性能优异、抗冲击和疲劳阻等优点,与单向纤维增强的复合材料层合板相比,平面编织复合材料层板对改进层间、层内强度和损伤容限等方面具有巨大的潜力。因此,编织复合材料在航空、航天、航海、汽车等领域都得到了广泛应用。实验手段直接测量平面编织复合材料残余热应力成本较高,且测试过程中易受到很多偶然因素的影响;有限元数值模拟方法需要建立复杂的有限元模型,计算复杂,计算效率低,计算精度难以保证;因此,本发明运用细观力学分析方法研究平面编织复合材料细观结构对其宏观性能的影响,获得其宏观力学模量的解析解,仅仅需要少量的组分材料性能参数就能快速准确地预测平面编织复合材料宏观力学模量,实现平面编织复合材料宏观性能的优化设计,可见本发明具有重要学术意义和工程应用价值。
发明内容
本发明提供了一种设计平面编织复合材料力学模量的方法,该方法具有计算简便,精度高等优点,其技术方案如下:
步骤一、根据纤维束的平面编织方式(如周期性和重复性等),选择最小的重复性单元作为代表性体积元,由此确定其胞体单元;图2为所选的代表性体积元(胞体单元),包含了2条正交的经向纱和纬向纱(纤维束),其中方向1定义为经向,方向2定义为纬向,将经向纱和纬向纱理想化成正弦的曲梁,可以得到经向纱和纬向纱的中心线Z坐标表达式:
根据图3a、图3b,可以得到截面面积A和截面惯性矩I的表达式为
式中a和b分别为经向纱和纬向纱截面的宽度和高度。
步骤二、根据外载荷施加方式以及步骤一中的胞体单元,对胞体单元内织布的纤维束进行受力分析,从而,建立平面编织复合材料胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能U*,并利用最小势能原理,求解胞体单元纤维织布的内力。
胞体单元的总应变余能U*的表达式为
式中,
其中,M、N和T分别为弯矩、轴力和扭矩;I、A和Ip分别为纤维束截面惯性矩、面积和极惯性矩;E为纤维束的弹性模量。根据最小势能原理,可确定胞体单元各个内力。
步骤三、根据卡式定理或单位载荷法求解胞体单元纤维织布的变形,再根据应力-应变本构方程,得到平面编织复合材料织布的力学模量。
由卡式定理确定的系统在外载荷(包括拉伸、压缩和剪切载荷)作用下平面编织复合材料织布的变形Δ为
其中P代表所受的外载荷。
根据应力和应变关系表达式,可以得到纤维织布的拉伸、压缩和剪切模量公式:
式中Et、Ec和G分别表示纤维织布的拉伸、压缩和剪切模量,σtx、σcx和τ分别表示拉伸、压缩和剪切应力,εtx、εcx和γ分别表示拉伸、压缩和剪切应变。
步骤四、根据混合定理,可以得到平面编织复合材料的弹性模量。
其中,在步骤一中所述的“体积元”是指具有三轴向编织特征的三维单元。
本发明是一种设计平面编织复合材料力学模量的方法,其特点是通过少量的纤维束和基体的性能参数便可方便快捷预测平面编织复合材料的宏观力学模量。
附图说明
图1为平面编织复合材料的编织方式。
图2为理想胞体单元。
图3a为胞体单元经向纱和纬向纱(纤维束)的理想截面。
图3b为胞体单元经向纱和纬向纱的简化截面。
图4a为胞体单元拉伸状态下受力图。
图4b为胞体单元压缩状态下受力图。
图4c为胞体单元剪切状态下受力图。
图5a为拉伸载荷下经向纤维束的内力分布图。
图5b为拉伸载荷下纬向纤维束的内力分布图。
图6a为压缩载荷下经向纤维束的内力分布图。
图6b为压缩载荷下纬向纤维束的内力分布图。
图7a和b为剪切载荷下经向纤维束的内力和合力分布图。
图7c和d为剪切载荷下纬向纤维束的内力和合力分布图。
图8是本发明所述方法的流程框图。
图中符号说明如下:
图2中的a1、a2分别为经向纱和纬向纱截面的宽度,b1、b2分别为经向纱和纬向纱截面的高度,L1、L2分别为经向纱和纬向纱纤维束的波长,g1、g2分别为经向纱和纬向纱的间距。
图3a中a2分别为纬向纱截面的宽度,b1、b2分别为经向纱和纬向纱截面的高度,L1为经向纱纤维束波长,g2分别为纬向纱的间距,h1为经向纱纤维束的高度。图3b中a和b分别为经向纱和纬向纱截面的宽度和高度。
图4a中的N1为拉伸载荷,N2为纤维束间的相互作用力。图4b中的N1为压缩载荷,N2为纤维束间的相互作用力。图4c中的N1为剪切载荷。
图5a、b中的L1、L2分别为经向纱和纬向纱的波长,N1为拉伸载荷,N2为纤维束间的相互作用力,M1、M2为由拉伸载荷产生的经向纱和纬向纱弯矩。
图6a、b中的N1为拉伸载荷,N2为纤维束间的相互作用力,M1、M2为由拉伸载荷产生的经向纱和纬向纱弯矩。
图7a-d中的N1为剪切载荷,N2、N3为由剪切载荷产生的垂向和横向内力,N4、N5为由N2、N3产生的合力,M1、M2为由剪切载荷产生的弯矩,T1、T2为由剪切载荷产生的扭矩。
具体实施方式
见图8,本发明一种设计平面编织复合材料力学模量的方法,该方法具体步骤如下:
步骤一、胞体单元的确定。根据图1所示的织布编织方式,考虑编织的周期性和重复性,选择代表性体积元模型,图2为所选的代表性体积元(胞体单元),包含了2条正交的经向纱和纬向纱(纤维束),其中方向1定义为经向,方向2定义为纬向,将经向纱和纬向纱理想化成正弦的曲梁,可以得到经向纱和纬向纱的中心线Z坐标表达式:
根据图3a、图3b,可以得到截面面积A和截面惯性矩I的表达式为
式中a和b分别为经向纱和纬向纱截面的宽度和高度。
步骤二、平面编织复合材料胞体单元纤维织布的细观力学模型。根据不同外载情况(包括拉伸、压缩和剪切载荷)(如图4a-图4c所示),分析纤维束之间的相互作用,可以得到不同外载下(包括拉伸、压缩和剪切载荷)经纬向纤维束的内力分布图(如图5a、图5b,图6a、图6b到图7a-图7d所示),根据能量原理,求解内力。
胞体单元的总应变余能U*的表达式为
式中,
其中,M、N和T分别为弯矩、轴力和扭矩;I、A和Ip分别为纤维束截面惯性矩、面积和极惯性矩;E为纤维束的弹性模量。根据最小势能原理,可确定胞体单元各个内力。
步骤三、力学模量(包括拉伸、压缩和剪切模量)计算公式的建立。由卡式定理确定的系统在外载荷(包括拉伸、压缩和剪切载荷)作用下平面编织复合材料织布的变形Δ为
其中P代表所受的外载荷。
根据应力和应变关系表达式,可以得到纤维织布的拉伸、压缩和剪切模量公式:
式中Et、Ec和G分别表示纤维织布的拉伸、压缩和剪切模量,σtx、σcx和τ分别表示拉伸、压缩和剪切应力,εtx、εcx和γ分别表示拉伸、压缩和剪切应变。
步骤四、根据混合定理,可以得到平面编织复合材料的弹性模量为
下面结合附图和具体实施方式对本发明作进一步详细说明。
实施例一拉伸模量的预测
选用如图2所示的编织方式和胞体单元,当单元受到拉伸外载时(如图4a所示),考虑纤维束间的挤压力N2,结合纤维束的受力分析(如图5a-b所示),可以得到相应的内力表达式:
式中,M1、M2为由拉伸载荷产生的经向纱和纬向纱弯矩,N1为拉伸载荷,N2为纤维束间的相互作用力,L1、L2分别为经向纱和纬向纱纤维束的波长,h1为经向纱纤维束的高度。
将式(10)代入式(6),可以得到经向和纬向纤维束的余应变能分别为
式中U1 *和U2 *分别为经向和纬向纤维束的余应变能,E1为纤维束的弹性模量,M1、M2为由拉伸载荷产生的经向纱和纬向纱弯矩,N1为拉伸载荷,N2为纤维束间的相互作用力,L1、L2分别为经向纱和纬向纱纤维束的波长,h1为经向纱纤维束的高度。
为方便计算,定义J1到J10如下:
式中L1、L2分别为经向纱和纬向纱纤维束的波长,h1、h2分别为经向纱和纬向纱纤维束的高度,I1、I2分别为经向纱和纬向纱纤维束的截面惯性矩。
将式(13)至式(21)分别代入式(11)、(12),再代入式(5)后化简,可以得到胞体单元的总余应变能U*为
根据最小余能原理得到
解方程组(23),可以得到
由卡氏定理,得到编织布的变形Δ:
将式(25)进一步化简,变为
其中J的表达式为
于是,受拉纤维布的拉伸模量为
将式(28)代入式(9),可以得到平面编织复合材料的拉伸模量为
同理,可以导出压缩和剪切模量公式。
Claims (2)
1.一种设计平面编织复合材料力学模量的方法,其特征在于:该方法具体步骤如下:
步骤一、根据纤维束的平面编织方式即周期性和重复性,选择最小的重复性单元作为代表性体积元,由此确定其胞体单元;该胞体单元包含了2条正交的经向纱和纬向纱,其中方向1定义为胞体单元的经向,指代经向纱长度方向;方向2定义为胞体单元的纬向,指代纬向纱长度方向;
将经向纱和纬向纱理想化成正弦的曲梁,得到经向纱和纬向纱的中心线Z坐标表达式:
从而得到截面面积A和截面惯性矩I的表达式为
式中a和b分别为经向纱和纬向纱截面的宽度和高度;式中h1为平面编织复合材料经向纤维束的厚度;h2为平面编织复合材料纬向纤维束的厚度;L1为平面编织复合材料胞体单元中经向纤维束的一个波动周期内长度;L2为平面编织复合材料胞体单元中纬向纤维束的一个波动周期内长度;
步骤二、根据外载荷施加方式以及步骤一中的胞体单元,对胞体单元内织布的纤维束进行受力分析,从而,建立平面编织复合材料胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能U*,并利用最小势能原理,求解胞体单元纤维织布的内力;
胞体单元的总应变余能U*的表达式为
式中,
其中,M、N和T分别为弯矩、轴力和扭矩;I、A和Ip分别为纤维束截面惯性矩、面积和极惯性矩;E为纤维束的弹性模量;根据最小势能原理,确定胞体单元各个内力;
步骤三、根据卡式定理或单位载荷法求解胞体单元纤维织布的变形,再根据应力-应变本构方程,得到平面编织复合材料织布的力学模量;
由卡式定理确定的系统在外载荷作用下平面编织复合材料织布的变形Δ为
其中P代表所受的外载荷,
根据应力和应变关系表达式,得到纤维织布的拉伸、压缩和剪切模量公式:
式中Et、Ec和G分别表示纤维织布的拉伸、压缩和剪切模量,σtx、σcx和τ分别表示拉伸、压缩和剪切应力,εtx、εcx和γ分别表示拉伸、压缩和剪切应变;
步骤四、根据混合定理,得到平面编织复合材料的弹性模量;
式中Etla为平面编织复合材料的拉伸弹性模量;Ecla为平面编织复合材料的压缩弹性模量;Gla为平面编织复合材料的剪切弹性模量;Em为基体的压缩弹性模量;Gm为基体的压缩剪切模量;Vf为纤维织布的体积分数。
2.根据权利要求1所述的一种设计平面编织复合材料力学模量的方法,其特征在于:步骤一中所述的“体积元”是指具有三轴向编织特征的三维单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310363526.1A CN103455713B (zh) | 2013-08-20 | 2013-08-20 | 一种设计平面编织复合材料力学模量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310363526.1A CN103455713B (zh) | 2013-08-20 | 2013-08-20 | 一种设计平面编织复合材料力学模量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103455713A CN103455713A (zh) | 2013-12-18 |
CN103455713B true CN103455713B (zh) | 2016-06-22 |
Family
ID=49738067
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310363526.1A Expired - Fee Related CN103455713B (zh) | 2013-08-20 | 2013-08-20 | 一种设计平面编织复合材料力学模量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103455713B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105354390B (zh) * | 2015-12-01 | 2018-11-06 | 北京航空航天大学 | 一种预测平面斜交编织复合材料双向拉伸模量与强度的方法 |
CN107633120B (zh) * | 2017-09-07 | 2018-07-17 | 东南大学 | 一种纤维增强复合材料动态剪切本构模型的构建方法 |
CN108268691A (zh) * | 2017-11-21 | 2018-07-10 | 北京航空航天大学 | 自动铺丝复合材料力学性能的模拟分析方法 |
CN111400947B (zh) * | 2020-03-16 | 2022-05-10 | 北京航空航天大学 | 一种预测平面正交编织复合材料压缩模量与强度的方法 |
-
2013
- 2013-08-20 CN CN201310363526.1A patent/CN103455713B/zh not_active Expired - Fee Related
Non-Patent Citations (6)
Title |
---|
三维四向编织复合材料有效弹性模量的预测;马忠辉等;《机械科学与技术》;20041130;第23卷(第11期);1334-1337 * |
三维编织复合材料细观结构与力学性能分析;郑锡涛;《万方学位论文》;20040427;全文 * |
三维编织复合材料细观结构平均刚度及其压缩性能的研究;杨朝坤;《万方学位论文》;20110818;全文 * |
基于微胞体划分的三维编织复合材料有效力学性能分析;洪平等;《宇航材料工艺》;20051117;第2005年卷(第5期);32-36 * |
夹芯结构复合材料弹性力学性能预测与试验研究;李海伟;《万方学位论文》;20130320;全文 * |
平纹织物复合材料的弹性模量预测;王瑞等;《复合材料学报》;20020228;第19卷(第1期);90-94 * |
Also Published As
Publication number | Publication date |
---|---|
CN103455713A (zh) | 2013-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rajesh et al. | Experimental investigation on buckling and free vibration behavior of woven natural fiber fabric composite under axial compression | |
Zhao et al. | Modeling the transverse tensile and compressive failure behavior of triaxially braided composites | |
Thanh et al. | A refined size-dependent couple stress theory for laminated composite micro-plates using isogeometric analysis | |
Li et al. | Repeated low-velocity impact response and damage mechanism of glass fiber aluminium laminates | |
Zhang et al. | A meso-scale finite element model for simulating free-edge effect in carbon/epoxy textile composite | |
CN104537259B (zh) | 一种基于细观力学模型的横观各向同性纤维材料性能的预测方法 | |
CN103455712B (zh) | 一种三轴向编织复合材料格栅力学模量的预测方法 | |
CN106066913A (zh) | 复杂复合材料结构等效材料性能多尺度计算方法 | |
CN111368389B (zh) | 一种预测复合材料层合板失效强度的方法 | |
Naveen et al. | Finite element analysis of natural fiber-reinforced polymer composites | |
CN103455713B (zh) | 一种设计平面编织复合材料力学模量的方法 | |
CN110688789B (zh) | 一种碳纤维复合单向板偏载冲击的多尺度响应分析方法 | |
CN107451309A (zh) | 一种多尺度计算复杂复合材料结构等效热膨胀系数的方法 | |
CN105354390A (zh) | 一种预测平面斜交编织复合材料双向拉伸模量与强度的新方法 | |
Bacarreza et al. | Micromechanical modeling of advanced composites | |
Mauget et al. | Large deformation nonlinear response of soft composite structures via laminate analogy | |
Sun et al. | Multiscale uncertainty propagation analysis and reliability optimization of the CFRP crossbeam of the twist beam axle | |
Ashir et al. | A statistical approach for the fabrication of adaptive pleated fiber reinforced plastics | |
Huang et al. | Computational analysis of torsional bulking behavior of 3D 4-directional braided composites shafts | |
Chen et al. | Numerical investigation of the compressive behavior of 2.5 D woven carbon-fiber (2.5 D-CFRP) honeycomb and experimental validation | |
CN103473440A (zh) | 一种预测平面编织复合材料残余热应力的新方法 | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
Boisse | Simulations of woven composite reinforcement forming | |
Liu et al. | Multiscale modeling of viscoelastic behaviors of textile composites using mechanics of structure genome | |
Ji et al. | Study of in-plane fatigue failure and life prediction of weave composites under constant and variable amplitude loading |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160622 Termination date: 20200820 |
|
CF01 | Termination of patent right due to non-payment of annual fee |