生长对超弹性球壳变形和稳定性的影响

应用连续介质力学有限变形理论建立受内压作用不可压超弹性球壳大变形问题的力学模型, 且运用基于变形梯度张量极分解的弹性体积生长理论分析生长对不可压超弹性球壳变形和稳定性的影响. 通过对球壳变形与内压关系式的数值计算得到不同生长条件下球壳的变形曲线和应力分布曲线及由生长引起的残余变形和残余应力分布. 计算结果表明生长对球壳变形特性有明显的影响, 生长影响球壳可产生不稳定变形的临界壁厚和临界内压, 且在某些情况下可改变球壳的稳定性. 生长在球壳中可产生一定的残余变形和残余应力, 对球壳中的应力分布有一定的影响. 另外当生长的程度足够大时, 即便没有外力作用,球壳仅在生长引起的残余应力作用下也可产生不稳定变形.      

双层球壳结构准静态线性化学弹性问题的位移势函数解

本文研究了各向同性固体的化学-力学耦合问题,在传统化学弹性理论描述的扩散-变形耦合关系基础上,进一步考虑了化学反应与固体变形的相互作用关系,发展了等温状态下固体-扩散-反应-变形耦合的线性化学弹性理论,拓展了化学弹性力学的应用范围．该理论能够同时描述固体内介质扩散和固体与介质之间化学反应两个不同时间尺度的化学过程,并给出由此引起的弹性范围内的应变和应力．为应用该模型求解具体化学弹性问题,本文通过构造扩散-反应位移势函数来获得位移特解形式,再与齐次Lamé方程通解叠加获得完整解;针对反应控制问题,引入化学弹性准静态假设,将反应-扩散-变形全耦合的瞬态过程分解为两个可解耦的相继过程,从而获得相应位移解．基于此解法,本文获得了反应控制的双层球壳结构化学弹性问题的解析解,并分析了化学反应、几何结构和弹性模量对应力分布的影响．     

一种近场动力学非普通状态理论零能模式控制方法

近场动力学非普通状态理论在采用节点积分时将引起零能模式,造成位移场、应力应变场的数值不稳定性,影响计算精度甚至会导致完全错误的结果,因此必须对其进行控制.目前国际上还没有十分有效的零能模式控制方法.本文针对零能模式问题,提出了一种通用的、高效的控制方法.根据近场动力学线性键理论,确定非均匀变形对应弹性张量的具体形式,考虑了微模量随不同作用键的变化.通过最小位能原理推导出非均匀变形引起的力状态,结合近场动力学力状态,得到稳定力状态表达式.从而建立起基于线性键理论的稳定关联材料模型,并应用于含圆孔平板、三点弯试件线弹性变形和损伤破坏过程模拟.数值结果表明,本文模型能有效抑制近场动力学非普通状态理论中的零能模式现象.与已有零能模式控制方法相比,其物理意义明确,不包含控制参数,避免了复杂的零能模式参数调节过程,提高了计算效率.      

考虑缠结效应的超弹性本构模型

经典熵弹性模型, 如 Neo-Hookean模型和Arruda-Boyce八链模型, 被广泛应用于预测橡胶等软材料的超弹性力学行为. 然而, 大量实验结果也显示仅采用一套模型参数, 这类模型不能同时准确地描述橡胶在多种加载模式下的应力响应. 为了克服上述模型的不足, 本文在熵弹性的模型基础上引入缠结约束效应. 微观上, 采用Langevin统计模型来表征熵弹性变形自由能, 通过管模型(tube model)引入缠结约束自由能, 并基于仿射假设, 建立微观变形与宏观变形之间的映射关系. 在宏观上, 所建立的超弹性模型的Helmholtz自由能同时包含熵弹性和缠结约束两部分, 其中熵弹性自由能与经典的Arruda-Boyce八链模型一致, 依赖于柯西-格林应变张量的第一不变量, 而缠结约束自由能依赖于柯西-格林应变张量的第二不变量. 与文献中的实验结果对比发现, 该三参数模型能准确地预测实验中所测得的橡胶材料在单轴拉伸、纯剪切和等双轴拉伸变形条件下的应力响应, 也能较好地描述不同预拉伸比条件下双轴拉伸实验结果. 最后, 本文比较了所建立的基于应变不变量的缠结约束模型与文献中相关的缠结约束模型在多种加载模式下自由能的异同. 总的来说, 本文所建立的本构理论能准确模拟橡胶等软材料的大变形力学行为, 对其工程应用有促进作用.      

欧拉描述的大变形固结理论

