收稿日期:2002-09-09;修改稿收到日期:2 00 3-03-31.
基金项目:国家自然科学基金(50075068);陕西省教育厅科
研基金(
OOJK
181)资助项目 .
作者简介:刘宏昭
*
(1954-),男,教授,博士生导师 .
第21卷第3期
2004 年 6 月
计算力学学报
Chinese Journal of Computational M echanics
Vol
.21,
No
.3
Ju n e
2004
文章编号:1007-4708(2004)03-0303-05
结构阻尼时域本构模型及其应用
刘宏昭
*
, 原大宁, 李冬平, 张忠明
(西安理工大学, 陕西 西安 710048)
摘 要:阻尼合金作为一种新型结构功能材料在不少领域已获应用,由于其阻尼值较大且随频率呈复杂变化关
系,传统的线性粘性阻尼理论或经典的非频变结构阻尼理论难以精确地描述其耗能行为。本文应用粘弹性阻尼
理论,根据阻尼合金储能模量和损耗因子在频域的实测数据,应用最优化方法拟合出标准线性体模型中的本构
参数;根据积分形式的三参量本构关系和变形体虚功原理,推导出了有限元形式的动力学方程;讨论了三参数初
值的选取;对包含卷积积分的有限元动力学方程通过数学推导将其化为三阶线性微分方程组,再转化为标准状
态变量方程,应用数值求解。数值计算实例证明了所提方法的正确和有效性。
关键词:结构阻尼;三参量;阻尼合金;时域本构;有限元法
中图分类号:
TU
31 文献标识码:
A
1 引 言
合理选择阻尼的数学模型,尽可能精确地描述
其力学行为对于正确预估结构的动力特性和合理
设计有重要意义。目前应用最为广泛的是粘性阻尼
和结构阻尼两类模型。根据粘性阻尼理论,振动一
周的耗能与频率成正比,这与许多试验观察结果不
符。人们又提出了结构阻尼摸型,即复刚度阻尼模
型,该模型与金属结构简谐振动时宏观阻尼力在相
当宽的频带内变化平稳这一事实相符,但它只适用
于有限频段内的简谐振动分析。人们在将结构阻尼
模型应用于自由振动和任意激励响应时域分析时
遇到了有悖于物理事实和数学方面的困难
[1, 2]
。为
了提供结构阻尼系统时域动响应实用分析方法,文
献[1]就单自由度结构阻尼系统提出了基于响应优
势谱成分的粘性阻尼近似模型和基于粘弹性理论
的三参数频变结构阻尼模型。文献[3]也对通过组
合标准线性体在时域逼近进行了研究。该方面存在
的主要问题是怎样合理地确定标准线性体模型的
参数,以及提供易于被工程技术人员掌握的有效的
多自由度分析方法。本文的工程应用背景是阻尼合
金结构的动力学分析与设计。阻尼合金也称减振合
金,其损耗因子可达普通钢材的数百倍至数千倍,
它兼具有常规金属材料的机械性能和阻尼材料的
耗能能力,集结构和减振功能于一体。实验表明阻
尼合金的耗能能力在频率域内呈非线性变化。无论
采用粘性阻尼等效理论还是将损耗因子视为常数
的经典结构阻尼模型都不能很好地描述阻尼合金
构件的耗能行为。有效的方法是在实验的基础上拟
合出反映其应力应变和耗能规律的本构关系式,进
而导出其动力学方程。本文应用粘弹性阻尼理论,
根据阻尼合金储能模量和损耗因子在频域的实测
数据,构造出加权优化目标函数,应用最优化方法
拟合出标准线性体模型中的本构参数;根据积分形
式的三参量本构关系和变形体虚功原理,推导出了
有限元形式的动力学方程;讨论了三参量初值的选
取;对包含卷积积分的有限元动力学方程通过数学
推导将其化为三阶线性微分方程组,再转化为标准
状态变量方程,应用数值方法求解。数值计算实例
证明了所提方法的正确和有效性。
2 结构阻尼三参量本构关系式
将减振合金阻尼材料视为一般粘弹性材料,其
本构关系可用粘弹性力学的三参量模型来描述,其
一般形式为
σ+ p
1
σ
·
= q
0
ε+ q
1
ε
·
(1)
对式(1) 两边取拉氏变换,并令 s = iω,则复模量
Y(i ω)=
σ
-
(iω)
ε
-
(iω)
= Y
1
(ω)+ iY
2
(ω)(2)
由上式得储能模量 Y
1
(ω),耗能模量 Y
2
(ω)和损耗
因子 Z (ω) 分别为