CN107330202A - 一种基于虚拟分层扩展试验的复合材料i型分层桥联法则确定方法 - Google Patents
一种基于虚拟分层扩展试验的复合材料i型分层桥联法则确定方法 Download PDFInfo
- Publication number
- CN107330202A CN107330202A CN201710540112.XA CN201710540112A CN107330202A CN 107330202 A CN107330202 A CN 107330202A CN 201710540112 A CN201710540112 A CN 201710540112A CN 107330202 A CN107330202 A CN 107330202A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- layering
- data
- bridging
- 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
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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
- G01N3/12—Pressure testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，包括以下步骤：(1)开展I型分层试验，测得分层扩展R曲线；(2)建立I型分层试验件有限元模型，采用双线性内聚力单元模拟分层界面；(3)将R曲线嵌入有限元模型，开展I型分层模拟；(4)从数值模拟结果中调出一组一一对应的载荷P、位移d、分层长度a和分层起始位置处张开位移δ*的数据，计算断裂韧度GIc(a)；(5)利用GIc(a)和δ*的数据，确定广义裂纹面桥联应力随δ*变化的σg(δ*)，并拟合得到桥联法则具体形式。本发明适用于复合材料I型分层桥联法则的确定，相比于纯试验测试手段具有显著优势，可方便实现张开位移与断裂韧度数据的对应，简化试验操作，并降低试验测试误差。
Description
技术领域
本发明涉及复合材料I型分层扩展过程中桥联法则的确定的技术领域，具体涉及一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，适用于飞行器结构中广泛使用的碳纤维增强树脂基复合材料层合板结构。
背景技术
在发生大范围纤维桥联的分层扩展过程中断裂韧性R曲线取决于试验件的几何尺寸，此时R曲线不能被继续当做反映材料内在属性的参量，取而代之的是桥联法则，因为研究表明桥联法则与试验件的几何尺寸无关。此外，在测得分层扩展过程中桥联准则的基础上，利用纤维桥联影响下R曲线的理论公式，或采用基于三线性本构的内聚力单元进行分层扩展的数值模拟，就能预测任意几何尺寸试验件的R曲线，因此桥联法则的确定一直是复合材料分层行为研究中的一个重要方面。
根据现有研究确定桥联法则的方法有以下五类：(1)基于微观分层阻力机制和相应的力学假设，提出包含不同微观增韧机制的桥联法则力学模型。(2)采用裂纹已经贯穿整个试验件，但试验件上、下两臂仍被裂纹面间桥联的纤维粘联在一起的试验件，进行横向加载，直接测量裂纹面张开位移和纤维桥联应力之间的关系。(3)根据试验测得的裂纹面张开位移与一个裂纹面上作用着桥联应力的力学模型，然后基于断裂力学原理确定桥联法则。(4)基于试验测量的断裂韧性数据反推桥联法则。(5)采用光纤光栅传感(FBG)测量应变，结合具有特定参数化桥联应力分布的有限元模型，通过迭代求解，获得了上述层合板I型分层的桥联法则。
纵观上述5种确定桥联法则的方法，第一种基于力学理论的桥联法则确定方法无法全面涵盖纤维桥联过程中发生的纤维断裂、纤维基体界面分离、纤维抽丝等复杂的断裂行为，这限制了它的应用。第二种方法操作难度较大，而第三种方法需要不断调试力学模型中的桥联应力分布，实施难度也很大。第五种方法需要特殊制备的嵌入了光纤光栅传感(FBG)的试验件。相比而言第四种方法操作简单，可实施性强，因此被广泛采用。但是采用这种方法存在一个关键的技术难点，就是测量预埋分层末端的张开位移δ*。这个预埋分层末端张开位移δ*的测量需与试验中实时测量的载荷、位移、分层长度数据一一对应起来，为此可在试验中设置一个相机对预埋分层末端位置进行实时拍照，最后通过图像处理确定δ*的数据，但随着分层的扩展这个张开位移δ*会逐渐增大甚至超过相机能捕捉的范围，此后将无法获得δ*的数据，此外通过图像处理获得δ*的数值可能存在精度低的问题。
发明内容
本发明要解决的技术问题为：克服现有技术的不足，改进现存的采用试验测得的断裂韧度数据反推桥联法则的方法，针对δ*数据测量难度大的问题，提出了一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法。
本发明解决上述技术问题采用的技术方案为：一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，包括以下步骤：
第一步，开展复合材料层板I型静力分层扩展试验，测得其分层扩展过程中的R曲线；
第二步，建立I型载荷作用下DCB试验件的三维有限元模型，采用双线性本构内聚力单元模拟分层界面；
第三步，将R曲线嵌入上述采用内聚力单元模拟分层界面的层合板DCB试验件的有限元模型，在软件的隐形分析模块，开展该DCB试验件的I型分层扩展模拟；
第四步，从数值模拟结果中调出一组一一对应的载荷P、位移d、分层长度a的数据、分层起始位置处张开位移δ*的数据。利用调出的P、d、a数据计算断裂韧度数据GIc(a)；
第五步，利用断裂韧度GIc(a)和分层起始位置处张开位移δ*的数据，确定广义裂纹面桥联应力随裂纹面张开位移变化规律的数据σg(δ*)，并采用曲线拟合确定桥联法则的具体形式σg(δ)。
进一步地，所述第一步测得复合材料层板I型分层扩展过程中R曲线的方法是：
A.开展静力载荷作用下的分层扩展试验，记录试验过程中载荷P、位移d、分层长度a的数据；
B.利用上述试验测得的P、d、a数据，根据基于修正梁理论的公式计算I型分层扩展过程中不同分层长度下的断裂韧度数据GIc：
其中，P是试验中的加载载荷，d是加载端的位移，B是试验件的宽度，Δ是分层长度修正量。采用基于实测柔度(C＝d/P)的分层长度修正方法获得，具体实施方法是：对试验测得的C1/3～a数据进行线性拟合，以拟合直线在横轴上截距的绝对值作为分层长度的修正量Δ；
C.采用如下式所示的R曲线理论公式，对分层扩展过程中获得的断裂韧度数据GIc进行拟合，获得该层合板试验件I型分层扩展的R曲线，
其中，<>表示一种特殊的运算法则：lbz表示纤维桥接区域的长度，GI,Prop是I型分层断裂韧度稳定值，GI,Init是I型分层断裂韧度初始值，△a表示分层扩展长度，为实测分层长度a与预埋分层长度a0之差，即△a＝a-a0。
进一步地，所述第三步将R曲线嵌入DCB试验件有限元模型，开展I型分层扩展模拟的方法是，编写子程序USDFLD，将R曲线嵌入内聚力单元本构定义模块。其中，在子程序USDFLD中以单元坐标为自变量，建立单元坐标和分层长度，以及分层长度与断裂韧度之间的映射关系，并将子程序USDFLD嵌入数值计算中。
进一步地，所述第四步根据数值模拟结果计算I型分层断裂韧度数据的方法是，利用数值模拟中调出的P、d、a数据，并采用基于修正梁理论的I型分层断裂韧度的计算公式，计算I型分层扩展过程中的断裂韧度数据GIc(a)。
进一步地，所述第五步利用断裂韧度GIc(a)和分层起始位置处张开位移δ*的数据，确定广义裂纹面桥联应力随裂纹面张开位移的变化规律σg(δ*)的方法是，将上述断裂韧度数据GIc(a)和分层起始位置处张开位移δ*的数据，代入如下式所示的差分形式的公式：
其中，GIc(ai)和GIc(aj)分别表示分层长度为ai和aj下的断裂韧度，δ*(ai)和δ*(aj)分别表示分层长度为ai和aj时分层起始位置处的张开位移。
根据上述确定的广义裂纹面桥联应力随裂纹面张开位移变化规律σg(δ*)的数据点，进一步利用曲线拟合确定桥联法则具体形式σg(δ)的方法是，采用下式所示的指数形式的桥联法则，对广义裂纹面桥联应力随裂纹面张开位移变化的数据点σg(δ*)进行拟合，得到DCB试验件I型分层过程中广义裂纹面桥联法则的具体形式：
其中，δ表示裂纹面张开位移，δ0是裂纹面张开位移的临界值，σg(δ)表示随张开位移变化的裂纹面桥联应力，表示裂纹面桥联应力的最大值，n是指数。
本发明与现有技术相比的优点在于：
(1)方便地实现分层起始位置处张开位移与断裂韧度数据的对应；
(2)在试验过程中只需测量R曲线，无需测量分层起始处的张开位移，简化了试验操作；
(3)有效避免了试验测试分层起始位置处张开位移可能引入的测量误差。
附图说明
图1是本发明方法的实现流程图；
图2是正交层合板I型分层过程中的断裂韧度数据和拟合得到的R曲线；
图3是正交层合板三维有限元模型；
图4是基于虚拟分层扩展试验获得的正交层合板I型分层的桥联法则。
具体实施方式
下面结合实施例对本发明作进一步详细说明。
实施例：采用碳/环氧树脂基复合材料正交层合板试验件I型分层桥联法则的确定
1.几何描述
正交层合板DCB试验件由上、下两臂构成。试验件的厚度为3.95mm，宽度为25mm，长度为200mm。
正交层合板的铺层顺序为[0°/90°]10，单向带的材料属性为：E1＝120.2GPa，E2＝E3＝7.3GPa，G12＝G13＝G23＝3.9GPa，ν12＝ν13＝0.28，ν23＝0.48。
采用聚四氟乙烯共聚物与乙烯制成的55mm长和20μm厚的薄膜(A6000)，作为插入层板中面的预埋分层，得到预定的0°/90°分层界面。
2.开展复合材料层板I型静力分层扩展试验，测得其分层扩展过程中的R曲线；
(A)将DCB试验件装载于MTS 880伺服液压试验机上，并选用1500N量程的载荷传感器。采用位移加载模式开展静力下的DCB试验，加载速率控制为0.1mm/min以保证分层缓慢、稳定地扩展。实时记录I型分层扩展试验中载荷P、位移d、分层长度a的数据；
(B)利用上述试验测得的P、d、a数据，根据基于修正梁理论的公式计算I型分层扩展过程中不同分层长度下的断裂韧度数据GIc，数据结果如图2所示。
其中，P是试验中的加载载荷，d是加载端的位移，B是试验件的宽度，Δ是分层长度修正量。采用基于实测柔度(C＝d/P)的分层长度修正方法获得，具体实施方法是：对试验测得的C1/3～a数据进行线性拟合，以拟合直线在横轴上截距的绝对值作为分层长度的修正量Δ；
(C)采用下式所示的R曲线理论公式，对分层扩展过程中获得的断裂韧度数据GIc进行拟合，获得该层合板试验件I型分层扩展的R曲线，如图2所示。
其中，<>表示一种特殊的运算法则：lbz表示纤维桥接区域的长度，GI,Prop是I型分层断裂韧度稳定值，GI,Init是I型分层断裂韧度初始值。△a表示分层扩展长度，为实测分层长度a与预埋分层长度a0之差，即△a＝a-a0。
3.建立I型载荷作用下DCB试验件的三维有限元模型，如图3所示。有限元模型中采用双线性本构内聚力单元模拟分层界面。
4.将拟合得到的R曲线嵌入上述采用内聚力单元模拟分层界面的层合板DCB试验件的有限元模型。为此，编写子程序USDFLD，将R曲线嵌入内聚力单元本构定义模块，在子程序USDFLD中以单元坐标为自变量，建立单元坐标和分层长度，以及分层长度与断裂韧度之间的映射关系，并将子程序USDFLD嵌入数值计算中。在软件的隐形分析模块，开展该DCB试验件的I型分层扩展模拟。
5.从数值模拟结果中调出一组一一对应的载荷P、位移d、分层长度a的数据、分层起始位置处张开位移δ*的数据。利用数值模拟中调出的P、d、a数据，并采用基于修正梁理论的I型分层断裂韧度的计算公式，计算I型分层扩展过程中的断裂韧度数据GIc(a)。
6.将上述断裂韧度数据GIc(a)和分层起始位置处张开位移δ*的数据，代入下式所示的差分形式的公式，获得一组广义裂纹面桥联应力随裂纹面张开位移变化的数据σg(δ*)，如图4中星点所示。
其中，GIc(ai)和GIc(aj)分别表示分层长度为ai和aj下的断裂韧度，δ*(ai)和δ*(aj)分别表示分层长度为ai和aj时分层起始位置处的张开位移。
7.根据上述确定的广义裂纹面桥联应力随裂纹面张开位移变化的数据点σg(δ*)，采用下式所示的指数形式的桥联法则，对广义裂纹面桥联应力随裂纹面张开位移变化的数据点σg(δ*)进行拟合，得到DCB试验件I型分层过程中广义裂纹面桥联法则的拟合曲线及其具体形式，如图4所示。
其中，δ表示裂纹面张开位移，δ0是裂纹面张开位移的临界值，σg(δ)表示随张开位移变化的裂纹面桥联应力，表示裂纹面桥联应力的最大值，n是指数。
本发明未详细阐述的部分属于本领域公知技术。
以上所述，仅为本发明部分具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本领域的人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。
Claims (6)
1.一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于包括以下步骤：
第一步，开展复合材料层板I型静力分层扩展试验，测得其分层扩展过程中的R曲线；
第二步，建立I型载荷作用下DCB试验件的三维有限元模型，采用双线性本构内聚力单元模拟分层界面；
第三步，将R曲线嵌入上述采用内聚力单元模拟分层界面的层合板DCB试验件的有限元模型，在软件的隐形分析模块，开展该DCB试验件的I型分层扩展模拟；
第四步，从数值模拟结果中调出一组一一对应的载荷P、位移d、分层长度a的数据、分层起始位置处张开位移δ*的数据，利用调出的P、d、a数据计算断裂韧度数据GIc(a)；
第五步，利用断裂韧度GIc(a)和分层起始位置处张开位移δ*的数据，确定广义裂纹面桥联应力随裂纹面张开位移变化规律的数据σg(δ*)，并采用曲线拟合确定桥联法则的具体形式σg(δ)。
2.根据权利要求1所述一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于：所述第一步测得复合材料层板I型分层扩展过程中R曲线的方法是：
A.开展静力载荷作用下的分层扩展试验，记录试验过程中载荷P、位移d、分层长度a的数据；
B.利用上述试验测得的P、d、a数据，根据基于修正梁理论的公式计算I型分层扩展过程中不同分层长度下的断裂韧度数据GIc：
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>P</mi>
<mi>d</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中，P是试验中的加载载荷，d是加载端的位移，B是试验件的宽度，Δ是分层长度修正量；采用基于实测柔度(C＝d/P)的分层长度修正方法获得，具体实施方法是：对试验测得的C1/3～a数据进行线性拟合，以拟合直线在横轴上截距的绝对值作为分层长度的修正量Δ；
C.采用如下式所示的R曲线理论公式，对分层扩展过程中获得的断裂韧度数据GIc进行拟合，获得该层合板试验件I型分层扩展的R曲线，
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mo>,</mo>
<mi>I</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mo>,</mo>
<mi>Pr</mi>
<mi>o</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mo>,</mo>
<mi>I</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo><</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>b</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>a</mi>
<mo>></mo>
</mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>b</mi>
<mi>z</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mo>,</mo>
<mi>Pr</mi>
<mi>o</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mo>,</mo>
<mi>I</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo><</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>b</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>a</mi>
<mo>></mo>
</mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>b</mi>
<mi>z</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中，<>表示一种特殊的运算法则：lbz表示纤维桥接区域的长度，GI,Prop是I型分层断裂韧度稳定值，GI,Init是I型分层断裂韧度初始值，△a表示分层扩展长度，为实测分层长度a与预埋分层长度a0之差，即△a＝a-a0。
3.根据权利要求1所述一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于：所述第三步将R曲线嵌入DCB试验件有限元模型，开展I型分层扩展模拟的方法是，编写子程序USDFLD，将R曲线嵌入内聚力单元本构定义模块，其中，在子程序USDFLD中以单元坐标为自变量，建立单元坐标和分层长度，以及分层长度与断裂韧度之间的映射关系，并将子程序USDFLD嵌入数值计算中。
4.根据权利要求1所述一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于：所述第四步根据数值模拟结果计算I型分层断裂韧度数据的方法是，利用数值模拟中调出的P、d、a数据，并采用基于修正梁理论的I型分层断裂韧度的计算公式，计算I型分层扩展过程中的断裂韧度数据GIc(a)。
5.根据权利要求1所述一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于：所述第五步利用断裂韧度GIc(a)和分层起始位置处张开位移δ*的数据，确定广义裂纹面桥联应力随裂纹面张开位移的变化规律σg(δ*)的方法是，将上述断裂韧度数据GIc(a)和分层起始位置处张开位移δ*的数据，代入如下式所示的差分形式的公式：
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&delta;</mi>
<mo>*</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>lim</mi>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&RightArrow;</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
</mrow>
</munder>
<mfrac>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>I</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&delta;</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&delta;</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中，GIc(ai)和GIc(aj)分别表示分层长度为ai和aj下的断裂韧度，δ*(ai)和δ*(aj)分别表示分层长度为ai和aj时分层起始位置处的张开位移。
6.根据权利要求1所述一种基于虚拟分层扩展试验的复合材料I型分层桥联法则确定方法，其特征在于：利用曲线拟合确定桥联法则具体形式σg(δ)的方法是，采用下式所示的指数形式的桥联法则，对广义裂纹面桥联应力随裂纹面张开位移变化的数据点进行拟合，得到DCB试验件I型分层过程中广义裂纹面桥联法则的具体形式：
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>g</mi>
<mn>0</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&delta;</mi>
<mo>/</mo>
<msub>
<mi>&delta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&delta;</mi>
<mo>&le;</mo>
<msub>
<mi>&delta;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>&delta;</mi>
<mo>></mo>
<msub>
<mi>&delta;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，δ表示裂纹面张开位移，δ0是裂纹面张开位移的临界值，σg(δ)表示随张开位移变化的裂纹面桥联应力，表示裂纹面桥联应力的最大值，n是指数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710540112.XA CN107330202B (zh) | 2017-07-05 | 2017-07-05 | 一种复合材料i型分层桥联法则确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710540112.XA CN107330202B (zh) | 2017-07-05 | 2017-07-05 | 一种复合材料i型分层桥联法则确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107330202A true CN107330202A (zh) | 2017-11-07 |
| CN107330202B CN107330202B (zh) | 2020-06-30 |
Family
ID=60196080
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710540112.XA Active CN107330202B (zh) | 2017-07-05 | 2017-07-05 | 一种复合材料i型分层桥联法则确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107330202B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109991077A (zh) * | 2019-03-18 | 2019-07-09 | 重庆大学 | 一种复合材料混合型分层阻力曲线的预测方法 |
| CN110376055A (zh) * | 2019-06-20 | 2019-10-25 | 重庆大学 | 一种基于新型内聚力本构关系的cfrp层板分层失效行为预测方法 |
| CN110470534A (zh) * | 2019-07-22 | 2019-11-19 | 中国航发北京航空材料研究院 | 一种无需追踪裂纹的复合材料层间i型断裂韧性测试方法 |
| CN110598293A (zh) * | 2019-09-03 | 2019-12-20 | 上海交通大学 | 微纳米纤维增强复合材料的断裂破坏行为预测方法 |
| CN111339697A (zh) * | 2020-01-14 | 2020-06-26 | 北京航空航天大学 | 一种基于修正梁理论的复合材料i型分层桥联律确定方法 |
| CN112857996A (zh) * | 2021-01-13 | 2021-05-28 | 北京航空航天大学 | 一种高效的复合材料层合板i型分层桥联法则确定方法 |
| CN114279844A (zh) * | 2021-12-29 | 2022-04-05 | 河北工业大学 | 纤维增强复合材料层合板i/iii混合型层间断裂韧度测试方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
| EP2769834A1 (en) * | 2013-02-26 | 2014-08-27 | Airbus Operations, S.L. | An artificially defective cured laminate. |
| CN105022898A (zh) * | 2015-08-17 | 2015-11-04 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料胶接结构的裂纹性能测定和优化求解方法 |
| CN105488310A (zh) * | 2016-01-22 | 2016-04-13 | 重庆大学 | 一种归一化cfrp多向层合板疲劳分层扩展速率预测方法 |
| CN106289975A (zh) * | 2016-08-12 | 2017-01-04 | 上海电气电站设备有限公司 | 材料微区断裂韧性的试验方法 |
| CN106599507A (zh) * | 2016-12-26 | 2017-04-26 | 北京航空航天大学 | 一种改进b‑k准则用于含纤维桥接影响复合材料多向层板分层预测的方法 |
-
2017
- 2017-07-05 CN CN201710540112.XA patent/CN107330202B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2769834A1 (en) * | 2013-02-26 | 2014-08-27 | Airbus Operations, S.L. | An artificially defective cured laminate. |
| CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
| CN105022898A (zh) * | 2015-08-17 | 2015-11-04 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料胶接结构的裂纹性能测定和优化求解方法 |
| CN105488310A (zh) * | 2016-01-22 | 2016-04-13 | 重庆大学 | 一种归一化cfrp多向层合板疲劳分层扩展速率预测方法 |
| CN106289975A (zh) * | 2016-08-12 | 2017-01-04 | 上海电气电站设备有限公司 | 材料微区断裂韧性的试验方法 |
| CN106599507A (zh) * | 2016-12-26 | 2017-04-26 | 北京航空航天大学 | 一种改进b‑k准则用于含纤维桥接影响复合材料多向层板分层预测的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 张肖肖 等: "复合材料分层模拟数值方法概述", 《北京力学会第17届学术年会》 * |
| 颜剑敏 等: "基于内聚力模型的复合材料分层扩展研究", 《北京力学会第21届学术年会暨北京振动工程学会第22届学术年会论文集》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109991077A (zh) * | 2019-03-18 | 2019-07-09 | 重庆大学 | 一种复合材料混合型分层阻力曲线的预测方法 |
| CN109991077B (zh) * | 2019-03-18 | 2021-10-19 | 重庆大学 | 一种复合材料混合型分层阻力曲线的预测方法 |
| CN110376055A (zh) * | 2019-06-20 | 2019-10-25 | 重庆大学 | 一种基于新型内聚力本构关系的cfrp层板分层失效行为预测方法 |
| CN110470534A (zh) * | 2019-07-22 | 2019-11-19 | 中国航发北京航空材料研究院 | 一种无需追踪裂纹的复合材料层间i型断裂韧性测试方法 |
| CN110470534B (zh) * | 2019-07-22 | 2021-10-15 | 中国航发北京航空材料研究院 | 一种无需追踪裂纹的复合材料层间i型断裂韧性测试方法 |
| CN110598293A (zh) * | 2019-09-03 | 2019-12-20 | 上海交通大学 | 微纳米纤维增强复合材料的断裂破坏行为预测方法 |
| CN111339697A (zh) * | 2020-01-14 | 2020-06-26 | 北京航空航天大学 | 一种基于修正梁理论的复合材料i型分层桥联律确定方法 |
| CN111339697B (zh) * | 2020-01-14 | 2022-03-25 | 北京航空航天大学 | 一种基于修正梁理论的复合材料i型分层桥联律确定方法 |
| CN112857996A (zh) * | 2021-01-13 | 2021-05-28 | 北京航空航天大学 | 一种高效的复合材料层合板i型分层桥联法则确定方法 |
| CN114279844A (zh) * | 2021-12-29 | 2022-04-05 | 河北工业大学 | 纤维增强复合材料层合板i/iii混合型层间断裂韧度测试方法 |
| CN114279844B (zh) * | 2021-12-29 | 2023-06-16 | 河北工业大学 | 纤维增强复合材料层合板i/iii混合型层间断裂韧度测试方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107330202B (zh) | 2020-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107330202A (zh) | 一种基于虚拟分层扩展试验的复合材料i型分层桥联法则确定方法 | |
| Sokolinsky et al. | Experimental and analytical study of nonlinear bending response of sandwich beams | |
| Bhandari et al. | Finite element analysis of thermoplastic polymer extrusion 3D printed material for mechanical property prediction | |
| Stöven et al. | Continuous monitoring of three-dimensional resin flow through a fibre preform | |
| Laffan et al. | Measurement of the in situ ply fracture toughness associated with mode I fibre tensile failure in FRP. Part I: Data reduction | |
| Jumbo et al. | Thermal residual stress analysis of epoxy bi-material laminates and bonded joints | |
| Gude et al. | Piezoelectrically driven morphing structures based on bistable unsymmetric laminates | |
| Rodini Jr et al. | An analytical and experimental investigation of edge delamination in composite laminates | |
| Davidson et al. | Influences of friction, geometric nonlinearities, and fixture compliance on experimentally observed toughnesses from three and four-point bend end-notched flexure tests | |
| Solis et al. | Influence of ply orientation on free-edge effects in laminates subjected to in-plane loads | |
| Melo et al. | Viscoelastic characterization of transversely isotropic composite laminae | |
| Bao et al. | Estimation of J-resistance curves for CT specimen based on unloading compliance method and CMOD data | |
| Dong | Dimension variation prediction and control for composites | |
| Knowles et al. | Thermal and volumetric property analysis of polymer networks and composites using elevated temperature digital image correlation | |
| Odegard et al. | Elasto-plastic analysis of the Iosipescu shear test | |
| Melo et al. | Determination of the elastic constants of a transversely isotropic lamina using laminate coefficients of thermal expansion | |
| Durelli | Experimental strain and stress analysis of solid propellant rocket motors | |
| Chen et al. | Integrated simulations of structural performance, molding process, and warpage for gas‐assisted injection‐molded parts. I. Analysis of part structural performance | |
| Brockman et al. | OMC processing simulation using an elastic-viscoplastic material model | |
| Cheng et al. | Simplified braiding through integration points model for triaxially braided composites | |
| Sessner et al. | Experimental and numerical characterisation of fibre‐metal‐elastomer laminates by using DMA regarding its damping behaviour | |
| Ghabezi et al. | Investigation and modeling of compaction behavior of plain fabrics | |
| El-Chiti | Experimental variability of e-glass reinforced vinyl ester composites fabricated by vartm/scrimp | |
| Montay et al. | Determination of the residual stresses in composite laminate using the compliance method | |
| Rodera Garcia | Damage onset modelling of curved composite laminates |
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 |