纤维排列方式对复合材料总体粘弹性常数的影响

对于金属基或高分子聚合物基复合材料，在特定情况下会表现出明显的粘弹性特性。本文采用Riemann—Liouville形式的分数阶导数模型描述基体的粘性特性，通过渐进均匀化方法给出了预测纤维加强复合材料整体本构关系的解析表达式，给出应用于基体具有Makris粘弹性关系的具体形式。最后，考察了圆截面纤维正方形排列和对角排列时的总体粘弹性弹性常数随纤维比的变化曲线。结果表明，这类复合材料仍具有粘弹性特性，其整体粘弹性本构关系的弹性部分综合了纤维弹性和基体弹性的贡献，粘性部分来自基体粘性的贡献，复合材料具有和基体相同的粘性系数和分数阶。为分析微结构特征对整体特性的贡献，须求解两类局部问题。在相同纤维体积比情况下，正方形排列的总体弹性系数大于正方形对角排列，而粘性常数相反。  

Novel numerical implementation of asymptotic homogenization method for periodic plate structures

The present paper develops and implements finite element formulation for the asymptotic homogenization theory for periodic composite plate and shell structures, earlier developed in  and , and thus adopts this analytical method for the analysis of periodic inhomogeneous plates and shells with more complicated periodic microstructures. It provides a benchmark test platform for evaluating various methods such as representative volume approaches to calculate effective properties. Furthermore, the new numerical implementation (Cheng et al., 2013) of asymptotic homogenization method of 2D and 3D materials with periodic microstructure is shown to be directly applicable to predict effective properties of periodic plates without any complicated mathematical derivation. The new numerical implementation is based on the rigorous mathematical foundation of the asymptotic homogenization method, and also simplicity similar to the representative volume method. It can be applied easily using commercial software as a black box. Different kinds of elements and modeling techniques available in commercial software can be used to discretize the unit cell. Several numerical examples are given to demonstrate the validity of the proposed methods.  

Effective elastic shear stiffness of a periodic fibrous composite with non-uniform imperfect contact between the matrix and the fibers

In this contribution, effective elastic moduli are obtained by means of the asymptotic homogenization method, for oblique two-phase fibrous periodic composites with non-uniform imperfect contact conditions at the interface. This work is an extension of previous reported results, where only the perfect contact for elastic or piezoelectric composites under imperfect spring model was considered. The constituents of the composites exhibit transversely isotropic properties. A doubly periodic parallelogram array of cylindrical inclusions under longitudinal shear is considered. The behavior of the shear elastic coefficient for different geometry arrays related to the angle of the cell is studied. As validation of the present method, some numerical examples and comparisons with theoretical results verified that the present model is efficient for the analysis of composites with presence of imperfect interface and parallelogram cell. The effect of the non uniform imperfection on the shear effective property is observed. The present method can provide benchmark results for other numerical and approximate methods.  

Homogenization of acoustic waves in strongly heterogeneous porous structures

We consider acoustic waves in fluid-saturated periodic media with dual porosity. At the mesoscopic level, the fluid motion is governed by the Darcy flow model extended by inertia terms and by the mass conservation equation. In this study, assuming the porous skeleton is rigid, the aim is to distinguish the effects of the strong heterogeneity in the permeability coefficients. Using the asymptotic homogenization method we derive macroscopic equations and obtain the dispersion relationship for harmonic waves. The double porosity gives rise to an extra homogenized coefficient of dynamic compressibility which is not obtained in the upscaled single porosity model. Both the single and double porosity models are compared using an example illustrating wave propagation in layered media.  

Estimation of very narrow bounds to the behavior of nonlinear incompressible elastic composites

