基于绝对节点坐标法的弹性线方法研究

相比于浮动坐标系法, 绝对节点坐标法(absolute nodal coordinateformulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势,ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标, 具有常数质量阵,不存在科氏离心力等优点, 然而弹性力阵为非线性项,其求解将比较耗时且占用资源. 据此, 在弹性力求解方法中,引入弹性线方法(elastic line method, ELM),该方法将格林--拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变, 转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程, 通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法(continuum mechanics method,CMM)和弹性线方法在求解弹性力上的不同点, 并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明:弹性线方法能在保证精度的前提下有效提高计算效率.      

基于绝对节点坐标法的弹性线方法研究

相比于浮动坐标系法,绝对节点坐标法(absolute nodal coordinate formulation, ANCF)在处理柔性体非线性大变形问题上具有显著优势, ANCF将单元节点坐标定义在全局坐标系下,采用斜率矢量代替节点转角坐标,具有常数质量阵,不存在科氏离心力等优点,然而弹性力阵为非线性项,其求解将比较耗时且占用资源.据此,在弹性力求解方法中,引入弹性线方法 (elastic line method, ELM),该方法将格林–拉格朗日应变张量定义在中心线上,采用曲率公式来定义弯曲应变,转角公式来定义扭转应变.同时采用有限元法对三维柔性梁位移场进行离散,求解梁单元常数质量阵、广义刚度阵、广义力阵,进而得到单元的动力学方程,通过转换矩阵得到三维梁的动力学方程.接着从理论上指出连续介质力学方法 (continuum mechanics method, CMM)和弹性线方法在求解弹性力上的不同点,并编制动力学仿真软件.最后分别采用连续介质力学方法和弹性线方法对柔性单摆以及履带式车辆的动力学问题进行仿真分析,结果表明:弹性线方法能在保证精度的前提下有效提高计算效率.     

热-动力学耦合多体系统建模与降阶

在轨运行的卫星天线受到太阳辐射热冲击后容易出现热致振动或指向不准确等问题, 严重时会导致航天器失效. 本文提出了一种基于改进模态综合法的刚柔热耦合多体系统建模与降阶方法. 采用绝对节点坐标法单元形函数对柔性天线的位移场与温度场进行统一离散插值, 避免了两种物理场网格不匹配带来的映射误差与效率问题, 并使用绝对节点坐标参考节点描述中心刚体. 在系统方程中考虑了热流输入和表面自热辐射. 针对绝对节点坐标法切线刚度阵高度非线性的特点, 利用一阶泰勒展开对系统动力学和传热学方程进行了分段线性化, 在线性化区间内切线刚度矩阵为常数矩阵, 避免了每个时间步上的弹性力及其雅克比矩阵的迭代计算, 并使得基于模态的降阶手段得以应用. 利用改进的模态综合方法划分子结构并缩减系统自由度. 相邻子结构之间通过约束方程保证连接精度和连续性. 通过纯导热半圆形薄板、薄板的热膨胀、柔性太阳能电池板和刚柔热耦合抛物线天线四个数值算例验证本文所提出方法的有效性. 结果表明, 本文提出的方法在保证仿真精度的前提下缩减了系统规模, 提高了仿真计算效率.      

弹性大变形问题中的应力状态描述、奇异性问题和余能原理

对弹性大变形理论中的3方面问题进行了综述.首先,对各种应变度量的共轭应力进行综述.大变形问题引起的应力状态描述的复杂性引起了许多学者的兴趣,对这个问题的研究也促进了大变形弹性理论的发展.在各种特定问题中,人们提出了不同的应力张量来描述应力状态,如Caucby应力张量、第一类和第一二类Piola-Kirchhoff应力张...     

含分数阻尼特性元件的多体系统动力学研究

在绝对节点坐标体系下研究了具有分数导数阻尼特性元件的多体系统动力学建模、求解问题. 采用基于绝对节点坐标的无闭锁效应剪变梁单元离散柔性构件,建立了含常数质量矩阵的系统动力学方程, 并采用数值耗散可控的广义a方法求解. 通过数值算例计算,对比研究了算法参数与阻尼项的分数指数对系统动力学响应的影响规律.该方法可以进一步扩展到众多工程实际问题研究中.      

正交曲线坐标系下弹性力学的应变张量——Jacobi矩阵的一个应用

借助于Jacobi矩阵的链式法则来导出曲线坐标系下的应变张量，导出曲线坐标系下矩阵形式的几何方程，作为例子，给出球坐标下的几何方程。     

中心刚体-柔性梁系统的耦合变形影响

本文对作大范围运动的中心刚体-柔性梁系统的耦合变形的影响进行研究.给出一种新的描述柔性梁耦合变形的有限元插值方法,该方法采用笛卡尔变形坐标对横向变形和纵向变形之间的耦合项进行描述,该耦合变形项只与本单元的节点变形坐标相关.分别讨论了耦合变形项对惯性力与弹性力的贡献,分析了它们对刚-柔耦合动力学行为的影响.通过研究指出当采用笛卡尔变形坐标描述时,如果在计算弹性力的时候考虑了耦合变形影响,无论在计算惯性力时是否考虑耦合变形影响,都可以得到稳定收敛的结果.反之,如果在计算弹性力时忽略了耦合变形影响,无论在计算惯性力时是否考虑耦合变形影响,当大范围运动的速度较高时,将会得到错误的发散的结果.因此,通过忽略耦合变形对质量分布的影响,只保留耦合变形对弹性力的影响,可实现对动力学方程的简化.     

求解大柔度梁/绳索的ANCF单元及隐式迭代格式

为模拟大柔度梁/绳索结构的变形和大范围运动,基于绝对节点坐标方法ANCF(Absolute nodal coordi-nate formulation)和HHT(Hilber-Hughes-Taylor)积分方法,建立了ANCF单元的隐式动力学迭代格式.得到了简洁的节点等效力向量,且进一步导出了切线刚度矩阵的全部公式,...     

空间展开折叠桁架结构动力学分析研究

本文以笛卡尔坐标系下节点自然坐标为未知量,建立了桁架结构系的基本运动力学方程,并首次推导出桁架结构中常用节点附加几何约束方程,相应约束Jacobi矩阵及其导数矩阵,采用奇异值分解法求约束Jacobi矩阵的零空间基和M-P广义逆,并由矩阵缩减法建立了带约束桁架体系的运动力学方程和求解方法。数值算例表明该方法适于可展折叠桁架结构运动力学分析。     

用于准脆性材料破坏过程数值模拟的虚拟连接单元法

将有限变形单元与虚拟连接单元相结合,用于模拟准脆性材料破坏过程.首先基于精确的有限变形理论,采用第二类Piola-Kirchhoff应力与Green-Lagrange应变作为能量共轭的应力、应变对,推导出虚拟连接单元的单元刚度矩阵;通过数值算例,验证该单元的正确性与合理性,给出虚拟连接单元高度的取值范围,并与无此单元时...     

Convergence of Peridynamics to Classical Elasticity Theory

The peridynamic model of solid mechanics is a nonlocal theory containing a length scale. It is based on direct interactions between points in a continuum separated from each other by a finite distance. The maximum interaction distance provides a length scale for the material model. This paper addresses the question of whether the peridynamic model for an elastic material reproduces the classical local model as this length scale goes to zero. We show that if the motion, constitutive model, and any nonhomogeneities are sufficiently smooth, then the peridynamic stress tensor converges in this limit to a Piola-Kirchhoff stress tensor that is a function only of the local deformation gradient tensor, as in the classical theory. This limiting Piola-Kirchhoff stress tensor field is differentiable, and its divergence represents the force density due to internal forces. The limiting, or collapsed, stress-strain model satisfies the conditions in the classical theory for angular momentum balance, isotropy, objectivity, and hyperelasticity, provided the original peridynamic constitutive model satisfies the appropriate conditions.      

Linear Constitutive Relations in Isotropic Finite Elasticity

For simple shearing and simple extension deformations of a homogeneous and isotropic elastic body, it is shown that a linear relation between the second Piola-Kirchhoff stress tensor and the Green-St. Venant strain tensor does not predict a physically reasonable response of the body. This constitutive relation implies that the slope of the curve between an appropriate component of the first Piola-Kirchhoff stress tensor and a deformation measure is an increasing functions of the deformation measure. This revised version was published online in August 2006 with corrections to the Cover Date.     

Objective constitutive relations in elasticity and viscoelasticity

Summary The compatibility between the objectivity principle and affine constitutive equations for the elastic Cauchy and Piola-Kirchhoff stress tensors with non-zero residual stress is examined. It is found that the Cauchy stress is allowed to be only a constant tensor, proportional to the identity tensor, while the Piola-Kirchhoff stress may be a linear function on the deformation gradient thus generalizing previous results by Fosdick and Serrin. The same conclusions are arrived at also by starting from viscoelasticity. Finally, in the case of Maxwell-like materials, the solutions to the objective evolution equations are shown to be objective functionals.

---

Sommario Si esamina la compatibilità tra il principio di obiettività ed equazioni costitutive affini per i tensori di stress elastici di Cauchy a Piola-Kirchhoff con stress residuo non nullo. Generalizzando risultati di Fosdick e Serrin si prova che il tensore di Cauchy può essere soltanto un tensore costante, proporzionale al tensore identità, mentre il tensore di Piola-Kirchhoff può essere una funzione lineare del gradiente di deformazione. Alle stesse conclusioni si perviene anche partendo dal funzionale della viscoelasticità. Infine si mostra che, per materiali tipo Maxwell, le soluzioni di equazioni di evoluzione obiettive sono funzionali obiettivi.

     

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.     

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.     

On Non-Symmetric Deformations of an Incompressible Nonlinearly Elastic Isotropic Sphere

A class of non-symmetric deformations of a neo-Hookean incompressible nonlinearly elastic sphere are investigated. It is found via the semi-inverse method that, to satisfy the governing three-dimensional equations of equilibrium and the incompressibility constraint, only three special cases of the class of deformation fields are possible. One of these is the trivial solution, one the solution describing radially symmetric deformation, and the other a (non-symmetric, non-homogeneous) deformation describing inflation and stretching. The implications of these results for cavitation phenomena are also discussed. In the course of this work, we also present explicitly the spherical polar coordinate form of the equilibrium equations for the nominal stress tensor for a general hyperelastic solid. These are more complicated than their counterparts for Cauchy stresses due to the mixed bases (both reference and deformed) associated with the nominal (or Piola-Kirchhoff) stress tensor, but more useful in considering general deformation fields. This revised version was published online in August 2006 with corrections to the Cover Date.     

Couple stress based strain gradient theory for elasticity

The deformation behavior of materials in the micron scale has been experimentally shown to be size dependent. In the absence of stretch and dilatation gradients, the size dependence can be explained using classical couple stress theory in which the full curvature tensor is used as deformation measures in addition to the conventional strain measures. In the couple stress theory formulation, only conventional equilibrium relations of forces and moments of forces are used. The couple's association with position is arbitrary. In this paper, an additional equilibrium relation is developed to govern the behavior of the couples. The relation constrained the couple stress tensor to be symmetric, and the symmetric curvature tensor became the only properly conjugated high order strain measures in the theory to have a real contribution to the total strain energy of the system. On the basis of this modification, a linear elastic model for isotropic materials is developed. The torsion of a cylindrical bar and the pure bending of a flat plate of infinite width are analyzed to illustrate the effect of the modification.     

A note on nonassociated plastic flow rules

The analytical properties of the constitutive equations in plasticity with a nonassociated flow rule are investigated. Under the assumption of small deformations the directional stiffness (and compliance) rule is considered and the relevant spectral properties of the tangent stiffness tensor are assessed. It is shown that the directional stiffness may be larger than the elastic. It may also be negative in the case of a formally perfectly plastic material and so the nonassociative flow rule represents (spurious) softening in terms of an associated flow rule. The issue of uniqueness at finite strains is briefly addressed, whereby use is made of the tangent stiffness tensor relating the velocity gradient to the first Piola-Kirchhoff stress rate. The relevant spectral properties, which generalise those from the small deformation case, are found explicit. A sufficient condition for uniqueness is given in terms of a critical (upper bound) value of the hardening modulus.     

拉格朗日几何非线性混合应变元

本文应用由Ｓｉｍｏ和Ｒｉｆａｉ建议的混合假定附加应变途径，采用第二Ｐｉｏｌａ－Ｋｉｒｃｈｈｏｆｆ应力张量和Ｇｒｅｅｎ－Ｌａｇｒａｎｇｅ应变张量作能量共轭的应力应变度量，导出了Ｌａｇｒａｎｇｅ几何非线性下的胡海昌－Ｗａｓｈｉｚｕ三变量变分原理的Ｇａｌｅｒｋｉｎ形式以及相应的混合假定应变元公式。