内容发布更新时间 : 2024/12/27 21:02:36星期一 下面是文章的全部内容请认真阅读。

为了确定对应给定阻尼比ξ的α和β值，通常假定α和β之和在某个频率范围内近似为恒定值（见图5）。这样，在给定阻尼比ξ和一个频率范围ωi～ωj后，解两个并列方程组便可求得α和β。

图5瑞利阻尼

Alpha阻尼在模型中引入任意大质量时会导致不理想的结果。一个常见的例子是在结构的基础上加一个任意大质量以方便施加加速度谱（用大质量可将加速度谱转化为力谱）。Alpha阻尼系数在乘上质量矩阵后会在这样的系统中产生非常大的阻尼力，这将导致谱输入的不精确，以及系统响应的不精确。

Beta阻尼和材料阻尼在非线性分析中会导致不理想的结果。这两种阻尼要和刚度矩阵相乘，而刚度矩阵在非线性分析中是不断变化的。由此所引起的阻尼变化有时会和物理结构的实际阻尼变化相反。例如，存在由塑性响应引起的软化的物理结构通常相应地会呈现出阻尼的增加，而存在Beta阻尼的ANSYS模型在出现塑性软化响应时则会呈现出阻尼的降低。

2. 和材料相关的阻尼

和材料相关的阻尼允许将Beta阻尼做为材料性质来指定[MP,DAMP]。但要注意在谱分析[ANTYPE,SPECTR]中的MP,DAMP是指定和材料相关的阻尼比ξ，而不是β。同样要注意对于多材料单元如SOLID46,SOLID65,SHELL91和SHELL99，只能对单元整体指定一个β值，而不能对单元中的每一种材料都指定。在这些情形下，β是由单元的材料指针（用MAT命令设置）决定的，而不是由单元实常数MAT指向的材料决定的。

3. 恒定阻尼比

恒定阻尼比是在结构中指定阻尼的最简单的方法。它表示实际阻尼和临界阻尼之比，是用DMPRAT命令指定的小数值。DMPRAT只可用于谱分析、谐响应分析和模态叠加法瞬态动力学分析。

4. 振型阻尼

振型阻尼可用于对不同的振动模态指定不同的阻尼比。它用MDAMP命令指定且只能用于谱分析和模态叠加法瞬态动力学分析、谐响应分析。

5. 单元阻尼

单元阻尼在用到有粘性阻尼特征的单元类型时会涉及到，如单元COMBIN7, COMBIN14,COMBIN37,COMBIN40等。

关于阻尼的更详细描述参见<>。

阻尼是动力分析的一大特点，也是动力分析中的一个易于引起困惑之处，而且由于它只是影响动力响应的衰减，出了错不容易觉察。阻尼的本质和表现是相当复杂的，相应的模型也很多。ANSYS提供了强大又丰富的阻尼输入，但也正以其强大和丰富使初学者容易发生迷惑这里介绍各种阻尼的数学模型在ANSYS中的实现，与在ANSYS中阻尼功能的使用。 1．比例阻尼

最常用也是比较简单的阻尼大概是Rayleigh阻尼，又称为比例阻尼。它是多数实用动力分析的首选，对许多实际工程应用也是足够的。在ANSYS里，它就是 阻尼与 阻尼之和，分别用ALPHD与BETAD命令输入。已知结构总阻尼比是 ，则用两个频率点上 阻尼与 阻尼产生的等效阻尼比之和与其相等，就可以求出近似的 阻尼与 阻尼系数来用作输入： （5.1.1）

求比例阻尼系数的拟合公式

用方程组（5.1.1）可以得到 阻尼与 阻尼系数值，然后用ALPHD与BETAD命令输入，这种阻尼输入既可以做full（完全）法的分析，也可以作减缩法与振型叠加法的分析，都是一样的有效。

