基于均匀化理论的多孔板弯曲问题新解

由于订型分布孔的存在,通常的有限元法不能有效地分析多孔板的弯曲问题。该文基于均匀化理论建立了该类问题的新解法。对含密集型分布的阶梯型圆孔板的分析结果,说明了该文方法是有效的。     

High-Velocity Laminar and Turbulent Flow in Porous Media

We model high-velocity flow in porous media with the multiple scale homogenization technique and basic fluid mechanics. Momentum and mechanical energy theorems are derived. In idealized porous media inviscid irrotational flow in the pores and wall boundary layers give a pressure loss with a power of 3/2 in average velocity. This model has support from flow in simple model media. In complex media the flow separates from the solid surface. Pressure loss effects of flow separation, wall and free shear layers, pressure drag, flow tube velocity and developing flow are discussed by using phenomenological arguments. We propose that the square pressure loss in the laminar Forchheimer equation is caused by development of strong localized dissipation zones around flow separation, that is, in the viscous boundary layer in triple decks. For turbulent flow, the resulting pressure loss due to average dissipation is a power 2 term in velocity.     

COMPUTATIONAL MULTISCALE METHODS FOR LINEAR HETEROGENEOUS POROELASTICITY

We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling.This poroelasticity problem suffers from rapidly oscillating material parameters,which calls for a thorough numerical treatment.In this paper,we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity.Therein,local corrector problems are constructed in line with the static equations,whereas we propose to consider the full system.This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure.We prove the optimal first-order convergence of this method and verify the result by numerical experiments.     

A boundary element method for homogenization of periodic structures

Homogenized coefficients of periodic structures are calculated via an auxiliary partial differential equation in the periodic cell. Typically, a volume finite element discretization is employed for the numerical solution. In this paper, we reformulate the problem as a boundary integral equation using Steklov–Poincaré operators. The resulting boundary element method only discretizes the boundary of the periodic cell and the interface between the materials within the cell. We prove that the homogenized coefficients converge super-linearly with the mesh size, and we support the theory with examples in two and three dimensions.     

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

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

An implicit-gradient-enhanced incremental-secant mean-field homogenization scheme for elasto-plastic composites with damage

This paper presents an incremental-secant mean-field homogenization (MFH) procedure for composites made of elasto-plastic constituents exhibiting damage. During the damaging process of one phase, the proposed method can account for the resulting unloading of the other phase, ensuring an accurate prediction of the scheme. When strain softening of materials is involved, classical finite element formulations lose solution uniqueness and face the strain localization problem. To avoid this issue the model is formulated in a so-called implicit gradient-enhanced approach, with a view toward macro-scale simulations. The method is then used to predict the behavior of composites whose matrix phases exhibit strain softening, and is shown to be accurate compared to unit cell simulations and experimental results. Then the convergence of the method upon strain softening, with respect to the mesh size, is demonstrated on a notched composite ply. Finally, applications consisting in a stacking plate, successively without and with a hole, are given as illustrations of the possibility of the method to be used in a multiscale framework.     

非均质Cosserat连续体细-宏观均匀化条件

基于平均场理论的多尺度模拟关键问题之一是给定恰当的表征元(RVE)边界条件,以使均匀化过程满足Hill-Mandel细宏观能量等价条件,也即Hill宏观均匀化条件。对于非均质Cosserat连续体,已有的研究工作只能得到合理的混合平动位移-偶应力表征元边界条件,常用的一致平动位移-转角以及周期边界条件等均不能使用,给计算均匀化算法推导和实施带来了困难,也阻碍了多尺度分析方法的进一步发展与应用。为此,本文在推导和建立一个新的Hill定理版本基础上,不仅成功地给定了多种强形式表征元边界条件,而且构造出了合理的弱形式周期边界条件,这些条件既满足细宏观能量等价也符合一阶平均场理论基本假定,可在均匀化方法中推广与应用。     

基于均匀材料微结构模型的热弹性结构与材料并发优化

研究由宏观上均匀多孔材料制成的结构的优化设计问题,待设计的结构受到给定的外力与温度载荷作用,优化设计旨在给定结构允许使用的材料体积约束下,设计宏观结构的拓扑及多孔材料的徼结构,使得结构柔度最小.本文提出了一种宏观结构与微观单胞构型并发优化设计的方法.在此方法中,引入宏观密度和微观密度两类设计变量,在微观层次上采用带惩罚的实心各向同性材料法SIMP(Solid Isotropic Material with Penalty),在宏观层次上采用带惩罚的多孔各向异性材料法PAMP(Porous Anisotropic Material with Pemlty),借助均匀化方法建立两个层次阃的联系,通过优化方法自动确定实体材料在结构与材料两个层次上的分配,得到优化设计;提供的数值算例检验了本文所提方法及计算模型,并讨论了温度变化、材料体积及计算参数对优化结果的影响.研究结果表明同时考虑热和机械载荷时,采用多孔材料可以降低结构柔顺性.     

PLASTIC LIMIT ANALYSIS OF DUCTILE COMPOSITE STRUCTURES FROM MICRO-TO MACRO-MECHANICAL ANALYSIS

The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （RVE） is selected to reflect the microstructures of the composite materials and the constituents are assumed to be elastic perfectly-plastic. Based on the homogenization theory and the static limit theorem, an optimization formulation to directly calculate the macroscopic strength domain of the RVE is obtained. The finite element modeling of the static limit analysis is formulated as a nonlinear mathematical programming and solved by the sequential quadratic programming method, where the temperature parameter method is used to construct the self-stress field. Secondly, Hill＇s yield criterion is adopted to connect the micromechanicai and macromechanical analyses. And the limit loads of composite structures are worked out on the macroscopic scale. Finally, some examples and comparisons are shown.     

Characterization of cuttlebone for a biomimetic design of cellular structures

Cuttlebone is a natural material possessing the multifunctional properties of high porosity, high flexural stiffness and compressive strength, making it a fine example of design optimization of cellular structures created by nature. Examination of cuttlebone using scanning electron micros- copy （SEM） reveals an approximately periodic microstruc- ture, appropriate for computational characterization using direct homogenization techniques. In this paper, volume fractions and stiffness tensors were determined based on two different unit cell models that were extracted from two different cuttlefish samples. These characterized results were then used as the target values in an inverse homogenization procedure aiming to re-generate microstructures with the same properties as cuttlebone. Unit cells with similar topologies to the original cuttlebone unit cells were achieved, attaining the same volume fraction （i.e. bulk density） and the same （or very close） stiffness tensor. In addition, a range of alternate unit cell topologies were achieved also attaining the target properties, revealing the non-unique nature of this inverse homogenization problem.