CN103593567B - 一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 - Google Patents
一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 Download PDFInfo
- Publication number
- CN103593567B CN103593567B CN201310571786.8A CN201310571786A CN103593567B CN 103593567 B CN103593567 B CN 103593567B CN 201310571786 A CN201310571786 A CN 201310571786A CN 103593567 B CN103593567 B CN 103593567B
- Authority
- CN
- China
- Prior art keywords
- unit
- length
- parallelepiped
- energy
- characteristic length
- 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
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,首先建立了平面单元特征长度的计算模型,在此基础上建立了三维平行六面体单元的特征长度计算方法,再通过对一般形状单元和楔形单元进行等体积转换,将其特征长度计算问题转化为平行六面体单元的特征长度计算问题。本发明用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法考虑了复合材料结构有限元模型中常用的典型单元构型,包括平行四边形单元,三角形单元,平行六面体单元,一般形状单元及楔形单元,建立了单元特征长度计算模型,提供了准确的单元特征长度计算公式,能够有效地计算单元特征长度,从而计算单元损伤的耗散能量,提高了复合材料结构失效有限元模拟的客观性和准确性。
Description
技术领域
本发明涉及复合材料结构强度和失效的有限元分析和模拟,适用于工程领域中各种复合材料结构有限元损伤模拟模型。
背景技术
随着现代航空航天飞行器对材料性能要求的提高,复合材料在航空航天飞行器结构中的用量越来越高。因此,发展复合材料结构分析技术日益重要。随着计算机和有限元技术的迅猛发展,以及复合材料在主结构上的逐步应用,复合材料结构的应力分析和强度预测越来越多地依赖于先进的有限元和数值分析技术。目前广泛应用的方法有裂纹闭合技术(VCCT),渐进损伤分析方法(PDM)及内聚力模型(CZM)等方法,这些方法或基于断裂力学、或基于损伤力学,或基于断裂-损伤力学,并结合先进的有限元分析技术,为深入理解复合材料结构的失效机理和有效预测复合材料结构破坏提供了有效的途径。
基于损伤力学方法应用有限元分析技术时,面临着损伤局部化问题,若处理不当,将导致预测结果失真。损伤局部化是指由于连续损伤力学模型中应力软化段的存在,在数值模拟中损伤会局限在一条窄带内,而由于损伤的数值模拟结果与网格的精细度密切相关,随着有限元模型网格尺寸的减小,数值计算的结构损伤耗散能量下降,模型预测的结构强度会下降。为克服损伤局部化,提出了裂纹带模型,假设材料损伤形成的宏观裂纹可以用一层单元的损伤区域来代替,将损伤裂纹的断裂韧度在整个单元体积内均匀化。裂纹带模型的本质是将材料损伤过程中耗散的总能量在整个单元体积(微元体)内均匀化,使得有限元模拟过程中材料损伤耗散的能量恒定。而单元特征长度是裂纹带模型中的重要参数,它不仅与单元的几何形状有关,而且与材料的属性和裂纹的方向有关,准确地计算单元特征长度是保证准确地计算耗散能量的前提条件。因此,有必要建立一般单元形状的特征长度的计算方法。
发明内容
本发明要解决的技术问题为:克服现有技术的不足,考虑了复合材料结构有限元模型中常用的典型单元构型,包括平行四边形单元,三角形单元,平行六面体单元,一般形状单元及楔形单元,建立了计算其特征长度的方法,该方法能够有效地计算单元特征长度,为裂纹带模型提供准确的单元特征长度值。
本发明解决上述技术问题采用的技术方案为:一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,该估计方法的步骤如下:
步骤一、裂纹带模型的本质是将材料损伤过程中耗散的总能量在整个单元体积即微元体内均匀化,使得有限元模拟过程中材料损伤耗散的能量恒定,而单元特征长度是裂纹带模型中的重要参数,它不仅与单元的几何形状有关,而且与材料的属性和裂纹的方向有关,准确地计算单元特征长度是保证准确地计算耗散能量的前提条件,复合材料结构失效有限元模拟中单元特征长度计算方法分为平面单元特征长度计算和三维单元特征长度计算,其中:
平面单元特征长度计算步骤如下:
设裂纹与单元一边的夹角为θcrc,裂纹长度为lcrc,单元相邻两边的夹角为边长分别为a和b,与AB边对应的距离为
步骤A1,首先假设一条裂纹穿过一层由n个平行四边形单元或三角形单元组成的区域;
步骤B1,根据单元特征长度的定义,可得:
步骤C1,当时,近似地认为n个平行四边形单元的面积为直角三角形AFE的2倍,则上式可简化为:
步骤D1,时,平行四边形单元或三角形单元的特征长度为
三维单元特征长度计算步骤如下:
步骤A2,首先将六面体与纤维平行的平面按照平行四边形特征长度等效转化方法转换为两边与纤维平行的矩形E'F'G'H',然后将六面体的侧面四条边在平面E'F'G'H'的垂线上进行投影,得到一立方体;
步骤B2,上述立方体A"B"C"D"E'F'G'H'各边的长度为:
式中,分别代表相应线段的长度,∠HEF、∠HEF代表角度,代表用AE表示的向量,代表单位向量;
步骤C2,由立方体A"B"C"D"E'F'G'H'易知,对应于纤维“裂纹”的单元特征长度为:
步骤D2,设基体裂纹面与等效立方体单元的E'F'G'H'面所成夹角为则对应于基体裂纹的单元特征长度为:
其中基体裂纹面与立方体单元边的相对位置可以通过裂纹面在模型整体坐标系中的方向坐标和平面E′F′H′G′的法向坐标求得;
步骤E2,对于形状更为一般的六面体单元,可以首先将其等体积转化为平行六面体单元,其中平行六面体单元的底面由四边形EFGH等面积转化来确定,进而采用平行六面体的计算公式求得单元的特征长度;
步骤F2,对于楔形单元,可以借鉴平面三角形单元特征长度的求解方法,将两个相同楔形单元组合、扩展为平行六面体后得到其单元特征长度。
步骤二、复合材料结构失效有限元模拟中单元耗散能量计算步骤如下:
步骤A3,在损伤生成阶段,微裂纹分散在整个单元体积内,随着载荷的增大,微裂纹生长、贯通,最终形成宏观裂纹。因此,单元在损伤(微裂纹)形成、扩展过程中耗散的总能量WD就是细观裂纹形成消耗的断裂能WC,即:
WD=WC;
步骤B3,假设材料的临界应变能释放率GC恒定,则在单元内材料生成面积为A的裂纹消耗的断裂能WC为:
WC=GC*A;
步骤C3,基于上述单元特征长度计算方法得到的特征长度,由WD=WC,得单元的损伤耗散能密度为:
其中,GC为材料的临界应变能释放率,V为单元的体积,A为与宏观裂纹平行的截面积,l*为单元的特征长度,由上述单元特征长度计算方法得到;
步骤D3,基于上面得到的单元损伤耗散能密度,最终得到从微裂纹产生、扩展到宏观裂纹形成过程中,单元损伤所耗散的能量WD为:
WD=Γ*V。
本发明与现有技术相比的优点在于:
(1)、本发明用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,首先建立了平面单元特征长度计算模型,进而得到其计算公式,在此基础上建立了三维平行六面体单元的特征长度近似计算方法,再通过对一般形状单元和楔形单元进行等体积转换,将其特征长度计算问题转化为平行六面体单元的特征长度计算问题。
(2)、本发明用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,考虑了复合材料结构有限元模型中常用的典型单元构型,包括平行四边形单元,三角形单元,平行六面体单元,一般形状单元及楔形单元,能够有效地计算单元特征长度,为裂纹带模型提供准确的单元特征长度值。
附图说明
图1是复合材料典型三维单元构型,其中图1(a)和(b)是六面体单元及材料方向,构型(b)更具有一般性,图1(c)和(d)是楔形单元及材料方向;
图2是单元特征长度的计算原理图;
图3是平面单元与裂纹示意图,其中图3(a)为平行四边形单元与裂纹示意图,图3(b)为三角形单元与裂纹示意图;
图4是平行六面体等体积转换示意图;
图5是一般形状六面体转换示意图;
图6是楔形单元扩展为六面体示意图。
具体实施方式
本发明用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,该估计方法的具体步骤如下:
步骤一、裂纹带模型的本质是将材料损伤过程中耗散的总能量在整个单元体积即微元体内均匀化,使得有限元模拟过程中材料损伤耗散的能量恒定,而单元特征长度是裂纹带模型中的重要参数,它不仅与单元的几何形状有关,而且与材料的属性和裂纹的方向有关,准确地计算单元特征长度是保证准确地计算耗散能量的前提条件,复合材料结构失效有限元模拟中单元特征长度计算方法分为平面单元特征长度计算和三维单元特征长度计算,其中:
平面单元特征长度计算步骤如下:
设裂纹与单元一边的夹角为θcrc,裂纹长度为lcrc,单元相邻两边的夹角为边长分别为a和b,与AB边对应的距离为首先假设一条裂纹穿过一层由n个平行四边形单元或三角形单元组成的区域,根据单元特征长度的定义,可得:
当时,近似地认为n个平行四边形单元的面积为直角三角形AFE的2倍,则上式可简化为:
当时,平行四边形单元或三角形单元的特征长度为:
三维单元特征长度计算步骤如下:
首先将六面体与纤维平行的平面按照平行四边形特征长度等效转化方法转换为两边与纤维平行的矩形E'F'G'H',然后将六面体的侧面四条边在平面E'F'G'H'的垂线上进行投影,得到如图4所示的立方体A"B"C"D"E'F'G'H',其各边的长度为:
式中,分别代表相应线段的长度,∠HEF、∠HEF代表角度,代表用AE表示的向量,代表单位向量。由立方体A"B"C"D"E'F'G'H'易知,对应于纤维“裂纹”的单元特征长度为:
设基体裂纹面与等效立方体单元的E'F'G'H'面所成夹角为则对应于基体裂纹的单元特征长度为:
其中基体裂纹面与立方体单元边的相对位置可以通过裂纹面在模型整体坐标系中的方向坐标和平面E′F′H′G′的法向坐标求得。对于形状更为一般的六面体单元,可以首先将其等体积转化为平行六面体单元,其中平行六面体单元的底面由四边形EFGH等面积转化来确定,进而采用平行六面体的计算公式求得单元的特征长度,其等体积转化过程如图5所示;对于楔形单元,可以借鉴平面三角形单元特征长度的求解方法,将两个相同楔形单元组合、扩展为平行六面体后得到其单元特征长度,楔形单元转化为平行六面体的过程如图6所示。
步骤二、复合材料结构失效有限元模拟中单元耗散能量计算步骤如下:
在损伤生成阶段,微裂纹分散在整个单元体积内,随着载荷的增大,微裂纹生长、贯通,最终形成宏观裂纹。因此,单元在损伤(微裂纹)形成、扩展过程中耗散的总能量WD就是细观裂纹形成消耗的断裂能WC,即:
WD=WC
假设材料的临界应变能释放率GC恒定,则在单元内材料生成面积为A的裂纹消耗的断裂能WC为:
WC=GC*A
基于上述单元特征长度计算方法得到的特征长度,由WD=WC,得单元的损伤耗散能密度为:
其中,GC为材料的临界应变能释放率,V为单元的体积,A为与宏观裂纹平行的截面积,l*为单元的特征长度,由上述单元特征长度计算方法得到。基于单元损伤耗散能密度,最终得到从微裂纹产生、扩展到宏观裂纹形成过程中,单元损伤所耗散的能量WD为:
WD=Γ*V。
本发明未详细公开的部分属于本领域的公知技术。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法,其特征在于,该估计方法的步骤如下:
步骤一、裂纹带模型的本质是将材料损伤过程中耗散的总能量在整个单元体积即微元体内均匀化,使得有限元模拟过程中材料损伤耗散的能量恒定,而单元特征长度是裂纹带模型中的重要参数,它不仅与单元的几何形状有关,而且与材料的属性和裂纹的方向有关,准确地计算单元特征长度是保证准确地计算耗散能量的前提条件,复合材料结构失效有限元模拟中单元特征长度计算方法分为平面单元特征长度计算和三维单元特征长度计算,其中:
平面单元特征长度计算步骤如下:
设裂纹与单元一边的夹角为θcrc,裂纹长度为lcrc,单元相邻两边的夹角为边长分别为a和b,与AB边对应的距离为
步骤A1,首先假设一条裂纹穿过一层由n个平行四边形单元或三角形单元组成的区域;
步骤B1,根据单元特征长度的定义,可得:
步骤C1,当时,近似地认为n个平行四边形单元的面积为直角三角形AFE的2倍, 为直角三角形AFE的边AF的长度;则上式可简化为:
步骤D1,当时,平行四边形单元或三角形单元的特征长度为:
三维单元特征长度计算步骤如下:
步骤A2,首先将六面体与纤维平行的平面EFGH按照平行四边形特征长度等效转化方法转换为两边与纤维平行的矩形E'F'G'H',然后将六面体的侧面四条边在平面E'F'G'H'的垂线上进行投影,得到一立方体;
步骤B2,上述立方体A"B"C"D"E'F'G'H'各边的长度为:
式中,分别代表相应线段的长度,为平行六面体EFGH面上EH边的长度,θf为平行六面体EFGH面上EF边与纤维方向的夹角,SEFGH为平行六面体EFGH面的面积,∠HEF、∠HEF代表角度,代表用AE表示的向量,代表单位向量;
步骤C2,由立方体A"B"C"D"E'F'G'H'易知,对应于纤维“裂纹”的单元特征长度为:
其中,为平行六面体等体积转化时平面E'F'G'H'上E'F'边的长度;
步骤D2,设基体裂纹面与等效立方体单元的E'F'G'H'面所成夹角为则对应于基体裂纹的单元特征长度为:
其中,为平行六面体单元等体积转化时平面A'B'E'F'上A'E'边的长度,基体裂纹面与立方体单元边的相对位置可以通过裂纹面在模型整体坐标系中的方向坐标和平面E′F′H′G′的法向坐标求得;
步骤E2,对于形状更为一般的六面体单元,可以首先将其等体积转化为平行六面体单元,其中平行六面体单元的底面由四边形EFGH等面积转化来确定,进而采用平行六面体的计算公式求得单元的特征长度;
步骤F2,对于楔形单元,可以借鉴平面三角形单元特征长度的求解方法,将两个相同楔形单元组合、扩展为平行六面体后得到其单元特征长度;
步骤二、复合材料结构失效有限元模拟中单元耗散能量计算步骤如下:
步骤A3,在损伤生成阶段,微裂纹分散在整个单元体积内,随着载荷的增大,微裂纹生长、贯通,最终形成宏观裂纹,因此,单元在损伤即微裂纹形成、扩展过程中耗散的总能量WD就是细观裂纹形成消耗的断裂能WC,即:
WD=WC;
步骤B3,假设材料的临界应变能释放率GC恒定,则在单元内材料生成面积为A的裂纹消耗的断裂能WC为:
WC=GC*A;
步骤C3,基于上述单元特征长度计算方法得到的特征长度,由WD=WC,得单元的损伤耗散能密度为:
其中,GC为材料的临界应变能释放率,V为单元的体积,A为与宏观裂纹平行的截面积,l*为单元的特征长度,由上述单元特征长度计算方法得到;
步骤D3,基于上面得到的单元损伤耗散能密度,最终得到从微裂纹产生、扩展到宏观裂纹形成过程中,单元损伤所耗散的能量WD为:
WD=Γ*V。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310571786.8A CN103593567B (zh) | 2013-11-13 | 2013-11-13 | 一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310571786.8A CN103593567B (zh) | 2013-11-13 | 2013-11-13 | 一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103593567A CN103593567A (zh) | 2014-02-19 |
CN103593567B true CN103593567B (zh) | 2016-08-17 |
Family
ID=50083704
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310571786.8A Expired - Fee Related CN103593567B (zh) | 2013-11-13 | 2013-11-13 | 一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103593567B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104517036A (zh) * | 2014-12-22 | 2015-04-15 | 华中科技大学 | 一种基于应变统计矩的简支件损伤识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484989A (zh) * | 2006-06-30 | 2009-07-15 | 英特尔公司 | 智能聚合物复合材料对集成电路封装的应用 |
CN101604357A (zh) * | 2009-06-24 | 2009-12-16 | 北京理工大学 | 一种基于材料显微组织图像的介观层次裂纹自动统计方法 |
CN101885596A (zh) * | 2010-06-11 | 2010-11-17 | 北京交通大学 | 一种植物纤维水泥基砌块材料及其制备方法 |
CN102819678A (zh) * | 2012-07-30 | 2012-12-12 | 西北工业大学 | 基于应变场的短纤维增强复合材料纤维取向预测方法 |
CN103020426A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种矩形板中心斜裂纹疲劳扩展寿命预测的简化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5375757B2 (ja) * | 2010-04-20 | 2013-12-25 | 富士電機株式会社 | 太陽電池直列接続構造セル特性解析・評価装置 |
-
2013
- 2013-11-13 CN CN201310571786.8A patent/CN103593567B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484989A (zh) * | 2006-06-30 | 2009-07-15 | 英特尔公司 | 智能聚合物复合材料对集成电路封装的应用 |
CN101604357A (zh) * | 2009-06-24 | 2009-12-16 | 北京理工大学 | 一种基于材料显微组织图像的介观层次裂纹自动统计方法 |
CN101885596A (zh) * | 2010-06-11 | 2010-11-17 | 北京交通大学 | 一种植物纤维水泥基砌块材料及其制备方法 |
CN102819678A (zh) * | 2012-07-30 | 2012-12-12 | 西北工业大学 | 基于应变场的短纤维增强复合材料纤维取向预测方法 |
CN103020426A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种矩形板中心斜裂纹疲劳扩展寿命预测的简化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103593567A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Viola et al. | On static analysis of composite plane state structures via GDQFEM and cell method | |
Li et al. | A layerwise/solid-element method for the composite stiffened laminated cylindrical shell structures | |
Tu et al. | Finite element modeling for bending and vibration analysis of laminated and sandwich composite plates based on higher-order theory | |
CN109885864A (zh) | 一种三维钢桥塔涡激振动计算方法 | |
Lin et al. | Stress intensity factors for corner cracks emanating from fastener holes under tension | |
Bhaskar et al. | An elasticity approach for simply-supported isotropic and orthotropic stiffened plates | |
CN103593567B (zh) | 一种用于复合材料结构失效有限元模拟中单元损伤耗散能量的估计方法 | |
Rong et al. | Static, buckling and free vibration analyses of composite sandwich plate with bi-directional trapezoidal cores using VAM-based reduced plate model | |
Yadav | A XFEM approach for the three-dimensional cracks in piezoelectric material using interaction integral | |
Hosseini-Toudeshky et al. | Finite element crack propagation of adhesively bonded repaired panels in general mixed-mode conditions | |
Zhu et al. | Modelling cohesive crack growth in concrete beams using scaled boundary finite element method based on super-element remeshing technique | |
CN105117516A (zh) | 一种含孔复合材料层合板孔周应力分布的计算方法 | |
Qin et al. | Special elements for composites containing hexagonal and circular fibers | |
Ismail | Multiple crack interactions in Bi-Material plates under mode I tension loading | |
Cao et al. | A VAM-based equivalent model for random vibration of composite sandwich plate with arrowhead-on cores | |
Zhao et al. | Simulation of the dynamic cracking of brittle materials using a nonlocal damage model with an effective strain rate effect | |
Prytz et al. | Numerical simulation of flow around the NREL S826 airfoil at moderate Reynolds number using delayed detached Eddy simulation (DDES) | |
Gao et al. | Aerodynamic performance of a blunt trailing-edge airfoil affected by vortex generators and a trailing-edge wedge | |
Lehmkuhl et al. | Large-eddy simulations of wind turbine dedicated airfoils at high Reynolds numbers | |
Senjanović et al. | Vibration analysis of thick plates: analytical and numerical approaches | |
RU2015156581A (ru) | Алгоритм формирования поверхности для дискретных решетчатых сеток | |
You et al. | A blended model for simulating massive flow separation and laminar-turbulence transition | |
Luangarpa et al. | Evaluation of intensity of singularity for three-materials joints with power-logarithmic singularities using an enriched finite element method | |
Amano et al. | Use of Algebraic-Stress Model for determination of near-wall Reynolds-Stresses in turbulent flow over a flat plate | |
Phan et al. | A Short Review on Numerical Modelling Approaches for Seismic Evaluation Performance of Nuclear Power Plant Structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20161113 |