但是尽管 阻尼与 阻尼概念简单明确，在使用中也要小心一些可能的误区。首先， 阻尼与质量有关，主要影响低阶振型，而 阻尼与刚度有关，主要影响高阶振型；如果要做的是非线性瞬态分析，同时刚度变化很大时，那么使用 阻尼很可能会造成收敛上的困难；一样的理由，有时在使用一些计算技巧时，比如行波效应分析的大质量法，加上了虚假的大人工质量，那么就不可以使用 阻尼。同样，在模型里加上了刚性连接时，也应该检查一下 阻尼会不会造成一些虚假的计算结果。 2．阻尼阵的计算

ANSYS中有多种办法可以输入阻尼特性。先概括几个在结构分析中常用的输入阻尼的命令： ALPHAD： 输入 阻尼参数 BETAD： 输入 阻尼参数 DMPRAT： 输入全结构的阻尼比

MDAMP： 输入与各频率的振型对应的模态阻尼比 MP，DAMP 输入对应于某种材料的材料阻尼??。

与以上几种命令的输入对应的ANSYS计算的总阻尼阵[C]是： （5.1.2）

ANSYS计算阻尼矩阵的公式

其中m是结构中有阻尼的材料种类数，n是具有特有阻尼的单元类型数。前两项是用 与 定义的Rayleigh阻尼，第三项是与全结构的阻尼比 对应的阻尼阵，第四项是材料阻尼，最后一项是一些单元特有的单元阻尼阵。 3．粘性阻尼比

粘性阻尼表现为类似物体在粘性流体中运动时的阻力，与速度成正比。 （5.1.3） 粘性阻尼力

对单自由度系统，c就是粘性阻尼系数，对多自由度系统，就是阻尼矩阵[C]。[C]是定义结构阻尼特性的最基本形式，然而对粘性阻尼，很少有直接定义阻尼阵[C]的，阻尼比才是定义粘性阻尼最简捷的方法。在ANSYS中，既可以定义在结构坐标系下的全结构阻尼比（DMPRAT命令），也可以在模态坐标下对各个模态定义各自的模态阻尼比（MDAMP命令）。ANSYS最终计算的各模态相应的模态阻尼比是MDAMP定义的模态阻尼比与DMPRAT定义的全结构阻尼比的叠加。 DMPRAT与MDAMP都是只对响应谱分析、谐分析及使用模态叠加法的瞬态分析有效，它们所

对应的阻尼阵[C]是随频率不同而变化的阻尼阵。已知模态阻尼比 后，则对应的阻尼阵[C]用下式求出： （5.1.4）

与输入的模态阻尼比对应的阻尼矩阵

其中 是第i个振型向量， 是对应的模态频率。

值得注意的是上述公式只有理论意义，在振型叠加中是直接使用定义的振型阻尼比与全结构阻尼比，没有哪个程序会用公式（3）去反求出阻尼阵来。（也许某些程序里可以反求出阻尼阵来，但至少ANSYS没有这么做）。所以在做Full（完全）积分法的瞬态分析时，用阻尼比定义的阻尼都被程序忽略掉了，那么许多时候我们需要用一个全结构的阻尼比去做full法的瞬态分析计算时间，（如一些规范上规定某些结构可以用0.005~0.05的阻尼比做分析），该怎么办呢？这时候一个简单的办法是用 阻尼与 阻尼来逼近一个常数阻尼比。

图5.1 用ALPHD与BETAD来拟合常数阻尼比

选定 与 ，就可以用公式（1）计算出做输入用的ALPHD与BETAD值来。 4．材料阻尼

与其它几种阻尼不同的是，材料阻尼是在材料参数里面进行定义的（命令：MP,DAMP），材料阻尼又叫滞回阻尼，其最显著的特点是与结构响应频率无关。

图5.2 两种阻尼与频率的关系

许多文献上常把它写成复数刚度的形式： 。其中k是结构刚度， ， 称做材料阻尼系数（又叫结构阻尼系数）。

在单自由度情况，质量m做简谐振动时， （c是对应的粘性阻尼系数），因此得到 对应的阻尼比为： （5.1.5）

材料阻尼系数与粘性阻尼比的关系式

