周期性蜂窝材料微极等效模型及其微观应力的一种快速算法

将周期性蜂窝材料等效为具有非局部本构的微极连续介质,以解释实验中出现的尺度效应和边界层效应.在评论相关的多种不同方法(能量法、体积平均的均匀化法等)之后,提出了一种基于位移连续和单胞力平衡的推导微极等效本构参数的新方法.以正方形单胞制成的结构为例,在不同的结构与单胞尺寸比下,考虑承受集中点载荷、均布轴力和均布剪力三种载荷工况,比较了离散完全计算、经典连续介质等效和不同微极连续体等效本构的计算结果,建议了较好的微极本构参数值.数值模拟表明,集中点载荷和剪切载荷作用时,在加载点附近和边界部分,微极等效可以显著提高计算精度.最后,给出了一种映射算法,可以根据微极等效连续体分析的结果,快速计算出对应微观单胞构件的应力,以开有圆孔的方板应力集中为例,验证并考察了所提快速算法的有效性和计算精度.     

Generalized continuum modeling of 2-D periodic cellular solids

A generalized continuum representation of two-dimensional periodic cellular solids is obtained by treating these materials as micropolar continua. Linear elastic micropolar constants are obtained using an energy approach for square, equilateral triangular, mixed triangle and diamond cell topologies. The constants are obtained by equating two different continuous approximations of the strain energy function. Furthermore, the effects of shear deformation of the cell walls on the micropolar elastic constants are also discussed. A continuum micropolar finite element approach is developed for numerical simulations of the cell structures. The solutions from the continuum representation are compared with the “exact” discrete simulations of these cell structures for a model problem of elastic indentation of a rectangular domain by a point force. The utility of the micropolar continuum representation is illustrated by comparing various cell structures with respect to the stress concentration factor at the root of a circular notch.     

Vibration analysis of periodic cellular solids based on an effective couple-stress continuum model

Cellular solids are usually treated as homogeneous continuums with effective properties. Nevertheless, these mechanical properties depend strongly on the ratio of the specimen size to the cell size. These size effects may be accounted for according to preliminary static analysis of effective continuums based on couple-stress theory. In this paper an effective dynamic continuum model, based on couple-stress theory, is proposed to analyze the behavior of free vibrations of periodic cellular solids. In this continuum model, the effective mechanical constants of the effective continuum are deduced by an equivalent energy method. The cellular solid structure is then replaced with the equivalent couple-stress continuum with same overall dimension and shape. Moreover, the finite element formulation of the couple-stress continuum for the generalized eigenvalue analysis is developed to implement the free vibration analysis. The eigenfrequencies of the effective continuum are then obtained via the shear beam theory or the finite element method. A conventional finite element analysis by discretizing each cell of the cellular solids is also carried out to serve as an exact solution. Several structural cases are calculated to demonstrate the accuracy and effectiveness of the proposed continuum model. Good agreement on structural eigenfrequencies between the effective continuum solutions and the exact solutions shows that the proposed continuum model can accurately simulate the dynamic behavior of the cellular solids.     

Effective couple-stress continuum model of cellular solids and size effects analysis

The purpose of this paper is to develop a homogeneous, orthotropic couple-stress continuum model to take the place of the periodic heterogeneous cellular solids. Through generalizing the definition of the characteristic length for isotropic couple-stress continuum, four characteristic lengths are introduced as material engineering constants for such kind of continuum. In order to determine the effective moduli and the characteristic lengths of the effective couple-stress continuum, a Representative Volume Element (RVE) method is constructed. The effective properties are obtained based on the response of the RVE under prescribed boundary conditions, and our results agree with the analytical solutions in literature. In addition, the influences of the relative density, the topology, the size, and the properties of the solid material of cellular materials on the effective moduli as well as the characteristic lengths are discussed, respectively. Furthermore, the size effects in cellular solid beams are investigated using our effective couple-stress continuum model. The results show that the developed continuum model in this paper can precisely capture the size effects in cellular solids.     

Continuum modeling of periodic brickwork

