ANSYS中的 阻尼 1
【拉布索思】谐响应中阻尼的设置及其工程应用

【拉布索思】谐响应中阻尼的设置及其工程应用讨论背景这里主要讨论的是ANSYS(Workbench)中模态叠加法(包括瞬态中用模态叠加)的阻尼比(包括α、β阻尼)设置问题,不考虑材料的阻尼比。
各阻尼比的解释这里会把阻尼都写成阻尼比的形式,因为阻尼比最直观,也可以在实验中得到。
那么,总阻尼比为:其中,是常值阻尼比,是作用于所有阶的模态,而且值是恒定的,在Constant Damping Ratio 设置;是第i阶模态的阻尼比,用来设定某些阶的阻尼比,要通过命令MDAMP设置(命令解释请看help,命令使用实例请看下面例子);,这是β阻尼相应的阻尼比,也是作用于所有阶的模态,但值随频率增大而线性增大。
HELP中说,在很多实际结构问题中,α会被忽略,所以上式中就没有α只有β,我是这么猜的,因为通常实际问题的频率都在几十到几千赫兹不等,那么就比较小,可以忽略。
反正,在Workbench界面中是没有α的设置项的,默认α=0,要设的话就要加命令,这里也不讨论α了。
那么β值就决定了这个阻尼比,而因为β值是恒定的,所以这个阻尼比会随频率增大而线性增大,就能起到抑制高频的作用。
AWB中有两种输入β阻尼的方法,一是直接输入β值(Direct Input),二是输入某个频率下的阻尼比(Dampingvs Frequency),系统就会根据来计算出β值,界面中的Beta Damping Measure就是。
通过两种输入方法设置了β值后,系统就会自动求出各阶的β阻尼比,个人喜欢用第二种方法,因为设置阻尼比更直观,结合使用另外两个阻尼比时会更方便;注意,三个阻尼比的效果是叠加的。
为什么要设置阻尼?首先,加阻尼对共振频率的影响很小,比如是固有频率乘以,所以阻尼对共振频率的影响可以忽略。
那我认为,阻尼的主要作用是压低共振处的幅值,使频响曲线变得更平缓。
而实际结构中一定有阻尼,所以分析中适当设置一定的阻尼会比较接近实际。
如果阻尼都取为0(默认值)的话,频呼曲线的峰值会相当大,理论上是无穷大的。
ansys中阻尼的设置

