多尺度方法在力学中的应用Word文档下载推荐.docx
- 文档编号:15814235
- 上传时间:2022-11-16
- 格式:DOCX
- 页数:11
- 大小:247.48KB
多尺度方法在力学中的应用Word文档下载推荐.docx
《多尺度方法在力学中的应用Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《多尺度方法在力学中的应用Word文档下载推荐.docx(11页珍藏版)》请在冰豆网上搜索。
首先讨论胞元模型。
胞元是材料的一个基本结构,它嵌含材料的细观几何和相结构的要素。
就复合材料来说,胞元应嵌含颗粒形状、颗粒百分比、颗粒分布几何、基本结构、界面状况等要素。
自洽方法是考虑宏观和细观交互作用的研究方法。
广义自洽方法则是将平均化的小尺度的胞元与大尺度的宏观等效介质进行自洽连接。
把宏观-细观平均化计算方法在多尺度思想上作一定的推广,即
并不要求达到细观尺度,而是相对于宏观大尺度来说,胞元尺寸构成相对的小尺度。
一个典型的例子就是复合材料等效模量计算中常用的复合圆柱模型。
下面就以复合圆柱模型(图1)为例,给出一个多尺度计算的具体算例。
图1•复合圆柱模型(坐标图)
我们的问题是计算横向剪切模量03。
在图2中我们为了视觉上的清晰起见,夸大了胞元的尺寸。
实际上,我们只取出宏观等效均匀介质中周期性分布的一个微小胞元,虽然胞元本身并没有达到细观尺度,但是可以肯定的是,胞元的尺寸与周围的宏观等效均匀介质相比起来,在尺度上已经相差很大的数量级。
而在计算的过程中,我们将
分别考虑较小尺度的胞元内的物理量和胞元周围较大尺度的等效均
匀介质中的物理量,然后再通过广义自洽方法将平均化的小尺度的胞
元与大尺度的宏观等效均匀介质进行自洽连接。
图2•复合圆柱模型(横截面图)
在大尺度的宏观等效均匀介质中,位移场是:
b2rbb3
Ure[('
1)-a3飞C3】COS2d
4匕3brr
b2rbb3.
口丁[('
-1)-a3=C3】sin2,
4^23brr
在平均化的小尺度的胞元中,基体和纤维内的位移场分别是:
其中
连接小尺度的胞元和大尺度的宏观等效均匀介质的条件在不同
的问题中可能不尽相同,在我们以上考虑的这个问题中,这一条件就是连续性条件,即二r,;
「中Ur,在界面r二a,b上连续。
这样就得到了一系列方程,另外再补充其他一些方程,例如应用能量原理得到的方程,可以想见,以上问题就归结为解一个非线性方程组,而我们所要求的鼻只是其中的一个未知量。
以上这个具体问题的求解比较特殊,消去不相关的未知元ai、a2等,只留下^23,我们得到一个二次方程:
A(子)2b(屮)C=0
mm
其中A,B,C是材料常数的函数,可以由给定的具体材料来决定,其形式不是本文所要讨论的重点,故不再赘述。
通过复合圆柱模型,我们可以总结出宏观-细观平均化算法的流程图,如图3所示
图3.宏观一细观平均化算法的流程图
2.2材料强度的统计计算方法
在材料强度的统计计算方法中,应用比较成熟的是有关带有强相互作用共线裂纹的脆性材料强度的统计计算。
这个算法与宏观-细观平均化算法有所不同,宏观-细观平均化算法连接了宏观尺度和细观尺度,换句话说,我们可以很明确地看到在尺度上的大幅度的跳跃。
而带有强相互作用共线裂纹的脆性材料强度的统计计算这一算法,它
注重的则是裂纹或者间距在尺度上小幅度的涨落,所以我们用统计理
论来处理这种涨落
这个算法的核心思想是:
用微裂纹长度和间距的统计分布来描述它们在尺度上的涨落,进而来确定材料的统计强度。
下面将通过一个具体的脆性材料强度统计计算的例子来介绍这个算法。
含有共线裂纹的无限大平板的破坏几率的统计分析算例:
问题描述:
无限大平板,包含N个共线裂纹,无穷远处作用有均匀拉应力(7^
a:
半裂纹长
c:
裂纹间距
a和c都是统计变量[见图4(a)],它们的统计分布用函数f(a)、p(c)来表示,都是正态分布,c-,c+,a-,a+是c和a的下届和上届。
求解应力强度因子非常繁琐,为简化问题,我们主要考虑相邻的两个微裂纹之间的强相互作用,这两个相邻的微裂纹的长度用a‘以及a'
表示,其它裂纹用周期性分布的裂纹代替(裂纹长度及间距分别是2ao以及co),如图4(b)所示。
(b)
图4(a)裂纹长度不同,间距不同
(b)只考虑临近的两裂纹之间的强相互作用,简化成远场裂
纹为周期性分布
A点的SIF(stressintensifyfacto):
皿F二岂旦0
laoa0a0a0丿
F是无量纲函数
为分析简便起见,下面具体计算一个特例:
N个长度相同的裂纹,间
距不同(即a'
a'
='
o)
则f(a)是Diracdelta函数S(a-比)当a'
=a0时,c的分布函数p(c)是
一个正态分布,c的取值范围为(c-,c+),平均间距为q=c=cp(c)dco
Lc
考虑到当a'
二'
o时,Ka是一个c/ao单调递减函数,见图5
^0.125
OOQ0.25050075100
图5应力强度因子Ka曲线(N=100,a'
时)
(T的临界值(Tth可被给定为
Gh=K^/F(=1,1単)
"
•:
ha0aoao
Kic是基体断裂刚度,对给定的。
二有一个与之对应的临界裂纹间距
Ccr1,满足
F(血,1,1単)二K—
aoao、-■■■a0
由以上两式可知:
1.如果C(Tthm,应力小于使基体断裂的最小应力,裂纹不会
扩展。
2.如果J》(7J,贝恫距小于Ccr1的相邻裂纹将连通。
裂纹连通后,裂纹长度和间距的分布函数p(c)和f(a)将改变
B(C)=¥
^2—&
)
1-a
fl(a)=5(a—a。
)+P(2a_4a°
)[H(2a—4ao—c」—H(2a—4a。
—&
)]
1_a1_a
上式中:
cCrp(c)dc代表连通概率c
H代表Heavisidestep方程
fi(a)的第一部分代表那些没有连通的裂纹,它们的长度仍然是2ao
第二部分代表连通的裂纹部分
第一次扩展后的微裂纹长度及间距的期望值分别是
2a0'
2ccrc-
2ai=2af1(a)daci1cp-i(c)dc
a0'
ccr
这是一个循环的过程
第一次扩展后的微裂纹右端点SIF为
Kright=严隔卩「二色,1,里
0a1a1/
如果Kright>
Kic,该裂纹将和临近的期望长度为2^1的裂纹连通。
与2a1相关联,我们又可以找到一个临界裂纹间距,用Ccr2表示,满足
_(2一、
Ke"
呱&
冬色,1,里
>
a1a1a1;
上式表明,如果裂纹间距c在(Ccr1,Ccr2)范围内,这个长度为2®
的裂
纹将和临近的期望长度为2;
1的裂纹连通,这一扩展过程的概率为
当这一步骤重复n次后裂纹长度变为
n
2an=4a02(n-1)aic('
ck(c1)
k=2
第k次扩展后的裂纹间距期望值为
-J:
cpi(c)dc
Sr
ck卫k=2,3,...
/Pi(c)dc
r
其中C0代表第k次扩展的临界裂纹长度,由下式给出
裂纹间的总连接次数M可用下式求得
1閃2
aM4:
:
—K©
/;
”-aM
JI
对于长度为2ai的裂纹,如果M=1,那么它的破坏概率为1;
否则等Mcn
于成功链接概率Pf(G):
丨丨Pi(c)dc
n=2&
进而该裂纹的存活概率为Ps(G)=1-Pf(cj
根据WLT,我们又可得到
存活概率FS(cJ-'
1-Pf(q)VP(c1)jexp—NPjGpgdc^
累积存活概率为
最终,我们可以得到对于含有N个微裂纹的脆性材料的破坏概率为
(0th时)
C1
Rail=1-exp「N「Pf(c)p(c)dc?
利用直接数值模拟进行校验
0.400.450.500.55060Q.550.70
nomorlizedstrength
图6统计预测与直接数值模拟的对比
(其中Ng意义为直接数值模拟中,在同样的s和N下,
生成Ng个不同裂纹分布状态进行计算,显然Ng越大越精确)
另外,Weibull提出用如下带三个参数(m,co)的分布函数描述脆性材料的强度
m
W(G=1-exp(-(-))
%
W(c)是应力为c时的破坏概率(横轴为KJ屁),c-表明累积破坏概率开始增长的位置,c0标示了破坏概率曲线的过渡区的尺度,
无量纲的参数m(称为Weibull模量)描述了脆性材料中的裂纹分布特性。
上式可化为
1
lnln(rWn)En(—u)"
*。
即上式在一个InIn-In的Weibull图中为一条斜率为m的直线
如果前面分析的累积破坏概率函数可以用Weibull分布近似,那么它
应该在InIn-In的Weibull图中呈直线。
我们可以把数据在Weibull图中标示出测定Weibull模量m,也可以估计m,cu,0这三个参数与s和N之间的关联。
这就是采用统计方法对材料强度进行多尺度分析的例子,例子中推算出了材料的破坏概率,并利用直接数值模拟进行校验,最后用分布进行拟和。
以上详细介绍的是宏观-细观平均化计算、材料强度的统计计算
这两种力学上的多尺度分析方法。
最后还需要强调的是,正如我们前
文所说,多尺度方法是迎合研究过程中的具体需要而产生的一种计算思想,它本身没有固定的计算格式,不论是在力学方面,还是在其他领域,多尺度方法的应用都必须结合其具体的研究模型来展开。
参考文献
1.“材料的多尺度力学与
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 尺度 方法 力学 中的 应用
![提示](https://static.bdocx.com/images/bang_tan.gif)