Continuum models of periodic masonry brickwork, viewed at a micro-level as a discrete system, are identified within the frame of linearized elasticity. The accuracy of various identification schemes is investigated for standard and micropolar continua, which are directly compared with the help of some numerical benchmarks, for different loading conditions that induce periodic and non-periodic deformation states. It is shown that periodic deformation states of brickwork are exactly reproduced by both continua, provided that a suitable identification scheme is adopted. For non-periodic states micropolar continuum is shown to better reproduce the discrete solutions, due to its capability to take scale effects into account. Both continua are asymptotically equivalent as the characteristic length of the discrete system tends to zero, while providing an upper and a lower bound of the discrete solution.     

An analytical procedure to determine constitutive coefficients of a mixture of two linear elastic solids

In this paper, we study the coefficients of constitutive equations of a binary mixture of elastic solids and give an analytical approach to determine them. Assuming a material of two-phase elastic composite with randomly distributed elastic spheres is equivalent to a mixture of two elastic solids, we find the values of unknown coefficients by making use of Boussinesq problem. Furthermore, a mean displacement vector definition is also given.     

基于多层模拟方法的双层梁断裂有限元分析

本文提出一种用于含分层的双层梁线弹性断裂分析的有限元方法．将上下子梁均模拟为多个子层,采用只有平动位移自由度的新型梁单元,假设单元内的位移沿纵向和横向均线性变化,推导了该单元的单元刚度矩阵．将开裂部分和未开裂部分的子梁进行单元刚度矩阵组装,施加相应的等效结点力,得到整体平衡方程,并结合边界条件进行求解．为验证该方法的有效性和精度,开展非对称双悬臂梁(Asymmetric Double Cantilever Beam, ADCB)和单臂弯曲梁(Single Leg Bending, SLB)试样的断裂分析,利用虚拟裂纹闭合技术(Virtual Crack Closure Technique, VCCT)得到了试样的能量释放率及其分量,并将求得的结果与解析解和二维有限元解进行对比．计算结果表明,相对于传统双层模拟方法,该多层模拟方法能够精确、高效地计算各类梁试样的能量释放率及其分量,并且无需引入界面连续条件．     

Electroelastic waves in dielectrics modeled as polarizable continua

On the basis of a dielectric microcontinuum model, we investigate the problem of bulk wave propagation in a dielectric crystal with hexagonal material symmetry. The present linear micropolar model allows to express electric polarization via mechanical macro and micro-strain measures so that the coupling between acoustic and polarization modes can be described in terms of intrinsic dipole and quadrupole densities. The governing differential systems for different coupled modes are equivalent to some previous results of the classical phenomenological approach to ferroelectrics but also hold for piezoelectric solids with null intrinsic polarization. Resonance couplings between polaritons and acoustic waves arise from the dispersion equations depending on suitable relations among the micropolar constitutive parameters. Exploiting the dynamical representation of polarization for the admitted modes, we obtain piezoelectric coefficients and electromechanical coupling factors as functions of the wavelength (or frequency). As an application, a numerical example is given for the hexagonal phase of zinc sulfide.     

Spatial Estimates for the Constrained Anisotropic Elastic Cylinder

This paper establishes spatial estimates in a prismatic (semi-infinite) cylinder occupied by an anisotropic homogeneous linear elastic material, whose elasticity tensor is strongly elliptic. The cylinder is maintained in equilibrium under zero body force, zero displacement on the lateral boundary and pointwise specified displacement over the base. The other plane end is subject to zero displacement (when the cylinder is finite, say). The limiting case of a semi-infinite cylinder is also considered and zero displacement on the remote end (at large distance) is not assumed in this case. A first approach is developed by considering two mean-square cross-sectional measures of the displacement vector whose spatial evolution with respect to the axial variable is studied by means of a technique based on a second-order differential inequality. Conditions on the elastic constants are derived that show the cross-sectional measures exhibit alternative behaviour and in particular for the semi-infinite cylinder that there is either at least exponential growth or at most exponential decay. A second approach considers cross-sectional integrals involving the displacement and its gradient and furnishes information upon the spatial evolution, without restricting the range of strongly elliptic elastic constants. Such models are principally based upon a first-order differential inequality as well as on one of second order. The general results are explicitly presented for transversely isotropic materials and graphically illustrated for a cortical bone.     

