Two-dimensional analyses of delamination buckling of symmetrically cross-ply rectangular laminates

The conventional approach to analysis the buckling of rectangular laminates containing an embedded delamination subjected to the in-plane loading is to simplify the laminate as a beam-plate from which the predicted buckling load decreases as the length of the laminate increases. Two-dimensional analyses are employed in this paper by extending the one-dimensional model to take into consideration of the influence of the delamination width on the buckling performance of the laminates. The laminate is simply supported containing a through width delamination. A new parameter β defined as the ratio of delamination length to delamination width is introduced with an emphasis on the influence of the delamination size. It is found that (i) when the ratio β is greater than one snap-through buckling prevails, the buckling load is determined by the delamination size and depth only; (ii) as the ratio β continues to increase, the buckling load will approach to a constant value. Solutions are verified with the well established results and are found in good agreement with the latter.     

Free-edge interlaminar stress analysis of piezo-bonded composite laminates under symmetric electric excitation

A stress function-based approach is proposed to analyze the free-edge interlaminar stresses of piezo-bonded symmetric laminates. The proposed method satisfies the traction free boundary conditions, as well as surface free conditions. The symmetric laminated structure was excited under electric fields that can generate induced strain, resulting in pure extension in the laminated plate. The governing equations were obtained by taking the principle of complementary virtual work. To verify the proposed method, cross-ply, angle-ply and quasi-isotropic laminates were analyzed. The stress concentrations predicted by the present method were compared with those analyzed by the finite element method. The results show that the stress function-based analysis of piezo-bonded laminated composite structures is an efficient and accurate method for the initial design stage of piezo-bonded composite structures.     

计及界面损伤的复合材料正交(混杂)叠层板的应力集中分析

对正交(混杂)叠层复合材料最终拉伸破坏过程中的细观应力集中问题,提出了一种修正的剪滞分析模型;研究了叠层中由于90°层的基体开裂、层间界面破坏、0°层中部分纤维断裂及纤维/基体界面损伤相互作用所导致的细观应力重新分布,获得了相应的应力集中因子和界面破坏区长度与界面剪切强度的定量关系。本文结果为进一步研究正交叠层复合材料的细观破坏机理、最终拉伸强度及协同效应等提供了重要的理论依据。     

Mean-field homogenization of elasto-viscoplastic composites based on a general incrementally affine linearization method

We propose a general formulation – which we believe to be new – for the mean-field homogenization of inclusion-reinforced elasto-viscoplastic composites assuming small strains. Our proposal is based on an interplay between constitutive equations and numerical algorithms, and the key ideas behind it are the following. The evolution equations for inelastic strain and internal variables at the beginning of each time interval are linearized around the ending time of the same interval. The linearized equations are then numerically integrated using a fully implicit backward Euler scheme. The obtained algebraic equations lead to an incrementally affine stress–strain relation which involves two important terms. The first one is the algorithmic tangent operator, obtained by consistent linearization of the time discretized constitutive equations. The second term is a new one and called an affine strain increment. The proposal leads to thermoelastic-like relations directly in the time domain, and not in the Laplace–Carson (L–C) one. There is no need for viscoelastic-type integral rewriting of the evolution equations, for L–C transforms, or for numerical inversion back from L–C to time domains. The proposed method can be readily applied to sophisticated elasto-viscoplastic models with an arbitrary set of scalar or tensor internal variables, and is valid for multi-axial, non-monotonic and non-proportional loading histories. The theory is applied in detail to a well-known constitutive model, and verified against finite element simulations of representative volume elements or unit cells, for a number of composite materials.     

Broadband acoustic cloaks based on the homogenization of layered materials

We present an analysis of acoustic cloaks based on the homogenization of two fluidlike materials, with an emphasis on periodically layered imperfect cloaks, by removing the singularities of the acoustic parameters required for ideal cloaks. The conditions that material layers should satisfy are systematically analyzed and critically discussed according to their feasibility.     

Material sensitivity analysis in homogenization of linear elastic composites

Summary  The main goal of the paper is to present theoretical aspects and the finite element method (FEM) implementation of the sensitivity analysis in homogenization of composite materials with linear elastic components, using effective modules approach. The deterministic sensitivity analysis of effective material properties is presented in a general form for an n-components periodic composite, and is illustrated by the examples of 1D as well as of 2D heterogeneous structures. The results of the sensitivity analysis presented in the paper confirm the usefulness of the homogenization method in computational analysis of composite materials the method may be applied to computational optimization of engineering composites, to the shape-sensitivity studies and, after some probabilistic extensions, to stochastic sensitivity analysis of random composites. Received 10 November 2000; accepted for publication 24 April 2001     