ANSYS/LS-DYNA中阻尼的设置ANSYS/LS-DYNA中阻尼的设置总结如下:EDDAMP, PART, LCID, VALDMP (用此命令流来定义阻尼)Defines mass weighted (Alpha) or stiffness weighted (Beta) damping for an explicit dynamics model.(默认有质量阻尼与刚度阻尼两种)图1 LS-DYNA中的阻尼选项PARTPART number [EDPART] identifying the group of elements to which damping should be applied. If PART = ALL (default), damping is applied to the entire model.(定义施加阻尼的Part 号,如果Part =ALL (or blank),整体质量阻尼将被应用于整个模型。
)LCIDLoad curve ID (previously defined with the EDCURVE command) identifying the damping coefficient versus time curve. If time-dependent damping is defined, an LCID is required.(用来指定相对于时间的质量阻尼)VALDMPConstant system damping coefficient or a scale factor applied to the curve defining damping coefficient versus time.(可以用来代替对时间的质量阻尼曲线,或者使用LCID命令中的比例系数)注意:Mass-weighted (Alpha) or stiffness-weighted (Beta) damping can be defined with the EDDAMP command. Generally, stiffness proportional or beta damping is effective foroscillatory motion at high frequencies. This type of damping is orthogonal to rigid body motion and so will not damp out rigid body motion. On the other hand, mass proportional or alpha damping is more effective for low frequencies and will damp out rigid body motion. The different possibilities are described below:1.Global DampingMass-weighted or Alpha damping (质量阻尼)When PART = (blank) or ALL (default), mass-weighted global damping can be defined in the following 2 ways. In this case, the same damping is applied for the entire structure.1.When the damping coefficient versus time curve (LCID) is specified usingthe EDCURVE command, VALDMP is ignored by LS-DYNA (although it iswritten in the LS-DYNA input file Jobname.K). The damping force appliedto each node in the model is given by f d= d(t)mv, where d(t) is thedamping coefficient as a function of time defined by the EDCURVEcommand, m is the mass, and v is the velocity.2.When the LCID is 0 or blank (default), a constant mass-weighted systemdamping coefficient can be specified using VALDMP.The constant and time-dependent damping, described above, cannot be defined simultaneously. The last defined global damping will overwrite any previously defined global damping.2.Damping defined for a PART(1)Mass-weighted or Alpha damping (质量阻尼)When both a valid PART number is specified and the damping coefficient versus time curve (LCID) is specified using the EDCURVE command, mass-weighted time-dependent damping will be defined for the particular PART. In this case, VALDMP will act as a scaling factor for the damping versus time curve (if VALDMP is not specified, it will default to 1). A valid PART number must be specified to define this type of damping. For example, use PART=1 (and not blank) when the entire model consists of only one PART. Issue the command repeatedly with different PART numbers in order to specify alpha damping for different PARTS.(2)Stiffness-weighted or Beta damping (刚度阻尼)When a valid PART number is specified with LCID= 0 or (blank) (default), a stiffness-weighted (Beta) constant damping coefficient for this particular PART canbe defined by VALDMP. The stiffness-weighted value corresponds to the percentage of damping in the high frequency domain. For example, 0.1 roughly corresponds to 10% damping in the high frequency domain. Recommended values range from 0.01 to 0.25. Values lower than 0.01 may have little effect. If a value larger than 0.25 is used, it may be necessary to lower the time step size significantly. Issue the command repeatedly with different PART numbers in order to specify beta damping for different PARTS. Time-dependent stiffness-weighted damping is not available in ANSYS LS-DYNA.The mass-weighted and stiffness-weighted damping, described above, cannot be defined simultaneously for a particular PART number. The last defined damping for the particular PART number will overwrite any previously defined mass-weighted or stiffness-weighted damping for this PART.In order to define the mass-weighted and stiffness-weighted damping simultaneously, you can use the MP,DAMP command to define stiffness-weighted (Beta) constant damping coefficient. However, do not use both of these commands together to define stiffness-weighted (Beta) constant damping coefficient for a particular PART. If you do, duplicate stiffness-weighted (Beta) constant damping coefficients for this PART will be written to the LS-DYNA input file Jobname.K. The last defined value will be used by LS-DYNA. Also, note that the MP,DAMP command is applied on the MAT number, and not on the PART number. Since a group of elements having the same MAT ID may belong to more than one PART (the opposite is not true), you need to issue the MP,DAMP command only once for this MAT ID and the stiffness-weighted (Beta) damping coefficients will be automatically defined for all the PART s with that MAT ID.Mass-weighted and stiffness-weighted damping can be defined simultaneously using the EDDAMP command only when mass-weighted damping (constant or time-dependent) is defined as global damping (EDDAMP, ALL, LCID, VALDMP) and stiffness-weighted damping is defined for all necessary PARTs (EDDAMP,PART, ,VALDMP).To remove defined global damping, reissue the EDDAMP, ALL command with LCID and VALDMP set to 0. To remove damping defined for a particular PART, reissue EDDAMP, PART, where PART is the PART number, with LCID and VALDMP set to 0. There is no default for the EDDAMP command, i.e., issuing the EDDAMP command with PART = LCID = VALDMP = 0 will result in an error. Stiffness-weighted damping defined by the MP,DAMP command can be deleted using MPDELE, DAMP, MAT.In an explicit dynamic small restart (EDSTART,2) or full restart analysis (EDSTART,3), you can only specify global alpha damping. This damping will overwrite any alpha damping input in the original analysis. If you do not input global alpha damping in the restart, the damping properties input in the original analysis will carry over to the restart.Damping specified by the EDDAMP command can be listed, along with other explicit dynamics specifications, by typing the command string EDSOLV$STAT into the ANSYS input window. Beta damping specified by the MP,DAMP command can be listed by MPLIST, MAT command.Menu PathsMain Menu>Preprocessor>Loads>Load Step Opts>Other>Change MatProps>DampingMain Menu>Preprocessor>Material Props>DampingMain Menu>Solution>Load Step Opts>Other>Change Mat Props>DampingK文件,Eg.1 $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ SYSTEM DAMPING $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $*DAMPING_PART_STIFFNESS2 1.000$K文件,Eg.2 $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ SYSTEM DAMPING $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $*DAMPING_GLOBAL00.5000E+02$。
ansys阻尼详解