Stability, bifurcation and chaos of a delayed oscillator with negative damping and delayed feedback control

This paper presents a detailed analysis on the dynamics of a delayed oscillator with negative damping and delayed feedback control. Firstly, a linear stability analysis for the trivial equilibrium is given. Then, the direction of Hopf bifurcation and stability of periodic solutions bifurcating from trivial equilibrium are determined by using the normal form theory and center manifold theorem. It shows that with properly chosen delay and gain in the delayed feedback path, this controlled delayed system may have stable equilibrium, or periodic solutions, or quasi-periodic solutions, or coexisting stable solutions. In addition, the controlled system may exhibit period-doubling bifurcation which eventually leads to chaos. Finally, some new interesting phenomena, such as the coexistence of periodic orbits and chaotic attractors, have been observed. The results indicate that delayed feedback control can make systems with state delay produce more complicated dynamics.     

复合材料扭转轴截面微结构拓扑优化设计

提出复合材料扭转轴截面微结构拓扑优化设计新模型，模型的优化目标是获得具有最大宏观剪切特性加权和的单胞形式．通过模型和均匀化方法及优化技术可以获得优化的微结构单胞，进而改善或者得到最优宏观弹性特性的复合材料．为了便于制造和应用，胞体材料用来获得复合材料的极值剪切模量．最后的优化结果表明，该模型连同数值处理技巧可以非常有效地实现微结构的拓扑优化设计．     

Constitutive theories for nonlocal micropolar continua with implicity and with multiple interactions

In this paper the constitutive theory for nonlocal micropolar continua which was proposed by A. C. Eringen is extended to the cases for nonlocal micropolar continua with implicity and with multiple interactions. Here nonlocal micropolar thermoelastic solids with implicity and with multiple interactions are cited as instances to illustrate the procedure for the establishment of their constitutive theories as well as two relevant theorems concerning the constitutive theories for those solids are given.     

Formulation of the high-fidelity generalized method of cells with arbitrary cell geometry for refined micromechanics and damage in composites

A new parametric formulation for high-fidelity generalized method of cells (HFGMC) is presented for the micromechanical analysis of multiphase periodic composites. To this end, a linear parametric and geometric mapping is employed to transform arbitrary quadrilateral cell shapes from the physical space to an auxiliary uniform square shapes. A complete quadratic displacement expansion is performed in the mapped space. Thus, a new bilinear term is added to the quadratic displacement expansion; unlike the original HFGMC for regular array of rectangular cells where this term in not required. The continuity of displacements, tractions, together with the periodicity and equilibrium conditions are imposed in the average sense, similar to the original HFGMC formulation, using both the physical and mapping variables. However, the addition of bilinear terms requires the introduction of the first averaged moments of the equilibrium equations. In order to demonstrate the ability the new HFGMC formulation, spatial stress fields are compared with analytical and numerical solutions of circular and elliptical fibers in an infinite medium. Furthermore, two progressive damage methodologies are coupled with the new HFGMC formulation in order to predict the strain softening and elastic degrading behaviors. The first methodology employs a cell extinction approach, while the second uses cohesive interfaces between the cells. Due to the strain softening, both damage methodologies require an iterative solution approach of the governing system nonlinear equations. Damage applications are presented for the transverse loading of composites with square and hexagonal repeating unit-cells (RUC).     

RESTUDY OF COUPLED FIELD THEORIES FOR MICROPOLAR CONTINUA (Ⅰ)──MICROPOLAR THERMOELASTICITY