A homogenization analysis of the field theoretic approach to the quasi-continuum method

Using the orbital-free density functional theory as a model theory, we present an analysis of the field theoretic approach to quasi-continuum method. In particular, by perturbation method and multiple scale analysis, we provide a formal justification for the validity of numerical coarse-graining of various fields in the quasi-continuum reduction of field theories by taking the homogenization limit. Further, we derive the homogenized equations that govern the behavior of electronic fields in regions of smooth deformations. Using Fourier analysis, we determine the far-field solutions for these fields in the presence of local defects, and subsequently estimate cell-size effects in computed defect energies.     

Three-dimensional stochastic analysis using a perturbation-based homogenization method for elastic properties of composite material considering microscopic uncertainty

This paper discusses evaluation of influence of microscopic uncertainty on a homogenized macroscopic elastic property of an inhomogeneous material. In order to analyze the influence, the perturbation-based homogenization method is used. A higher order perturbation-based analysis method for investigating stochastic characteristics of a homogenized elastic tensor and an equivalent elastic property of a composite material is formulated.As a numerical example, macroscopic stochastic characteristics such as an expected value or variance, which is caused by microscopic uncertainty in material properties, of a homogenized elastic tensor and homogenized equivalent elastic property of unidirectional fiber reinforced plastic are investigated. The macroscopic stochastic variation caused by microscopic uncertainty in component materials such as Young’s modulus or Poisson’s ratio variation is evaluated using the perturbation-based homogenization method. The numerical results are compared with the results of the Monte-Carlo simulation, validity, effectiveness and a limitation of the perturbation-based homogenization method is investigated. With comparing the results using the first-order perturbation-based method, effectiveness of a higher order perturbation is also investigated.     

A consistent nonlocal scheme based on filters for the homogenization of heterogeneous linear materials with non-separated scales

In this work, the question of homogenizing linear elastic, heterogeneous materials with periodic microstructures in the case of non-separated scales is addressed. A framework if proposed, where the notion of mesoscopic strain and stress fields are defined by appropriate integral operators which act as low-pass filters on the fine scale fluctuations. The present theory extends the classical linear homogenization by substituting averaging operators by integral operators, and localization tensors by nonlocal operators involving appropriate Green functions. As a result, the obtained constitutive relationship at the mesoscale appears to be nonlocal. Compared to nonlocal elastic models introduced from a phenomenological point of view, the nonlocal behavior has been fully derived from the study of the microstructure. A discrete version of the theory is presented, where the mesoscopic strain field is approximated as a linear combination of basis functions. It allows computing the mesoscopic nonlocal operator by means of a finite number of transformation tensors, which can be computed numerically on the unit cell.     

基于均匀化方法的单向纤维增强体渗透率预报

针对具有周期性分布细观结构的纤维增强体，从Stokes方程出发，用均匀化理论建立了预报纤维预制体渗透率的数学模型. 将Stokes方程与线弹性力学中的Lame方程进行类比，给出了用线弹性平面应变问题的有限元分析程序求解Stokes方程的方法. 据该方法编写了FORTRAN程序HAPS求解控制方程，并以此预报单向纤维增强体渗透率，与有关文献的结果进行比较证明了该方法的合理性.     

含孔复合材料层合板静拉伸三维逐渐损伤分析

针对面内静拉伸纤维增强复合材料含中孔层合板，发展了参数化三维逐渐损伤模型. 该模型 可以模拟含中孔层合板损伤起始、发展及最终结构破坏整个过程，并能较好地预测含中孔层 合板的破坏模式和破坏强度. 采用所发展的模型和有限元三维逐渐损伤分析技术即应力分 析、失效判定准则及损伤过程中材料性能退化等，对其他文献所提供的9种不同类型含中孔层合板进行了损伤扩展分析及强度预测，同时对层合板的损伤基本机理、类型及其相互关联作用进行了探讨，计算结果与文献实验结果非常吻合.     

Microscopic symmetric bifurcation condition of cellular solids based on a homogenization theory of finite deformation