antype,modal !模态分析
modopt,lanb,1
! important
mxpand,1,,,yes
lumpm,on
,,,,
sole
/solu
antype,harmic !谐分析
1.比例阻尼
最常用也是比较简单的阻尼大概是Rayleigh阻尼,又称为比例阻尼。它是多数实用动力分析的首选,对许多实际工程应用也是足够的。
在ANSYS里,它 就是 阻尼与 阻尼之和,分别用ALPHD与BETAD命令输入。已知结构总阻尼比是 ,则用两个频率点上 阻尼与 阻尼产生的等效阻尼比之和与其相等,就可以求出近似的 阻尼与 阻尼系数来用作输入:
ALPHAD: 输入 阻尼参数
BETAD: 输入 阻尼参数
DMPRAT: 输入全结构的阻尼比
MDAMP: 输入与各频率的振型对应的模态阻尼比
MP,DAMP 输入对应于某种材料的材料阻尼??。
与以上几种命令的输入对应的ANSYS计算的总阻尼阵[C]是:
(5.1.2)
ANSYS计算阻尼矩阵的公式
其中m是结构中有阻尼的材料种类数,n是具有特有阻尼的单元类型数。前两项是用 与 定义的Rayleigh阻尼,第三项是与全结构的阻尼比 对应的阻尼阵,第四项是材料阻尼,最后一项是一些单元特有的单元阻尼阵。
3.粘性阻尼比
粘性阻尼表现为类似物体在粘性流体中运动时的阻力,与速度成正比。
(5.1.3)
5.模态阻尼比的计算
当采用模态叠加法时,ANSYS对模态阻尼比与结构阻尼比是直接使用的,对其它阻尼则是计算多种阻尼产生的模态阻尼比来计算各模态的响应。在各种阻尼输入下,ANSYS程序计算出的第i个模态的总模态阻尼比是
ansys动力学瞬态分析详解

§3.1瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。
可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。
载荷和时间的相关性使得惯性力和阻尼作用比较重要。
如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。
瞬态动力学的基本运动方程是:其中:[M] =质量矩阵[C] =阻尼矩阵[K] =刚度矩阵{}=节点加速度向量{}=节点速度向量{u} =节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力([M]{})和阻尼力([C]{})的静力学平衡方程。
ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。
两个连续时间点间的时间增量称为积分时间步长(integration time step)。
§3.2学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。
可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。
例如,可以做以下预备工作:1.首先分析一个较简单模型。
创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。
在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。
通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。
同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。
<<高级技术分指南>>中将讲述子结构。
§3.3三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。
ANSYS中阻尼的定义和使用方式