以往大变形固结理论主要基于一般的固体力学模型，其控制方程忽视了固结过程中排水引起 的质量变化. 提出饱和土的连续介质模型，并基于连续介质力学的公理体系推导了反映质量 变化的欧拉描述的大变形固结控制方程. 发现传统固结理论中：(1)忽视了渗流速度对土体平衡条件的影响；(2)决定土体平衡的总应力张量只有在土体变形速度和渗流速度方向相同时才具有对称性等. 在忽略变质量效应等条件下，传统理论成为本文理论的特例. 通过算例 的有限元分析，比较了欧拉描述与两种物质描述方法的差别，得到初步结论：(1)欧拉描述 方法计算的地基沉降量要小于物质描述方法的结果；(2)欧拉描述方法计算的侧向位移偏大 于两种物质描述结果.     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

一种近场动力学非普通状态理论零能模式控制方法

近场动力学非普通状态理论在采用节点积分时将引起零能模式,造成位移场、应力应变场的数值不稳定性,影响计算精度甚至会导致完全错误的结果,因此必须对其进行控制.目前国际上还没有十分有效的零能模式控制方法.本文针对零能模式问题,提出了一种通用的、高效的控制方法.根据近场动力学线性键理论,确定非均匀变形对应弹性张量的具体形式,考虑了微模量随不同作用键的变化.通过最小位能原理推导出非均匀变形引起的力状态,结合近场动力学力状态,得到稳定力状态表达式.从而建立起基于线性键理论的稳定关联材料模型,并应用于含圆孔平板、三点弯试件线弹性变形和损伤破坏过程模拟.数值结果表明,本文模型能有效抑制近场动力学非普通状态理论中的零能模式现象.与已有零能模式控制方法相比,其物理意义明确,不包含控制参数,避免了复杂的零能模式参数调节过程,提高了计算效率.     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性–塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.     

不同应力分量下广义开尔文模型粘性系数探讨

对不同应力分量下的广义开尔文模型应力应变关系进行了研究,推导了在不同应力分量下的广义开尔文模型的粘性应变增量计算式;通过对这些粘性应变增量计算式的比较分析,得到结论:对于线性粘弹性模型,当应力张量引起粘性变形的规律与应力偏量和球应力分别引起粘性变形的规律相同时,它们的系数满足关系式Ek/ηk=Gsk/ηsk=Kmk/ηmk;否则,这个关系式不成立.现有文献采用应力张量表示的粘性变形有限元计算式隐含假定了球应力与应力偏量产生的粘性变形规律相同.对于复杂的工程材料而言,这种假定并不总是合适的.这在工程问题粘性分析时值得注意.     

非局部平面应变和平面应力问题界定及其精确性讨论

在非局部弹性理论框架下对平面应变和平面应力状态重新界定.首先,分别在其相应简化假设下推导控制方程,并与经典局部情况进行比较.然后,引入变形协调条件对两类非局部平面问题的精确性进行讨论.其中,对于非局部平面应力状态,通过应变协调方程的Fourier变换形式来进行研究,使问题得以简化.通过以上分析,最终得到一些有价值的结论.     

Compatibility conditions and equations of motion in the linear micropolar theory of elasticity

An analog of Cesàro’s formula and several compatibility conditions are given for the three-dimensional and two-dimensional linear micropolar theory of elasticity in the form different from that used in the literature. A number of formulas are obtained to determine the antisymmetric part of the strain (stress) tensor in terms of the symmetric part of the strain tensor and the symmetric part of the bending-torsion (stress and couple-stress) tensor and to determine the antisymmetric part of the bending-torsion (couple-stress) tensor in terms of the symmetric part of the bending-torsion (couplestress) tensor. Some integro-differential equations of motion expressed in terms of the symmetric parts of the stress and couple-stress tensors are proposed for the micropolar theory of elasticity.     

On non-linear Beltrami-Mitchell equations

Beltrami-Mitchell equations for non-linear elasticity theory are derived using the first Piola-Kirchhoff stress and the deformation gradient tensors as field variables so as to yield linear equilibrium and compatibility equations, respectively. In the derivation it is assumed that a strain energy density and, correspondingly, a complementary strain energy density exist, and satisfy the axiom of objectivity. Substitution for the deformation gradient in the compatibility equations yields non-linear differential equations in terms of the first Piola-Kirchhoff stress tensor which may be regarded as the Beltrami-Mitchell equations of non-linear elasticity. The equations are also derived for “semi-linear” isotropic elastic materials and the theory is illustrated by three simple examples.     

DECOMPOSITION OF SYMMETRIC TENSOR AND ITS APPLICATION

In this paper any symmetric tensor is decomposed into the sum of two tensors. One of them is a "type of stress" tensor, and another is a "type of strain" tensor. The inner product space of symmetric tensor is decomposed into the sum of two orthogonal subspaces. The geometric meaning of several principles in the theory of elasticity is given.     

Dislocations in second strain gradient elasticity

