CN108256211A - 基于abaqus的木材本构关系数值模拟方法 - Google Patents
基于abaqus的木材本构关系数值模拟方法 Download PDFInfo
- Publication number
- CN108256211A CN108256211A CN201810039597.9A CN201810039597A CN108256211A CN 108256211 A CN108256211 A CN 108256211A CN 201810039597 A CN201810039597 A CN 201810039597A CN 108256211 A CN108256211 A CN 108256211A
- Authority
- CN
- China
- Prior art keywords
- formula
- timber
- strain
- plastic
- stress
- 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
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供一种基于ABAQUS的木材本构关系数值模拟方法,属于木材本构研究技术领域。本发明的模拟方法通过嵌入自定义的木材本构子程序VUMAT,能准确描述木材的正交各向异性,并预测木材在复杂受力状态下的破坏模式。
Description
技术领域
本发明属于木材本构研究技术领域,提供了木材在复杂应力状态下应力-应变关系的数值模拟方法。
背景技术
木材是各向异性材料,其本构的复杂主要表现为在受压作用下发生塑性变形,而在拉、剪作用下发生脆性破坏,同时各向拉、压强度均不相等。在木结构设计中,除极少数情况外(如单轴受拉、受压或受剪),木构件受力是复杂的,一个能准确描述木材各向异性本构关系的模型是木结构有限元模型的重要组成部分。
在有限元分析中,木材通常可以被简化为正交各向异性材料。目前,ABAQUS内嵌的材料本构中,针对正交各向异性材料发生受压塑性变形的只有Hill屈服准则,但该本构模型不能同时描述木材顺纹理想塑性和横纹塑性硬化。
此外,ABAQUS内嵌的渐进损伤和失效模型尽管可以描述拉、剪作用下发生的脆性破坏,但会引起顺纹和横纹方向同时失效,这与木材的真实受力情况不符。综上所述,在目前的ABAQUS软件中,其内嵌的材料本构还不能准确描述木材在不同方向受压的塑性变形和拉、剪的脆性破坏。
发明内容
本发明提供了一种基于ABAQUS的木材本构关系数值模拟方法,通过嵌入自定义的木材本构子程序VUMAT,能准确描述木材的正交各向异性,并预测木材在复杂受力状态下的破坏模式。
本发明的技术方案:
基于ABAQUS的木材本构关系数值模拟方法,步骤如下:
步骤一:建立弹性阶段本构方程,表达弹性阶段正交各向异性材料的应力-应变关系,如式(1)所示:
σ=Dε (1)
式中:σ为应力矩阵、ε为应变矩阵、D为刚度矩阵;
展开式(1)得到式(2):
式中,D11=E1(1-v23v32)γ,D22=E2(1-v13v31)γ,D33=E3(1-v21v12)γ,D12=E1(v21+ν23ν31)γ,D13=E1(ν21+ν23ν31)γ,D23=E2(ν32+ν12ν31)γ,D44=G12,D55=G23,D66=G13,γ=(1-ν21ν12-ν23ν32-v13v31-2v21v32v13)-1;
E1,E2,E3分别为L、R、T三个方向的弹性模量;G12,G13,G23分别为L-R、L-T和R-T平面内的剪切模量;vij为泊松比;
步骤二:建立强度准则
采用Yamada-Sun强度准则建立顺纹方向L、横纹径向R和切向T三向独立的强度准则,描述木材的正交各向异性,当满足式(3)时,材料进入非弹性阶段。
式中,X,Y,Z分别为L、R、T方向的抗拉或抗压强度,当σ11,σ22和σ33为拉或压应力时选用相应的抗拉或抗压强度;Sij(i,j=x,y,z)为i-j平面内的抗剪强度;
步骤三:当木材达到强度准则后,进行非弹性阶段木材的本构关系的模拟,具体如下:
(一)建立损伤模型
木材在拉、剪作用下时发生脆性破坏,通过引入连续损伤因子di(i=L,R,T),对应力进行折减来模拟材料的应变软化;当材料出现应变软化时,则会出现应变局域化;此时随着网格的细化,能量的耗散将减少,有限元计算结果将严重依赖网格尺寸的大小,所以采用能量释放率来减小网格的依赖性,损伤因子计算如式(4)所示
式中,为i方向抗拉强度;LC为网格特征长度,Gi(i=L,R,T)为i方向断裂能,LC与Gi的比值即为断裂能释放率;
应力更新得到式(5):
式中,为无损应力;
(二)塑性发展
设定木材横纹受压时,发生应变硬化,设定木材顺纹受压时,发生少许硬化后进入理想弹塑性;
(1)建立硬化法则:
采用Ziegler随动硬化模型,控制屈服面由初始屈服面向最终破坏面的转移来实现硬化,屈服面表达式如式(6)所示:
式中,ni(i=L,R,T)为屈服强度与极限强度的比值,0.9≤nL≤1;bij为背应力,增量dbij的表达式如式(7)所示:
式中,Cα,i(i=L,R,T)为硬化参数;为屈服面转移约束方程,保证最终的屈服面不超过破坏面;
其中,为最终屈服应力,由式(3)可以求得;计算式如式(8)所示:
为等效塑性应变,其表达式如式(9)所示:
(2)流动法则:
当应力达到屈服面后,通过弹塑性径向返回算法,将应力点约束在屈服面上,如式(10)所示:
式中,dεij、和分别为总应变增量、弹性应变增量和塑性应变增量;dλ为塑性比例因子;g为塑性势能函数,由于在本发明中采用关联塑性,令g为屈服函数f;
(A)硬化阶段塑性因子求解:
类推各向同性材料塑性因子推导方法,推出正交各向异性材料的塑性因子表达式为:
式中,
其中,hL、hR、hT为塑性模量,由一致性条件求得
(B)理想弹塑性阶段塑性因子求解:
基于式(3),屈服函数f的表达式为:
由塑性理论的一致性条件,可推得木材顺纹(L)、横纹径向(R)以及切向(T)的塑性因子,三者的表达式如式(15)所示:
式中为试探应力,
(3)塑性应变的求解:
塑性应变求解的公式如下:
屈服面交线上点的塑性应变增量取在该点相交各面法线方向所确定增量的线性组合,即其中Ak为线性组合系数。
本发明的有益效果:本发明通过嵌入自定义的木材本构子程序VUMAT,能准确描述木材的正交各向异性,并预测木材在复杂受力状态下的破坏模式。
附图说明
图1是木材L,R和T三个方向的定义图示。
图2是本发明的本构关系模型。图中:1顺纹抗拉;2横纹抗拉;3顺纹抗压;4横纹抗压;5木材抗剪;6弹性阶段;7应变软化;8理想弹塑性;9应变硬化;10材料断裂;XC,YC和ZC分别为L,R和T三个方向的抗压强度;XT,YT和ZT分别为L,R和T三个方向的抗拉强度;S为木材平均抗剪强度;ni(i=L,R,T)为三个方向的屈服强度与极限强度的比值。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
1.弹性阶段本构方程
弹性阶段正交各向异性材料应力-应变关系式如式(1)所示:
σ=Dε (1)
式中:σ为应力矩阵、ε为应变矩阵、D为刚度矩阵。
展开表达式如式(2):
式中,D11=E1(1-v23v32)γ,D22=E2(1-ν13v31)γ,D33=E3(1-v21v12)γ,D12=E1(v21+v23v31)γ,D13=E1(v21+v23v31)γ,D23=E2(v32+v12v31)γ,D44=G12,D55=G23,D66=G13,γ=(1-v21v12-v23v32-v13v31-2v21v32v13)-1。
E1,E2,E3分别为L、R、T三个方向的弹性模量;G12,G13,G23分别为L-R、L-T和R-T平面内的剪切模量;vij为泊松比;
2.强度准则
采用Yamada-Sun强度准则描述木材的正交各向异性,具体表达式如下式(3)所示:
式中X,Y,Z分别为L,R,T方向的抗拉或抗压强度,当σ11,σ22和σ33为拉或压应力时选用相应的抗拉或抗压强度;Sij(i,j=x,y,z)为i-j平面内的抗剪强度。
3.非弹性阶段木材本构关系描述
木材达到强度准则后(如图2所示),对于拉、剪作用下的脆性破坏通过引入连续损伤因子以及断裂能释放率来模拟。对于压力作用下的塑性变形,其流动法则由Yamada-Sun屈服函数推导得到,而硬化法则采用Ziegler随动硬化模型。顺纹受压作用时的应力-应变关系为理想弹塑性;横纹受压作用时发生应变硬化,当压应变达到极限压应变时,材料在该方向发生断裂,同时该方向失去承载力。
3.1损伤演化
木材在拉、剪作用下时发生脆性破坏,通过引入连续损伤因子di(i=L,R,T),对应力进行折减来模拟材料的应变软化。当材料出现应变软化时,会出现应变局域化。这时,随着网格的细化,能量的耗散将减少,有限元计算结果将严重依赖网格尺寸的大小。基于此,本文中采用了能量释放率来减小网格的依赖性,损伤因子计算如式(4)所示
式中,为i方向抗拉强度;
LC为网格特征长度,Gi(i=L,R,T)为i方向断裂能,两者的比值即为断裂能释放率;
应力更新如式(5)所示:
式中为无损应力。
3.2塑性发展
木材在受压时,发生塑性变形。当顺纹受压时为理想弹塑性;横纹受压时发生硬化。
3.2.1硬化法则
木材横纹受压时,发生应变硬化,顺纹受压为理想弹塑性,本发明为了便于求解背应力,假定顺纹受压也是有少许硬化后进入理想弹塑性的。本发明通过采用Ziegler随动硬化模型,控制屈服面由初始屈服面向最终破坏面的转移来实现硬化,屈服面表达式如式(6)所示:
式中,ni(i=L,R,T)为屈服强度与极限强度的比值,0.9≤nL≤1;bij为背应力,增量dbij表达式如式(7)所示:
式中,Cα,i(i=L,R,T)为硬化参数;
为屈服面转移约束方程,保证最终的屈服面不超过破坏面;其中,为最终屈服应力,由式(3)可以求得。计算式如式(8)所示:
为等效塑性应变,其表达式如式(9)所示:
3.2.2流动法则
当应力达到屈服面后,通过弹塑性径向返回算法,将应力点约束在屈服面上,如式(10)所示:
式中,dεij、和分别为总应变增量、弹性应变增量和塑性应变增量;dλ为塑性比例因子;g为塑性势能函数,由于在本发明中采用关联塑性,可令g为屈服函数f。
3.2.2.1硬化阶段塑性因子求解
类推各向同性材料塑性因子推导方法可推出正交各向异性材料的塑性因子表达式为:
式中,
其中,hL、hR、hT为塑性模量,由一致性条件求得
3.2.2.2理想弹塑性阶段塑性因子求解
基于式(3),屈服函数f的表达式为:
由塑性理论的一致性条件,可推得木材顺纹(L)、横纹径向(R)以及切向(T)的塑性因子,三者的表达式如式(15)所示:
式中为试探应力,
3.2.3塑性应变的求解
塑性应变求解的公式如下:
屈服面交线上点的塑性应变增量取在改点相交各面法线方向所确定增量的线性组合,即其中Ak为线性组合系数。
4.编写材料本构子程序
利用Fortran编写VUMAT材料子程序,通过ABAQUS提供的二次开发接口将VUMAT嵌入主程序中进行木材构件的数值模拟。
Claims (1)
1.一种基于ABAQUS的木材本构关系数值模拟方法,其特征在于,步骤如下:
步骤一:建立弹性阶段本构方程,表达弹性阶段正交各向异性材料的应力-应变关系,如式(1)所示:
σ=Dε (1)
式中:σ为应力矩阵、ε为应变矩阵、D为刚度矩阵;
展开式(1)得到式(2):
式中,D11=E1(1-ν23ν32)γ,D22=E2(1-ν13ν31)γ,D33=E3(1-ν21ν12)γ,D12=E1(ν21+ν23ν31)γ,D13=E1(ν21+ν23v31)γ,D23=E2(ν32+v12v31)γ,D44=G12,D55=G23,D66=G13,γ=(1-v21v12-v23v32-v13v31-2v21ν32v13)-1;
E1,E2,E3分别为L、R、T三个方向的弹性模量;G12,G13,G23分别为L-R、L-T和R-T平面内的剪切模量;νij为泊松比;
步骤二:建立强度准则
采用Yamada-Sun强度准则建立顺纹方向L、横纹径向R和切向T三向独立的强度准则,描述木材的正交各向异性,当满足式(3)时,材料进入非弹性阶段;
式中,X,Y,Z分别为L、R、T方向的抗拉或抗压强度,当σ11,σ22和σ33为拉或压应力时选用相应的抗拉或抗压强度;Sij(i,j=x,y,z)为i-j平面内的抗剪强度;
步骤三:当木材达到强度准则后,进行非弹性阶段木材的本构关系的模拟,具体如下:
(一)建立损伤模型
木材在拉、剪作用下时发生脆性破坏,通过引入连续损伤因子di(i=L,R,T),对应力进行折减来模拟材料的应变软化;当材料出现应变软化时,则会出现应变局域化;此时随着网格的细化,能量的耗散将减少,有限元计算结果将严重依赖网格尺寸的大小,所以采用能量释放率来减小网格的依赖性,损伤因子计算如式(4)所示
式中,为i方向抗拉强度;LC为网格特征长度,Gi(i=L,R,T)为i方向断裂能,LC与Gi的比值即为断裂能释放率;
应力更新得到式(5):
式中,为无损应力;
(二)塑性发展
设定木材横纹受压时,发生应变硬化,设定木材顺纹受压时,发生少许硬化后进入理想弹塑性;
(1)建立硬化法则:
采用Ziegler随动硬化模型,控制屈服面由初始屈服面向最终破坏面的转移来实现硬化,屈服面表达式如式(6)所示:
式中,ni(i=L,R,T)为屈服强度与极限强度的比值,0.9≤nL≤1;bij为背应力,增量dbij的表达式如式(7)所示:
式中,Cα,i(i=L,R,T)为硬化参数;为屈服面转移约束方程,保证最终的屈服面不超过破坏面;
其中,为最终屈服应力,由式(3)求得;计算式如式(8)所示:
为等效塑性应变,其表达式如式(9)所示:
(2)流动法则:
当应力达到屈服面后,通过弹塑性径向返回算法,将应力点约束在屈服面上,如式(10)所示:
式中,dεij、和分别为总应变增量、弹性应变增量和塑性应变增量;dλ为塑性比例因子;g为塑性势能函数,由于在本发明中采用关联塑性,令g为屈服函数f;
(A)硬化阶段塑性因子求解:
类推各向同性材料塑性因子推导方法,推出正交各向异性材料的塑性因子表达式为:
式中,
其中,hL、hR、hT为塑性模量,由一致性条件求得
(B)理想弹塑性阶段塑性因子求解:
基于式(3),屈服函数f的表达式为:
由塑性理论的一致性条件,推得木材顺纹L、横纹径向R以及切向T的塑性因子,三者的表达式如式(15)所示:
式中:为试探应力,
(3)塑性应变的求解:
塑性应变求解的公式如下:
屈服面交线上点的塑性应变增量取在该点相交各面法线方向所确定增量的线性组合,即其中Ak为线性组合系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810039597.9A CN108256211B (zh) | 2018-01-08 | 2018-01-08 | 基于abaqus的木材本构关系数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810039597.9A CN108256211B (zh) | 2018-01-08 | 2018-01-08 | 基于abaqus的木材本构关系数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108256211A true CN108256211A (zh) | 2018-07-06 |
CN108256211B CN108256211B (zh) | 2020-06-16 |
Family
ID=62741091
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810039597.9A Active CN108256211B (zh) | 2018-01-08 | 2018-01-08 | 基于abaqus的木材本构关系数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108256211B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117573A (zh) * | 2018-08-28 | 2019-01-01 | 中国海洋大学 | 一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法 |
CN109858108A (zh) * | 2019-01-14 | 2019-06-07 | 北京市轨道交通建设管理有限公司 | 一种异形结构板极限承载力及破坏模式获取方法及系统 |
CN111950206A (zh) * | 2020-07-21 | 2020-11-17 | 西安理工大学 | 一种微型动压气体止推轴承表面抛物线凹槽的设计方法 |
CN112100820A (zh) * | 2020-08-25 | 2020-12-18 | 华南理工大学 | 一种基于Excel的各向异性屈服函数参数标定方法 |
CN114112676A (zh) * | 2021-12-03 | 2022-03-01 | 中国林业科学研究院木材工业研究所 | 一种木材横纹抗压全时程本构关系的构建方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
CN106227928A (zh) * | 2016-07-20 | 2016-12-14 | 福州大学 | 木材各向异性塑性屈服本构模型的数值模拟方法 |
-
2018
- 2018-01-08 CN CN201810039597.9A patent/CN108256211B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
CN106227928A (zh) * | 2016-07-20 | 2016-12-14 | 福州大学 | 木材各向异性塑性屈服本构模型的数值模拟方法 |
Non-Patent Citations (3)
Title |
---|
M.OUDJENE等: "Elasto-plastic constitutive law for wood behaviour under compressive loadings", 《CONSTRUCTION AND BUILDING MATERIALS》 * |
陈志勇等: "复杂应力状态下木材力学性能的数值模拟", 《计算力学学报》 * |
韩建龙: "梁柱和斗栱构架的振动分析及材性编程", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117573A (zh) * | 2018-08-28 | 2019-01-01 | 中国海洋大学 | 一种考虑复合材料参数三维各向异性非线性的有限元数值模拟方法 |
CN109858108A (zh) * | 2019-01-14 | 2019-06-07 | 北京市轨道交通建设管理有限公司 | 一种异形结构板极限承载力及破坏模式获取方法及系统 |
CN111950206A (zh) * | 2020-07-21 | 2020-11-17 | 西安理工大学 | 一种微型动压气体止推轴承表面抛物线凹槽的设计方法 |
CN111950206B (zh) * | 2020-07-21 | 2024-04-26 | 西安理工大学 | 一种微型动压气体止推轴承表面抛物线凹槽的设计方法 |
CN112100820A (zh) * | 2020-08-25 | 2020-12-18 | 华南理工大学 | 一种基于Excel的各向异性屈服函数参数标定方法 |
CN114112676A (zh) * | 2021-12-03 | 2022-03-01 | 中国林业科学研究院木材工业研究所 | 一种木材横纹抗压全时程本构关系的构建方法 |
CN114112676B (zh) * | 2021-12-03 | 2023-12-19 | 中国林业科学研究院木材工业研究所 | 一种木材横纹抗压全时程本构关系的构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108256211B (zh) | 2020-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108256211A (zh) | 基于abaqus的木材本构关系数值模拟方法 | |
Kong et al. | Modified K&C model for cratering and scabbing of concrete slabs under projectile impact | |
Peerlings et al. | Gradient‐enhanced damage modelling of concrete fracture | |
Kang et al. | Uniaxial ratcheting and fatigue failure of tempered 42CrMo steel: Damage evolution and damage-coupled visco-plastic constitutive model | |
Damjanac et al. | Three-dimensional numerical model of hydraulic fracturing in fractured rock masses | |
Kim et al. | Failure simulation of RC structures under highly dynamic conditions using random lattice models | |
Sirumbal-Zapata et al. | A three-dimensional plasticity-damage constitutive model for timber under cyclic loads | |
CN109902362B (zh) | 考虑纤维增强复合材料复杂非线性行为的弹粘塑性本构模型构建方法 | |
CN102364489A (zh) | 木材复杂各向异性本构关系模型的数值模拟方法 | |
Aráoz et al. | Modeling concrete like materials under sever dynamic pressures | |
Jia et al. | Numerical study on the mechanical behavior of a polyurethane adhesive under high strain rate | |
Kontou et al. | Viscoplastic response and creep failure time prediction of polymers based on the transient network model | |
Chen et al. | Three-dimensional modelling on the impact fracture of glass using a GPGPU-parallelised FDEM | |
Shin et al. | Experimental, theoretical and numerical studies on plain concrete fracture in the low-strain rate regime—A state-of-the-art review | |
Sugihardjo et al. | FE model of low grade rubber for modeling housing’s low-cost rubber base isolators | |
Xu et al. | A comparative study of modelling RC slab response to blast loading with two typical concrete material models | |
Vaz et al. | A note on the thermal effects upon a Gurson-type material model | |
Xu et al. | Numerical modelling for reinforced concrete response to blast load: understanding the demands on material models | |
Kang et al. | Failure analysis of R/C columns using a triaxial concrete model | |
Zhang et al. | Progressive failure of bamboo-steel adhesive bonding interface subjected low-energy impact and tension in sequence | |
Abdia et al. | Rock failure analysis based on a coupled elastoplastic-logarithmic damage model | |
Wang et al. | Plasticity tool for predicting shear nonlinearity of unidirectional laminates under multiaxial loading | |
Ruggiero et al. | ADI 1050-6 mechanical behavior at different strain rates and temperatures | |
Tsyss et al. | Analysis of the applicability in the ansys interface of finite element models for modeling composite rubber-cord designs of vibration-isolating pipes | |
Nikolic et al. | Modelling of internal fluid flow in cracks with embedded strong discontinuities |
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 |