Ansys中的阻尼

Ansys中的阻尼ANSYS动力学分析中提供了各种的阻尼形式，这些阻尼在分析中是如何计算，并对分析有什么影响呢？本文将就此做一些说明何介绍.一．首先要清楚，在完全方法和模态叠加法中定义的阻尼是不同。因为前者使用节点坐标，而后者使用总体坐标.1．在完全的模态分析、谐相应分析和瞬态分析中，振动方程为：阻尼矩阵为下面的各阻尼形式之和：α为常值质量阻尼（α阻尼）（ALPHAD命令）β为常值刚度阻尼（β阻尼）（BETA命令）ξ为常值阻尼比，f为当前的频率（DMPRAT命令）βj为第j种材料的常值刚度矩阵系数（MP,DAMP命令）[C]为单元阻尼矩阵（支持该形式阻尼的单元）where:[C]=structuredampingmatrixα=massmatrixmultiplier(inputonALPHADcommand)[M]=structuremassmatrixβ=stiffnessmatrixmultiplier(inputonBETADcommand)βc=variablestiffnessmatrixmultiplier(seeEquation15–23)[K]=structurestiffnessmatrixNm=numberofmaterialswithDAMPorDMPRinput=stiffnessmatrixmultiplierformaterialj(inputasDAMPonMPcommand)=constant(frequency-independent)stiffnessmatrixcoefficientformaterialj(inputasDMPRonMPcommand)Ω=circularexcitationfrequencyKj=portionofstructurestiffnessmatrixbasedonmaterialjNe=numberofelementswithspecifieddampingCk=elementdampingmatrixCξ=frequency-dependentdampingmatrix(seeEquation15–21)2．对模态叠加方法进行的谐相应分析、瞬态分析何谱分析，动力学求解方程为：每个模态产生有效阻尼比ξid而不是创建阻尼矩阵α为常值质量阻尼β为常值刚度阻尼ξ为常值阻尼比ξmi为第i个模态的常值阻尼比ξj为第j个材料的阻尼系数Ejs为第j个材料的应变能，ANSYS由{f}T[K]{f}计算得到。二．对谱分析，阻尼仅仅包含在模态组合里，而在计算模态系数的时候并没有考虑。当使用模态叠加法时，材料阻尼被添加到扩展的模态中，因此，用户必须在进行模态分析之前，就包括材料阻尼（MP,DAMP）并进行单元应力的计算（MXPAND命令）。三．模态叠加法支持使用QR阻尼，但是用户必须知道尽管是模态组合方法，阻尼在模态分析中已包含了，所以应该使用上面的完全阻尼矩阵[C]来计算阻尼。如果使用QR阻尼的的模态提取方法（MOPT,QRDAMP），并且在前处理或模态分析中指定了任何形式的阻尼，那么ANSYS将在进行模态叠加时忽略阻尼。四．了解MP,DAMP在不同的情况下有不同的作用非常重要。在完全分析中，材料阻尼代表了该材料的一个刚度矩阵乘子，与粘性阻尼（与频率成线性关系，但针对所有的材料）类似。因此，在这种情况下，对单自由度结构，材料阻尼值等于ξ/πf或c/k。如果有多种材料，那么，阻尼矩阵就是简单得将材料的βj乘以相应的材料的刚度矩阵：但是，在使用模态叠加法时材料阻尼值类似于结构阻尼（与频率无关），也即材料阻尼值会等于单自由度体系的ξ。如果有多种材料，使用模态应变能方法（MSE）来计算系统的有效阻尼比：也就是说，一个有效的常值材料阻尼将被用于所有模态的计算。五．下面的表格列出了不同分析类型中可以用到的阻尼。§3.8.3阻尼大多数系统中存在阻尼，而且在动力学分析中应当指定阻尼。在ANSYS程序可以指定五种形式的阻尼：·Aplha和Beta阻尼（Rayleigh阻尼）·和材料相关的阻尼·恒定的阻尼比·振型阻尼·单元阻尼在ANSYS/Professional程序中只有恒定阻尼比和振型阻尼可用。可以在模型中指定多种形式的阻尼，程序按所指定的阻尼之和形成阻尼矩阵[C]。下表列出了在不同结构分析中可用的阻尼类型。不同分析类型可用的阻尼分析类型α,β阻尼[ALPHAD,BETAD]材料相关阻尼[MP,DAMP]恒定阻尼比[DMPRAT]振型阻尼[MDAMP]单元阻尼[3](COMBIN7等)静力学分析N/AN/AN/AN/AN/A模态分析无阻尼No[5]No[5]No[5]NoNo有阻尼YesYesNoNoYes谐响应分析完全法YesYesYesNoYes缩减法YesYesYesNoYes模态叠加法Yes[6]Yes[4，6]Yes[7]Yes[7]Yes[6]瞬态分析完全法YesYesNoNoYes缩减法YesYesNoNoYes模态叠加法Yes[6]Yes[4，6]Yes[7]Yes[7]Yes[6]谱分析SPRS,MPRS[2]Yes[1]YesYesYesNoDDAM[2]Yes[1]YesYesYesNoPSD[2]YesNoYesYesNo屈曲分析N/AN/AN/AN/AN/A子结构YesYesNoNoYesN/A表示不能使用[1]表示只可用β阻尼，不可用α阻尼[2]表示阻尼只用于模态合并，不用于计算模态系数[3]表示包括超单元阻尼矩阵[4]表示如果经模态扩展转换成了振型阻尼[5]表示如果指定了，程序会计算出一个用于随后的谱分析的有效阻尼比[6]表示如果使用QR阻尼模态提取方法[MODOPT,QRDAMP]，在前处理或模态分析过程中指定任何阻尼，但ANSYS在执行模态叠加分析时将忽略任何阻尼。[7]如果你使用QR阻尼模态提取方法[MODOPT,QRDAMP]，DMPART和MDAMP不能使用。1.Alpha阻尼和Beta阻尼Alpha阻尼和Beta阻尼用于定义瑞利（Rayleigh）阻尼常数α和β。阻尼矩阵是在用这些常数乘以质量矩阵[M]和刚度矩阵[K]后计算出来的。命令ALPHAD和BETAD分别用于确定瑞利（Rayleigh）阻尼常数α和β。通常α和β的值不是直接得到的，而是用振型阻尼比计算出来的。是某个振型i的实际阻尼和临界阻尼之比。如果是模态i的固有角频率，则α和β满足下列关系：在许多实际问题中，Alpha阻尼（或称质量阻尼）可以忽略（α=0）。这种情形下，可以由已知的和计算出β：由于在一个载荷步中只能输入一个β值，因此应该选取该载荷步中最主要的被激活频率来计算β值。为了确定对应给定阻尼比ξ的α和β值，通常假定α和β之和在某个频率范围内近似为恒定值（见图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中的实现，与在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个振型向量，是对应的模态频率。值得注意的是上述公式只有理论意义，在振型叠