ANSYS的轴承座结构分析教程一、实体模型的建立建立实体模型可以通过自上而下和自下而上两个途径:1、自上而下建模,首先要建立体(或面),对这些体或面按一定规则组合得到最终需要的形状。
2、自下而上建模,首先要建立关键点,由这些点建立线、由线连成面等一般建模原则是充分利用对称性,合理考虑细节。
根据题中的轴承座,由于轴承座具有对称性,只需建立轴承座的半个实体对称模型,在进行镜像操作即可。
采用自下而上的建模方法得到如下图1所示的三维实体模型:(1)生成长方体Main Menu:Preprocessor>Modeling->Create>Volumes->Block>By Dimensions输入x1=0,x2=60,y1=0,y2=20,z1=0,z2=60平移并旋转工作平面Utility Menu>WorkPlane>Offset WP by IncrementsX,Y,Z Offsets 输入45,25,15 点击ApplyXY,YZ,ZX Angles输入0,-90,0点击OK。
创建圆柱体Main Menu:Preprocessor>Create>Cylinder> Solid CylinderRadius输入15/2, Depth输入-30,点击OK。
拷贝生成另一个圆柱体Main Menu:Preprocessor>Copy>Volume拾取圆柱体,点击Apply, DZ输入30然后点击OK从长方体中减去两个圆柱体Main Menu:Preprocessor>Operate>Subtract Volumes首先拾取被减的长方体,点击Apply,然后拾取减去的两个圆柱体,点击OK。
使工作平面与总体笛卡尔坐标系一致Utility Menu>WorkPlane>Align WP with> Global Cartesian(2)创建支撑部分Main Menu: Preprocessor -> -Modeling-Create -> -Volumes-Block -> By 2 corners & Z在创建实体块的参数表中输入下列数值:WP X = 0WP Y = 20Width = 30Height = 35Depth = 15OKToolbar: SAVE_DB(3)偏移工作平面到轴瓦支架的前表面Utility Menu: WorkPlane -> Offset WP to -> Keypoints +1. 在刚刚创建的实体块的左上角拾取关键点2. OKToolbar: SAVE_DB(4)创建轴瓦支架的上部Main Menu: Preprocessor -> Modeling-Create -> Volumes-Cylinder -> Partial Cylinder +1). 在创建圆柱的参数表中输入下列参数:WP X = 0WP Y = 0Rad-1 = 0Theta-1 = 0Rad-2 = 30Theta-2 = 90Depth = -15或者在by dimensions 下建立圆柱体,输入相应的参数,其余圆柱的创建方式相同。
Ansys中的阻尼

