CN111709174B - 一种基于失效面理论的复合材料层合板强度分析方法 - Google Patents
一种基于失效面理论的复合材料层合板强度分析方法 Download PDFInfo
- Publication number
- CN111709174B CN111709174B CN202010559632.7A CN202010559632A CN111709174B CN 111709174 B CN111709174 B CN 111709174B CN 202010559632 A CN202010559632 A CN 202010559632A CN 111709174 B CN111709174 B CN 111709174B
- Authority
- CN
- China
- Prior art keywords
- failure
- composite material
- coordinate system
- laminated plate
- damage
- 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.)
- Active
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 127
- 238000004458 analytical method Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 32
- 239000000463 material Substances 0.000 claims abstract description 27
- 230000015556 catabolic process Effects 0.000 claims abstract description 9
- 238000006731 degradation reaction Methods 0.000 claims abstract description 9
- 238000011161 development Methods 0.000 claims abstract description 6
- 230000000977 initiatory effect Effects 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 239000000835 fiber Substances 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 7
- 230000000750 progressive effect Effects 0.000 claims description 6
- 206010010214 Compression fracture Diseases 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims 1
- 230000003902 lesion Effects 0.000 claims 1
- 238000013461 design Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000011165 3D composite Substances 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000010959 steel 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于失效面理论的复合材料层合板强度分析方法,涉及复合材料结构失效分析领域,该方法包括:基于三维失效面理论确定材料损伤柔度矩阵的计算公式,从而建立起复合材料力学本构关系;采用失效面上的应力构造失效准则以判断复合材料是否失效及其失效模式,并基于失效面坐标系进行复合材料的刚度退化。通过对有限元软件进行二次开发建立了复合材料层合板的三维有限元模型,仿真模拟了复合材料层合板损伤起始、损伤演化和最终破坏的完整过程。该方法物理机制明确,能较好地预测复合材料层合板的失效载荷与失效模式,在很大程度上提高了复合材料层合板强度分析的预测精度,从而避免进行大量耗时长、成本高的试验测试。
Description
技术领域
本发明涉及复合材料结构失效分析领域,尤其是一种基于失效面理论的复合材料层合板强度分析方法。
背景技术
复合材料凭借其高比强度、高比刚度、可设计、抗疲劳和耐腐蚀等突出特点,在航空航天、风电叶片、轨道交通等领域得到了广泛应用。然而在工程实际中,复合材料结构的设计仍是相当保守的。工程人员在设计时一般会通过选取较大的安全系数,并进行大量试验以最大程度地保证结构的安全性。但采用该方法会导致设计的结构效率偏低,也会造成设计周期变长,成本增加。因此,有必要建立一种准确预测复合材料层合板强度的分析方法,这有助于提高结构效率,减轻结构重量,降低试验成本并缩短研发周期。
复合材料层合板的渐进失效分析(Progressive Failure Analysis,PFA)方法目前在学术界得到了广泛应用,并取得了相对较好的预测效果。但因为复合材料具有高度各向异性且损伤模式多样化,导致其失效机理异常复杂,所以不少学者提出的复合材料层合板渐进失效分析模型都带有一定经验性。例如,广泛用于预测复合材料层合板失效的Tsai-Wu准则中系数F12的取值没有公认的确定方法,基于应力不变量理论的Hashin准则中的一些参数因缺少实验数据而被人为设置成0。此外,在失效发生后复合材料刚度性能的折减方式也不尽相同。这些问题给建立一种基于物理机制且预测精度较高的复合材料层合板强度分析方法带来了挑战。
发明内容
本发明人针对上述问题及技术需求,提出了一种基于失效面理论的复合材料层合板强度分析方法,该方法采用失效面上的应力构造失效准则以判断复合材料是否失效及其失效模式,并基于失效面坐标系进行复合材料刚度退化。该方法物理机制明确,能较好地预测复合材料层合板的失效全过程,在很大程度上提高了复合材料强度分析的预测精度。
本发明的技术方案如下:
一种基于失效面理论的复合材料层合板强度分析方法,包括如下步骤:
建立基于三维失效面理论的复合材料力学本构关系,包括:
复合材料在未损伤和损伤这两种状态下的应力—应变本构方程分别为:ε=S0σ和ε=Sdσ,其中ε为自然坐标系下的应变,σ为自然坐标系下的应力向量,S0为未损伤状态下的柔度矩阵,S0通过试验测得的复合材料工程弹性常数确定;Sd为损伤状态下的柔度矩阵,Sd与复合材料损伤程度有关,并且Sd由三维失效面理论推导得到,包括:
首先采用基于三维失效面理论构造的失效准则来判断复合材料是否发生损伤;根据三维失效面理论,材料失效是由失效面上的应力决定的,因此先通过坐标转换矩阵T将应力向量σ=[σ1 σ2 σ3 τ23 τ13 τ21]T由自然坐标系转换到三维失效面坐标系:σfp=T-1σ,其中σfp=[σl σn σt τnt τlt τnl]T,代表失效面坐标系下的真实应力向量,然后将失效面坐标系下的应力分量σl,σn,τnt,τnl代入基于三维失效面理论构造的失效准则以判断复合材料是否失效及其失效模式;
如果复合材料满足失效准则,则发生损伤,需进行刚度性能退化;根据三维失效面理论,损伤变量矩阵Dfp同样应该定义在失效面坐标系下,而非传统的自然坐标系;由损伤力学知识可知,失效面上的有效应力和真实应力σfp满足关系式:/>基于能量等效假设,在失效面坐标系中未损伤和损伤这两种状态下的柔度矩阵满足关系式/>最后通过坐标转换矩阵T将失效面坐标系中的两柔度矩阵由失效面坐标系转换回自然坐标系,得到Sd=(TT)-1(Dfp)TTTS0TDfpT-1;在复合材料发生损伤后,应力—应变本构关系由ε=Sdσ确定;
如果复合材料不满足失效准则,则未发生损伤,材料刚度性能不变,复合材料的应力—应变本构关系由ε=S0σ确定;
基于有限元软件二次开发的复合材料层合板三维建模,包括:
根据复合材料结构的几何特征,采用参数化方法建立复合材料层合板的三维有限元模型,设置材料属性,再划分网格,并施加载荷和边界条件;
复合材料层合板的渐进失效分析,包括:
对复合材料层合板的三维有限元模型进行应力分析,模拟复合材料层合板的损伤起始、损伤演化和最终破坏的全过程,在每个载荷增量步,采用基于三维失效面理论构造的失效准则判断材料是否损伤,对失效的材料基于失效面坐标系进行刚度退化,然后继续施加载荷增量,直至确定复合材料层合板最终破坏;
根据有限元数值分析结果得到复合材料层合板的失效载荷与失效模式。
其进一步的技术方案为,首先采用基于三维失效面理论构造的失效准则来判断复合材料是否发生损伤,包括:
失效准则及其失效模式为:
①纤维拉伸失效:
②纤维压缩失效:
③基体拉伸失效:
④基体压缩失效:
式中XT是复合材料的纤维纵向拉伸,XC是复合材料的纤维压缩强度,ΥT是复合材料的基体横向拉伸,γC是复合材料的基体压缩强度,S21是纵向剪切强度,为复合材料横向压缩断裂角度。
其进一步的技术方案为,对失效的材料基于失效面坐标系进行刚度退化,包括:
在失效面坐标系中定义损伤变量矩阵Dfp:
式中di(i=l,n,t)和dij(i,j=n,t,l)是对应不同失效模式的损伤变量,若材料未发生失效,则损伤变量的值都为0;基于能量等效假设,最终得到损伤状态下的柔度矩阵Sd=(TT)-1(Dfp)TTTS0TDfpT-1,式中S0为复合材料未损伤状态下的柔度矩阵,T为坐标转换矩阵,对应的计算公式为:
其中θ表示失效面坐标系t方向与自然坐标系对应的方向之间的夹角。
本发明的有益技术效果是:
1、按照本申请所提供的方法,基于三维失效面理论建立复合材料的力学本构关系,其物理机制更加明确。
2、本方法采用失效面上的应力构造失效准则以判断复合材料是否失效及其具体失效模式,并基于失效面坐标系进行复合材料的刚度退化。与现有方法相比,能够更为准确地预测复合材料层合板的损伤起始和损伤后的力学性能。
3、按照本申请所提供的方法,建立复合材料层合板的三维有限元模型,较二维模型预测精度更高。
4、本申请能模拟复合材料层合板损伤起始、演化和最终失效的全过程,可较为准确地预测复合材料层合板的失效载荷与失效模式。
附图说明
图1是本申请公开的基于三维失效面理论的复合材料层合板强度分析方法流程图。
图2是复合材料失效面的示意图。
图3是一实施例提供的复合材料层合板的几何尺寸示意图。
图4是一实施例提供的复合材料层合板的三维有限元模型。
图5是一实施例提供的复合材料层合板的应力—应变曲线。
图6是一实施例提供的复合材料层合板的基体失效模式。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。
本申请公开了一种基于三维失效面理论的复合材料层合板强度分析方法,其方法流程图如图1所示,该方法包括如下步骤:
步骤S1:建立基于三维失效面理论的复合材料力学本构关系。
复合材料在未损伤和损伤这两种状态下的应力—应变本构方程分别为:ε=S0σ和ε=Sdσ,ε为自然坐标系下的应变,σ为自然坐标系下的应力向量,S0为未损伤状态下的柔度矩阵,S0通过试验测得的复合材料工程弹性常数确定;Sd为损伤状态下的柔度矩阵,Sd与复合材料损伤程度有关,在本方法中需基于三维失效面理论推导得到,包括:
首先采用基于三维失效面理论构造的失效准则来判断复合材料是否发生损伤。如图2所示,根据失效面理论,材料失效是由失效面(Failure Plane)上的应力决定的。因此先通过坐标转换矩阵T将应力向量σ=[σ1 σ2 σ3 τ23 τ13 τ21]T由自然坐标系(1-2-3坐标系)转换到失效面坐标系(l-n-t坐标系):σfp=T-1σ,其中失效面坐标系中的应力向量为σfp=[σlσn σt τnt τlt τnl]T,σfp代表失效面坐标系下的真实应力向量,
θ表示失效面坐标系t方向与自然坐标系3方向之间的夹角。
然后将失效面坐标系下的应力分量σl,σn,τnt,τnl代入基于三维失效面构造的失效准则中以判断复合材料是否失效及其失效模式:
①纤维拉伸失效:
②纤维压缩失效:
③基体拉伸失效:
④基体压缩失效:式中XT是复合材料的纤维纵向拉伸,XC是复合材料的纤维压缩强度,YT是复合材料的基体横向拉伸,YC是复合材料的基体压缩强度,S21是纵向剪切强度,/>复合材料横向压缩断裂角度,根据试验测量结果通常取53°。
如果复合材料满足失效准则,则发生损伤,需进行刚度性能退化,包括:
根据三维失效面理论,损伤变量矩阵Dfp同样应该定义在失效面坐标系下:
式中di(i=l,n,t)和dij(i,j=n,t,l)是对应不同失效模式的损伤变量,若材料未发生失效,则损伤变量的值都为0。
由损伤力学知识可知,失效面上的有效应力和真实应力σfp满足关系式:基于能量等效假设,在失效面坐标系中未损伤和损伤这两种状态下的柔度矩阵满足关系式/>最后通过坐标转换矩阵T将失效面坐标系中的两柔度矩阵由失效面坐标系转换回自然坐标系,得到Sd=(TT)-1(Dfp)TTTS0TDfpT-1。在复合材料发生损伤后,应力—应变本构关系由ε=Sdσ确定。
如果复合材料不满足失效准则,则未发生损伤,材料刚度不变。复合材料的应力—应变本构关系由ε=S0σ确定。
上述基于三维失效面理论的复合材料力学本构关系通过在有限元软件Abaqus中编写相应的材料子程序VUMAT实现。
步骤S2:基于有限元软件二次开发的复合材料层合板三维建模。
在有限元软件Abaqus中进行二次开发,使用Python语言编写自动化建模脚本程序。在实施例中,受单向拉伸载荷的复合材料层合板由IM7/8552复合材料制成,材料工程常数为:E1=171.42GPa,E2=9.08GPa,G12=5.29GPa,υ12=0.32,XT=2326.2MPa,XC=1200.1MPa,YT=62.3MPa,YC=199.8MPa,S21=92.3MPa。复合材料层合板单层厚度为0.131mm,铺层顺序为[90/0/±45]3s。如图3所示,复合材料层合板的尺寸为150mm×48mm×3.144mm,圆孔直径为8mm。将上述所有尺寸参数定义为几何模型的特征参数,采用参数化建模方法生成的复合材料层合板三维有限元模型如图4所示。然后设置材料属性,再划分网格,并施加载荷和边界条件。本实施例中的复合材料层合板一端为固支边界条件,在另一端施加拉伸载荷。
步骤S3:复合材料层合板的渐进失效分析。
对复合材料层合板的三维有限元模型进行应力分析,模拟复合材料层合板损伤起始、演化和最终失效的全过程。如图1所示,在每个载荷增量步,自然坐标系下的应力向量σ=[σ1 σ2 σ3 τ23 τ13 τ21]T是输入参数,应力向量σ通过坐标转换矩阵T可求得失效面上的真实应力向量σfp=[σl σn σt τnt τlt τnl]T。然后将应力分量σl,σn,τnt,τnl和材料的强度参数XT、XC、YT、YC、S21及代入基于三维失效面理论构造的失效准则来判断复合材料是否发生损伤,若材料发生损伤,则对失效的材料基于失效面坐标系进行刚度退化,并判断结构是否完全失效,若没有完全失效,则继续施加载荷增量并重新执行对复合材料层合板的三维有限元模型进行应力分析的步骤,若完全失效,则根据有限元分析结果得到复合材料层合板的失效载荷与失效模式。若材料未发生损伤,则执行继续施加载荷增量并重新执行对复合材料层合板的三维有限元模型进行应力分析的步骤,直至确定结构最终失效。
按照相关经验,采用显示算法进行准静态有限元分析,可以满足预测精度要求,并避免隐式算法收敛性差的问题。
图5所示为复合材料层合板的应力—应变曲线,图中最大应力即对应失效载荷。从计算结果可知,采用本方法能够较为准确地预测复合材料层合板的失效载荷。图6所示为复合材料层合板的最终破坏情况,中间黑色竖线表示该处材料的刚度已完全退化,说明复合材料沿该方向发生断裂,失效模式为基体拉伸失效,仿真结果与实验实际情况也相当吻合。
以上所述的仅是本申请的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。
Claims (2)
1.一种基于失效面理论的复合材料层合板强度分析方法,其特征在于,所述方法包括:
建立基于三维失效面理论的复合材料力学本构关系,包括:
复合材料在未损伤和损伤这两种状态下的应力—应变本构方程分别为:和/>,其中/>为自然坐标系下的应变,/>为自然坐标系下的应力向量,/>为未损伤状态下的柔度矩阵,/>通过试验测得的复合材料工程弹性常数确定;/>为损伤状态下的柔度矩阵,/>与复合材料损伤程度有关,并且/>由所述三维失效面理论推导得到,包括:
首先采用基于所述三维失效面理论构造的失效准则来判断所述复合材料是否发生损伤;根据所述三维失效面理论,材料失效是由失效面上的应力决定的,因此先通过坐标转换矩阵T将所述应力向量由自然坐标系转换到三维失效面坐标系:/>,其中/>,代表失效面坐标系下的真实应力向量,然后将所述失效面坐标系下的应力分量/>代入基于所述三维失效面理论构造的失效准则以判断所述复合材料是否失效及其失效模式;
如果所述复合材料满足所述失效准则,则发生损伤,需进行刚度性能退化;根据所述三维失效面理论,损伤变量矩阵同样应该定义在所述失效面坐标系下,而非传统的自然坐标系;由损伤力学知识可知,失效面上的有效应力/>和真实应力/>满足关系式:/>;基于能量等效假设,在失效面坐标系中未损伤和损伤这两种状态下的柔度矩阵满足关系式;最后通过所述坐标转换矩阵T将失效面坐标系中的两柔度矩阵由失效面坐标系转换回自然坐标系,得到/>;在所述复合材料发生损伤后,应力—应变本构关系由/>确定;
如果所述复合材料不满足所述失效准则,则未发生损伤,材料刚度性能不变,所述复合材料的应力—应变本构关系由确定;
基于有限元软件二次开发的复合材料层合板三维建模,包括:
根据复合材料结构的几何特征,采用参数化方法建立所述复合材料层合板的三维有限元模型,设置材料属性,再划分网格,并施加载荷和边界条件;
所述复合材料层合板的渐进失效分析,包括:
对所述复合材料层合板的三维有限元模型进行应力分析,模拟所述复合材料层合板的损伤起始、损伤演化和最终破坏的全过程,在每个载荷增量步,采用基于三维失效面理论构造的失效准则判断材料是否损伤,对失效的材料基于失效面坐标系进行刚度退化,然后继续施加载荷增量,直至确定所述复合材料层合板最终破坏;
根据有限元数值分析结果得到所述复合材料层合板的失效载荷与失效模式;
所述对失效的材料基于失效面坐标系进行刚度退化,包括:
在所述失效面坐标系中定义所述损伤变量矩阵:
式中d l、d n、d t和d nt、d lt、d nl是对应不同失效模式的损伤变量,若材料未发生失效,则损伤变量的值都为0;基于能量等效假设,最终得到损伤状态下的柔度矩阵,式中/>为复合材料未损伤状态下的柔度矩阵,T为所述坐标转换矩阵,对应的计算公式为:,
其中表示所述失效面坐标系t方向与自然坐标系对应的方向之间的夹角。
2.根据权利要求1所述的基于失效面理论的复合材料层合板强度分析方法,其特征在于,所述首先采用基于所述三维失效面理论构造的失效准则来判断所述复合材料是否发生损伤,包括:
所述失效准则及其失效模式为:
① 纤维拉伸失效:,
② 纤维压缩失效:,
③ 基体拉伸失效:,
④ 基体压缩失效:,
式中X T是复合材料的纤维纵向拉伸,X C是复合材料的纤维压缩强度,是复合材料的基体横向拉伸,/>是复合材料的基体压缩强度,S 21是纵向剪切强度,/>为复合材料横向压缩断裂角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010559632.7A CN111709174B (zh) | 2020-06-18 | 2020-06-18 | 一种基于失效面理论的复合材料层合板强度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010559632.7A CN111709174B (zh) | 2020-06-18 | 2020-06-18 | 一种基于失效面理论的复合材料层合板强度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111709174A CN111709174A (zh) | 2020-09-25 |
CN111709174B true CN111709174B (zh) | 2024-04-09 |
Family
ID=72541473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010559632.7A Active CN111709174B (zh) | 2020-06-18 | 2020-06-18 | 一种基于失效面理论的复合材料层合板强度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111709174B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112231858B (zh) * | 2020-10-09 | 2024-04-19 | 北京航天发射技术研究所 | 一种复合材料筒体多排螺接结构的失效预测方法和装置 |
CN113139313B (zh) * | 2021-04-25 | 2023-03-14 | 重庆大学 | 基于有限元理论的等效轻量化结构分析方法 |
CN114611389A (zh) * | 2022-03-04 | 2022-06-10 | 北京航空航天大学 | 一种基于人工智能的复合材料失效高效模拟方法 |
CN117275633B (zh) * | 2023-11-20 | 2024-01-23 | 北京航空航天大学 | 一种航空复合材料结构破坏过程分析方法及计算机设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698199A (zh) * | 2013-09-26 | 2014-04-02 | 北京航空航天大学 | 一种基于细观力学退化模型的复合材料结构失效分析方法 |
CN106503292A (zh) * | 2016-09-20 | 2017-03-15 | 浙江大学 | 预测低速冲击下复合材料层合板渐进失效的有限元方法 |
CN106777769A (zh) * | 2017-01-08 | 2017-05-31 | 浙江大学 | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 |
CN107832560A (zh) * | 2017-11-29 | 2018-03-23 | 北京航空航天大学 | 一种基于渐近损伤模型的全SiC复合材料多钉连接结构失效分析方法 |
-
2020
- 2020-06-18 CN CN202010559632.7A patent/CN111709174B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698199A (zh) * | 2013-09-26 | 2014-04-02 | 北京航空航天大学 | 一种基于细观力学退化模型的复合材料结构失效分析方法 |
CN106503292A (zh) * | 2016-09-20 | 2017-03-15 | 浙江大学 | 预测低速冲击下复合材料层合板渐进失效的有限元方法 |
CN106777769A (zh) * | 2017-01-08 | 2017-05-31 | 浙江大学 | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 |
CN107832560A (zh) * | 2017-11-29 | 2018-03-23 | 北京航空航天大学 | 一种基于渐近损伤模型的全SiC复合材料多钉连接结构失效分析方法 |
Non-Patent Citations (2)
Title |
---|
基于Puck理论的复合材料层合板横向剪切失效分析;贾利勇等;复合材料学报;第1-8页 * |
飞机结构材料的失效准则研究;顾杰斐;知网;第6章 * |
Also Published As
Publication number | Publication date |
---|---|
CN111709174A (zh) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111709174B (zh) | 一种基于失效面理论的复合材料层合板强度分析方法 | |
Maziz et al. | A progressive damage model for pressurized filament-wound hybrid composite pipe under low-velocity impact | |
CN110793853B (zh) | 基于基本力学参量的拉扭稳态循环应力应变建模方法 | |
Asaee et al. | A practical analytical model for predicting the low-velocity impact response of 3D-fiber metal laminates | |
Patel et al. | N-Layer concentric cylinder model (NCYL): an extended micromechanics-based multiscale model for nonlinear composites | |
Cheng et al. | Modelling pre-fatigue, low-velocity impact and post-impact fatigue behaviours of composite helicopter tail structures under multipoint coordinated loading spectrum | |
Branner et al. | Anisotropic beam model for analysis and design of passive controlled wind turbine blades | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
CN116882232A (zh) | 低速冲击下的碳纤维复合材料层合板损伤状态的预测方法 | |
CN115391947A (zh) | 轨道交通车辆复合材料结构虚拟分析方法及系统 | |
Alhussein et al. | Simulation of the effect laminate sequence on delamination mode-I with elastic couplings | |
Lou et al. | Vibration fatigue probabilistic life prediction model and method for blade | |
Li et al. | Simulations and tests of composite marine structures under low-velocity impact | |
Yi et al. | An improved energy prediction method to predict the fatigue life of laminated rubber-alloy spherical thrust elastomeric bearing under multiaxial loads | |
Pham et al. | Experimental and numerical failure analyses of composite cruciform under biaxial static and cyclic loading | |
Sadeghi et al. | Local and microscopic damage indices applicable to RC structures and concretes subjected to cyclic loading | |
Widanto et al. | Numerical simulation investigation of split element modeling with variation impactor diameter effect on low-velocity impact response of unidirectional CFRP | |
Borkowski et al. | Rapid Analysis Method for Composite Compression after Impact Strength Prediction | |
Ye et al. | Test and simulation of connection and failure performance of high-strength bolts in steel structure | |
Li et al. | Reliability analysis of laminated composite plates based on progressive failure method and universal gray system approach | |
Cheng et al. | Modelling Pre-fatigue, Low-velocity Impact and Fatigue behaviours of Composite Helicopter Tail Structures under Multipoint Coordinated Loading Spectrum | |
Goldberg et al. | Verification and validation of a three-dimensional composite impact model with tabulated input | |
Hu | The Simulation of Different Types of Failure | |
Goldberg et al. | Simulation of Impact Response of Polymer Composites Utilizing an Orthotropic Three-Dimensional Model with Tabulated Input | |
Sobrinho et al. | Modelling of the ferritic stainless steel grade 430 bolted connections with damage implementation |
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 |