Variational bounds for the effective behavior of nonlinear composites are improved by incorporating more-detailed morphological information. Such bounds, which are obtained from the generalized Hashin–Shtrikman variational principles, make use of a reference material with the same microstructure as the nonlinear composite. The geometrical information is contained in the effective properties of the reference material, which are explicitly present in the analytical formulae of the nonlinear bounds. In this paper, the variational approach is combined with estimates for the effective properties of the reference composite via the asymptotic homogenization method (AHM), and applied to a hexagonally periodic fiber-reinforced incompressible nonlinear elastic composite, significantly improving some recent results.  

Material forces in microfractured bodies

In the present work, we study the overall behavior of a microfractured elastic body within the configurational mechanics framework. Micro and macro scales are considered and scale changes are carried out by asymptotic developments homogenization. The homogenized equations of material momentum and scalar moment of material momentum are obtained. In these equations the microcrack length appears as an internal variable, describing the damage evolution. Both quasistatic and dynamic formulations are presented.  

MULTI-SCALE METHOD OF PLAIN WOVEN COMPOSITES SUBJECTED TO LOW VELOCITY IMPACT BASED ON ASYMPTOTIC HOMOGENIZATION1)

针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明：宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.  

STUDY ON NONLINEAR ELASTIC HOMOGENIZATION WITH ITERATIVE METHOD

微观结构对复合材料的宏观力学性能具有至关重要的影响, 通过合理设计复合材料微观结构可以得到期望的宏观性能. 均质化方法作为一种有效的设计方法, 它从微观结构的角度出发, 利用均匀化的概念, 实现了对复合材料宏观力学性能的预测和设计. 而当考虑非线性因素, 均质化的实现就非常困难. 本文利用双渐近展开方法, 将位移按照宏观位移和微观位移展开, 推导了非线性弹性均质化方程. 通过直接迭代法, 对非线性弹性均质化方程进行了求解, 并给出了具体的迭代方法和实现步骤. 本文基于迭代步骤和非线性弹性均质化方程编写MATLAB 程序, 对3种典型本构关系的周期性多孔材料平面问题进行了计算, 对比细致模型的应变能、最大位移和等效泊松比, 对程序及迭代方法的准确性进行了验证. 之后对一种三元橡胶基复合材料进行多尺度均质化, 将其分为芯丝尺度和层间尺度. 用线弹性的均质化方法得到了芯丝尺度的等效弹性参数, 并将其作为层间尺度的材料参数. 在层间尺度应用非线性弹性均质化方法对结构进行计算, 得到材料的宏观等效性能, 并以实验结果为基准进行评价.  

一维周期性梁结构等效性能计算方法讨论

具有轴向周期微结构的复合梁结构,通常在宏观上简化为一维欧拉-伯努利梁。由于缺乏基于严格数学理论、同时考虑降维及均匀化的等效性能计算方法,已有研究或采用基于平截面假定的弯曲能量近似方法,或采用基于三维周期性介质等效性质的方法。本文首先介绍了基于一维周期性梁的渐近均匀化理论求解新方法,并在此基础上与上述两种方法进行比较。结果表明,基于平截面假定的近似方法忽视了这类梁结构内的三维应力状态,过高地估计了梁的等效性质。  

纤维加强复合材料的总体粘弹性响应

在某些纤维增强复合材料(FRC)中使用金属或高分子聚合物作为基体材料。在高温等情况下，这类材料具有明显的粘弹性特性。本文采用Riemann—Liouville形式的分数阶导数模型描述基体的粘弹性特性。通过渐近均匀化方法给出了预测FRC整体三维本构关系的解析表达式。给出了应用于基体具有Makris粘弹性关系的具体形式。以圆截面纤维正方形排列的情形为例，给出了等效模量随纤维体积比的变化曲线。结果说明，这类复合材料仍具有粘弹性特性，其整体粘弹性本构关系的弹性部分综合了纤维弹性和基体弹性的贡献，粘性部分来自基体粘性的贡献，复合材料具有和基体相同的粘性系数和分数阶。为分析微结构特征对整体特性的贡献，须求解两类局部问题。可以看出，在整体的等效模量中包含了局部变形的贡献，局部变形增加了复合材料的耦合刚度。