CN108388691A - 复材分层缺陷在横刃挤压阶段扩展行为计算方法 - Google Patents
复材分层缺陷在横刃挤压阶段扩展行为计算方法 Download PDFInfo
- Publication number
- CN108388691A CN108388691A CN201810039038.8A CN201810039038A CN108388691A CN 108388691 A CN108388691 A CN 108388691A CN 201810039038 A CN201810039038 A CN 201810039038A CN 108388691 A CN108388691 A CN 108388691A
- Authority
- CN
- China
- Prior art keywords
- fiber
- defect
- point
- load
- lamination defect
- 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.)
- Granted
Links
- 230000007547 defect Effects 0.000 title claims abstract description 108
- 238000003475 lamination Methods 0.000 title claims abstract description 51
- 239000000463 material Substances 0.000 title claims abstract description 24
- 238000000205 computational method Methods 0.000 title claims abstract description 17
- 230000006835 compression Effects 0.000 title claims abstract description 13
- 238000007906 compression Methods 0.000 title claims abstract description 13
- 239000000835 fiber Substances 0.000 claims abstract description 90
- 238000000034 method Methods 0.000 claims abstract description 36
- 239000011159 matrix material Substances 0.000 claims abstract description 30
- 238000000465 moulding Methods 0.000 claims abstract description 27
- 238000005553 drilling Methods 0.000 claims abstract description 19
- 230000000694 effects Effects 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000005452 bending Methods 0.000 claims abstract description 6
- 238000001125 extrusion Methods 0.000 claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 19
- 238000004458 analytical method Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 101100532451 Rattus norvegicus Slc22a17 gene Proteins 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims description 2
- 238000010008 shearing Methods 0.000 claims description 2
- 239000011347 resin Substances 0.000 claims 1
- 229920005989 resin Polymers 0.000 claims 1
- 239000004918 carbon fiber reinforced polymer Substances 0.000 abstract description 3
- 239000002131 composite material Substances 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000032798 delamination Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 229920000049 Carbon (fiber) Polymers 0.000 description 2
- 239000004917 carbon fiber Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000002990 reinforced plastic Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Moulding By Coating Moulds (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明复材分层缺陷在横刃挤压阶段扩展行为计算方法属于机械加工技术领域,涉及一种碳纤维复合材料构件内部的成型分层缺陷在钻削横刃挤压阶段扩展行为的计算方法。该计算方法将横刃对缺陷层的挤压过程简化为集中载荷作用下的弹性地基梁弯曲过程,并求得梁的弯曲变形控制方程;确定刀具作用位置以及成型分层缺陷大小,根据约束情况及受力特点对缺陷层中的纤维进行区域划分,利用有限元计算方法及叠加原理,得到整根纤维的整体刚度矩阵;最后,采用向刀具作用点逐步加载方式对成型分层缺陷扩展的临界条件及扩展长度进行计算。该方法适用于不同加工位置、不同大小形状成型分层缺陷扩展行为的计算,抑制分层缺陷的扩展,保证最终制孔质量。
Description
技术领域
本发明属于机械加工技术领域,涉及一种碳纤维复合材料构件内部的成型分层缺陷在钻削横刃挤压阶段扩展行为的计算方法,是一种抑制构件内部缺陷在机械加工过程中进一步扩展的方法。
背景技术
碳纤维复合材料承力构件常采用铺放、固化、加工一体化设计制造,相比于传统材料构件的制造方式,有效缩短了工艺流程,大幅提高了生产效率。然而,由于碳纤维复合材料的各向异性、层叠成型特征,以及热压固化过程中的热、化学收缩非均匀性,其在成型过程中易产生分层缺陷。碳纤维复合材料构件成型分层缺陷分布具有很大的随机性,对构件性能有很大影响,若成型分层缺陷在静态载荷作用下或钻削过程中发生进一步扩展,构件性能及制孔质量将难以保证。
目前,部分学者们针对碳纤维复合材料成型分层缺陷在静态载荷下的扩展行为进行了研究,如Pekbey和Sayman等人发表的《A Numerical and ExperimentalInvestigation of Critical Buckling Load of Rectangular Laminated CompositePlates with Strip Delamination》一文在《Journal of Reinforced Plastics&Composites》2006年第25期第685到697页中,对含带状分层损伤的复合材料层合板的临界屈曲载荷进行了研究,指出了试样尺寸、铺层角度及铺层顺序对临界屈曲载荷均有影响。然而,实际的钻削过程是在刀具连续进给下的动态过程,其钻削力大小非恒定且作用点位置随刀刃的回转和进给运动发生动态变化,此外由于成型缺陷分布的随机性,刀具与缺陷的实际作用位置不可控。因此,静态载荷下对分层扩展行为的研究无法适用于碳纤维复合材料的制孔过程,需对碳纤维复合材料成型缺陷在钻削加工中的扩展行为进行研究。
发明内容
本发明为解决现有技术的缺陷,有效抑制成型分层缺陷在钻削过程中的扩展,提高制孔质量,降低复材构件生产成本,针对在钻削过程中轴向力最大的横刃挤压阶段缺陷层成型分层缺陷的扩展行为进行分析,发明了复材分层缺陷在横刃挤压阶段扩展行为计算方法。该计算方法将刀具横刃对缺陷层的挤压过程简化为集中载荷作用下的弹性地基梁弯曲过程,求梁的弯曲变形控制方程;通过确定刀具作用位置以及成型分层缺陷大小,根据约束情况及受力特点对缺陷层中的纤维进行区域划分,利用有限元计算方法及叠加原理,得到整根纤维的整体刚度矩阵;并采用向刀具作用点逐步加载方式对成型分层缺陷扩展的临界条件及扩展长度进行计算,能够适用于不同加工位置下的,不同大小形状成型分层缺陷扩展行为的计算,抑制分层缺陷的扩展,保证最终制孔质量。
本发明的技术方案是一种复材分层缺陷在横刃挤压阶段扩展行为计算方法,其特征是,该计算方法将刀具横刃对缺陷层的挤压过程简化为集中载荷作用下的弹性地基梁弯曲过程,并求得梁的弯曲变形控制方程;确定刀具作用位置以及成型分层缺陷大小,根据约束情况及受力特点对缺陷层中的纤维进行区域划分,利用有限元计算方法及叠加原理,得到整根纤维的整体刚度矩阵;最后,采用向刀具作用点逐步加载方式对成型分层缺陷扩展的临界条件及扩展长度进行计算。该方法的具体步骤如下:
1)首先确定刀具作用位置及成型分层缺陷大小参数,被加工材料的性能参数;
2)简化模型
将钻削刀具的横刃与工件相互作用的过程简化为集中载荷作用下的弹性地基梁模型,选取长度为dx的微小单元体,分析其受力情况,设,pm为单位长度上周围材料对纤维的反作用力,qb为界面对纤维的粘结作用力,
pm=kmw (1)
qb=kbw (2)
其中,km为周围材料等效刚度,kb为界面的等效模量,w为微小单元体等效挠度;km通过Boit’s方程求解:
其中,为周围材料的等效模量,Ef为纤维杨氏模量,If为纤维截面惯性矩,ν为周围材料泊松比,df为纤维直径;
其次,由微小单元体的平衡方程,可得单纤维温克尔地基梁的变形控制方程为:
当界面处存在缺陷时,则缺陷的下一层,记作层S,缺陷区域为AB,载荷作用于点O,在缺陷区域纤维将不再受到界面的粘结作用,此时上述单纤维变形控制方程简化为:
3)求解缺陷层纤维整体刚度矩阵
输入刀具作用位置及分层缺陷大小,根据约束情况和受力特点,将缺陷层S层中的纤维划分为FB、BO、OA、AE四段。考虑到O点左右两侧结构相似,以O点左侧为例,采用有限元计算方法,对纤维的弯曲断裂及界面缺陷的扩展情况进行求解。对于OA段,根据简化后的变形控制方程求得挠度解为:
w(x)=eαx(B1 cos βx+B2 sin βx)+e-αx(B3 cos βx+B4 sin βx) (6)
式中,B1~B4为积分常数,且有:
梁单元内部的剪力Q和扭矩M分别表示为:
设OA段位移矢量dOA、载荷矢量rOA分别为:
dOA={wO,θO,wA,θA}T (10)
rOA={QO,MO,QA,MA}T (11)
则OA段单元刚度矩阵kOA为:
对于AE段,由于点E为无穷远处,其挠度为0,即满足:
wE=w|x=+∞=0 (13)
则AE段挠度解简化为:
w(x)=e-αx(B3 cos βx+B4 sin βx) (14)
其中,B3~B4为积分常数,且有:
设AE段位移矢量dAE、载荷矢量rAE分别为:
dAE={wA,θA}T (16)
rAE={QA,MA}T (17)
则AE段单元刚度矩阵kAE为:
综上,得到O点左侧纤维的整体载荷、位移、刚度之间的关系式:
Rleft=Kleftdleft (19)
其中,Kleft为单元刚度矩阵经叠加得到的O点左侧纤维整体刚度矩阵,dleft为O点左侧纤维整体位移矢量,Rleft t为O点左侧纤维整体载荷矢量。
采用类似的计算方法,对O点右侧纤维进行分析求解,最终确定整体纤维的载荷、位移、刚度之间的关系式:
R=Kd (20)
式中,整体载荷矢量R为:
R={QF,QB,MB,QO,MO,QA,MA,QE}T (21)
整体位移矢量d为:
矩阵K为根据FB、BO、OA、AE四段的单元刚度矩阵经叠加原理求得的整根纤维的整体刚度矩阵。
4)判断成型分层缺陷是否扩展
当求解出K之后,通过向刀具作用点O逐步施加载荷,计算出各节点位移。随着载荷不断增加,纤维的弯曲变形不断增加,界面应力若满足界面强度则点A(或点B)处界面将发生开裂,然而,在切削过程中,纤维的弯曲变形存在极限,当纤维内的截面的最大拉应力σmax达到其拉伸强度极限σt时,纤维将发生断裂。一旦纤维发生断裂,纤维的弯曲变形将不再增加,则点A(或点B)的位移也不再发生变化。因此,点A(或点B)处的界面若发生开裂,其一定发生在纤维弯曲断裂之前。结合以上分析,得到在横刃挤压阶段预埋缺陷扩展的临界条件为:
其中,设纤维半径为r,则纤维各截面最大拉应力为:
其中kb为层间界面的等效模量,σb为树脂粘结强度,wA为A点弯曲挠度,wB为B点弯曲挠度算法流程;若缺陷扩展则转至步骤5,若缺陷不扩展停止计算。
5)计算成型分层缺陷扩展长度
若满足缺陷扩展的临界条件,预埋缺陷将发生扩展,则迭代过程继续。扩展后每个增量步将逐渐增大BO或OA段长度,以及载荷矢量R的值,依据新的约束条件重新对梁内的各部分单元进行求解,更新缺陷扩展后纤维整体刚度矩阵,并求解wA、wB以及σmax的值。当σmax达到σt时,即纤维发生弯曲断裂,迭代停止,此时AB段的长度与迭代之前的长度之差即为缺陷扩展长度。
本发明的有益效果是该方法通过对成型分层缺陷在钻削横刃挤压阶段的扩展行为进行分析,求得分层缺陷扩展的临界条件及其扩展长度的计算方法,可有效抑制分层缺陷在钻削加工中的进一步扩展,保证制孔质量,降低实际复材构件的生产成本。该方法包括简化加工模型,计算纤维刚度矩阵,判断分层缺陷扩展临界条件,计算分层缺陷扩展长度四个部分,能够适用于不同加工位置下的,不同大小形状成型分层缺陷扩展行为的计算,指导实际生产中钻削工艺参数的选择,抑制分层缺陷的扩展,保证最终制孔质量。
附图说明
图1是横刃挤压阶段微小单元体的受力分析图。其中,pm为单位长度上周围材料对纤维的反作用力(N);qb为界面对纤维的粘结作用力(N);M为梁单元内部扭矩(N·m);Q为梁单元内部剪力(N·m);dM为梁单元扭矩增量(N·m)dQ为梁单元剪力增量(N·m);
图2是横刃挤压成型分层缺陷时的弹性地基梁模型。其中,FZ为刀具横刃轴向力(N),缺陷层S层中的纤维划分为FB、BO、OA、AE四段。
图3是成型分层缺陷扩展临界条件计算方法的流程图,图4是成型分层缺陷扩展长度计算方法的流程图。其中,σb(pa·m)为树脂粘结强度,σt(Gpa)为纤维拉伸强度极限,σmax(Gpa)为纤维截面最大正应力,r(μm)为纤维半径,Qo(N)为初始加载载荷,Qoi(N)为第i步加载载荷,q(N)为载荷增量步步长,i为载荷迭代次数,wA(μm)为A点弯曲挠度,wB(μm)为B点弯曲挠度,j为位移迭代次数,s(μm)为位移增量步步长,m为A点位移扩展步数,n为B点位移扩展步数,lAj(μm)为迭代第j次后A点与O点的距离,lBj(μm)为迭代第j次后B点与O点的距离,lAj-1(μm)为迭代第j-1次后A点与O点的距离,lBj-1(μm)为迭代第j-1次后B点与O点的距离,le(μm)为分层缺陷扩展长度,Kj为缺陷扩展后纤维整体刚度矩阵。
图5是成型分层缺陷扩展前边界处弯曲挠度图。
其中,横坐标-加载载荷(N),纵坐标-纤维最大弯曲正应力(MPa),1-A点挠度曲线,2-B点挠度曲线。
图6是成型分层缺陷扩展后A点扩展位移图,其中,横坐标-加载载荷(N),纵坐标-纤维最大弯曲正应力(MPa),1-弯曲正应力曲线,2-A点扩展位移曲线。
具体实施方式
下面结合技术方案和附图详细说明本发明的具体实施。
1)输入被加工材料的性能参数、刀具作用位置及成型分层缺陷大小等参数。本实例中所采用的是T800/977-2预浸料,预浸料中纤维的杨氏模量Ef为295GPa,纤维拉伸强度σt为4.5GPa,纤维半径r为3.25μm,纤维等效模量kb为115GPa,树脂粘结强度σb为50pa·m,周围材料的等效模量为9.65Gpa,周围材料泊松比υ为0.3。刀具中心位置距离缺陷左边缘0.1mm,构件中所含成型分层缺陷大小及形状为直径7mm的圆形。
2)将横刃对缺陷层的挤压过程简化为集中载荷作用下弹性地基梁弯曲过程,如图1所示,由微小单元体的平衡方程,根据公式(4)、(5)可求得单纤维的温克尔地基梁的变形控制方程和简化缺陷区域纤维的变形控制方程。其中km可通过Boit’s方程求解为:
3)每个增量步q逐渐增加Qo的值以增加载荷矢量R,计算不同载荷下纤维整体刚度矩阵K。以刀具作用位置为原点,纤维变形前所在的轴线为X轴,将缺陷左边无穷远处到缺陷左边缘区域定义为AE段,将缺陷左边缘与刀具作用位置之间0.1mm区域定义为OA段,将刀具作用点与缺陷右边缘6.9mm区域定义为OB段,将缺陷右边缘到右侧无穷远区域定义为BF段。如图2所示。图3是成型分层缺陷扩展临界条件计算方法的流程图,图4是成型分层缺陷扩展长度计算方法的流程图。利用有限元计算方法求得各段纤维、载荷、刚度之间的关系,最终确定整体纤维的载荷、位移、刚度之间的关系式为公式(20),因此,整体载荷矢量R为:
R={0,0,0,QOi,0,0,0,0}T
其中QOi=i·q,增量步q=1N,i为迭代次数。
整体位移矢量d为:
矩阵K为根据FB、BO、OA、AE四段的单元刚度矩阵经叠加原理求得的整根纤维的整体刚度矩阵。
4)求解每个增量步中缺陷边缘处的弯曲挠度wA、wB以及最大弯曲正应力σmax的值,如图5所示。根据材料参数可计算出成型缺陷边缘处产生开裂的临界失效挠度为结合所求得的刚度矩阵K,可计算出A点弯曲挠度达到260.9μm时,所需施加的载荷为53.37N,而此时纤维截面最大弯曲正应力小于纤维拉伸强度极限4500Mpa,满足缺陷扩展临界条件:
因此,A点成型分层缺陷扩展,转至步骤5求取成型分层缺陷扩展长度。
5)随着加载载荷的持续增加,缺陷边缘将不断远离刀具作用点,依据新的约束条件重新对梁内的各部分单元进行求解,更新缺陷扩展后纤维整体刚度矩阵Kj,并重新求解wA、wB以及σmax的值,如图6所示。当加载载荷增加至69.86N时,纤维截面最大正应力σmax达到纤维拉伸强度极限σt,纤维发生断裂,此时,缺陷边缘扩展位移为84.5μm,计算结束。
Claims (1)
1.一种复材分层缺陷在横刃挤压阶段扩展行为计算方法,其特征是,在机械加工时先确定刀具作用位置以及成型分层缺陷大小,根据约束情况及受力特点对缺陷层中的纤维进行区域划分,用该计算方法将刀具横刃对缺陷层的挤压过程简化为集中载荷作用下的弹性地基梁弯曲过程,求出梁的弯曲变形控制方程;再利用有限元计算方法及叠加原理,得到整根纤维的整体刚度矩阵;最后,采用向刀具作用点逐步加载方式对成型分层缺陷扩展的临界条件及扩展长度进行计算;方法的具体步骤如下:
1)首先确定刀具作用位置,输入被加工材料的性能参数和分层缺陷相关参数;
2)简化模型
将钻削刀具的横刃与工件相互作用的过程简化为集中载荷作用下的弹性地基梁模型,选取长度为dx的微小单元体,分析其受力情况,设pm为单位长度上周围材料对纤维的反作用力,qb为界面对纤维的粘结作用,有公式:
pm=kmw (1)
其中,w为挠度,km通过Boit’s方程求解,
其中,为周围材料的等效模量,与纤维、树脂的力学性能及含量相关;
qb=kbw (3)
式中,kb为界面的等效模量
其次,由微小单元体的平衡方程,得到单纤维温克尔地基梁的变形控制方程为:
当界面处存在缺陷时,则缺陷的下一层记作层S,缺陷区域为AB,载荷作用于点O,在缺陷区域纤维将不再受到界面的粘结作用,此时上述单纤维变形控制方程简化为:
3)求解缺陷层纤维整体刚度矩阵
输入刀具作用位置及分层缺陷大小,根据约束情况和受力特点,将缺陷层S层中的纤维划分为FB、BO、OA、AE四段;考虑到O点左右两侧结构相似,以O点左侧为例,采用有限元计算方法,对纤维的弯曲断裂及界面缺陷的扩展情况进行求解;对于OA段,根据简化后的变形控制方程求得挠度解为:
w(x)=eαx(B1cosβx+B2sinβx)+e-αx(B3cosβx+B4sinβx) (6)
式中,B1-B4为积分常数,且有:
梁单元内部的剪力和扭矩可分别表示为:
设OA段位移矢量dOA、载荷矢量rOA分别为:
dOA={wO,θO,wA,θA}T (10)
rOA={QO,MO,QA,MA}T (11)
则单元刚度矩阵kOA为:
对于AE段,由于点E为无穷远处,其挠度为0,即满足:
wE=w|x=+∞=0 (13)
则AE段挠度解简化为:
w(x)=e-αx(B3cosβx+B4sinβx) (14)
其中,B3-B4为积分常数,且有:
设AE段位移矢量dAE、载荷矢量rAE分别为:
dAE={wA,θA}T (16)
rAE={QA,MA}T (17)
则单元刚度矩阵kAE为:
综上,可得O点左侧纤维的整体载荷、位移、刚度之间的关系式:
Rleft=Kleftdleft (19)
其中,Kleft为单元刚度矩阵经叠加得到的O点左侧纤维整体刚度矩阵;采用类似的计算方法,对O点右侧纤维进行分析求解,最终确定整体纤维的载荷、位移、刚度之间的关系式:
R=Kd (20)
式中,整体载荷矢量R为:
R={QF,QB,MB,QO,MO,QA,MA,QE}T (21)
整体位移矢量d为:
矩阵K为根据FB、BO、OA、AE四段的单元刚度矩阵经叠加原理求得的整根纤维的整体刚度矩阵;
4)判断成型分层缺陷是否扩展
当求解出K之后,通过向刀具作用点O逐步施加载荷,计算出各节点位移;随着载荷不断增加,纤维的弯曲变形不断增加,界面应力若满足界面强度则点A(或点B)处界面将发生开裂,然而,在切削过程中,纤维的弯曲变形存在极限,当纤维内的截面的最大拉应力σmax达到其拉伸强度极限σt时,纤维将发生断裂;一旦纤维发生断裂,纤维的弯曲变形将不再增加,则点A(或点B)的位移也不再发生变化;因此,点A(或点B)处的界面若发生开裂,其一定发生在纤维弯曲断裂之前;结合以上分析,得到在横刃挤压阶段预埋缺陷扩展的临界条件为:
其中,设纤维半径为r,则纤维各截面最大拉应力为:
其中kb为层间界面的等效模量,σb为树脂粘结强度,wA为A点弯曲挠度,wB为B点弯曲挠度算法流程;若缺陷扩展则转至步骤5,若缺陷不扩展停止计算;
5)计算成型分层缺陷扩展长度
若满足缺陷扩展的临界条件,预埋缺陷将发生扩展,则迭代过程继续;扩展后每个增量步将逐渐增大BO或OA段长度,以及载荷矢量R的值,依据新的约束条件重新对梁内的各部分单元进行求解,更新缺陷扩展后纤维整体刚度矩阵,并求解wA、wB以及σmax的值;当σmax达到σt时,即纤维发生弯曲断裂,迭代停止,此时AB段的长度与迭代之前的长度之差即为缺陷扩展长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810039038.8A CN108388691B (zh) | 2018-01-16 | 2018-01-16 | 复材分层缺陷在横刃挤压阶段扩展行为计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810039038.8A CN108388691B (zh) | 2018-01-16 | 2018-01-16 | 复材分层缺陷在横刃挤压阶段扩展行为计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108388691A true CN108388691A (zh) | 2018-08-10 |
CN108388691B CN108388691B (zh) | 2020-04-28 |
Family
ID=63076953
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810039038.8A Expired - Fee Related CN108388691B (zh) | 2018-01-16 | 2018-01-16 | 复材分层缺陷在横刃挤压阶段扩展行为计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108388691B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110188446A (zh) * | 2019-05-24 | 2019-08-30 | 大连理工大学 | 一种考虑变形的复合材料板钻削分层临界轴向力计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014065069A (ja) * | 2012-09-27 | 2014-04-17 | Yasuda Kogyo Kk | 金属丸断面線材の製造方法、および金属丸断面線材、ニッケルフリーオーステナイトステンレス鋼丸断面線材、溶接線用・抵抗用・電熱線用丸断面線材 |
CN105891102A (zh) * | 2016-04-01 | 2016-08-24 | 江苏藤仓亨通光电有限公司 | 一种铜或铜合金接触线分层缺陷的检测方法 |
CN106503384A (zh) * | 2016-10-28 | 2017-03-15 | 大连理工大学 | 一种碳纤维复合材料表层切削损伤的综合抑制方法 |
CN106650119A (zh) * | 2016-12-27 | 2017-05-10 | 沈阳航空航天大学 | Cfrp与钛合金叠层结构钻孔刀具寿命的预测方法 |
CN106769853A (zh) * | 2016-12-27 | 2017-05-31 | 沈阳航空航天大学 | Cfrp与钛合金叠层结构钻孔分层缺陷的预测方法 |
-
2018
- 2018-01-16 CN CN201810039038.8A patent/CN108388691B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014065069A (ja) * | 2012-09-27 | 2014-04-17 | Yasuda Kogyo Kk | 金属丸断面線材の製造方法、および金属丸断面線材、ニッケルフリーオーステナイトステンレス鋼丸断面線材、溶接線用・抵抗用・電熱線用丸断面線材 |
CN105891102A (zh) * | 2016-04-01 | 2016-08-24 | 江苏藤仓亨通光电有限公司 | 一种铜或铜合金接触线分层缺陷的检测方法 |
CN106503384A (zh) * | 2016-10-28 | 2017-03-15 | 大连理工大学 | 一种碳纤维复合材料表层切削损伤的综合抑制方法 |
CN106650119A (zh) * | 2016-12-27 | 2017-05-10 | 沈阳航空航天大学 | Cfrp与钛合金叠层结构钻孔刀具寿命的预测方法 |
CN106769853A (zh) * | 2016-12-27 | 2017-05-31 | 沈阳航空航天大学 | Cfrp与钛合金叠层结构钻孔分层缺陷的预测方法 |
Non-Patent Citations (3)
Title |
---|
AM AMARO等: "Delamination Effect on Bending Behaviour in Carbon–Epoxy Composites", 《WILEY ONLINE LIBARARY》 * |
张璐: "含分层缺陷复合材料层合板分层扩展行为与数值模拟研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
朱炜垚: "含低速冲击损伤复合材料层板剩余强度及疲劳性能研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110188446A (zh) * | 2019-05-24 | 2019-08-30 | 大连理工大学 | 一种考虑变形的复合材料板钻削分层临界轴向力计算方法 |
CN110188446B (zh) * | 2019-05-24 | 2021-07-16 | 大连理工大学 | 一种考虑变形的复合材料板钻削分层临界轴向力计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108388691B (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Deng et al. | Isogeometric analysis of continuum damage in rotation-free composite shells | |
Dandekar et al. | Multiphase finite element modeling of machining unidirectional composites: prediction of debonding and fiber damage | |
Chedad et al. | A high-order closed-form solution for interfacial stresses in externally sandwich FGM plated RC beams | |
Xu et al. | On design of carbon fiber reinforced plastic (CFRP) laminated structure with different failure criteria | |
Jia et al. | Nonlinear viscoelastic multi-scale repetitive unit cell model of 3D woven composites with damage evolution | |
Alessi et al. | Damage evolution and debonding in hybrid laminates with a cohesive interfacial law | |
Wang et al. | An extended phantom node method study of crack propagation of composites under fatigue loading | |
CN108388691A (zh) | 复材分层缺陷在横刃挤压阶段扩展行为计算方法 | |
Maung et al. | Automated manufacture of optimised shape-adaptive composite hydrofoils with curvilinear fibre paths for improved bend-twist performance | |
Zhao et al. | Experimental study on tensile mechanics of arrow combination structure with carbon fiber–epoxy resin composite | |
Liu et al. | Micro model of carbon fiber/cyanate ester composites and analysis of machining damage mechanism | |
Xie et al. | Tunable properties and responses of architected lattice-reinforced cementitious composite components induced by versatile cell topology and distributions | |
Naqi et al. | Nonlinear finite element analysis of RCMD beams with large circular opening strengthened with cfrp material | |
Zhao et al. | Mechanical performance of aluminum reinforced wood plastic composites under axial tension: an experimental and numerical investigation | |
Cai et al. | A new interfacial model for transverse mechanical properties of unidirectional fiber reinforced composites | |
Prabhakaran et al. | Modelling and simulation of natural fibre/epoxy composites-prediction of stress state and deformations | |
Heise et al. | Investigation of Damage Mechanisms in Composites due to the Interaction of Ply and Fiber Waviness Defects | |
Wang et al. | Numerical investigation of two-dimensional Mode-II delamination in composite laminates | |
Vermes et al. | Thermally induced mechanical work and warpage compensation of asymmetric laminates | |
Dronamraju et al. | Studies of fiber-matrix debonding | |
Lei et al. | Macroscopic mechanical characterization of SMAs fiber-reinforced hybrid composite under uniaxial loading | |
Chinta et al. | Investigation of fracture toughness of bidirectional jute/epoxy composite and analysis by using FEA | |
Schaefer | Matrix dominated failure of fiber-reinforced composite laminates under static and dynamic loading | |
Juliyana et al. | A Methodology for Enhancing the Shear Response of Sandwich Composite Panel: Sandwich Composite Panel with Stair Keys | |
Sahu | Design and analysis of thin walled composite beam |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200428 |