拉压性能不同材料全量型本构关系及厚壁筒的应力分析

将经典全量理论作了推广，考虑了应力状态及塑性体积变形对拉压性能不同材料的塑性行为的影响。应用该本构模型分别计算了厚壁筒在内压和外压作用下的应力分布。给出了径向应力、环向应力和轴向应力沿壁厚的分布图。将本文的计算解与拉压性能相同(不考虑体积变形、强化曲线唯一)的幂函数强化材料的厚壁筒的理论解进行了比较。结果表明，材料的拉压性能不同对厚壁筒的环向应力和轴向应力影响较大。因此，对于拉压性能不同材料，考虑到其对应力状态及塑性体积变形敏感时，是不能将其简化成拉压性能相同、体积不可压缩、强化曲线唯一的理想材料。     

粘弹性大变形动力响应的有限元分析

本文基于Ｔｏｔａｌ Ｌａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

粘弹性大变形动力响应的有限元分析

本文基于ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量叠加方法，采用Ｋｉｒｃｈｈｏｆｆ应力增量和Ｇｒｅｅｎ应变增量表示的动力虚功方程和Ｋｉｒｃｈｈｏｆｆ应力－Ｇｒｅｅｎ应变的单积分型本构关系，导出粘弹性大变形的动力变分方程。依此采用Ｎｅｗｍａｒｋ法和八节点轴对称等参数元与二十节点三维等参数元编制了轴对称及三维问题的动力响应计算程序，典型例题的计算结果表明分析符合结构的物理性质。     

基于Hellinger-Reissner变分原理的应变梯度杂交元设计

从一般的偶应力理论出发，基于Hellinger-Reissner变分原理，通过对有限元 离散体系的位移试解引入非协调位移函数，得到了偶应力理论下有限元离散系统的能量相容 条件，并由此建立了应变梯度杂交元的应力函数优化条件. 根据该优化条件，构造了一 个C0类的平面4节点梯度杂交元，数值结果表明，该单元对可压缩和不可压缩状态的 梯度材料均可给出合理的数值结果，再现材料的尺度效应.     

以Kirchhoff应力张量-Green应变张量表示本构关系的粘弹性大变形平面问题的有限元法

本文基于Updated Lagrangian增量叠加方法,采用以现时Kirchhoff应力增量和现时Green应变增量表示的虚功方程和Kirchhoff应力-Green应变的积分本构关系,导出粘弹性大变形的变分方程.依此采用八节点二维等参数元编制相应的平面问题的计算程序,算例的计算结果与以往的工作完全符合.     

基于有理插值方法模拟可压缩类橡胶材料多轴变形行为

提出一个多轴应变能函数来模拟可压缩类橡胶材料在载荷作用下的大变形行为。首先,基于对数应变构造三个不变量,使得应变能函数分别具备可压缩功能、多轴有效功能以及适用于多个变形模式的功能;其次,利用有理插值方法构造单轴形函数,并结合前面引入的三个不变量,通过哈密顿插值方法构建统一的应变能函数;最后,在单轴情况下,应变能函数通过对Hencky应变求导,可以得到对应变形模式下的应力-应变关系。本文通过调节泊松比的大小来控制体积变化从而实现模型的可压缩性;此应变能函数可以模拟单轴拉伸压缩、平面应变以及等双轴拉伸等至少三个基准试验。新方法得到的伸长比-应力关系、伸长比-体积比关系与经典试验数据结果作比较,进而证明方法的准确性和有效性。     

棒材通过锥形模静液挤压成形的无网格法分析

建立了一种基于初始构形及有限变形的粘塑性弹性本构关系，并由空间描述的Galerkin能量弱变分原理，经一致转换得一整体拉格朗日方程描述下的动量平衡方程．同时经线性化处理给出了显示中心差分法求解格式，以棒材通过锥形模的静液挤压成形为例进行了全面的EFG法数值模拟，从而证明了有限变形粘塑性EFG法对实际成形工艺分析、优化及设计的有效性。     

纯钼烧结体压缩变形本构模型建立

基于可压缩连续材料的屈服准则和塑性本构关系,推导了粉末烧结材料的压缩变形屈服条件,采用真空烧结技术制备了不同密度的钼,通过在MTS810实验机上进行单向压缩试验,研究了在不同的初始密度,变形温度和应变速率的条件下,纯钼烧结休压缩变形流变应力的变化规律.建立了其流变应力的经验回归模型,进而确定了纯钼烧结材料压缩变形的屈服...     

有限变形计算中的体积不可压缩问题

在金属材料塑性成形的数值模拟中 ,由于变形通常很大 ,有限变形理论被广泛采用。但由于有限变形下 ,塑性本构理论的研究相对欠缺 ,还未建立起较为公认的理论。塑性体积不可压缩是公认的结果 ,但在有限变形下如何描述体积不可压缩则必须为人们所关注。本文重点讨论了这一问题 ,认为必须完全抛弃以εii=0作为体积不可压缩条件 ,否则不仅将会给计算带来很大误差 ,而且还会导致计算不收敛     

刚塑性有限元体积可压缩法的简化形式

基于传统的刚塑性有限元体积可压缩法中一直忽视了等效应变速率ε的矩阵表达式可进一步简化，本文运用矩阵运算方法简化了等效应变速率的矩阵表达式。并结合刚塑性可压缩材料的变分原理，导出了刚塑性有限元体积可压缩法的简化计算格式     

Large deformation of incompressible rubber cylinder under plane strain

The large deformation of incompressible rubber cylinder under inner pressure is analyzed by a kind of new rubber materials strain energy function. The theory formulation for the displacement and stress is presented. The penalty finite element formulation is established in order to analyze nonlinear rubber materials, and the results of finite element method agree well with theoretical ones. A new method for controlling the calculating stability and convergence rates is put forward. The selection of the appropriate penalty factor and its influence on calculated results are discussed.     

