微极介质混合物的等效本构方程

对微极介质混合物引入代表性体积的等效均匀体,用代表性体积边界上的面力和面力偶定义等效应力和等效偶应力,提出了建立微极介质混合物的等效本构方程的一般原理和方法.讨论了以十字形框架为胞元的多胞材料面内变形问题的等效本构方程,给出合理的解析结果.     

一类多孔固体的等效偶应力动力学梁模型

一维多孔固体结构可采用等效连续介质梁模型来研究其动力学行为. 当类梁结构的高度尺寸和多孔固体单胞结构尺寸相近时,等效模型的力学行为会产生尺寸效应现象. 等效经典模型由于不包含尺度参数而无法描述尺寸相关特点,而广义连续介质力学模型则可以准确地考虑尺寸效应的影响. 基于偶应力理论,对一类单胞含有圆形孔洞的周期性多孔固体类梁结构,给出了分析其横向自由振动的等效连续介质铁木辛柯梁模型. 通过对单胞分析,在应变能等价和几何平均的意义下,定义了等效偶应力介质的材料常数. 利用已有的材料常数,推导了等效铁木辛柯梁的动力学微分方程. 将实际多孔固体结构进行完全的动力学有限元离散计算,所获得的解作为精确解以检验等效梁模型所获得的频率和振型的精度. 振型的比较借助于模态置信准则矩阵方法. 大量算例表明,等效偶应力铁木辛柯梁模型在频率和振型两方面均具有较高的计算精度. 重点研究了单胞孔径的相对大小、类梁结构高度与单胞尺寸比以及类梁结构长高比对等效梁模型精度的影响. 在此基础上,偏保守地建议了多孔固体类梁结构自振分析方法.      

类石墨烯二维原子晶体的微态理论模型

基于微态(Micromorphic) 连续介质理论,提出了针对类石墨烯二维原子晶体的新力学模型. 该模型以有限大小的布拉维单胞为基元体,考虑基元粒子的宏观位移和微观变形,依据微态理论基本方程,推导了全局坐标系下模型的主导方程. 然后针对布拉维单胞中含有两个原子的类石墨烯晶体,通过分析单胞中声子振动模式与基元体自由度的关系,获得了微态形式下声子色散关系的久期方程,并根据二维晶体声子色散特性对久期方程进行了简化,进而确定了类石墨烯晶体模型的本构方程. 最后,以石墨烯和单层六方氮化硼为例,利用简化的表达式拟合了它们面内声子色散关系数据,计算了模型材料的常数,石墨烯模型的等效杨氏模量、泊松比分别为1.05 TPa 和0.197,氮化硼分别为0.766 TPa 和0.225,均与已有的实验值相符合.      

点阵材料微极连续介质模型的应力优化设计

在将二维周期性点阵类材料等效为具有非局部化本构的微极连续介质的基础上，运用优化技术，探讨了基于材料相对密度和微单胞特征尺度两类变量的优化结构应力的方法，给出了针对最小化结构关键部位应力、结构最大应力最小化、最大化结构关键部位安全储备三类特定目标的结构与材料一体化协同优化结果．利用圆板小孔应力集中的数值算例验证了方法的有效性．     

基于应变梯度理论的微尺度蜂窝等效模量计算方法

本文提出了一种新的能够计及尺度效应的微纳米蜂窝等效模量的计算方法。将一种单参数应变梯度理论引入到本构方程当中,并基于能量等效原理推导了蜂窝面内等效模量地计算公式。算例分析表明,本文方法能够有效地计及尺度效应对蜂窝等效模量的影响。尺度效应与胞壁厚度和长度的值都有关,当胞壁厚度较小时,尺度效应显著,本文方法预测的模量会明显高于传统方法;而当胞壁厚度较大时,尺度效应变得微弱乃至可以忽略不计。但如果胞壁的长度/厚度比很大,则面内等效模量会趋近于0,此时是否考虑尺度效应意义不大。     

剪切和扭转工况下微纳米薄壁蜂窝的等效剪切模量计算

分别针对剪切和扭转两种工况给出了微纳米薄壁蜂窝等效剪切模量的解析计算方法.该方法综合考虑了由面板对芯层的约束导致的高度效应和当蜂窝胞壁厚度进入微纳米量级时引起的尺度效应.首先对蜂窝各胞壁选取了可反映面板约束以及受力状态的三角级数位移场,然后在本构关系中引入修正偶应力理论以描述尺度效应,最后应用能量均匀化方法求得蜂窝的等效剪切模量.以典型六边形蜂窝为例,给出了完整的计算过程和结果.与文献中的等效剪切模量结果进行对比,讨论了不同工况下等效剪切模量随芯层高度和胞壁厚度的变化趋势,以及高度效应和尺度效应之间的相互影响.     

