多尺度方法在力学中的应用.docx

1、多尺度方法在力学中的应用多尺度方法在力学中的应用1.背景概述多尺度科学是一门研究各种不同长度或者时间尺度相互耦合现 象的一门科学。 多尺度科学的研究领域十分宽广， 涵盖的学科之多难 以一一罗列。在诸如流体动力学、复合材料力学、生物力学、环境科 学、化学、地质学、气象学和高能物理之类的各门科学中，多尺度科 学及其相应的方法发挥着相当重要的作用。 正如同随机现象和非线性 科学受到了广泛的重视一样， 多尺度科学因其处于当代科学的许多极 富挑战性问题的核心地位，未来的发展前途不可限量。在材料科学领域中， 材料的动态特性就是多尺度的问题。 金属的 塑性变形问题是从位错流动着手研究的， 但是位错理论本身并

2、不能预 测塑性流动率和屈服强度位错与晶界、 点缺陷以及原子振动之间 的相互作用才是导致诸如应变强化和材料强度特性动态变化等现象 的主导因素。所以将固体的微观结构与原子层次的组成成分相结合来 预测固体材料的宏观特性， 就是材料科学的宏伟理想， 并可期望达到 人工设计材料的终极目标。在气象学领域中，在大气环流模拟中计算尺度的典型数量级为100km，但是局部降水量、水汽含量以及某些风暴系统的数量级则要 小得多，因而必须在较小尺度层次上进行模拟， 这也是典型的多尺度 问题，应该用多尺度方法来处理。必须说明的是， 正是因为多尺度科学广泛的应用背景， 多尺度方 法作为一种研究的手段和方法， 在各种截然不同

3、的研究领域的应用过 程中，往往与该研究领域的具体背景相结合，具有一定的特殊性。从 算法的角度来说， 与线性方程组的解法等常规算法不同的是， 目前多 尺度方法本身没有固定的算法格式， 它所体现的更多的是一种研究的 需求和应用的思想， 在程序上的实现必须结合具体的研究模型， 这将 在下文中得到充分的体现。2.多尺度的力学分析方法在多尺度的分析方法中已经发展了若干力学分析的方法， 目前比 较典型算法有： 宏观细观平均化计算方法、 材料强度的统计计算方 法等。下面将详细介绍这两种方法。2 . 1宏观一细观平均化计算方法典型的宏观细观平均化算法是： 利用材料的细观周期性的胞元 模型和强调宏观与细观之间相

4、连接的广义自洽模型相结合所进行的 计算。首先讨论胞元模型。 胞元是材料的一个基本结构，它嵌含材料 的细观几何和相结构的要素。 就复合材料来说， 胞元应嵌含颗粒形状、 颗粒百分比、颗粒分布几何、基本结构、界面状况等要素。自洽方法 是考虑宏观和细观交互作用的研究方法。 广义自洽方法则是将平均化 的小尺度的胞元与大尺度的宏观等效介质进行自洽连接。把宏观-细观平均化计算方法在多尺度思想上作一定的推广， 即并不要求达到细观尺度，而是相对于宏观大尺度来说，胞元尺寸构成 相对的小尺度。一个典型的例子就是复合材料等效模量计算中常用的 复合圆柱模型。下面就以复合圆柱模型（图1）为例，给出一个多尺 度计算的具体算

