一种梯度材料宏观等效弹性模量的计算方法

标题：一种梯度材料宏观等效弹性模量的计算方法

摘要：本发明涉及梯度材料力学分析技术领域，具体涉及一种梯度材料宏观等效弹性模量的计算方法，以解决现有梯度材料的宏观等效弹性模量获取效率低的问题。计算方法包括如下步骤：得到均匀介质层的平均应力与各相材料内部应力的关系式；得到夹杂相的扰动应变差值ε′与等效本征应变间ε*的关系式；得到扰动应力差σ′、基体相的扰动应力以及基体相的扰动应变；确定各均匀介质层平均应力和平均应变的关系式；得到各均匀介质层的等效模量L；得到各均匀介质层的等效体积模量K和等效剪切模量G。本发明能够获取梯度区宏观等效弹性模型及泊松比数值，从而得到梯度区处材料性能分布情况，以便于进行细节应力分析，能够极大地缩短对梯度区进行细节分析的建模工作量。

申请号：CN201510654378.8

申请日：2015/10/10

申请人：中国飞机强度研究所

首项权利要求：一种梯度材料宏观等效弹性模量的计算方法，其中，梯度材料形成于两种预定的各向同性材料的过渡区中，所述过渡区材料的力学性能沿一预定方向呈梯度变化，且所述过渡区是具有颗粒夹杂微结构形式的多个均匀介质层的集合体，其特征在于，计算方法包括如下步骤：步骤一、根据场量平均理论得到均匀介质层的平均应力与各相材料内部应力的关系式为：σ ‾ =(1-V)σ 0+Vσ 1—(1), 其中，V表示所述梯度材料中夹杂相的体积分数，各相材料内部应力包括基体相的平均应力σ0以及所述夹杂相的平均应力σ1；所述基体相中平均应力σ0与所述基体相中平均应变ε0关系式为：σ 0=σ ‾ +σ ~=L0(ϵ 0+ϵ ~)—(2), 其中，表示所述基体相的扰动应力，L0表示所述基体相的刚度张量，表示所述基体相的扰动应变；所述夹杂相中平均应力σ1与所述基体相中平均应变ε0关系式为：σ 1=σ ‾ +σ ~+σ ′ =L1(ϵ 0+ϵ ~+ϵ ′ )—(3), 其中，σ′表示所述夹杂相的扰动应力相对于所述基体相的扰动应力的扰动应力差值，ε′表示所述夹杂相的扰动应变相对于所述基体相的扰动应变的扰动应变差值，L1表示所述夹杂相的刚度张量；步骤二、根据Eshelby等效夹杂理论引入等效本征应变ε*，进一步得到所述夹杂相中平均应力σ1与所述基体相中平均应变ε0关系式为：σ 1=L1(ϵ 0+ϵ ~+ϵ ′ )L0(ϵ 0+ϵ ~+ϵ ′ -ϵ *)—(4); 步骤三、得到所述夹杂相的扰动应变差值ε′与等效本征应变间ε*的关系式为：ε′＝Sε*???(5)，其中，S为Eshelby张量；步骤四、得到所述扰动应力差σ′、所述基体相的扰动应力以及所述基体相的扰动应变分别为：σ′＝L0(ε′-ε*)＝L0(S-I)ε*???(6)，σ ~=-Vσ ′ =-VL0(S-I)ϵ *—(7), ϵ ~=-V(S-I)ϵ *—(8), 其中，I表示一个与S同阶的单位矩阵；步骤五、得到所述等效本征应变ε*为：ε*＝{L0+(L1-L0)[VI+(1-V)S]}-1(L0-L1)ε0???(9)；步骤六、确定各均匀介质层平均应力和平均应变的关系式为：ϵ ‾ =(1-V)ϵ 0+Vϵ 1=ϵ 0+Vϵ *={I+V[ L0+(L1-L0)(VI+(1-V)S)] -1(L0-L1)}ϵ 0={I+V[ L0+(L1-L0)(VI+(1-V)S)] -1(L0-L1)}L0-1σ ‾ —(10); 步骤七、根据均匀介质层平均应力和平均应变的关系式，得到各均匀介质层的等效模量L为：L＝{I+V[L0+(L1-L0)(VI+(1-V)S)]-1(L0-L1)}-1L0???(11)；步骤八、各均匀介质层的等效体积模量K和等效剪切模量G分别为：K=K0+(K1-K0)(V+9(1-V)K0Kp)1+9(1-V)(K1-K0)Kp—(12), G=G0+(V+4(1-V)G0Gp)(G1-G0)1+4(1-V)(G1-G0)Gp—(13), 其中，K0表示所述基体相的等效体积模量，G0表示所述基体相等效剪切模量，K1表示所述夹杂相的等效体积模量、G0表示所述夹杂相等效剪切模量，Kp、Gp为中间量，Kp=13(4G0+3K0), Gp=3(2G0+K0)10G0(4G0+3K0).

专利类型：发明申请