结构阻尼时域本构模型及其应用

阻尼合金作为一种新型结构功能材料在不少领域已获应用，由于其阻尼值较大且随频率呈复杂交化关系，传统的线性粘性阻尼理论或经典的非频变结构阻尼理论难以精确地描述其耗能行为。本文应用粘弹性阻尼理论，根据阻尼合金储能模量和损耗因子在频域的实测数据．应用最优化方法拟合出标准线性体模型中的本构参数；根据积分形式的三参量本构关系和变形体虚功原理，推导出了有限元形式的动力学方程；讨论了三参数初值的选取；对包含卷积积分的有限元动力学方程通过数学推导将其化为三阶线性微分方程组，再转化为标准状态变量方程，应用数值求解。数值计算实例证明了所提方法的正确和有效性。     

非粘滞阻尼系统时程响应分析的精细积分方法

考虑一个具有非粘滞阻尼特性的多自由度系统响应的时程分析问题.该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式体现为阻尼力等于质点速度与某一核函数的卷积.在利用状态空间方法将系统运动方程转换成一阶的状态方程的基础上,采用精细积分方法对状态方程进行数值求解,得到一种求解该阻尼系统时程响应的精确、高效的计算方法.通过两个数值算例表明,采用该方法得到几乎精确的数值计算结果,而且计算效率有成数量级的提高.     

热粘塑性体的积分—微分型本构关系

应用ИЛЪЮШИН关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史泛；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分型本构方程这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价。对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实     

热粘塑性体的积分－微分型本构关系

应用　　　关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内＊时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史的泛函；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分到本构方程．这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价．对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实验结果符合较好．     

率相关本构方程积分新算法

提出一种积分率相关本构方程的隐式积分新算法,引入0～1范围内的缩放因子λ对本构方程进行间接求解,可以避免直接求解等效塑性应变或等效塑性应变率时,由于其数值过大或过小而造成的收敛困难或收敛失败,实现对率相关本构方程的快速准确求解.以B-P统一本构方程及双曲正弦本构方程为例,验证了算法的可行性.结果表明,新算法对于准静态变形条件下的无硬化本构方程也可以得出准确的解.     

粘弹性结构动力稳定性分析的谐波平衡法

本文分析了线粘弹性结构的长期动力稳定特性。设材料具积分型本构关系，且其松弛模量能用Prony级数描述，将微分-积分型控制方程化成微分型方程，应用谐波平衡法确定动力稳定区域，着重讨论了材料参数及系统振动频率对动力稳定区域的影响发现该类粘弹性结构具有与一般阻尼系统不同的动力稳定特性。文中也将系统平衡法直接应用于微分-积分型控制方程，忽略卷积积分运算后产生的随时间衰减的非谐波项来得到决定动力稳定边界的特征方程，并对两种应用所得结果进行了比较。     

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

固体推进剂非线性拉伸行为的实验研究

固体推进剂非线性拉伸行为的实验研究王晓明，沈亚鹏（西安交通大学，７１００４９）关键词本构关系，粘弹性，非线性，推进剂１引言Ｈａｄｄａｄ（１）和杨挺青（２）对粘弹性非线性本构方程作了比较全面的综述，其中最简单，而又有利于实验确定材料常数的本构方程是其中...     

粘弹性板的非线性动力稳定特性分析

采用Boltzman积分型本构关系,分析了线粘弹性薄板在考虑几何线性与非线性时的长期动力稳定特性,设材料为标准线性固体,将系统的微分一积分型控制方程转化成微分型控制方程,由增量谐波平衡法确定主要动力不稳定区域的边界,发现粘弹性结构具有与一般阻尼系统不同的动力稳定特性,由于材料的粘性阻尼与松弛效应的综合影响,动力不稳定区域有不同程度的缩小与偏移,且在考虑几何线性与非线性情形下,其影响程度又不一样。     

Dynamical behavior of nonlinear viscoelastic beams

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＧａｌｅｒｋｉｎｔｒｕｎｃａｔｉｏｎｉｓｗｉｄｅｌｙｕｓｅｄｔｏｓｔｕｄｙｔｈｅｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒ(ｅｓｐｅｃｉａｌｌｙｔｈｅｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒ)ｏｆｓｔｒｕｃｔｕｒｅｓ[1].Ｈｏｗｅｖｅｒ,ａｓｆａｒｔｈｅｒｅｉｓｎｏｄｉｒｅｃｔｅｖｉｄｅｎｃｅｔｏｐｒｏｖｅｔｈｅｐｌａｕｓｉｂｉｌｉｔｙｏｆｔｈｅｌｏｗｏｒｄｅｒＧａｌｅｒｋｉｎｔｒｕｎｃａｔｉｏｎ,ａｌｔｈｏｕｇｈｉｔｃａｎｂｅｉｎｆｅｒｒｅｄｆｒｏｍｃｅｒｔａｉｎｉ…     