5、例。图1 复合圆柱模型（坐标图）我们的问题是计算横向剪切模量 03。在图2中我们为了视觉上 的清晰起见，夸大了胞元的尺寸。实际上，我们只取出宏观等效均匀 介质中周期性分布的一个微小胞元，虽然胞元本身并没有达到细观尺 度，但是可以肯定的是，胞元的尺寸与周围的宏观等效均匀介质相比 起来，在尺度上已经相差很大的数量级。而在计算的过程中，我们将分别考虑较小尺度的胞元内的物理量和胞元周围较大尺度的等效均匀介质中的物理量，然后再通过广义自洽方法将平均化的小尺度的胞元与大尺度的宏观等效均匀介质进行自洽连接。图2 复合圆柱模型（横截面图）在大尺度的宏观等效均匀介质中，位移场是:b 2r b b3Ure ( 1

6、)-a3 飞 C3】COS2d4匕3 b r rb 2r b b3 .口丁 ( -1)-a3 =C3】sin2,423 b r r在平均化的小尺度的胞元中，基体和纤维内的位移场分别是:其中连接小尺度的胞元和大尺度的宏观等效均匀介质的条件在不同的问题中可能不尽相同，在我们以上考虑的这个问题中，这一条件就 是连续性条件，即二r, ；中Ur, 在界面r二a, b上连续。这样就得到 了一系列方程，另外再补充其他一些方程，例如应用能量原理得到的 方程，可以想见，以上问题就归结为解一个非线性方程组，而我们所 要求的 鼻只是其中的一个未知量。以上这个具体问题的求解比较特 殊，消去不相关的未知元ai、a2等，

7、只留下23，我们得到一个二次 方程:A(子)2 b(屮)C = 0m m其中A，B，C是材料常数的函数，可以由给定的具体材料来决定， 其形式不是本文所要讨论的重点，故不再赘述。通过复合圆柱模型，我们可以总结出宏观-细观平均化算法的流 程图，如图3所示图3.宏观一细观平均化算法的流程图2. 2材料强度的统计计算方法在材料强度的统计计算方法中，应用比较成熟的是有关带有强相 互作用共线裂纹的脆性材料强度的统计计算。 这个算法与宏观-细观 平均化算法有所不同，宏观-细观平均化算法连接了宏观尺度和细观 尺度，换句话说，我们可以很明确地看到在尺度上的大幅度的跳跃。 而带有强相互作用共线裂纹的脆性材料强度的