Pressure Dependent Properties of a Compressible Polymer

This paper deals with the investigations of a porous carbon black-filled rubber, tested with regard to its pressure and tension behaviour. In the tension range only uniaxial tests are performed while in the pressure range uniaxial as well as hydrostatic tests are performed. The uniaxial experiments are carried out in a custom-made uniaxial device and the hydrostatic tests in a pressure chamber which is specially developed for this application. The construction and use of the pressure chamber is clearly described in this paper. All experiments are related to the basic elasticity of the material. The viscoelastic behaviour is completely disregarded at this point. Not only the experiments are discussed, also the modelling of the material is looked at. The tested cellular rubber is composed of an incompressible solid phase and a compressible gas phase. For that reason a so-called structural compressibility is observed. The compressible behaviour of cellular rubber is an important property. So the main focus of the paper is on the pressure tests and the simulation of these. The existing material models for rubber like materials only deal with incompressible rubber structures. To represent the compressible behaviour, the Theory of Porous Media is used. The constitutive model is based on a polynomial approach for an incompressible material. This is complemented by a volumetric expansion term with a point of compaction to model the structural compressibility.     

Determining Stress–Strain in Rubber with In-Plane Biaxial Tensile Tester

The nominal stress–strain relationships of industrial rubber materials under multiaxial deformation are essential for precisely determining the constitutive laws of those materials. This paper proposes a new method for precisely estimating the nominal stress–stretch relationships of carbon-black-filled styrene butadiene rubbers (SBRs) under uniaxial tension, pure shear, and equibiaxial tension by using an in-plane biaxial tensile tester. The proposed method employs sheet-shaped rubber samples with notches for the pure-shear and equibiaxial tension tests to mitigate the influence of non-uniform deformation around the clamps. Finite element analysis and biaxial tensile tests were performed to verify the effectiveness of the proposed method. Performance evaluations based on both numerical calculations and experiments revealed that the proposed method enabled the precise calculation of the nominal stress–stretch relationship for uniform deformation from a tensile load and deformation of the reference area defined at the center of the samples.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

一种适合橡胶类材料的非线性粘弹性本构模型

借助非线性流变模型建立大变形情况非线性粘弹材料的本构关系,考虑到大多数橡胶类材料具有的几乎不可压缩性,以及体积响应和剪切响应的流变性能不同,将变形梯度乘法分解为等容部分和体积变形部分,给出了一种适合橡胶类材料的非线性粘弹性本构模型,并模拟了粘滞效应。对于极快或极慢的过程,该模型退化为橡胶弹性理论;在小变形情况下退化为经典的广义Maxwell粘弹性材料。模型与热力学第二定律相容,适合于大规模数值分析。     

非连续变形计算力学模型的粘弹性分析方法Ⅰ:基本理论

基于粘弹性广义有限单元和接触力元，发展了适用于多体相互作用系统非连续变形分析的粘弹性数值分析方法，通过虚功原理，给出了其分区参变量最小势能原理，从而阐明了其理论基础。粘弹性广义有限单元的本构关系可由粘弹性退化为弹性或刚性，因此本文所提出的方法可对由刚体、弹性体和粘弹性体所构成的复杂多体系统在外荷载作用下的力学行为进行数值模拟，同时能够比本文精确地直接得到多体之间的接触应力。     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

平面应变不可压缩橡胶圆柱的大变形

针对一类新的橡胶材料应变能函数，推导了受内压作用下不可压缩橡胶圆柱的大变形公式，给出了位移、应力的解析表达式（用积分形式表示）。建立了适用于分析非线性不可压缩橡胶材料的罚有限元列式，算例表明：位移与应力能很好地与理论解吻合，并提出了控制计算稳定的方法，特别是详细地讨论罚因子的选取及对计算结果的影响。     

Fractional Derivative Viscoelasticity at Large Deformations

A time domain viscoelastic model for large three-dimensional responses underisothermal conditions is presented. Internal variables with fractional orderevolution equations are used to model the time dependent part of the response. By using fractional order rate laws, the characteristics of the timedependency of many polymeric materials can be described using relatively fewparameters. Moreover, here we take into account that polymeric materials are often used in applications where the small deformations approximation does nothold (e.g., suspensions, vibration isolators and rubber bushings). A numerical algorithm for the constitutive response is developed and implemented into a finite element code forstructural dynamics. The algorithm calculates the fractional derivatives by means of the Grünwald–Lubich approach.Analytical and numerical calculations of the constitutive response in the nonlinearregime are presented and compared. The dynamicstructural response of a viscoelastic bar as well as the quasi-static response of athick walled tube are computed, including both geometrically and materiallynonlinear effects. Moreover, it isshown that by applying relatively small load magnitudes, the responses ofthe linear viscoelastic model are recovered.     

Thermo-Chemo-Electro-Mechanical Formulation of Saturated Charged Porous Solids

A thermo-chemo-electro-mechanical formulation of quasi-static finite deformation of swelling incompressible porous media is derived from a mixture theory including the volume fraction concept. The model consists of an electrically charged porous solid saturated with an ionic solution. Incompressible deformation is assumed. The mixture as a whole is assumed locally electroneutral. Different constituents following different kinematic paths are defined: solid, fluid, anions, cations and neutral solutes. Balance laws are derived for each constituent and for the mixture as a whole. A Lagrangian form of the second law of thermodynamics for incompressible porous media is used to derive the constitutive restrictions of the medium. The material properties are shown to be contained in one strain energy function and a matrix of frictional tensors. A principle of reversibility results from the constitutive restrictions. Existing theories of swelling media should be evaluated with respect to this principle.