两相混合流体的非线性本构理论

本文从连续介质力学的基本原理出发,建立了微极流体与经典流体两相流动的非线性扩散理论。给出了混合流体本构方程的一般形式。对单相流体、单相微极流体及稀悬浮体三种特殊情形,得到了具体形式的二阶非线性本构方程,并同已有的理论进行了比较。     

基于分子动力学模拟的多晶结构微观热-力耦合行为的计算

基于分子动力学模拟,建立了一套可用于表征微观下多晶结构热-力耦合行为的算法框架。该算法的要点是将连续模型和分子模拟耦合起来,并使守恒定律在微观连续模型和原子层次上都得到满足,与利用传统的连续介质力学建立晶界与晶粒的本构方程相比,本模型中的连续流是通过原子模型准确计算得到的,从而避免了使用经验的本构方程。     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.     

双向粘弹性节理岩体时域等效分析

以双向粘弹性节理岩体为研究对象,基于一种简单等效假定和自适应算法,并结合有限元方法,得到递推格式的等效本构模型和等效粘弹性场计算模型.通过算例,将等效解与基于ANSYS的非均质解进行比较.从计算精度和效率平衡的角度而言,等效模型的结果是令人满意的.     

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.     

Finite element study for conical indentation of elastoplastic micropolar material

Numerous experiments have repetitively shown that the material behavior presents effective size dependent mechanical properties at scales of microns or submicrons. In this paper, the size dependent behavior of micropolar theory under conical indentation is studied for different indentation depths and micropolar material parameters. To illustrate the effectiveness of the micropolar theory in predicting the indentation size effect (ISE), an axisymmetric finite element model has been developed for elastoplastic contact analysis of the micropolar materials based on the parametric virtual principle. It is shown that the micropolar parameters contribute to describe the characteristic of ISE at different scales, where the material length scale regulates the rate of hardness change at large indentation depth and the value of micropolar shear module restrains the upper limit of hardness at low indentation depth. The simulation results showed that the indentation loads increase as the result of increased material length scale at any indentation depth, however, the rate of increase is higher for lower indentation depth, relative to conventional continuum. The numerical results are presented for perfectly sharp and rounded tip conical indentations of magnesium oxide and compared with the experimental data for hardness coming from the open literature. It is shown that the satisfactory agreement between the experimental data and the numerical results is obtained, and the better correlation is achieved for the rounded tip indentation compared to the sharp indentation.     

金属材料的强度与应力-应变关系的球压入测试方法

压入法获取材料单轴应力-应变关系和抗拉强度对服役结构完整性评价有重要的基础意义.假定材料均匀连续、各向同性、应力应变关系符合Hollomon律,基于能量等效假定,即代表性体积单元(representativevolume element, RVE)的vonMises等效和有效变形域内能量中值等效假定,本文提出了关联材料载荷、深度、球压头直径和Hollomon律的四参数半解析球压入(semi-analyticalspherical indentation,SSI)模型.通过球压入载荷-深度试验关系获得材料的应力-应变关系和抗拉强度.考虑压入过程中的损伤效应,针对金属材料提出了用于球压入测试的材料弹性模量修正模型.对11种延性金属材料完成了球压入试验,采用本文提出的球压入试验方法测到的弹性模量、应力-应变关系和抗拉强度与单轴拉伸试验结果吻合良好.      

Macrohomogeneity condition in dynamics of micropolar media

We determine the macrohomogeneity (Hill-Mandel type) condition in the dynamic response of inhomogeneous micropolar (Cosserat) materials. The setting calls for small deformation gradients and curvatures, but without restrictions on the constitutive behavior and without any requirements of spatial periodicity. The condition gives admissible boundary loadings, along with extra terms representing kinetic energy contributions of both classical type and micropolar type. The said loadings involve various combinations of average stresses and strains, along with couple-stresses and curvature-torsion tensors. If applied to a specific microstructure in a computational mechanics approach, these boundary loadings will allow one to determine scale-dependent homogenization toward a representative volume element (RVE) of an equivalent homogeneous micropolar medium in either elastic or inelastic settings. By restricting the continuum model to an inhomogeneous Cauchy continuum and/or a quasi-static setting, the macrohomogeneity condition simplifies to conventional versions.     

Non-linear micropolar and classical continua for anisotropic discontinuous materials