8、统计计算这一算法， 它注重的则是裂纹或者间距在尺度上小幅度的涨落， 所以我们用统计理论来处理这种涨落这个算法的核心思想是： 用微裂纹长度和间距的统计分布来描述 它们在尺度上的涨落， 进而来确定材料的统计强度。 下面将通过一个 具体的脆性材料强度统计计算的例子来介绍这个算法。含有共线裂纹的无限大平板的破坏几率的统计分析算例： 问题描述：无限大平板，包含N个共线裂纹，无穷远处作用有均匀拉应力 (7a：半裂纹长c：裂纹间距a和c都是统计变量见图4(a),它们的统计分布用函数f(a)、p(c)来 表示，都是正态分布，c-,c+,a-,a+是 c和a的下届和上届。求解应力强度因子非常繁琐, 为简化问题,

9、 我们主要考虑相邻的两个 微裂纹之间的强相互作用，这两个相邻的微裂纹的长度用 a以及a表示,其它裂纹用周期性分布的裂纹代替(裂纹长度及间距分别是 2ao以及co),如图4 (b)所示。(b)图4 (a)裂纹长度不同，间距不同(b)只考虑临近的两裂纹之间的强相互作用，简化成远场裂纹为周期性分布A 点的 SIF(stress intensify facto):皿F二岂旦0lao a0 a0 a0 丿F是无量纲函数为分析简便起见，下面具体计算一个特例：N个长度相同的裂纹，间距不同（即aa=o）则f（a）是Dirac delta函数S （a-比）当aa =a0时，c的分布函数p（c）是一个正态分布，c

10、的取值范围为（c-,c+）,平均间距为q =c= cp（c）dc oLc考虑到当aa二o时，Ka是一个c/ao单调递减函数，见图50.125OOQ 0.25 0 50 075 100图5应力强度因子Ka曲线 （N = 100, aa=时）（T的临界值（T th可被给定为Gh = K/F（=1,1単）:h a0 ao aoKic是基体断裂刚度，对给定的。二有一个与之对应的临界裂纹间距Ccr1,满足F （血,1,1単）二 K ao ao 、- a0由以上两式可知: 1 . 如果 C（T thm，应力小于使基体断裂的最小应力，裂纹不会扩展。2. 如果 J(7 J，贝恫距小于Ccr1的相邻裂纹将连通。

11、裂纹连通后，裂纹长度和间距的分布函数 p(c)和f(a)将改变B(C)=2&)1 -afl (a) = 5(a a。)+ P(2a_4a)H(2a 4ao cH(2a 4a。 &)1_a 1_a上式中：cCr p(c)dc代表连通概率 cH 代表 Heaviside step方程fi(a )的第一部分代表那些没有连通的裂纹，它们的长度仍然是 2ao第二部分代表连通的裂纹部分第一次扩展后的微裂纹长度及间距的期望值分别是2a0 2ccr c -2ai =2 af1 (a)da ci 1 cp-i (c)dca0 ccr这是一个循环的过程第一次扩展后的微裂纹右端点 SIF为Kright =严隔卩二色

12、,1,里0 a1 a1 /如果Kright Kic，该裂纹将和临近的期望长度为21的裂纹连通。与2a1相关联，我们又可以找到一个临界裂纹间距，用 Ccr2表示，满足,_ ( 2 一 、Ke 呱&冬色,1,里a1 a1 a1 ；上式表明，如果裂纹间距c在(Ccr1, Ccr2)范围内，这个长度为2的裂纹将和临近的期望长度为 2；1的裂纹连通，这一扩展过程的概率为当这一步骤重复n次后裂纹长度变为n2an =4a0 2(n -1)ai c( ck(c1)k=2第k次扩展后的裂纹间距期望值为- J：cpi(c)dcSrck 卫 k = 2,3,./ Pi(c)dcr其中C0代表第k次扩展的临界裂纹长度

13、，由下式给出裂纹间的总连接次数 M可用下式求得1 閃2aM 4 ： K / ；” - aMJI对于长度为2ai的裂纹，如果M = 1,那么它的破坏概率为1;否则等 M cn于成功链接概率Pf(G)：丨丨Pi(c)dcn=2 &进而该裂纹的存活概率为Ps(G)=1 -Pf (cj根据WLT,我们又可得到存活概率 FS(cJ - 1 - Pf(q)VP(c1) j expNPjGpgdc累积存活概率为最终，我们可以得到对于含有 N个微裂纹的脆性材料的破坏概率为(0 th 时)C1Rail =1-expNPf(c)p(c)dc?利用直接数值模拟进行校验0.40 0.45 0.50 0.55 0 60

14、 Q.55 0.70nomorlized strength图6统计预测与直接数值模拟的对比（其中Ng意义为直接数值模拟中，在同样的s和N下，生成Ng个不同裂纹分布状态进行计算，显然 Ng越大越精确）另外，Weibull提出用如下带三个参数（m, co）的分布函数描述 脆性材料的强度mW（G =1 - exp（-（ -）%W（ c ）是应力为c时的破坏概率（横轴为KJ 屁）,c -表明累积破 坏概率开始增长的位置，c 0标示了破坏概率曲线的过渡区的尺度，无量纲的参数m （称为Weibull模量）描述了脆性材料中的裂纹分布 特性。上式可化为1lnln（rWn）En（ u）*。）即上式在一个InIn

15、-In的Weibull图中为一条斜率为 m的直线如果前面分析的累积破坏概率函数可以用 Weibull分布近似，那么它应该在InIn-In的Weibull图中呈直线。我们可以把数据在 Weibull图 中标示出测定 Weibull模量m,也可以估计m,c u, 0这三个参数与s 和N之间的关联。这就是采用统计方法对材料强度进行多尺度分析的例子， 例子中推算 出了材料的破坏概率，并利用直接数值模拟进行校验，最后用分布进 行拟和。以上详细介绍的是宏观-细观平均化计算、 材料强度的统计计算这两种力学上的多尺度分析方法。 最后还需要强调的是，正如我们前文所说，多尺度方法是迎合研究过程中的具体需要而产生的一种计算 思想，它本身没有固定的计算格式，不论是在力学方面，还是在其他 领域，多尺度方法的应用都必须结合其具体的研究模型来展开。参考文献1.“材料的多尺度力学与