论拖带坐标系中应力的客观速率

应力客观速率是有限变形力学中一个十分重要的问题,本文利用非线性几何场论方法,推导出拖带坐标系中应力客观速率公式,并应用新公式计算了拉伸、拉伸与转动复合、简单剪切大变形问题.通过将计算结果与用Jaumann等其它应力速率计算的结果进行比较,说明了本文所得结果是合理的.     

简单剪切振荡现象及弹塑性本构的限制条件

基于在真应力空间刻划弹塑性物质的强化、软化和理想塑性特性,本文以刚塑性随动强化模型为例说明产生简单剪切振荡现象,是与模型在简单剪切变形下,其强化和软化特性发生交替变化的现象有关。为使弹塑性本构模型更符合实际,要求它们必须满足如下条件:即对任意弹塑性加载变形过程,本构模型所给出的应力应该是非振荡的,所描述的强化和软化特性,不存在从软化阶段至强化阶段的过渡。     

形状记忆合金相变过程三维大变形有限元模拟

形状记忆合金（SMA）一直被作为智能材料开发,并被用于阻尼器、促动器和智能传感器元件．形状记忆合金（SMA）的一项重要特性,是它具有恢复在机械加卸载周期下产生的大变形而不表现出永久变形的能力．该文旨在介绍一种由应力产生的相变且可以描述马氏体和奥氏体之间的超弹性滞回环现象本构方程．形状记忆合金的马氏体系数假设为应力偏张量的函数,因此形状记忆合金在相变过程中锁定体积．本构模型是在大变形有限元的基础上执行的,采用了现时构型Lagrange大变形算法．为了方便地使用Cauchy应力和线性应变本构关系,使用了与旋转无关的Jaumann应力增率计算应力．数值分析结果表明,相变引起的超弹性滞回环可以有效地通过该文提出的本构方程和大变形有限元模拟．     

粘弹性流体在渐缩毛管中渗流规律研究

驱油用聚合物溶液是粘弹性流体,当其在多孔介质中流动时由于拉伸变形会产生额外耗散,因此拉伸压降是不可忽略的.将Binding等人通过机械能平衡原理和最小能量原理建立的非牛顿流体入口收敛流动理论推广应用于渐缩毛管,并以幂指数关系表征剪切粘度与剪切速率关系、第一法向应力差与剪切速率关系、拉伸粘度与拉伸速率关系,从而建立了适用于描述粘弹性流体流动特性的渗流模型,分析了粘弹性流体在渐缩毛管中的渗流规律,利用建立的模型可以准确反映粘弹流体的渗流特性.     

弹塑性有限变形的广义Prandtl—Reuss本构方程和应力共旋率研究

本文通过一种新的途径研究弹塑性有限变形的广义Ｐｒａｎｄｔｌ＿Ｒｅｕｓ本构方程·研究表明对于广义Ｐｒａｎｄｔｌ＿Ｒｅｕｓ本构方程，变形率弹塑性和分解的假设并非必须·研究了采用物质共旋率的广义Ｐｒａｎｄｔｌ＿Ｒｅｕｓ本构方程，从理论上分析了简单剪切应力振荡的原因·提出一种用于构造广义Ｐｒａｎｄｔｌ＿Ｒｅｕｓ本构方程中应力和背应力共旋率的修正相对旋率·最后，对简单剪切变形进行应力计算·     

非线性弹性大变形问题的率型解法

本文分别采用Ｊａｕｍａｎｎ应力率、Ｔｒｕｅｓｄｅｌｌ应力率和Ｇｒｅｅｎ－Ｎａｇｈｄｉ应力率导出了非线性各向同性弹性体的率型本构表达形式，通过对Ｍｏｏｎｅｙ－Ｒｉｖｌｉｎ材料的简单剪切大变形分析表明，三种率型的本构关系均与全量本构关系相等价。文中还给出了相应的率型变分原理，并采应Ｒｉｔｚ法，数值求解了受单轴拉伸的矩形橡皮薄膜的大变形问题．     

弹性动脉血管中血液流动特性的模拟和分析

