Dissipation consistent fabric tensor definition from DEM to continuum for granular media

In elastoplastic soil models aimed at capturing the impact of fabric anisotropy, a necessary ingredient is a measure of anisotropic fabric in the form of an evolving tensor. While it is possible to formulate such a fabric tensor based on indirect phenomenological observations at the continuum level, it is more effective and insightful to have the tensor defined first based on direct particle level microstructural observations and subsequently deduce a corresponding continuum definition. A practical means able to provide such observations, at least in the context of fabric evolution mechanisms, is the discrete element method (DEM). Some DEM defined fabric tensors such as the one based on the statistics of interparticle contact normals have already gained widespread acceptance as a quantitative measure of fabric anisotropy among researchers of granular material behavior. On the other hand, a fabric tensor in continuum elastoplastic modeling has been treated as a tensor-valued internal variable whose evolution must be properly linked to physical dissipation. Accordingly, the adaptation of a DEM fabric tensor definition to a continuum constitutive modeling theory must be thermodynamically consistent in regards to dissipation mechanisms. The present paper addresses this issue in detail, brings up possible pitfalls if such consistency is violated and proposes remedies and guidelines for such adaptation within a recently developed Anisotropic Critical State Theory (ACST) for granular materials.     

厚板的高阶剪切变形理论研究

已有多种厚板理论和高阶剪切变形模型,但仍需要进一步研究以更加完善.首先根据平均转角及上下表面剪应力自由这两个条件,提出了具有统一高阶剪切变形模型的中面位移模式,并将之表示为正交分解形式.根据正交特性,定义了板的广义应力;运用板问题应变能密度表示的等价性,提出了与广义应力功共轭的广义应变表示形式,建立了板的本构关系.证明了不同转角定义时虚功原理板理论表示的客观性,以及与三维弹性理论表示的等价性.运用虚功原理,建立了变分自洽的高阶厚板理论和变分渐近的低阶厚板理论,推导了相应的平衡方程及边界条件,分析了与已有板理论的异同.以广义应力形式建立了厚板理论的平衡方程,厘清了不同转角表示时板理论间的关系、低阶厚板理论与高阶厚板理论间的关系以及剪切系数计算等若干基本问题.对圣维南扭转问题的求解证明了该理论的正确性.      

裂隙岩体损伤局部化破坏分岔模型及其应用

采用概率统计方法分析节理裂隙岩体的几何特征，定义了反映裂隙岩体几何特征的组构张量．基于不可逆热力学理论，通过裂纹扩展的细观分析，得出了损伤的发展机理和演化方程，把损伤演化和裂隙的几何特征的变化联系起来，建立了弹塑性损伤本构关系．为分析含有节理裂隙岩体在发生局部化破坏时的特征，通过对发生局部化时的裂隙岩体的分析。构造了适用于节理裂隙岩体局部化分析的不连续分岔模型．利用非线性规划数值解法，可以得出局部化破坏的方向特征．在有限元方法中，根据该模型给出了节理裂隙岩体相关的算例，分析表明该模型用于分析裂隙岩体的局部化破坏是有效的．     

An improved dynamic model for a silicone material beam with large deformation

Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate formulation (ANCF) beam element employs cubic interpolation in the longitudinal direction and linear interpolation in the transverse direction, whereas it cannot accurately describe the large bending deformation. On this account, a novel modeling method for studying the dynamic behavior of nonlinear materials is proposed in this paper. In this formulation, a higher-order beam element characterized by quadratic interpolation in the transverse directions is used in this investigation. Based on the Yeoh model and volumetric energy penalty function, the nonlinear elastic force matrices are derived within the ANCF framework. The feasibility and availability of the Yeoh model are verified through static experiment of nonlinear incompressible materials. Furthermore, dynamic simulation of a silicone cantilever beam under the gravity force is implemented to validate the superiority of the higher-order beam element. The simulation results obtained based on the Yeoh model by employing three different ANCF beam elements are compared with the result achieved from a commercial finite element package as the reference result. It is found that the results acquired utilizing a higher-order beam element are in good agreement with the reference results, while the results obtained using a lower-order beam element are different from the reference results. In addition, the stiffening problem caused by volumetric locking can be resolved effectively by applying a higher-order beam element. It is concluded that the proposed higher-order beam element formulation has satisfying accuracy in simulating dynamic motion process of the silicone beam.     

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

利用松质骨理想化模型对骨小梁刚度的研究

利用结构张量的概念建立了考虑松质骨结构各向异性的理想模型，将杆型与板型两种传统模型统一起来．采用均匀化理论，通过有限单元的数值计算建立了松质骨弹性常数与结构张量、固体体积比之间的数值关系，并预测了骨小梁的弹性模量．结果表明，本文模型较好地体现了松质骨的各向异性力学性质．对骨小梁的计算结果表明，骨小梁的弹性模量在６７６ＧＰａ～１０９ＧＰａ范围内，平均意义上的结果是９２１ＧＰａ．     