Ansys中的阻尼ANSYS动力学分析中提供了各种的阻尼形式,这些阻尼在分析中是如何计算,并对分析有什么影响呢?本文将就此做一些说明何介绍.一.首先要清楚,在完全方法和模态叠加法中定义的阻尼是不同。
因为前者使用节点坐标,而后者使用总体坐标.1.在完全的模态分析、谐相应分析和瞬态分析中,振动方程为:阻尼矩阵为下面的各阻尼形式之和:α为常值质量阻尼(α阻尼)(ALPHAD命令)β为常值刚度阻尼(β阻尼)(BETA命令)ξ为常值阻尼比,f为当前的频率(DMPRAT命令)βj为第j种材料的常值刚度矩阵系数(MP,DAMP命令)[C]为单元阻尼矩阵(支持该形式阻尼的单元)where: [C] = structure damping matrixα = mass matrix multiplier (input on ALPHAD command)[M] = structure mass matrixβ = stiffness matrix multiplier (input on BETAD command)βc = variable stiffness matrix multiplier (see Equation 15–23)[K] = structure stiffness matrixNm = number of materials with DAMP or DMPR input= stiffness matrix multiplier for material j (input as DAMP on MP command)= constant (frequency-independent) stiffness matrix coefficient for material j (input as DMPR on MP command)Ω = circular excitation frequencyKj = portion of structure stiffness matrix based on material jNe = number of elements with specified dampingCk = element damping matrixCξ = frequency-dependent damping matrix (see Equation 15–21)2.对模态叠加方法进行的谐相应分析、瞬态分析何谱分析,动力学求解方程为:每个模态产生有效阻尼比ξid而不是创建阻尼矩阵α为常值质量阻尼β为常值刚度阻尼ξ为常值阻尼比ξmi为第i个模态的常值阻尼比ξj为第j个材料的阻尼系数Ejs为第j个材料的应变能,ANSYS由{f}T[K]{f}计算得到。
(完整版)ansys动力学瞬态分析详解
!Continue with normal transient solution procedures
...
参见<<ANSYS命令参考手册>>中关于命令ACEL、TIME、NSUBST、KBC、LSWRITE、DDELE和KBC的论述。
2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。在某些场合,动力学分析中是没必要包括非线性特性的。
3.掌握结构动力学特性。通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。同时,固有频率对计算正确的积分时间步长十分有用。
4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。<<高级技术分指南>>中将讲述子结构。
·允许在实体模型上施加的载荷。
完全法的主要缺点是它比其它方法开销大。
§
模态叠加法通过对模态分析得到的振型(特征值)乘上因子并求和来计算结构的响应。此法是ANSYS/Professional程序中唯一可用的瞬态动力学分析法。
模态叠加法的优点是:
·对于许多问题,它比缩减法或完全法更快开销更小;
·只要模态分析不采用PowerDynamics方法,通过 LVSCALE 命令将模态分析中施加的单元载荷引入到瞬态分析中;
TIME,...! Realistic time interval
DDELE,ALL,UY! Remove displacement constraints
KBC,0! Ramped loads (if appropriate)
!Continue with normal transient solution procedures
ansysworkbench常见材料设置
Ansys workbench常用材料属性1.isotropic secant coefficient of expansion 各向同性的热胀系数需要输入基准温度、热膨胀系数。
Propertea of Outline Row 4:8Ftoperiy Unit j 审3 _ Iso tropic Secant CceFfideni of Thermal闩篦1©0(pSFEbl2囱Coefficient of Erpiansion ReferenceTenizierahjre基准温度,默认22度热膨胀系数2.orthotropic secant coefficient of expansion 各向异性的热胀系数需要输入基准温度、三个方向的热膨胀系数。
3.isotropic instantaneous coefficient of expansion 各向同性的热胀系数(随温度变化)需要输入基准温度、热膨胀系数。
(随温度变化)Table of Properties Row 2: Isotropc Instartaneous Coefficient of Thermal Expansion1 ATemperature (C ;___________________ ? _____________________ Coefficient of Thermal E^ansion (C A -l) ▼ 2 20 2.3E-O5 3 100 3E-O5 42004E-O5*■ ।Chart of Roperties Row 2: Isotropic IrEtantaneous Coefficent of Thermal Expansion4. orthotropic instantaneous coefficient of expansion 各向异性的热胀系数需要输入基准温度、三个方向的热膨胀系数。
ANSYS中的阻尼(转载)
阻尼是动力分析的一大特点,也是动力分析中的一个易于引起困惑之处,而且由于它只是影响动力响应的衰减,出了错不容易觉察。
阻尼的本质和表现是相当复杂的,相应的模型也很多。
ANSYS提供了强大又丰富的阻尼输入,但也正以其强大和丰富使初学者容易发生迷惑这里介绍各种阻尼的数学模型在ANSYS中的实现,与在ANSYS中阻尼功能的使用。
1.比例阻尼最常用也是比较简单的阻尼大概是Rayleigh阻尼,又称为比例阻尼。
它是多数实用动力分析的首选,对许多实际工程应用也是足够的。
在ANSYS里,它就是阻尼与阻尼之和,分别用ALPHD与BETAD命令输入。
已知结构总阻尼比是,则用两个频率点上阻尼与阻尼产生的等效阻尼比之和与其相等,就可以求出近似的阻尼与阻尼系数来用作输入:(5.1.1)求比例阻尼系数的拟合公式用方程组(5.1.1)可以得到阻尼与阻尼系数值,然后用ALPHD与BETAD命令输入,这种阻尼输入既可以做full(完全)法的分析,也可以作减缩法与振型叠加法的分析,都是一样的有效。
但是尽管阻尼与阻尼概念简单明确,在使用中也要小心一些可能的误区。
首先,阻尼与质量有关,主要影响低阶振型,而阻尼与刚度有关,主要影响高阶振型;如果要做的是非线性瞬态分析,同时刚度变化很大时,那么使用阻尼很可能会造成收敛上的困难;一样的理由,有时在使用一些计算技巧时,比如行波效应分析的大质量法,加上了虚假的大人工质量,那么就不可以使用阻尼。
同样,在模型里加上了刚性连接时,也应该检查一下阻尼会不会造成一些虚假的计算结果。
2.阻尼阵的计算ANSYS中有多种办法可以输入阻尼特性。
先概括几个在结构分析中常用的输入阻尼的命令:ALPHAD:输入阻尼参数BETAD:输入阻尼参数DMPRA T:输入全结构的阻尼比MDAMP:输入与各频率的振型对应的模态阻尼比MP,DAMP 输入对应于某种材料的材料阻尼??。
与以上几种命令的输入对应的ANSYS计算的总阻尼阵[C]是:(5.1.2)ANSYS计算阻尼矩阵的公式其中m是结构中有阻尼的材料种类数,n是具有特有阻尼的单元类型数。
阻尼的问题
阻尼的问题wudingyi积分 49 帖子 37#12007-1-21 19:51我在论坛上看到可以用mp,damp 命令来输入材料的材料阻尼系数;但是在ansys 的帮助文件中是这样说明的: MP, Lab, MAT, C0, C1, C2, C3, C4 LabDAMP ― K matrix multiplier for damping.Note: If used in an explicit dynamic analysis, the value corresponds to the percentage of damping in the high frequency domain. For example, 0.1 roughly corresponds to 10% damping in the high frequency domain.DMPR ― Constant material damping coefficient .我对上面的理解是:输入材料的阻尼系数应该是用mp,dmpr 而mp,damp 到底输入的是什么,从上面我感觉和瑞利阻尼中的刚度阻尼系数重复,或者他是用来将材料阻尼系数等效为刚度阻尼系数来输入,如果是,该如何等效?请同行明示,谢谢!wms328积分 220 帖子 153#22007-1-23 20:21我也想问下阻尼的问题:一般做什么分析要考虑阻尼?是不是实际工程的计算分析都要考虑阻尼,从而在材料属性中输入阻尼系数! 希望有人能帮忙解开疑惑!谢谢mysuper积分 3 帖子 3#32007-1-24 12:38楼主说的材料属性中阻尼项我也很想知道在用combin14单元模拟阻尼单元时要求材料属性中mp,damp 一项不知道具体表示什么,在combin14单元实常数中定义了阻尼系数后,随意改变mp,damp 项似乎对结构的动力响应没有影响,不只是何原因wudingyi#42007-1-26 21:39你应该做的是完全瞬态分析,在ansys 中对其有两种阻尼输入模式,一积分 49 帖子 37 个是mp,damp 输入材料阻尼,另一个是单元阻尼,而对于其他一切的以阻尼比的形式输入的阻尼参数,程序一律忽略。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Ansys中的阻尼(2008-11-02 16:46:49)转载分类:工程技术标签:教育AnsysANSYS动力学分析中提供了各种的阻尼形式,这些阻尼在分析中是如何计算,并对分析有什么影响呢?本文将就此做一些说明何介绍.一.首先要清楚,在完全方法和模态叠加法中定义的阻尼是不同。
因为前者使用节点坐标,而后者使用总体坐标.1.在完全的模态分析、谐相应分析和瞬态分析中,振动方程为:阻尼矩阵为下面的各阻尼形式之和:α为常值质量阻尼(α阻尼)(ALPHAD命令)β为常值刚度阻尼(β阻尼)(BETA命令)ξ为常值阻尼比,f为当前的频率(DMPRAT命令)βj为第j种材料的常值刚度矩阵系数(MP,DAMP命令)[C]为单元阻尼矩阵(支持该形式阻尼的单元)where: [C] = structure damping matrixα = mass matrix multiplier (input on ALPHAD command)[M] = structure mass matrixβ = stiffness matrix multiplier (input on BETAD command)βc = variable stiffness matrix multiplier (see Equation 15–23)[K] = structure stiffness matrixNm = number of materials with DAMP or DMPR input= stiffness matrix multiplier for material j (input as DAMP on MP command)= constant (frequency-independent) stiffness matrix coefficient for material j (input as DMPR on MP command)Ω = circular excitation frequencyKj = portion of structure stiffness matrix based on material jNe = number of elements with specified dampingCk = element damping matrixCξ = frequency-dependent damping matrix (see Equation 15–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)来计算系统的有效阻尼比:也就是说,一个有效的常值材料阻尼将被用于所有模态的计算。
五.下面的表格列出了不同分析类型中可以用到的阻尼。
/document/200707/article1154_7.htm§3.8.3阻尼大多数系统中存在阻尼,而且在动力学分析中应当指定阻尼。
在ANSYS程序可以指定五种形式的阻尼:·Aplha和Beta阻尼(Rayleigh阻尼)·和材料相关的阻尼·恒定的阻尼比·振型阻尼·单元阻尼在ANSYS/Professional程序中只有恒定阻尼比和振型阻尼可用。
可以在模型中指定多种形式的阻尼,程序按所指定的阻尼之和形成阻尼矩阵[C]。
下表列出了在不同结构分析中可用的阻尼类型。
不同分析类型可用的阻尼N/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理论参考手册>>。
在Transient 分析的Full 方法中,可用的阻尼就只有Rayleigh 阻尼(Alphad, betad)和材料相关的阻尼(MP,DAMP),以及单元阻尼。
而和材料相关的阻尼在Full 方法中其实还是和材料相关的Beta 阻尼。
现在有一结构,包括主体钢结构和柔索两部分组成,要求风荷载下的动力响应。
但两部分结构的阻尼比和频率相差许多,分别取1%和0.5%,遇到问题如下:1.有什么方法可是两部分的阻尼取值不一样?能够分别定义吗?我曾通过对柔索ω1=0.04Hz处ξ1=0.5%和钢结构ω1=0.4Hz处ξ1=1%来定义α,β,但是所得结果在柔索和钢结构的二阶频率处ω2处的阻尼与实际相差实在太远。
又通过和材料相关的阻尼来分别定义两部分的Beta 值,仍无法达到理想的ξ 曲线。
不知有谁处理过类似的问题,请指教2. 在ANSYS 手册中有:“Beta阻尼和材料阻尼在非线性分析中会导致不理想的结果。
这两种阻尼要和刚度矩阵相乘,而刚度矩阵在非线性分析中是不断变化的。
由此引起的阻尼变化有时会和物理结构的实际阻尼变化相反。
”而现在的情况是:若Beta 阻尼取为零,结构可正常计算;若Beta 阻尼取到0.002 以上,则算了几步马上就不收敛了?两部分结构的阻尼比和频率相差许多,我的意见是用材料阻尼输入,即MP,DAMP命令。
此命令虽然是材料阻尼,实际上也可以用来输入粘滞阻尼等,只是需要先作一个换算。
各种阻尼输入方式都是可以互相换算的,所以,您说:“和材料相关的阻尼在Full 方法中其实还是和材料相关的Beta 阻尼。
”我认为是不确的,尽管形式可能相同,但实际上物理本质是不一样的,所以互相要有一定的换算才可以代入。
Rayleigh 阻尼(Alphad, betad)及全结构阻尼比还有振型阻尼比等都是对正个结构而言,所以要对不同部分输入不同的阻尼只有用MP,DAMP,及单元阻尼了。
无论是粘滞阻尼还是材料阻尼都可以用MP,DAPM输入,但是相应的换算关系是不一样的,不能直接代入。
您的阻尼建议您应该先考虑一下其本质,再找到最合适的阻尼本构关系,再找到相应的数值,在换算为MP,DAMP的输入数值。
模态分析,谱分析,瞬态分析。
根据ANSYS的帮助文件,针对分析的不同,阻尼的输入方式也是有分别的,现归纳如下:1、模态分析:只谈有阻尼模态分析。