Non-linear Cosserat and Cauchy anisotropic continua equivalent to masonry-like materials, like brick/block masonry, jointed rocks, granular materials or matrix/particle composites, are presented.An integral procedure of equivalence in terms of mechanical power has been adopted to identify the effective elastic moduli of the two continuous models starting from a Lagrangian system of interacting rigid elements. Non-linear constitutive functions for the interactions in the Lagrangian system are defined in order to take into account both the low capability to carry tension and the friction at the interfaces between elements. The non-linear problem is solved through a finite element procedure based on the iterative adjustment of the continuum constitutive tensor due to the occurrence of some limit situation involving the contact actions of the discrete model.Differences between the classical and the micropolar model are investigated with the aid of numerical analyses carried out on masonry walls made of blocks of different size. The capability of the micropolar continuum to discern, unlike the classical continuum, the behaviour of systems made of elements of different size is pointed out. It is also shown that for anisotropic materials, even in the elastic case, the micropolar solution in general does not tend to the classical solution when the size of the elements vanishes.     

A new micromechanical approach of micropolar continuum modeling for 2-D periodic cellular material

In this paper, we present a new united approach to formulate the equivalent micropolar constitutive relation of two-dimensional (2-D) periodic cellular material to capture its non-local properties and to explain the size effects in its structural analysis. The new united approach takes both the displacement compatibility and the equilibrium of forces and moments into consideration, where Taylor series expansion of the displacement and rotation fields and the extended aver-aging procedure with an explicit enforcement of equilibrium are adopted in the micromechanical analysis of a unit cell. In numerical examples, the effective micropolar constants obtained in this paper and others derived in the literature are used for the equivalent micropolar continuum simulation of cellular solids. The solutions from the equivalent analysis are compared with the discrete simulation solutions of the cellu-lar solids. It is found that the micropolar constants developed in this paper give satisfying results of equivalent analysis for the periodic cellular material.     

A continuum micromechanical theory of overall plasticity for particulate composites including particle size effect

Classical continuum micromechanics cannot predict the particle size dependence of the overall plasticity for composite materials, a simple analytical micromechanical method is proposed in this paper to investigate this size dependence. The matrix material is idealized as a micropolar continuum, an average equivalent inclusion method is advanced and the Mori–Tanaka's method is extended to a micropolar medium to evaluate the effective elastic modulus tensor. The overall plasticity of composites is predicted by a new secant moduli method based on the second order moment of strain and torsion of the matrix in a framework of micropolar theory. The computed results show that the size dependence is more pronounced when the particle's size approaches to the matrix characteristic length, and for large particle sizes, the prediction coincides with that predicted by classical micromechanical models. The method is analytical in nature, and it can capture the particle size dependence on the overall plastic behavior for particulate composites, and the prediction agrees well with the experimental results presented in literature. The proposed model can be considered as a natural extension of the widely used secant moduli method from a heterogeneous Cauchy medium to a micropolar composite.     

Modeling elasto-plastic indentation on layered materials using the equivalent inclusion method

This paper develops a fast semi-analytical model for solving the three-dimensional elasto-plastic contact problems involving layered materials using the Equivalent Inclusion Method (EIM). The analytical elastic solutions of a half-space subjected to a unit surface pressure and a unit subsurface eigenstrain are employed in this model; the topmost layer is simulated by an equivalent inclusion with fictitious eigenstrain. Accumulative plastic deformation is determined by a procedure involving an iterative plasticity loop and an incremental loading process. Algorithms of the fast Fourier transform (FFT) and the Conjugate Gradient Method (CGM) are utilized to improve the computation efficiency. An analytical elastic solution of layered body contact (O’Sullivan and King, 1988) and an indentation experiment result involving a layered substrate (Michler et al., 1999) are used to examine the accuracy of this model. Comparisons between numerical results from the present model and a commercial FEM software (Abaqus) are also presented. Case studies of a rigid ball loaded against a layered elasto-plastic half-space are conducted to explore the effects of the modulus, yield strength, and thickness of the coating on the hardness, stiffness, and plastic deformation of the composite body.     

Size effects in two-dimensional Voronoi foams: A comparison between generalized continua and discrete models

In view of size effects in cellular solids, we critically compare the analytical results of generalized continuum theories with the computational results of discrete models. Representatives are studied from two classes of generalized continuum theories: the strain divergence theory from the class of higher-grade continua and the micropolar theory from the class of higher-order continua. In the former, the divergence of strain is proposed as an additional deformation measure, while in the latter the microrotation gradients act as the source for extra internal energy. We analytically solve a range of basic boundary value problems (simple shear, pure bending and the strain concentration around a hole) and compare the results with discrete, numerical calculations that are based on a Voronoi representation of the cellular microstructure. By comparing both the local deformation fields and the overall elastic response, we critically assess the capabilities of both theories in capturing size effects in cellular solids.