CN112926244A - 一种复合材料层合板开孔件极限载荷确定方法 - Google Patents
一种复合材料层合板开孔件极限载荷确定方法 Download PDFInfo
- Publication number
- CN112926244A CN112926244A CN202110181638.XA CN202110181638A CN112926244A CN 112926244 A CN112926244 A CN 112926244A CN 202110181638 A CN202110181638 A CN 202110181638A CN 112926244 A CN112926244 A CN 112926244A
- Authority
- CN
- China
- Prior art keywords
- damage
- plane
- composite material
- laminated plate
- fiber
- 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.)
- Pending
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000000463 material Substances 0.000 claims abstract description 10
- 230000015556 catabolic process Effects 0.000 claims abstract description 7
- 238000006731 degradation reaction Methods 0.000 claims abstract description 7
- 238000010008 shearing Methods 0.000 claims abstract description 7
- 239000000835 fiber Substances 0.000 claims description 47
- 238000006073 displacement reaction Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 20
- 230000006835 compression Effects 0.000 claims description 11
- 238000007906 compression Methods 0.000 claims description 11
- 239000011229 interlayer Substances 0.000 claims description 8
- 239000010410 layer Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000032798 delamination Effects 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 6
- 230000000977 initiatory effect Effects 0.000 claims description 5
- 102100036790 Tubulin beta-3 chain Human genes 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000008602 contraction Effects 0.000 claims description 2
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 2
- 239000000758 substrate Substances 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 abstract description 2
- 230000009022 nonlinear effect Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 102100021503 ATP-binding cassette sub-family B member 6 Human genes 0.000 description 1
- 101100000375 Homo sapiens ABCB6 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
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/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Laminated Bodies (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种复合材料层合板开孔件极限载荷确定方法,包括如下步骤:建立复合材料层合板开孔件有限元模型;建立复合材料损伤本构模型;基于ABAQUS‑UMAT有限元用户动态子程序模块,使用FORTRAN语言编写用户自定义子程序实现步骤SS2提出的损伤本构模型,求解应力、应变和损伤;对所述有限元模型进行计算,预测复合材料层合板开孔件极限载荷。本发明利用ABAQUS‑UMAT用户自定义子程序来数值实现所建立的三维损伤本构模型,该模型同时考虑了剪切非线性和损伤累积导致材料性能退化的影响,能准则预测复合材料的极限载荷。
Description
技术领域
本发明涉及一种复合材料层合板开孔件极限载荷确定方法,属于复合材料力学性能分析技术领域。
背景技术
在航空领域,飞行器上的结构关于复合材料的应用已十分普遍。纤维增强树脂基复合材料具有高比模量和比强度、优良的能量吸收性能,特别是各向刚度与强度可设计性等特点,广泛应用在航空航天、军事、海洋、土木和机械等工程领域。近年来针对复合材料层合板提出的损伤本构模型未能合理考虑复合材料存在的非线性力学行为、面外应力的影响、合适的失效准则及损伤演化。现有研究表明,如果模型未考虑非线性效应可能会低估复合材料结构的能量吸收能力;ABAQUS软件使用二维Hashin准则来判断层内损伤的起始萌生,二维Hashin准则虽然可以区分纤维与基体的不同损伤模式,但无法说明损伤产生的物理机制,对于适当的横向压缩能抑制剪切破坏发生的现象无法给出较为合理的解释。因此,如何准确预测试验件的极限载荷,需要运用较有优势的三维Hashin失效准则和Ye分层失效准则,对有非线性力学行为的复合材料进行极限载荷的预测。
发明内容
本发明的目的在于,克服现有技术存在的技术缺陷,解决上述技术问题,提出一种复合材料层合板开孔件极限载荷确定方法。
本发明具体采用如下技术方案:一种复合材料层合板开孔件极限载荷确定方法,包括如下步骤:
步骤SS1:建立复合材料层合板开孔件有限元模型;
步骤SS2:建立复合材料损伤本构模型;
步骤SS3:基于ABAQUS-UMAT有限元用户动态子程序模块,使用FORTRAN语言编写用户自定义子程序实现步骤SS2提出的损伤本构模型,求解应力、应变和损伤;
步骤SS4:对所述有限元模型进行计算,预测复合材料层合板开孔件极限载荷。
作为一种较佳的实施例,所述步骤SS1具体包括:
复合材料层合板的铺层角度沿厚度方向中平面对称布置,并且每一层厚度方向只划分一个单元;
网格类型为C3D8R,对孔边周围区域进行网格细化;
建立参考点与自由端面之间加载方向位移一致性约束条件:受拉荷载均采用位移加载方式,左加载面施加固支约束,右端自由端面外设置一个参考点,然后把参考点和端面进行绑定,在Abaqus/CAE模块中,采用creat constraint方法建立coupling耦合约束方程,此时,将位移荷载施加在参考点上,同时只要输出参考点上的位移和反力即U和RF1,就能够获得加载端面上的位移与反力。
作为一种较佳的实施例,所述步骤SS2具体包括:
步骤SS21:建立含损伤的复合材料层合板本构关系;
步骤SS22:建立三维Hashin强度失效准则来判断纤维和基体损伤,Ye分层失效准则来判断分层损伤;
步骤SS23:建立剪切非线性模型;
步骤SS24:建立连续损伤退化模型。
作为一种较佳的实施例,所述步骤SS21具体包括:
其中,符号“:”表示对两个张量指标的缩并计算;σ是有效应力张量;是名义应力张量;εe是弹性应变张量;e表示弹性;C(d)是含损伤单向复合材料层合板的四阶刚度张量;C是未损伤单向复合材料层合板的四阶线弹性刚度张量;d是一维向量,(d1,d2,d3,d23,d13,d12),其中d1、d2、d3分别为纤维方向纤维损伤的损伤变量、平面内垂直于纤维方向基体损伤的损伤变量、层间平面外方向分层损伤的损伤变量;d12、d23、d13分别为12、23、13平面内的剪切损伤变量;定义坐标系x1-x2-x3为单向板的自然坐标系,x1-xn-x1为断裂面的局部坐标系,两个坐标系下的x1轴重合;12、23、13平面分别对应坐标系下的x1x2平面、x2x3平面、x1x3平面。
作为一种较佳的实施例,所述步骤SS21具体还包括:
把损伤变量引入刚度矩阵,使刚度随着损伤的发展而逐渐变弱,即:
C(d)=M-1(d):C:MT,-1(d);
其中,M-1(d)为M(d)的逆矩阵,MT,-1(d)为M(d)转置矩阵的逆矩阵;M(d)为损伤因子张量,其损伤主轴系下矩阵形式可以表示如下:
复合材料主坐标系中单层板的三维正交各项异性损伤本构模型如下:
所述复合材料主坐标系为单向板的自然坐标系x1-x2-x3;
其中:
其中,σ1、σ2和σ3分别为纤维方向、垂直于纤维方向、层间平面外方向的名义正应力;τ23、τ12和τ13分别为x1-x2-x3坐标系下x1x2平面、x2x3平面、x1x3平面内的剪应力;ε1、ε2和ε3分别为纤维方向、垂直于纤维方向、层间平面外方向的工程正应变;γ23、γ13和γ12为分别为x1-x2-x3坐标系下x1x2平面、x2x3平面、x1x3平面内的工程剪应变;E1、E2、E3、分别为纤维方向、垂直于纤维方向、层间平面外方向的未损伤单向复合材料单层的弹性模量,G23、G13、G12分别为x1x2平面、x2x3平面、x1x3平面内未损伤单向复合材料单层的剪切模量,v12、v13、v23分别为纤维方向与垂直于纤维方向的、纤维方向与层间平面外方向的、垂直于纤维方向与层间平面外方向的泊松比v21、v31、v32分别为垂直于纤维方向与纤维方向的、层间平面外方向与纤维方向的、层间平面外方向与垂直于纤维方向的泊松比,满足关系式,
作为一种较佳的实施例,所述步骤SS22具体包括:
(a)对于纤维拉伸和压缩,损伤初始判据为:
纤维拉伸失效(ε11≥0):
纤维压缩失效(ε11<0):
(b)对于基体拉伸和压缩,损伤初始判据为:
基体拉伸失效(ε22+ε33≥0):
基体压缩失效(ε22+ε33≥0):
(c)分层损伤初始判据:
拉伸引起分层失效(ε33≥0):
压缩引起分层失效(ε33<0):
式中:fi(i=1,2,3)分别代表纤维、基体、分层的损伤状态;Cii表示材料的刚度系数;分别代表纤维在i方向上拉伸强度和压缩强度对应的正应变;γ12、γ13、γ23分别表示各平面剪切强度对应的剪应变;XT、XC分别为单向板沿纤维方向拉伸及压缩强度;YT、YC分别为单向板横向拉伸及压缩强度;ZT为法向拉伸强度;S12、S13、S23为相应平面的剪切强度。
作为一种较佳的实施例,所述步骤SS23具体包括:
所述剪切非线性模型具体建立方式为:
考虑剪切非线性后的剪切模量G的表达式为:
式中的τ和γ分别为剪应变和剪应力,G0为初始剪切模量,τ0为极限剪切强度,n是定义了剪切非线性关系曲线形状的参数。
作为一种较佳的实施例,所述步骤SS24具体包括:
所述连续损伤退化模型的具体建立方式为:
作为一种较佳的实施例,所述步骤SS3具体包括:
开始当前增量步,读取前一时刻收敛状态量及当前增量步中应变增量,更新应变和有效应力;
根据有效应力代入到步骤SS2的步骤SS22、步骤SS23判断是否有损伤出现,如有损伤出现再通过步骤SS2的步骤SS24更新损伤变量,再通过有效应力和损伤变量计算名义应力。
作为一种较佳的实施例,所述步骤SS4具体包括:
将步骤SS1建立的复合材料层合板有限元模型文件和步骤SS3建立的ABAQUS-UMAT用户子程序联合,完成对复合材料层合板失效强度的预测;先在ABAQUS软件里面建立复合材料层合板有限元模型,然后调用编写好的子程序进行应力应变分析,最后得到的载荷位移曲线就是该模型的力学行为反应,得到的最大值为极限载荷。
本发明所达到的有益效果:第一,对于ABAQUS内嵌的基于二维Hashin失效准则的二维弹性损伤本构,本发明基于更准确的三维Hashin失效准则和Ye分层失效准则连续退化模型以及包含剪切非线性效应的三维损伤本构,更符合工程实际。第二,本发明利用ABAQUS-UMAT用户自定义子程序来数值实现所建立的三维损伤本构模型,具有更高的计算效率和计算精度。
附图说明
图1是本发明的一种复合材料层合板开孔件极限载荷确定方法的流程图;
图2是ABAQUS用户材料的子程序流程图;
图3是试件几何尺寸图;
图4是利用本文确定方法计算得到的层合板的仿真与试验载荷位移曲线对比图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
在ABAQUS中建立有限元模型。所试验的试件为含缺口的复合材料层合板,层合板的几何尺寸和纤维铺层顺序见表1。材料弹性常数和强度参数见表2和表3。
表1复合材料层合板试件几何尺寸和铺层顺序表
表2复合材料材料弹性系数表
表3复合材料强度参数表
对图3何模型进行有限元网格划分,网格类型为C3D8R,属性中选择了增强沙漏控制。
模型的边界条件为:一端固支,另一端采用参考点与面耦合的方式。对参考点施加拉伸载荷,并最终输出参考点的载荷位移曲线图。
载荷采用施加位移载荷的形式。位移载荷须适中,若载荷过小,则无法达到材料的失效强度极限,无法判定材料的失效起始点和最终失效载荷:若载荷过大,则可能导致计算不收敛。本例中拟对模型采用2mm的位移载荷。
利用ABAQUS/STANDRAD计算模拟复合材料层合板的拉伸失效过程,利用用户子程序UMAT读取当前应变增量,更新应变和有效应力,根据有效应力判别单元是否进入损伤,进入损伤阶段则根据损伤演化模型计算损伤变量,从而得到名义应力,最后可得到模型的载荷-位移曲线。
图4为利用本文计算方法得到的层合板的仿真与试验载荷位移曲线对比图。从图4中可以发现,仿真的载荷-位移曲线中的非线性关系拟合较好。
本发明是在ABAQUS软件的基础上进行的用户子程序的开发,提出的三维损伤本构模型同时考虑了剪切非线性效应和损伤累积导致的材料性能退化的影响,能理想的预测复合材料层合板的失效强度,为深入阐明复合材料结构的损伤失效特性,提升轻量化强度设计水平提供了技术支撑。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (10)
1.一种复合材料层合板开孔件极限载荷确定方法,其特征在于,包括如下步骤:
步骤SS1:建立复合材料层合板开孔件有限元模型;
步骤SS2:建立复合材料损伤本构模型;
步骤SS3:基于ABAQUS-UMAT有限元用户动态子程序模块,使用FORTRAN语言编写用户自定义子程序实现步骤SS2提出的损伤本构模型,求解应力、应变和损伤;
步骤SS4:对所述有限元模型进行计算,预测复合材料层合板开孔件极限载荷。
2.根据权利要求1所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS1具体包括:
复合材料层合板的铺层角度沿厚度方向中平面对称布置,并且每一层厚度方向只划分一个单元;
网格类型为C3D8R,对孔边周围区域进行网格细化;
建立参考点与自由端面之间加载方向位移一致性约束条件:受拉荷载均采用位移加载方式,左加载面施加固支约束,右端自由端面外设置一个参考点,然后把参考点和端面进行绑定,在Abaqus/CAE模块中,采用creat constraint方法建立coupling耦合约束方程,此时,将位移荷载施加在参考点上,同时只要输出参考点上的位移U和反力RF1,就能够获得加载端面上的位移与反力。
3.根据权利要求1所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS2具体包括:
步骤SS21:建立含损伤的复合材料层合板本构关系;
步骤SS22:建立三维Hashin强度失效准则来判断纤维和基体损伤,Ye分层失效准则来判断分层损伤;
步骤SS23:建立剪切非线性模型;
步骤SS24:建立连续损伤退化模型。
4.根据权利要求3所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS21具体包括:
其中,符号“:”表示对两个张量指标的缩并计算;σ是有效应力张量;是名义应力张量;εe是弹性应变张量;e表示弹性;C(d)是含损伤单向复合材料层合板的四阶刚度张量;C是未损伤单向复合材料层合板的四阶线弹性刚度张量;d是一维向量(d1,d2,d3,d23,d13,d12),其中d1、d2、d3分别为纤维方向纤维损伤的损伤变量、平面内垂直于纤维方向基体损伤的损伤变量、层间平面外方向分层损伤的损伤变量;d12、d23、d13分别为12、23、13平面内的剪切损伤变量;定义坐标系x1-x2-x3为单向板的自然坐标系,x1-xn-x1为断裂面的局部坐标系,两个坐标系下的x1轴重合;12、23、13平面分别对应坐标系下的x1x2平面、x2x3平面、x1x3平面。
5.根据权利要求4所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS21具体还包括:
把损伤变量引入刚度矩阵,使刚度随着损伤的发展而逐渐变弱,即:
C(d)=M-1(d):C:MT,-1(d);
其中,M-1(d)为M(d)的逆矩阵,MT,-1(d)为M(d)转置矩阵的逆矩阵;M(d)为损伤因子张量,其损伤主轴系下矩阵形式可以表示如下:
复合材料主坐标系中单层板的三维正交各项异性损伤本构模型如下:
所述复合材料主坐标系为单向板的自然坐标系x1-x2-x3;
其中:
其中,σ1、σ2和σ3分别为纤维方向、垂直于纤维方向、层间平面外方向的名义正应力;τ23、τ12和τ13分别为x1-x2-x3坐标系下x1x2平面、x2x3平面、x1x3平面内的剪应力;ε1、ε2和ε3分别为纤维方向、垂直于纤维方向、层间平面外方向的工程正应变;γ23、γ13和γ12为分别为x1-x2-x3坐标系下x1x2平面、x2x3平面、x1x3平面内的工程剪应变;E1、E2、E3、分别为纤维方向、垂直于纤维方向、层间平面外方向的未损伤单向复合材料单层的弹性模量,G23、G13、G12分别为x1x2平面、x2x3平面、x1x3平面内未损伤单向复合材料单层的剪切模量,v12、v13、v23分别为纤维方向与垂直于纤维方向的、纤维方向与层间平面外方向的、垂直于纤维方向与层间平面外方向的泊松比v21、v31、v32分别为垂直于纤维方向与纤维方向的、层间平面外方向与纤维方向的、层间平面外方向与垂直于纤维方向的泊松比,满足关系式:
6.根据权利要求3所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS22具体包括:
(a)对于纤维拉伸和压缩,损伤初始判据为:
纤维拉伸失效(ε11≥0):
纤维压缩失效(ε11<0):
(b)对于基体拉伸和压缩,损伤初始判据为:
基体拉伸失效(ε22+ε33≥0):
基体压缩失效(ε22+ε33≥0):
(c)分层损伤初始判据:
拉伸引起分层失效(ε33≥0):
压缩引起分层失效(ε33<0):
9.根据权利要求3所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS3具体包括:
开始当前增量步,读取前一时刻收敛状态量及当前增量步中应变增量,更新应变和有效应力;
根据有效应力代入到步骤SS2的步骤SS22、步骤SS23判断是否有损伤出现,如有损伤出现再通过步骤SS2的步骤SS24更新损伤变量,再通过有效应力和损伤变量计算名义应力。
10.根据权利要求1所述的一种复合材料层合板开孔件极限载荷确定方法,其特征在于,所述步骤SS4具体包括:
将步骤SS1建立的复合材料层合板有限元模型文件和步骤SS3建立的ABAQUS-UMAT用户子程序联合,完成对复合材料层合板失效强度的预测;先在ABAQUS软件里面建立复合材料层合板有限元模型,然后调用编写好的子程序进行应力应变分析,最后得到的载荷位移曲线就是该模型的力学行为反应,得到的最大值为极限载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110181638.XA CN112926244A (zh) | 2021-02-10 | 2021-02-10 | 一种复合材料层合板开孔件极限载荷确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110181638.XA CN112926244A (zh) | 2021-02-10 | 2021-02-10 | 一种复合材料层合板开孔件极限载荷确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112926244A true CN112926244A (zh) | 2021-06-08 |
Family
ID=76169672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110181638.XA Pending CN112926244A (zh) | 2021-02-10 | 2021-02-10 | 一种复合材料层合板开孔件极限载荷确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112926244A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705047A (zh) * | 2021-08-25 | 2021-11-26 | 南京航空航天大学 | 一种拉伸载荷下cfrtp干涉连接孔周应力分布计算方法 |
CN115906570A (zh) * | 2022-11-22 | 2023-04-04 | 大连理工大学 | 一种基于数据驱动的复合材料开孔层合板渐进损伤预测方法 |
-
2021
- 2021-02-10 CN CN202110181638.XA patent/CN112926244A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113705047A (zh) * | 2021-08-25 | 2021-11-26 | 南京航空航天大学 | 一种拉伸载荷下cfrtp干涉连接孔周应力分布计算方法 |
CN113705047B (zh) * | 2021-08-25 | 2024-04-05 | 南京航空航天大学 | 一种拉伸载荷下cfrtp干涉连接孔周应力分布计算方法 |
CN115906570A (zh) * | 2022-11-22 | 2023-04-04 | 大连理工大学 | 一种基于数据驱动的复合材料开孔层合板渐进损伤预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106503292B (zh) | 预测低速冲击下复合材料层合板渐进失效的有限元方法 | |
CN111368389B (zh) | 一种预测复合材料层合板失效强度的方法 | |
CN106777769B (zh) | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 | |
Reddy et al. | Bending analysis of laminated composite plates using finite element method | |
Mines et al. | Numerical simulation of the progressive collapse of polymer composite sandwich beams under static loading | |
Liao et al. | Multi-scale modelling of dynamic progressive failure in composite laminates subjected to low velocity impact | |
Dadej et al. | On the effect of glass and carbon fiber hybridization in fiber metal laminates: Analytical, numerical and experimental investigation | |
Hirwani et al. | Nonlinear finite element solutions of thermoelastic deflection and stress responses of internally damaged curved panel structure | |
Juhász et al. | The effect of delamination on the critical buckling force of composite plates: Experiment and simulation | |
CN109101692B (zh) | 基于最大应力准则的复合材料层合板极限载荷计算方法 | |
Grover et al. | An efficient C0 finite element modeling of an inverse hyperbolic shear deformation theory for the flexural and stability analysis of laminated composite and sandwich plates | |
Hirwani et al. | Nonlinear transient analysis of delaminated curved composite structure under blast/pulse load | |
McElroy | Use of an enriched shell finite element to simulate delamination-migration in a composite laminate | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
McElroy et al. | Simulation of delamination–migration and core crushing in a CFRP sandwich structure | |
Xiaohui et al. | A new zig-zag theory and C 0 plate bending element for composite and sandwich plates | |
CN107688686A (zh) | 一种基于弹脆性损伤模式的层压板结构破坏模拟方法 | |
CN114139300A (zh) | 一种混杂纤维复合材料夹芯板抗冲击性能评价方法 | |
Shao et al. | Numerical comparison between Hashin and Chang-Chang failure criteria in terms of inter-laminar damage behavior of laminated composite | |
Elangovan et al. | Dynamic characterization of tapered composite sandwich plate with honeycomb core: Numerical and experimental investigations | |
Geng et al. | An interface damage model for high-cycle fatigue | |
Salami et al. | An advanced high-order theory for bending analysis of moderately thick faced sandwich beams | |
Dinulović et al. | 3D random fiber composites as a repair material for damaged honeycomb cores | |
CN113158508A (zh) | 一种复合材料层合板贴补修复后极限载荷确定方法 | |
Ginot et al. | Local buckling on large sandwich panels applied to light aviation: Experimental and computation dialogue |
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 |