In this paper, we establish a homogenization framework to analyze the microscopic symmetric bifurcation buckling of cellular solids subjected to macroscopically uniform compression. To this end, describing the principle of virtual work for infinite periodic materials in the updated Lagrangian form, we build a homogenization theory of finite deformation, which satisfies the principle of material objectivity. Then, we state a postulate that at the onset of microscopic symmetric bifurcation, microscopic velocity becomes spontaneous, yet changing the sign of such spontaneous velocity has no influence on the variation in macroscopic states. By applying this postulate to the homogenization theory, we derive the conditions to be satisfied at the onset of microscopic symmetric bifurcation. The resulting conditions are verified by analyzing numerically the in-plane biaxial buckling of an elastic hexagonal honeycomb. It is thus shown that three kinds of experimentally observed buckling modes of honeycombs i.e., uniaxial, biaxial and flower-like modes, are attained and classified as microscopic symmetric bifurcation. It is also shown that the multiplicity of bifurcation gives rise to the complex cell-patterns in the biaxial and flower-like modes.     

3-D numerical study on the bending of symmetric composite laminates

３－ＤＮＵＭＥＲＩＣＡＬＳＴＵＤＹＯＮＴＨＥＢＥＮＤＩＮＧＯＦＳＹＭＭＥＴＲＩＣＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＳＣｈｉｅｎＷｅｉ－ｚａｎｇ（钱伟长）;ＨｕａｎｇＱｉａｎ（黄默）;ＦｅｎｇＷｅｉ（冯伟）（ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｔｙｏｆＴｅｃ...     

Homogenization of microsystem interconnects based on micropolar theory and discontinuous kinematics

In the present paper, the homogenized mechanical response of an interface in a microsystem interconnection is established on the basis of micropolar theory. The interface is treated as a finite RVE (representative volume element), across which macroscopic discontinuities occur as expressed in terms of the regularized discontinuous displacement and rotation fields. For the microstructure within the interfacial RVE, the micro-macro kinematical coupling is introduced as a second-order Taylor series expansion, along with a fluctuation term representing the microscopic displacement variation. In the second-order term of the expansion a restriction for the curvature is made, which motivates the adopted micropolar kinematics. Explicit expressions for the homogenized traction vector and the couple stress traction, associated with displacement and rotational jumps across the interface surface, are derived. A planar elastic interface is subjected to three basic deformation modes, i.e. the standard modes I, II and a non-conventional rotation mode, which are considered in the numerical examples representing a typical interconnect. A comparison to the results from the Taylor assumption is made, which shows that the Taylor assumption method produces an overstiffening of the interface.     

基于Bézier提取的三维等几何分析

等几何分析(IGA)将非均匀有理B样条(NURBS)函数作为有限元形函数,具有几何精确、高阶连续和精度高等优点。与常规有限元法C0连续的形函数不同,高阶IGA基函数不是定义在一个单元上,而是跨越由几个单元组成的参数空间,因而编程复杂且无法嵌入现有的有限元法计算框架及相应算法。本文建立了基于Bézier提取的三维IGA,将NURBS函数分解成伯恩斯坦多项式的线性组合,从而实现把NURBS单元分解为C0连续Bézier单元,这些单元与Lagrange单元相似,使IGA的实现和常规有限元一样,以便将IGA分析嵌入现有的有限元软件中。两个三维算例结果表明,基于Bézier提取的IGA和传统IGA的收敛性和精度相同。     

电子封装元件的高阶层离散均匀化分析

将均匀化理论应用于具有非完全(单层内)周期性微结构的倒装焊底充胶电子封装元件,建立了高阶逐层离散层板模型,用解析法分析热载荷下结构的温度应力. 计算结果与有限元解的比较表明,该分析模型和方法是有效的,而且比较简便. 算例分析结果显示,胶层厚度、焊点密度、胶与焊点材料的模量比和体积比,对于焊点温度应力有明显影响.     

Empirical test of a microscopic three-phase traffic theory