研究了两个不同的非牛顿血液流动模型:低粘性剪切简单幂律模型和低粘性剪切及粘弹性振荡流的广义Maxwell模型．同时利用这两个非牛顿模型和牛顿模型,研究了磁场中刚性和弹性直血管中血液的正弦型脉动．在生理学条件下,大动脉中血液的弹性对其流动性态似乎并不产生影响,单纯低粘性剪切模型可以逼真地模拟这种血液流动．利用高剪切幂律模型模拟弹性血管中的正弦型脉动流,发现在同一压力梯度下,与牛顿流体相比较,幂律流体的平均流率和流率变化幅度都更小．控制方程用Crank-Niclson方法求解．弹性动脉中血液受磁场作用是产生此结果的直观原因．在主动脉生物流的模拟中,与牛顿流体模型比较,发现在匹配流率曲线上,幂律模型的平均壁面剪切应力增大,峰值壁面剪切应力减小．讨论了弹性血管横切磁场时的血液流动,评估了血管形状和表面不规则等因素的影响．     

基于第二粘性用局部微分求积法计算激波

考虑第二粘性效应,采用局部微分求积法数值求解激波问题．首先解释了在激波计算时,有必要考虑第二粘性,然后基于粘性模型,对一维和二维激波进行了数值模拟,还分别考察了剪切粘性力和第二粘性力对数值结果的影响．结果表明,采用粘性模型加上局部微分求积法能够模拟出激波特征,具有客观、简单的优点．     

有限变形下的混凝土动态本构关系研究

讨论有限变形和小变形假设下本构关系的区别,并将其运用于混凝土的弹-粘塑性本构关系研究,提出了一个应变率相关的动态力学模型．模型基于Ottosen的4参数屈服准则,分别考虑混凝土在硬化阶段和软化阶段加载面的不同变化规律,建立冲击荷载下的混凝土本构关系．该模型可以应用于冲击载荷下混凝土材料响应的模拟．引进Green-Naghdi客观率建立有限变形的混凝土模型．根据大量实验结果对应变率和材料强度的关系提出合理假设,使模型可以反映混凝土大变形的动态力学行为,为相关工程问题的研究提供有益的思路和有效的工具．     

剪切变形下非晶态高聚物的力学行为

基于非平衡态热力学理论，提出了一个适用于不可压材料的新的热粘弹性本构模型．该模型将橡胶弹性理论中的非高斯分子网络模型推广到计及粘性和热效应的情形．通过引入一组二阶张量形式的内变量，建议了一个新的Helmholtz自由能表达式，从而可以用来合理描述内变量的演化规律．根据以上模型，重点研究了热粘弹性材料在简单剪切变形下的力学行为，考察了由于分子链取向分布的变化而产生的“粘性耗散诱导”各向异性，讨论了应变率效应和由于粘性耗散而导致的热软化效应对剪应力的影响．理论预测结果与G’Sell等人的实验数据的定性比较表明了新的本构模型的有效性。     

Nonlinear,finite deformation,finite element analysis

The roles of the consistent Jacobian matrix and the material tangent moduli, which are used in nonlinear incremental finite deformation mechanics problems solved using the finite element method, are emphasized in this paper, and demonstrated using the commercial software ABAQUS standard. In doing so, the necessity for correctly employing user material subroutines to solve nonlinear problems involving large deformation and/or large rotation is clarified. Starting with the rate form of the principle of virtual work, the derivations of the material tangent moduli, the consistent Jacobian matrix, the stress/strain measures, and the objective stress rates are discussed and clarified. The difference between the consistent Jacobian matrix (which, in the ABAQUS UMAT user material subroutine is referred to as DDSDDE) and the material tangent moduli (C^e) needed for the stress update is pointed out and emphasized in this paper. While the former is derived based on the Jaumann rate of the Kirchhoff stress, the latter is derived using the Jaumann rate of the Cauchy stress. Understanding the difference between these two objective stress rates is crucial for correctly implementing a constitutive model, especially a rate form constitutive relation, and for ensuring fast convergence. Specifically, the implementation requires the stresses to be updated correctly. For this, the strains must be computed directly from the deformation gradient and corresponding strain measure (for a total form model). Alternatively, the material tangent moduli derived from the corresponding Jaumann rate of the Cauchy stress of the constitutive relation (for a rate form model) should be used. Given that this requirement is satisfied, the consistent Jacobian matrix only influences the rate of convergence. Its derivation should be based on the Jaumann rate of the Kirchhoff stress to ensure fast convergence; however, the use of a different objective stress rate may also be possible. The error associated with energy conservation and work-conjugacy due to the use of the Jaumann objective stress rate in ABAQUS nonlinear incremental analysis is viewed as a consequence of the implementation of a constitutive model that violates these requirements.     