极分解的一类近似计算与比较

基于级数展开给出了极分解中右伸长张量 的级数表示,通过对级数的项的选取得到右伸长张量的不同近似表达式。针对不同级数展开表示,得到表达式最小误差的级数展开形式。进而结合一些简单实例,验证误差了近似公式的有效性。最后与黄模佳等关于计算右伸长张量 和转动张量 的近似表达式进行了比较,本文的级数展开方式得到的右伸长张量 和转动张量 的近似表达式不但简洁,而且计算精度更高、适用范围更广。     

Reducing the Order of the Boundary Condition in a Theory of Plane Strain Elasticity with Boundary Reinforcement

We consider an elastic solid whose boundary is partly reinforced in such a way that the coated part of the boundary can be represented by the union of a finite number of open curves. We show how the (nonstandard) boundary condition on the reinforcement can be integrated resulting in an equivalent (lower-order) boundary condition which, unlike its higher-order counterpart, permits analysis of the corresponding boundary value problems using the boundary integral equation method. Mathematics Subject Classifications (2000) 74B05, 45E05, 35J25     

基于Oda裂隙结构张量的流纹岩各向异性试验研究

岩石内天然存在长度、倾角和形态不同的裂隙,造成岩石的各向异性特征。为揭示岩石内天然随机裂隙发育特征对岩石物理力学特性的影响规律,以泥巴山隧址区采集裂隙性流纹岩为研究对象,首先对试样裂隙进行素描统计分析;然后基于Oda裂隙结构张量,获得天然随机分布裂隙的几何统计参数;最后对裂隙性流纹岩试样分别进行单轴和常规三轴压缩试验,得到不同应力路径下流纹岩的应力-应变曲线及物理力学参数。分析Oda裂隙结构张量定义的各向异性参数与试验获得的力学参数之间的规律,研究结果表明:(1)Oda裂隙结构张量适用于天然随机分布裂隙的几何统计分析,各向异性参数A(F)越大,裂隙优势方向越明显;(2)单轴压缩下,随着各向异性参数I1和A(F)的增大,流纹岩各向异性程度增大,弹性模量减小,泊松比增大;(3)常规三轴压缩下,流纹岩弹性模量和泊松比随各向异性参数改变的规律较不明显,Oda裂隙结构张量不再适用。     

Some results for approximate strain and rotation tensor formulations in geometrically non-linear Reissner-Mindlin plate theory

Finite element deflection and stress results are presented for four flat plate configurations and are computed using kinematically approximate (rotation tensor, strain tensor or both) non-linear Reissner-Mindlin plate models. The finite element model is based on a mixed variational principle and has both displacement and force field variables. High order interpolation of the field variables is possible through p-type discretization. Results for some of the higher order approximate models are given for what appears to be the first time. It is found that for the class of example problems examined, exact strain tensor but approximate rotation tensor theories can significantly improve the solution over approximate strain tensor models such as the von Kármán and moderate rotation models when moderate deflections/rotations are present. However, for each of the problems examined (with the exception of a postbuckling problem) the von Kármán and moderate rotation model results compared favorably with the higher order models for deflection magnitudes which could be reasonably expected in typical aeroelastic configurations.     

Simulation of fabric drape using a thin plate element with finite rotation

The draping behavior of fabric is simulated by using four node quadrilateral thin plate elements with finite rotation. The finite element formulation is based on the total Lagrangian approach. An exact representation of finite rotation is introduced. The strain energy function accounting for the material symmetry is obtained by the tensor representation theory. To avoid shear locking, the assumed strain technique for transverse shear is adopted. The conjugate gradient method with a proposed line search algorithm is employed to minimize energy and reach the final shape of fabric. The draping behavior of a rectangular piece of fabric over a rectangular table is simulated.     

INCREMENTAL MICRO-MECHANICAL MODEL OF PLAIN WOVEN FABRIC

I. INTRODUCTION Previous research on woven fabric ignored the micro-weaving structures in fabric and modeled thewoven fabric as a ?exible orthotropic plate[1??3]. However, such a model is unable to predict certainbuckling phenomena of fabric commonly ob…     

基于微面有效应力矢量的各向异性屈服准则

基于微面模型，定义损伤变量为微面上有效承载面积的减少. 将Kachanov的一维有效 应力概念推广到三维，提出微面有效应力矢量的概念. 根据微面的有效应力矢量，将无损材 料的宏观应力张量及不变量与微面应力矢量的积分关系拓展到有损材料，得到了有损材料的 宏观有效应力张量及其不变量与宏观名义应力张量、微面面积损伤组构张量之间的关系. 将 无损材料的以应力张量不变量表示的Drucker-Prager准则推广到有损材料，建立了含缺陷 材料的各向异性屈服准则. 对有损材料，宏观有效应力张量与Murakami的有效应力张量具 有相同的形式，各向异性强度准则与Liu等提出的扩展Hill准则有相同的形式，当不考虑 静水应力对屈服的影响时，它与Hill准则具有相同的形式.     