A review of dynamic nonlinear features of spatiotemporal congested patterns in freeway traffic is presented. The basis of the review is a comparison of theoretical features of the congested patterns that are shown by a microscopic traffic flow model in the context of the Kerner's three-phase traffic theory and empirical microscopic and macroscopic pattern characteristics measured on different freeways over various days and years. In this test of the microscopic three-phase traffic flow theory, a model of an "open" road is applied: Empirical time-dependence of traffic demand and drivers' destinations are used at the upstream model boundaries. At downstream model boundary conditions for vehicle freely leaving a modeling freeway section(s) are given. Spatiotemporal congested patterns emerge, develop, and dissolve in this open freeway model with the same types of bottlenecks as those in empirical observations. It is found that microscopic three-phase traffic models can explain all known macroscopic and microscopic empirical congested pattern features (e.g., probabilistic breakdown phenomenon as a first-order phase transition from free flow to synchronized flow, moving jam emergence in synchronized flow rather than in free flow, spatiotemporal features of synchronized flow and general congested patterns at freeway bottlenecks, intensification of downstream congestion due to upstream congestion at adjacent bottlenecks). It turns out that microscopic optimal velocity (OV) functions and time headway distributions are not necessarily qualitatively different, even if local congested traffic behavior is qualitatively different. Model performance with respect to spatiotemporal pattern emergence and evolution cannot be tested using these traffic characteristics. The reason for this is that important spatiotemporal features of congested traffic patterns are lost in these and many other macroscopic and microscopic traffic characteristics, which are widely used as the empirical basis for a test of traffic flow models. PACS: 89.40. + k, 47.54. + r, 64.60.Cn, 64.60.Lx     

带攻角平头弹侵彻不同厚度芳纶层合板的数值模拟

为研究攻角对不同厚度芳纶层合板抗平头弹侵彻性能的影响,构建了三维有限元计算模型,首先通过对比实验结果验证了其可靠性,然后基于该数值模型,进一步计算了0°～30°攻角范围内,4、8和16 mm靶板的弹道响应,从子弹剩余速度、靶板能量吸收率、极限弹道速度与穿孔能量阈值4个方面,综合评估了芳纶层合板的抗侵彻性能。结果表明：攻角的影响与靶板厚度及子弹入射速度有关,随着攻角的增大,靶板的极限速度和穿孔能量阈值均有所降低,降低的程度随厚度的增加而减小;入射速度接近芳纶层合板弹道极限速度时,子弹剩余速度随着攻角增大而增大,但速度远高于弹道极限速度时,子弹剩余速度随着攻角增大而减小;攻角对芳纶层合板弹道性能的影响机理随靶板的破坏模式不同而改变。

     

基于降阶均匀化理论的混凝土双尺度力学性能分析

混凝土材料细观特性对宏观力学性能有着重要影响。为进一步分析混凝土细观特性对宏观力学行为的影响规律,将混凝土材料简化为由骨料、砂浆和界面三相组成,编制了随机凸多面体骨料生成、投放和网格剖分算法,建立可用于有限元计算的满足级配要求的随机细观模型。针对直接使用细观力学模型计算量较大的问题,采用降阶均匀化理论,对混凝土细观胞元模型进行预处理并编制了相应的双尺度计算程序。对不同强度混凝土进行了单轴静态压缩双尺度计算,与实验数据和细观力学模拟结果符合较好。研究表明,降阶均匀化理论在加快求解速度的同时具有较高的精度,可以用于混凝土的多尺度力学性能分析。     

Homogenization of inelastic solid materials at finite strains based on incremental minimization principles. Application to the texture analysis of polycrystals

The paper presents new continuous and discrete variational formulations for the homogenization analysis of inelastic solid materials undergoing finite strains. The point of departure is a general internal variable formulation that determines the inelastic response of the constituents of a typical micro-structure as a generalized standard medium in terms of an energy storage and a dissipation function. Consistent with this type of finite inelasticity we develop a new incremental variational formulation of the local constitutive response, where a quasi-hyperelastic micro-stress potential is obtained from a local minimization problem with respect to the internal variables. It is shown that this local minimization problem determines the internal state of the material for finite increments of time. We specify the local variational formulation for a distinct setting of multi-surface inelasticity and develop a numerical solution technique based on a time discretization of the internal variables. The existence of the quasi-hyperelastic stress potential allows the extension of homogenization approaches of finite elasticity to the incremental setting of finite inelasticity. Focussing on macro-deformation-driven micro-structures, we develop a new incremental variational formulation of the global homogenization problem for generalized standard materials at finite strains, where a quasi-hyperelastic macro-stress potential is obtained from a global minimization problem with respect to the fine-scale displacement fluctuation field. It is shown that this global minimization problem determines the state of the micro-structure for finite increments of time. We consider three different settings of the global variational problem for prescribed displacements, non-trivial periodic displacements and prescribed stresses on the boundary of the micro-structure and develop numerical solution methods based on a spatial discretization of the fine-scale displacement fluctuation field. Representative applications of the proposed minimization principles are demonstrated for a constitutive model of crystal plasticity and the homogenization problem of texture analysis in polycrystalline aggregates.