ＩｎｔｒｏｄｕｃｔｉｏｎＳｏｍｅａｕｔｈｏｒｓｓｔｕｄｉｅｄｔｈｅｃｏｕｐｌｅｄｆｉｅｌｄｐｒｏｂｌｅｍｓｆｏｒｍｉｃｒｏｐｏｌａｒｃｏｎｔｉｎｕａ .Ｅｓｐｅｃｉａｌｌｙ ,Ｗ .Ｎｏｗａｃｋｉｐｕｂｌｉｓｈｅｄａｓｅｒｉｅｓｏｆａｂｏｕｔ 4 0ｓｃｉｅｎｔｉｆｉｃｐａｐｅｒｓｄｅａｌｉｎｇｗｉｔｈｔｈｅｍｉｃｒｏｐｏｌａｒｔｈｅｒｍｏｅｌａｓｔｉｃｉｔｙａｓｗｅｌｌａｓｔｈｅｐｒｏｂｌｅｍｓｏｆｄｉｓｔｏｒｔｉｏｎ ,ｔｈｅｒｍｏｄｉｆｆｕｓｉｏｎ ,ｔｈｅｒｍｏｐｉｅｚｏｅｌｅｃｔｒｉｃｉｔｙａｎｄｍ…     

Finite Anticlastic Bending of Hyperelastic Solids and Beams

This paper deals with the equilibrium problem in nonlinear elasticity of hyperelastic solids under anticlastic bending. A three-dimensional kinematic model, where the longitudinal bending is accompanied by the transversal deformation of cross sections, is formulated. Following a semi-inverse approach, the displacement field prescribed by the above kinematic model contains three unknown parameters. A Lagrangian analysis is performed and the compressible Mooney-Rivlin law is assumed for the stored energy function. Once evaluated the Piola-Kirchhoff stresses, the free parameters of the kinematic model are determined by using the equilibrium equations and the boundary conditions. An Eulerian analysis is then accomplished to evaluating stretches and stresses in the deformed configuration. Cauchy stress distributions are investigated and it is shown how, for wide ranges of constitutive parameters, the obtained solution is quite accurate. The whole formulation proposed for the finite anticlastic bending of hyperelastic solids is linearized by introducing the hypothesis of smallness of the displacement and strain fields. With this linearization procedure, the classical solution for the infinitesimal bending of beams is fully recovered.     

Micropolar modeling of planar orthotropic rectangular chiral lattices

Rectangular chiral lattices possess a two-fold symmetry; in order to characterize the overall behavior of such lattices, a two-dimensional orthotropic chiral micropolar theory is proposed. Eight additional material constants are necessary to represent the anisotropy in comparison with triangular ones, four of which are devoted to chirality. Homogenization procedures are also developed for the chiral lattice with rigid or deformable circles, all material constants in the developed micropolar theory are derived analytically for the case of the rigid circles and numerically for the case of the deformable circles. The dependences of these material constants and of wave propagation on the microstructural parameters are also examined.     

Invariant characteristic representations for classical and micropolar anisotropic elasticity tensors

By extending and developing the characteristic notion of the classical linear elasticity initiated by Lord Kelvin, a new type of representation for classical and micropolar anisotropic elasticity tensors is introduced. The new representation provides general expressions for characteristic forms of the two kinds of elasticity tensors under the material symmetry restriction and has many properties of physical and mathematical significance. For all types of material symmetries of interest, such new representations are constructed explicitly in terms of certain invariant constants and unit vectors in directions of material symmetry axes and hence they furnish invariants which can completely characterize the classical and micropolar linear elasticities. The results given are shown to be useful. In the case of classical elasticity, the spectral properties disclosed by our results are consistent with those given by similar established results.     

Chiral effect in plane isotropic micropolar elasticity and its application to chiral lattices

In continuum mechanics, the non-centrosymmetric micropolar theory is usually used to capture the chirality inherent in materials. However, when reduced to a two dimensional (2D) isotropic problem, the resulting model becomes non-chiral. Therefore, influence of the chiral effect cannot be properly characterized by existing theories for 2D chiral solids. To circumvent this difficulty, based on reinterpretation of isotropic tensors in the 2D case, we propose a continuum theory to model the chiral effect for 2D isotropic chiral solids. A single material parameter related to chirality is introduced to characterize the coupling between the bulk deformation and the internal rotation, which is a fundamental feature of 2D chiral solids. Coherently, the proposed continuum theory is applied for the homogenization of a triangular chiral lattice, from which the effective material constants of the lattice are analytically determined. The unique behavior in the chiral lattice is demonstrated through the analyses of a static tension problem and a plane wave propagation problem. The results, which cannot be predicted by the non-chiral model, are verified by the exact solution of the discrete model.