A second strain gradient elasticity theory is proposed based on first and second gradients of the strain tensor. Such a theory is an extension of first strain gradient elasticity with double stresses. In particular, the strain energy depends on the strain tensor and on the first and second gradient terms of it. Using a simplified but straightforward version of this gradient theory, we can connect it with a static version of Eringen’s nonlocal elasticity. For the first time, it is used to study a screw dislocation and an edge dislocation in second strain gradient elasticity. By means of this second gradient theory it is possible to eliminate both strain and stress singularities. Another important result is that we obtain nonsingular expressions for the force stresses, double stresses and triple stresses produced by a straight screw dislocation and a straight edge dislocation. The components of the force stresses and of the triple stresses have maximum values near the dislocation line and are zero there. On the other hand, the double stresses have maximum values at the dislocation line. The main feature is that it is possible to eliminate all unphysical singularities of physical fields, e.g., dislocation density tensor and elastic bend-twist tensor which are still singular in the first strain gradient elasticity.     

Invariants for Rari- and Multi-Constant Theories with Generalization to Anisotropy in Biological Tissues

The rari-constant theory of linear elasticity is based on the assumption that elasticity in solids is caused by only pair potentials with coaxial forces acting between atoms. The strain energy of each pair potential depends on the square of the strain between the atoms in the pair. This strain can be determined by taking the inner product of the strain tensor with a structural tensor that is the tensor product of a unit vector with itself. It is shown that the 15 independent constants in the rari-constant theory can be generated by a complete set of 15 structural tensors. It is also shown that the 6 additional independent constants in the multi-constant theory can be generated by taking the inner product of 6 of these structural tensors with the square of the strain tensor. A generalization of these results for nonlinear elasticity is discussed with reference to recent work which compares the structural and generalized structural tensor approaches to modeling fibrous tissues.     

Relation between the stress and couple-stress tensors in the microcontinuum theory of elasticity

The stress tensor is expressed in terms of an arbitrary symmetric tensor field of second rank and the couple-stress tensor. The stress and couple-stress tensors are represented by arbitrary tensor fields satisfying the homogeneous equilibrium equations. These tensors are also given in the form of the expressions satisfying the inhomogeneous equilibrium equations used in the microcontinuum theory of elasticity. The stress tensor functions are considered.     

Dissymétrie de comportement élastique anisotrope couplé ou non à l'endommagement

A thermodynamic framework is proposed to model anisotropic elasticity different in tension and in compression. Based on Kelvin decomposition of the compliance tensor, it applies to any 3D loading. Coupling with damage is made considering fourth- and second-order damage tensors. The proposed formulation automatically satisfies the continuity of the stress tensor and of the energy release rate tensor. The particular case of initially isotropic materials is exposed.     

Wedge and twist disclinations in second strain gradient elasticity

The aim of this paper is to study disclinations in the framework of a second strain gradient elasticity theory. This second strain gradient elasticity has been proposed based on the first and second gradients of the strain tensor by Lazar et al. [Lazar, M., Maugin, G.A., Aifantis, E.C., 2006. Dislocations in second strain gradient elasticity. Int. J. Solids Struct. 43, 1787–1817]. Such a theory is an extension of the first strain gradient elasticity [Lazar, M., Maugin, G.A., 2005. Nonsingular stress and strain fields of dislocations and disclinations in first strain gradient elasticity. Int. J. Eng. Sci. 43, 1157–1184] with triple stress. By means of the stress function method, the exact analytical solutions for stress and strain fields of straight disclinations in an infinitely extended linear isotropic medium have been found. An important result is that the force stress, double stress and triple stress produced by wedge and twist disclinations are nonsingular. Meanwhile, the corresponding elastic strain and its gradients are also nonsingular. Analytical results indicate that the second strain gradient theory has the capacity of eliminating all unphysical singularities of physical fields.     

A variational differential quadrature solution to finite deformation problems of hyperelastic shell-type structures: a two-point formulation in Cartesian coordinates

A new numerical approach is presented to compute the large deformations of shell-type structures made of the Saint Venant-Kirchhoff and Neo-Hookean materials based on the seven-parameter shell theory. A work conjugate pair of the first Piola Kirchhoff stress tensor and deformation gradient tensor is considered for the stress and strain measures in the paper. Through introducing the displacement vector, the deformation gradient, and the stress tensor in the Cartesian coordinate system and by means of the chain rule for taking derivative of tensors, the difficulties in using the curvilinear coordinate system are bypassed. The variational differential quadrature (VDQ) method as a pointwise numerical method is also used to discretize the weak form of the governing equations. Being locking-free, the simple implementation, computational efficiency, and fast convergence rate are the main features of the proposed numerical approach. Some well-known benchmark problems are solved to assess the approach. The results indicate that it is capable of addressing the large deformation problems of elastic and hyperelastic shell-type structures efficiently.