Macro–micro relations in granular mechanics

In granular mechanics, macroscopic approaches treat a granular material as an equivalent continuum at macro-scale, and study its constitutive relationship between macro-quantities, such as stresses and strains. On the other hand, microscopic approaches consider a granular material as an assembly of individual particles interacting with each other at micro-scale (i.e., particle-scale), and the physical quantities under study are forces and displacements. This paper aims at linking up the findings from these two scales and to establish the macro–micro relations in granular mechanics.Three aspects of the macro–micro relations are investigated. They are about the internal structure, the stress tensor and the strain tensor. The internal structure is described with geometrical systems at the particle scale. Micro-structural definitions of the stress and strain tensors are derived, which link the macro-stress tensor with the contact forces, and the macro-strain tensor with the relative displacements at contact. In addition to a brief review of the past research work on these topics, further generalizations are made in this paper. In particular, the two cell systems proposed by Li and Li (2009), namely the solid cell system and the void cell system, are introduced and used for the derivation of the macro-structural expressions. The stress expression is derived based on Newton’s second law of motion. The result is valid for both static and dynamic cases. The strain expression is derived based on the compatibility requirement. And the expression is valid for any tessellation subdividing the granular assembly into polyhedral elements.The homogenization for deriving a macroscopic constitutive relationship from microscopic behaviour is discussed. Attention is placed on the macroscopic quantification of the internal structure in terms of a second rank tensor, known as the fabric tensor. Existing definitions of the fabric tensors have been reviewed. The correlations among different fabric tensors and their relations with the stress–strain behaviours have been investigated.     

Gradient elasticity in statics and dynamics: An overview of formulations,length scale identification procedures,finite element implementations and new results

In this paper, we discuss various formats of gradient elasticity and their performance in static and dynamic applications. Gradient elasticity theories provide extensions of the classical equations of elasticity with additional higher-order spatial derivatives of strains, stresses and/or accelerations. We focus on the versatile class of gradient elasticity theories whereby the higher-order terms are the Laplacian of the corresponding lower-order terms. One of the challenges of formulating gradient elasticity theories is to keep the number of additional constitutive parameters to a minimum. We start with discussing the general Mindlin theory, that in its most general form has 903 constitutive elastic parameters but which were reduced by Mindlin to three independent material length scales. Further simplifications are often possible. In particular, the Aifantis theory has only one additional parameter in statics and opens up a whole new field of analytical and numerical solution procedures. We also address how this can be extended to dynamics. An overview of length scale identification and quantification procedures is given. Finite element implementations of the most commonly used versions of gradient elasticity are discussed together with the variationally consistent boundary conditions. Details are provided for particular formats of gradient elasticity that can be implemented with simple, linear finite element shape functions. New numerical results show the removal of singularities in statics and dynamics, as well as the size-dependent mechanical response predicted by gradient elasticity.     

Evolution of fabric tensors in damage mechanics of solids with micro-cracks: Part I – Theory and fundamental concepts

The evolution of fabric tensors based on micro-crack distributions is formulated based on sound thermodynamic principles. In Part I of this work, the exact definition of fabric tensors based on micro-crack distributions is presented. This definition is seen to incorporate both the orientation and length of a micro-crack. In this regard, the micro-crack distribution is assumed to be radially symmetric, i.e. symmetric about a line through the origin.The equations of thermodynamics are employed in order to derive the exact evolution equations of the fabric tensors defined in the first part. In this regard, a thermodynamic force that is associated with the fabric tensor is defined and utilized in the derivation of the evolution equations. The application of the theory to the case of uniaxial tension is derived in Part II (companion paper) of this work.     

基于静动态混合重构的DG/FV混合格式

通过比较紧致格式和间断Galerkin(DG)格式, 提出了``静态重构'和``动态重构'的概念,对有限体积方法和DG有限元方法进行统一的表述. 借鉴有限体积的思想, 发展了基于``混合重构'技术的一类新的DG格式, 称之为间断Galerkin有限元/有限体积混合格式(DG/FV格式). 该类混合格式通过适当地扩展模板(拓展至紧邻单元)重构单元内的高阶多项式分布, 在提高精度的同时, 减少了传统DG格式的计算量和存储量. 通过典型一维和二维标量方程的计算发现新的混合格式在有些情况下具有超收敛(superconvergence)性质.      

Elastic Rods with Moderate Roration

Based on the theory of a directed or Cosserat curve of Green, Naghdi and several of their co-workers, several properly invariant approximate theories of an elastic rod are developed. The approximate theories are applicable when a rotation tensor associated with the motion of the rod is moderate. In one of these theories, some of the strains are moderate while others are small. A second theory where all of the strains are small is also developed. This revised version was published online in August 2006 with corrections to the Cover Date.