均匀化理论和多尺度方法课件.ppt

第六章均匀化理论和多尺度方法,高等复合材料力学,Advanced Mechanics of Composite Materials,陈玉丽 航空科学与工程学院,-,第六章均匀化理论和多尺度方法高等复合材料力学Advanc,6.1 引言,在考察实际复合材料微结构状态变量和材料系数的演化时，由于热载荷和机械载荷都是施加在宏观结构层面，所以研究采用的细观力学模型必须能够把细观响应和宏观行为联系起来。单胞模型通过在非均匀结构中提取出一个代表性体积单元(RVE)从而可以求得有效的材料响应和演化过程。这里假设微结构是周期性重复排列的单胞，与复合材料的宏观尺寸相比，它的不均匀性是很小的，此种类型的材料被称作具有周期性微观结构的复合材料（第三章 ） 。但是，单胞法还是存在许多不足。周期性假设用于预测最优材料性能非常有效，然而实际的非均匀材料很少具有完全的周期性微结构，宏观结构上不同的点可能具有不同的微结构形态。这种假设在处理复杂载荷条件下非线性非均匀结构变形问题时也存在不足。为了解决上述问题，单胞模型应该包含大的区域，采用大的模型。,2,-,6.1 引言在考察实际复合材料微结构状态变量和材料系数的演化,6.1 引言,20世纪70年代，学者们在研究非均匀材料时引入了一种替代的数学方法，Benssousan和Sanchez-Palencia等称之为均匀化理论。这种方法用于分析具有两个或者多个尺度的物质系统，它可以把含有第二相空间的细观尺度和整体结构上的宏观尺度联系起来。通过对位移和应力场进行渐进展开以及适当的变分原理，均匀化方法不仅可以求出等效的（均匀化）材料常数，而且可以得到两个尺度上的应力和应变分布。相对于单胞法，这种方法的优点在于不必作全局的周期性假设，在宏观结构的不同点可以有不同的微结构。然而，这种分析通过引入空间重复排列单胞作了局部周期性假设。Toledano和Murakami，Guedes和Kikuchi以及Devries等成功地把有限元方法和均匀化方法结合起来用于分析复合材料的线弹性问题。在这些研究当中，通过计算机模拟宏观结构的平均应力和应变场得到了全局的响应，同时借助局部应力和应变场的描述得到了微结构的行为。,3,-,6.1 引言20世纪70年代，学者们在研究非均匀材料时引入了,6.2 多尺度模型,一具有周期性结构的复合材料弹性体 ，受体力f，边界 t 上受表面力t，边界 u 上给定位移边界条件。宏观某点 x 处的细观结构可以看成是非均匀单胞在空间中周期性重复堆积而成。单胞的尺度 y 相对于宏观几何尺度为小量。,x,f,y,4,-,6.2 多尺度模型一具有周期性结构的复合材料弹性体 ，受体,6.2 多尺度模型,x,y=x/,宏观尺度：,微观尺度：,例如：宏观尺度为 m，微观尺度为 nm， = 10-9,实际为 1m 的尺寸，即 x=1 (m)， 在微观尺度下 y=x/= 109 (nm)实际为1nm的尺寸，即 y=1 (nm)，在宏观尺度下 x=y= 10-9 (m),y,5,-,6.2 多尺度模型0432104321xy=x/宏观尺度：,6.2 多尺度模型,对于非均匀的复合材料，当宏观结构受外部作用时，位移和应力等结构场变量将随宏观位置的改变而不同。同时由于细观结构的高度非均匀性，使得这些结构场变量在宏观位置 x 非常小的邻域 内也会有很大变化。因此所有变量都假设依赖于宏观与细观两种尺度，即：,上标 表示该函数具有两尺度的特征。,Y-周期性：微观单胞的周期为Y,6,-,6.2 多尺度模型对于非均匀的复合材料，当宏观结构受外部作用,6.2 多尺度模型,在 中，弹性张量 和柔度张量 分别为 假设应力场和位移场都满足平衡方程、几何方程和本构方程，有 其中 是细观坐标系 y 中的具有 Y-周期的位移场。同时，在给定力边界和给定位移边界分别满足,7,-,6.2 多尺度模型在 中，弹性张量,均匀化方法,3）以傅里叶变换为基础的多尺度方法,2）泰勒级数近似法（Taylor Series Approximation）,1）渐进展开法（Asymptotic expansion）,8,-,均匀化方法8-,6.3 渐进展开法,在均匀化理论中， Y-周期位移场可以近似为宏观坐标 x 的展开式，渐进展开是其中比较常用的一种展开方法中，其展开形式为：,注意到任意一个依赖于两个尺度的函数 对宏观坐标 x 的偏微分为,应变张量,Asymptotic expansion,9,-,6.3 渐进展开法在均匀化理论中， Y-周期位移场可以近似为,6.3 渐进展开法,代入本构方程，可得应力场的渐进展开式：,其中,将应力的渐进展开式代入平衡方程，有,10,-,6.3 渐进展开法代入本构方程，可得应力场的渐进展开式：其中,6.3 渐进展开法,令 i (i=-2,-1,0,1) 的系数为零，得到一系列控制方程：,（1）,（2）,（3）,11,-,6.3 渐进展开法令 i (i=-2,-1,0,1),6.3 渐进展开法,根据Y-周期性，可以证明（Devries et al. 1989）,可以得到,（2）式,（1）式,其中,细观平衡方程,细观本构方程,12,-,6.3 渐进展开法根据Y-周期性，可以证明（Devries,6.3 渐进展开法,在Y 内积分，有,均匀化弹性常数,（3）式,均匀化的宏观平衡方程,令,宏观弹性问题的解,13,-,6.3 渐进展开法在Y 内积分，有均匀化弹性常数（3）式均匀,6.3 渐进展开法,x,y=x/,宏观,微观,尺 度,z=x/2,对位移渐进展开,代入平衡方程,得到控制方程,得到均匀化方程,动态问题怎么办？,14,-,6.3 渐进展开法xy=x/宏观 微观尺 度z=x/2,6.4 含时间的渐进展开（1）,弹性动力学问题:,参考文献：Fish, J. and Chen, W. (2001). Higher-Order Homogenization of Initial/ Boundary-Value Problem. J. Eng. Mech., 127(12), 12231230.,x,y=x/,15,-,6.4 含时间的渐进展开（1）弹性动力学问题:参考文献：Fi,6.4 含时间的渐进展开（1）,令 i (i=-2,-1,0,1) 的系数为零，得到一系列控制方程：,不同阶的应力为：,16,-,6.4 含时间的渐进展开（1）令 i (i=-2,-1,6.4 含时间的渐进展开（1）,17,-,6.4 含时间的渐进展开（1）17-,6.4 含时间的渐进展开（1）,线性问题通解：,代入原式得：,如何求解 L(y) ？,提示：,1. 周期性（Periodicity）：,2. 连续性（Continuity）：,3. 正交性（Normalization）：,请求出L(y)（分段表达），进而求出,18,-,6.4 含时间的渐进展开（1） 线性问题通解： 代入原式得：,6.4 含时间的渐进展开（1）,均匀化后的材料性质与静态问题是一致的。因此，0阶问题是无色散的。为了反映波的色散效应（dispersion effect），必须考虑更高阶的项。,请证明,19,-,6.4 含时间的渐进展开（1）均匀化后的材料性质与静态问题是,6.4 含时间的渐进展开（1）,其中，Ed 表征了非均匀对宏观行为的影响。 Ed具有如下特性： 1）正比于单元尺寸的平方； 2）均匀材料（=0 或 =1） Ed =0。 右端力项正比于宏观应变梯度，梯度越小，右端项越小。,20,-,6.4 含时间的渐进展开（1） 其中，Ed 表征了非均匀对宏,6.4 含时间的渐进展开（1）,其中，Eg表征了微观结构对宏观行为的影响。 Eg具有如下特性： 1）强依赖于单元尺寸； 2）均匀材料（=0 或 =1） Eg =0。 如果材料的非均匀尺度很小，则色散效应可以忽略。 若 ，界面没有反射，则波不发生色散。（物理角度）,21,-,6.4 含时间的渐进展开（1） 其中，Eg表征了微观结构对宏,6.4 含时间的渐进展开（2）,空间尺度,x,y=x/,宏观尺度：,微观尺度：,时间尺度,=2t, = t,慢尺度：,快尺度：,弹性动力学问题:,22,-,6.4 含时间的渐进展开（2）空间尺度 0432104,6.4 含时间的渐进展开（2）,弹性动力学问题:,参考文献：Fish, J. et. al. (2002). Non-local dispersive model for wave propagation in heterogeneous media: one-dimensional case. Int. J. Numer. Meth. Engng, 54, 331346.,x,y=x/,23,-,6.4 含时间的渐进展开（2）弹性动力学问题:参考文献：Fi,6.4 含时间的渐进展开（2）,令 i (i=-2,-1,0,1) 的系数为零，得到一系列控制方程：,24,-,6.4 含时间的渐进展开（2）令 i (i=-2,-1,6.4 含时间的渐进展开（2）,25,-,6.4 含时间的渐进展开（2）25-,本 章 结 束,高等复合材料力学,Advanced Mechanics of Composite Materials,-,本 章 结 束高等复合材料力学Advanced Mechan,高等复合材料力学,Advanced Mechanics of Composite Materials,第一章 绪论+张量基础 复合材料力学的三个重要特征、各向异性本构 张量的基本概念、爱因斯坦求和约定 符号ij与erst 坐标与坐标转换 张量的分量转换规律，张量方程 张量代数，商法则 常用特殊张量，主方向与主分量 张量函数及其微积分、高斯公式（散度定理）,-,高等复合材料力学Advanced Mechanics of,高等复合材料力学,Advanced Mechanics of Composite Materials,第二章 复合材料的有效性质和均质化方法 尺度和代表单元(RVE) 周期单元的选取 局部化和均匀化 等效柔度和等效刚度 Reuss近似和Voigt近似 等效性质的能量定义,-,高等复合材料力学Advanced Mechanics of,高等复合材料力学,Advanced Mechanics of Composite Materials,第三章 单层复合材料刚度的细观力学分析 长纤维复合材料刚度的细观力学分析 E2 、E2 、21和12 、G12 材料力学和弹性理论分析 周期性边界条件 最小势能原理和最小余能原理（Voigt上限和Reuss下限） 长纤维复合材料强度的细观力学分析 Xt 、Xc 、Yt 、Yc 、S 热膨胀的力学分析（ a1 、a2 ） 短纤维复合材料的细观力学分析 应力传递理论(剪滞法)、刚度、强度,-,高等复合材料力学Advanced Mechanics of,高等复合材料力学,Advanced Mechanics of Composite Materials,第四章 复合材料的单夹杂问题 Green函数 弹性力学的一般解 椭球型夹杂问题 Eshelby夹杂理论 夹杂的能量,-,高等复合材料力学Advanced Mechanics of,高等复合材料力学,Advanced Mechanics of Composite Materials,第五章 复合材料线性有效模量预测的近似方法 宏观整体坐标系和局部坐标系 稀疏方法 Mori-Tanaka方法 自洽方法 广义自洽方法 有效自洽方法 微分法 Hashin-Shtrikman界限 几种方法的比较 复合材料有效热膨胀系数,-,高等复合材料力学Advanced Mechanics of,高等复合材料力学,Advanced Mechanics of Composite Materials,第六章 均匀化理论和多尺度方法 多尺度模型 渐进展开法 空间多尺度问题（弹性静力学） 含时间的空间多尺度问题（弹性动力学） 空间和时间的双重多尺度问题（弹性动力学）,-,高等复合材料力学Advanced Mechanics of,谢 谢 ！,高等复合材料力学,Advanced Mechanics of Composite Materials,-,谢 谢 ！高等复合材料力学Advanced Mechan,