AN INTERNAL STATE VARIABLE VISCOELASTIC-VISCOPLASTIC CONSTITUTIVE EQUATION WITH DAMAGE

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.     

Practical comparison of differential viscoelastic constitutive equations in finite element analysis of planar 4:1 contraction flow

Finite element modeling of planar 4:1 contraction flow (isothermal incompressible and creeping) around a sharp entrance corner is performed for favored differential constitutive equations such as the Maxwell, Leonov, Giesekus, FENE-P, Larson, White-Metzner models and the Phan Thien-Tanner model of exponential and linear types. We have implemented the discrete elastic viscous stress splitting and streamline upwinding algorithms in the basic computational scheme in order to augment stability at high flow rate. For each constitutive model, we have obtained the upper limit of the Deborah number under which numerical convergence is guaranteed. All the computational results are analyzed according to consequences of mathematical analyses for constitutive equations from the viewpoint of stability. It is verified that in general the constitutive equations proven globally stable yield convergent numerical solutions for higher Deborah number flows. Therefore one can get solutions for relatively high Deborah number flows when the Leonov, the Phan Thien-Tanner, or the Giesekus constitutive equation is employed as the viscoelastic field equation. The close relationship of numerical convergence with mathematical stability of the model equations is also clearly demonstrated.     

A viscoelastic constitutive model of rubber under small oscillatory load superimposed on large static deformation

Summary  A viscoelastic constitutive equation of rubber that is under small oscillatory load superimposed on large static deformation is proposed. The model is derived through linearization of Simo's nonlinear viscoelastic constitutive model and reference configuration transformation. Most importantly, in this model, static deformation correction factor is introduced to consider the influence of pre-strain on the relaxation function. Natural statically pre-deformed state is served as reference configuration. The proposed constitutive equation is extended to a generalized viscoelastic constitutive equation that includes widely used Morman's model as a special case using objective stress increment. The proposed constitutive model is tested for dynamic behavior of rubber specimens with different carbon black content. It is concluded from the test that the assumption that the effects of static deformation can be separated from time effects, which is the basis of Morman's model, is only applicable to unfilled rubber. The viscoelastic constitutive equation for filled rubber must include, therefore, the influence of the static deformation on the time effects. The suggested constitutive equation with static deformation correction factor shows good agreement with test values. Received 4 January 2001; accepted for publication 13 June 2001     

A large deformation framework for compressible viscoelastic materials: Constitutive equations and finite element implementation

For modeling the constitutive properties of viscoelastic solids in the context of small deformations, the so-called three-parameter solid is often used. The differential equation governing the model response may be derived in a thermodynamically consistent way considering linear spring-dashpot elements. The main problem in generalizing constitutive models from small to finite deformations is to extend the theory in a thermodynamically consistent way, so that the second law of thermodynamics remains satisfied in every admissible process. This paper concerns with the formulation and constitutive equations of finite strain viscoelastic material using multiplicative decomposition in a thermodynamically consistent manner. Based on the proposed constitutive equations, a finite element (FE) procedure is developed and implemented in an FE code. Subsequently, the code is used to predict the response of elastomer bushings. The finite element analysis predicts displacements and rotations at the relaxed state reasonably well. The response to coupled radial and torsional deformations is also simulated.     

粘弹性地基上粘弹性输流管道的稳定性分析

从Winkler假设和单轴线性粘弹性本构方程出发，推导了Kelvin-Voigt粘弹性地基上三参量固体模型输流管道的运动微分方程，采用改进的有限差分法，分析了管道和地基的粘弹性参数对输流管道无量纲复频率和无量纲流速之间的变化关系的影响。     

非线性粘弹性大挠度梁的动力学模型及其简化

本文建立了描述几何非线性均匀梁动力学行为的偏微分－－积分方程,梁的材料满足Ｌｅａｄｅｍａｎ非线性本构关系,对于两端简支的情形用Ｇａｌｅｒｋｉｎ方法进行了截断简化为常微分－－积分方程,然后引进附加变量的方法进一步简化为常微分方程组。     

NONLINEAR DYNAMIC RESPONSE AND ACTIVE VIBRATION CONTROL OF THE VISCOELASTIC CABLE WITH SMALL SAG

IntroductionCablesareveryefficientstructuralmembersandhencehavebeenwidelyusedinmanylong_spanstructures,includingcable_supportbridges,guyedtowersandcable_supportroofs.Sincecablesarelight,veryflexibleandlightlydamped ,structuresutilizingcables,i.e .,cable_structuresystems,usuallyhavevariousdynamicproblems.Theirmodelsarethereforeverimportantinpredictingandcontrollingtheirresponses.Inthelastdecade,thenonlineardynamicvibrationandstabilitybehaviorofcablesandcable_structureshavedrawntheattentionofman…