（在日本的结构减震规范中，用来定义阻尼的减衰系数就是此材料阻尼系数 。） 在ANSYS里，它是刚度矩阵的乘子，产生的阻尼阵是各材料对应刚度的加权和。 （5.1.6）

ANSYS计算材料阻尼对应阻尼矩阵的公式

很明显，它对应的阻尼阵[C]是可以对角化的，所以既能在full（完全）法瞬态分析中使用，也可以在振型叠加法分析中使用。上一小节里介绍了：ANSYS在做Full积分的瞬态分析时，用阻尼比定义的阻尼都被程序忽略掉，在许多时候，已知的是粘性阻尼的阻尼比，又要做full法的瞬态分析，那怎么办？此时一种办法是把粘性阻尼比换算为材料阻尼系数再用MP，DAMP输入。材料阻尼系数与粘性阻尼比的换算关系是： ，在单自由度情况下： （c是粘性阻尼系数）。

表5.1 常见材料的材料阻尼系数 纯铝 钢 铅 铸铁

0.00002~0.002 0.001~0.008 0.008~0.014 0.003~0.03

天然橡胶 硬橡胶 玻璃 混凝土

0.1~0.3 1.0 0.0006~0.002 0.01~0.06

以上材料来自：《结构振动分析》, C.F.比尔茨（作者对其使用不负任何责任）

金属的阻尼是比较低的，不知道这算不算是钢结构的一个缺点。一般来说高阻尼的金属其强度延性硬度均低。但是也有例外，如锰铜合金其强度硬度延性阻尼都高，但是相应价格也很高。 5．模态阻尼比的计算

当采用模态叠加法时，ANSYS对模态阻尼比与结构阻尼比是直接使用的，对其它阻尼则是计算多种阻尼产生的模态阻尼比来计算各模态的响应。在各种阻尼输入下，ANSYS程序计算出的第i个模态的总模态阻尼比是 （5.1.7）

ANSYS计算模态阻尼比的公式

其中前两项是 阻尼与 阻尼对应的模态阻尼比，第三项是输入的全结构阻尼比，第四项是输入的模态阻尼比，最后一项是M种材料的材料阻尼系数 产生的模态阻尼比。其中 是第j种材料对应的模态应变能，在日本减震规范中，就是采用此此应变能公式来计算结构阻尼比的。 ?注意：

如前所述，在做Full积分法的瞬态分析时，用阻尼比定义的阻尼都被ANSYS程序忽略掉了，所以同一个模型采用full法和模态叠加法的瞬态分析，ANSYS计算采用的阻尼可能不一样，造成结果也有差别。

以下是结构分析中常用的几种阻尼输入的ANSYS命令流演示。 1）用MP,damp来输入粘滞阻尼

DAMPRATO=0.025 ! 已知粘滞阻尼的阻尼比

LOSSMODM=2*DAMPRATO ! 粘滞阻尼的阻尼比乘以2是等价的材料阻尼系数(日 !本规范的“减衰系数”)

CRITFREQ=2.6 ! 此为粘性阻尼等效为材料阻尼时的换算频率

MP_BETAD=DAMPRATO/(acos(-1)*CRITFREQ) ! 粘滞阻尼与频率有关 /prep7

mp,damp,1,MP_BETAD !定义iscous damping，与频率有关 /solu

antype,modal modopt,lanb,1

! 要使模态计算考虑阻尼的影响，必须用材料阻尼，材料阻尼必须在求解前指定 ! mxpand,,,,yes,选项！阻尼比输入只在对求出的振型求反应再叠加中有用, ! ansys不会把阻尼比还原计算为阻尼阵[C]的 mxpand,1,,,yes ,,, sole,

2）用MP,Damp输入材料阻尼 DAMPRATO=0.025

LOSSMODM=2*DAMPRATO ! 材料阻尼系数，书上给的一般是LOSSMODM /prep7

mp,damp,1,DAMPRATO !常数，如果已知的是材料阻尼系数LOSSMODM，就要除以2 /solu

antype,modal ! 使用模态叠加法 modopt,lanb,1 ! important