受轴向压缩圆柱壳塑性屈曲的内时分析

采用内时塑性本构方程的增量和全量表达式分析了受轴向压缩圆柱壳的塑性屈曲，得到了塑性屈曲临界应力与圆柱壳特征尺寸间的关系。对ＡＭＦ和铝合金圆柱壳塑性屈曲进行了分析，与实验结果的比较表明：除对于ＡＭＦ圆柱壳由内时塑性本构方程的全量表达式给出了较经典塑性理论全量分析略为保守的结果外，在其它杨合下，内时分析均给出了较经典塑性理论更符合实验数据的结果。     

Computational implementation of stress integration in FE analysis of elasto-plastic large deformation problems

In this paper, we compare different numerical implementation algorithms for the rate type constitutive equation and present an integration scheme based on the physical meaning of the stress. Numerical implementation of various schemes is investigated in conjunction with the return mapping algorithm and the conditions to maintain plastic consistency. Jaumann and Truesdell rates are taken as the objective stress rates in the constitutive equation. An alternative numerical treatment for rate of deformation tensor D_ij is presented and is shown to maintain incremental objectivity. Numerical examples included a single element under rigid body rotation, a necking bifurcation of a bar in tension and a punch indentation process. It is shown that the use of Truesdell stress rate with specific numerical integration procedure gives more accurate results than other procedures presented.     

Development of a 2D Isoparametric Finite-Element Model Based on Reddy’s Third-Order Theory for the Bending Behavior Analysis of Composite Laminated Plates

Mechanics of Composite Materials - The objective of this work is to propose a simple efficient finite element, based on Reddy’s third-order shear deformation theory, for analyzing the bending...     

A time integration procedure for plastic deformation in elastic-viscoplastic metals

A simple unconditionally stable numerical procedure for time integration of the flow rule for large plastic deformation of an elastic-viscoplastic metal is developed. Specific attention is focused on a unified set of constitutive equations which represents a generalization (for large deformation and thermomechanical response) of the Bodner-Partom model [6, 7]. An analytical solution is obtained for large deformation simple shear at constant shear rate. Numerical examples of simple shear, a corner test exhibiting the transition from uniaxial compression to shear, and simple tension are considered which demonstrate the stability and accuracy of the procedure. It is shown that the same procedure can be used for a rate insensitive metal characterized by a yield function as well as for a rate sensitive metal characterized by an overstress model. Finally, an appendix is provided which records the basic equations associated with the small deformation theory.     

Finite element dynamic analysis of unsymmetric composite laminated beams with shear effect and rotary inertia under the action of moving loads

The finite element dynamic response of an unsymmetric composite laminated orthotropic beam, subjected to moving loads, has been studied. One-dimensional finite element based on classical lamination theory, first-order shear deformation theory, and higher-order shear deformation theory having 16, 20 and 24 degrees of freedom, respectively, are developed to study the effects of extension, bending, and transverse shear deformation. The theories also account for the Poisson effect, thus, the lateral strains and curvatures can be expressed in terms of the axial and transverse strains and curvatures and the characteristic couplings (bend–stretch, shear–stretch and bend–twist couplings) are not lost. The dynamic response of symmetric cross-ply and unsymmetric angle-ply laminated beams under the action of a moving load have been compared to the results of an isotropic simple beam. The formulation also has been applied to the static and free vibration analysis.