A novel implementation algorithm of asymptotic homogenization for predicting the effective coefficient of thermal expansion of periodic composite materials

Asymptotic homogenization (AH) is a general method for predicting the effective coefficient of thermal expansion (CTE) of periodic composites. It has a rigorous mathematical foundation and can give an accurate solution if the macrostructure is large enough to comprise an infinite number of unit cells. In this paper, a novel implementation algorithm of asymptotic homogenization (NIAH) is devel-oped to calculate the effective CTE of periodic composite materials. Compared with the previous implementation of AH, there are two obvious advantages. One is its implemen-tation as simple as representative volume element (RVE). The new algorithm can be executed easily using commercial finite element analysis (FEA) software as a black box. The detailed process of the new implementation of AH has been provided. The other is that NIAH can simultaneously use more than one element type to discretize a unit cell, which can save much computational cost in predicting the CTE of a complex structure. Several examples are carried out to demonstrate the effectiveness of the new implementation. This work is expected to greatly promote the widespread use of AH in predicting the CTE of periodic composite materials.     

On the implementation of symmetric and antisymmetric periodic boundary conditions for incompressible flow

In this paper we consider symmetric and antisymmetric periodic boundary conditions for flows governed by the incompressible Navier-Stokes equations. Classical periodic boundary conditions are studied as well as symmetric and antisymmetric periodic boundary conditions in which there is a pressure difference between inlet and outlet. The implementation of this type of boundary conditions in a finite element code using the penalty function formulation is treated and also the implementation in a finite volume code based on pressure correction. The methods are demonstrated by computation of a flow through a staggered tube bundle.     

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.     

Novel implementation of homogenization method to predict effective properties of periodic materials

Representative volume element (RVE) method and asymptotic homogenization (AH) method are two widely used methods in predicting effective properties of periodic materials. This paper develops a novel implementation of the AH method, which has rigorous mathematical foundation of the AH method, and also simplicity as the RVE method. This implementation can be easily realized using commercial software as a black box, and can use all kinds of elements available in commercial software to model unit cells with rather complicated microstructures, so the model may remain a fairly small scale. Several examples were carried out to demonstrate the simplicity and effectiveness of the new implementation.     

Modelling of functionally graded materials by numerical homogenization

Summary  In this contribution, the mechanical behaviour of different ZrO₂/NiCr 80 20 compositions is analysed and compared with experimental findings. The microwave-sintered material is found to possess a slightly dominant ceramic matrix for intermediate volume fractions. Its thermal expansion coefficient deviates from the rule of mixture. The modulus and the stress strain behaviour can be simulated by a numerical homogenization procedure, and the influence of residual stresses is found to be negligible. A newly introduced parameter (matricity) describes the mutual circumvention of the phases and is found to strongly control the stress level of the composite, globally as well as locally. Finally, a graded component and a metal/ceramic bi-material are compared for thermal as well as mechanical loading. Received 23 November 1999; accepted for publication 26 May 2000     

AN EXACT ELEMENT METHOD FOR THE BENDING OF NONHOMOGENEOUS REISSNER’S PLATE

In this paper,based on the step reduction method and exact analytic method,a new method,the exact element method for constructing finite element,is presented.Since the new method doesn’t need variational principle,it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficients.By this method,a triangle noncompatible element with15 degrees of freedom is derived to solve the bending of nonhomogenous Reissner’s plate.Because the displacement parameters at the nodal point only contain deflection and rotation angle.it is convenient to deal with arbitrary boundary conditions.In this paper,the convergence of displacement and stress resultants is proved.The element obtained by the present method can be used for thin and thick plates as well,Four numerical examples are given at the end of this paper,which indicates that we can obtain satisfactory results and have higher numerical precision.     

Zienkiewicz薄板单元的数值流形法实现

Zienkiewicz薄板单元是最早出现的板单元,因节点参数简单,所以在70年代之前曾起过重要作用;但同时它又是非协调的,仅对由三组间距相等的平行直线所形成的网格才能满足分片检验,这限制了这一单元的应用.因此,先后有不少学者对其进行过改进,使其能对任意形状的网格收敛.考虑到数值流形元法(NMM)总是用最佳质量的数学网格来进行逼近,本文尝试用NMM来解决Zienkiewicz单元的网格依赖问题.基于板的最小势能原理,提出了适用于NMM的混合变分提法,然后利用有解析解的椭圆板问题对方法进行了验证.本文所提出的方法,还可被用于其他类型的非协调单元的改造.     

A computational homogenization approach for the yield design of periodic thin plates. Part II : Upper bound yield design calculation of the homogenized structure

In the first part of this work (Bleyer and de Buhan, 2014), the determination of the macroscopic strength criterion of periodic thin plates has been addressed by means of the yield design homogenization theory and its associated numerical procedures. The present paper aims at using such numerically computed homogenized strength criteria in order to evaluate limit load estimates of global plate structures. The yield line method being a common kinematic approach for the yield design of plates, which enables to obtain upper bound estimates quite efficiently, it is first shown that its extension to the case of complex strength criteria as those calculated from the homogenization method, necessitates the computation of a function depending on one single parameter. A simple analytical example on a reinforced rectangular plate illustrates the simplicity of the method. The case of numerical yield line method being also rapidly mentioned, a more refined finite element-based upper bound approach is also proposed, taking dissipation through curvature as well as angular jumps into account. In this case, an approximation procedure is proposed to treat the curvature term, based upon an algorithm approximating the original macroscopic strength criterion by a convex hull of ellipsoids. Numerical examples are presented to assess the efficiency of the different methods.     

Stochastic boundary element method for reliability analysis of plate and beams composite structures

I.IntroductionTheengineeringstructuresareoftells,'1>:'rectcdtotheactionofthestochasticloadingthatvarieswiththetime,forexample,theengineeringstructuresactedonbytheearthquake,theoceanstructuresactedonbydynamicpressure,andthevehiclesofthetransportationinflue…     

An asymptotic iteration method for the numerical analysis of near-critical free-surface flows

Satisfying the boundary conditions at the free surface may impose severe difficulties to the computation of turbulent open-channel flows with finite-volume or finite-element methods, in particular, when the flow conditions are nearly critical. It is proposed to apply an iteration procedure that is based on an asymptotic expansion for large Reynolds numbers and Froude numbers close to the critical value 1.The iteration procedure starts by prescribing a first approximation for the free surface as it is obtained from solving an ODE that has been derived previously by means of an asymptotic expansion (Grillhofer and Schneider, 2003). The numerical solution of the full equations of motion then gives a surface pressure distribution that differs from the constant value required by the dynamic boundary condition. To determine a correction to the elevation of the free surface we next solve an ODE that is obtained from the asymptotic analysis of the flow with a prescribed pressure disturbance at the free surface. The full equations of motion are then solved for the corrected surface, and the procedure is repeated until criteria of accuracy for surface elevation and surface pressure, respectively, are satisfied.The method is applied to an undular hydraulic jump as a test case.