ANSYS动力学分析.docx
- 文档编号:5683230
- 上传时间:2022-12-31
- 格式:DOCX
- 页数:49
- 大小:55.41KB
ANSYS动力学分析.docx
《ANSYS动力学分析.docx》由会员分享,可在线阅读,更多相关《ANSYS动力学分析.docx(49页珍藏版)》请在冰豆网上搜索。
ANSYS动力学分析
标准化管理处编码[BBX968T-XBB8968-NNJ668-MM9N]
ANSYS动力学分析
第5章动力学分析
结构动力学研究的是结构在随时间变化载荷下的响应问题,它与静力分析的主要区别是动力分析需要考虑惯性力以及运动阻力的影响。
动力分析主要包括以下5个部分:
模态分析:
用于计算结构的固有频率和模态。
谐波分析(谐响应分析):
用于确定结构在随时间正弦变化的载荷作用下的响应。
瞬态动力分析:
用于计算结构在随时间任意变化的载荷作用下的响应,并且可涉及上述提到的静力分析中所有的非线性性质。
谱分析:
是模态分析的应用拓广,用于计算由于响应谱或PSD输入(随机振动)引起的应力和应变。
显式动力分析:
ANSYS/LS-DYNA可用于计算高度非线性动力学和复杂的接触问题。
本章重点介绍前三种。
【本章重点】
区分各种动力学问题;
各种动力学问题ANSYS分析步骤与特点。
.1动力学分析的过程与步骤
模态分析与谐波分析两者密切相关,求解简谐力作用下的响应时要用到结构的模态和振型。
瞬态动力分析可以通过施加载荷步模拟各种何载,进而求解结构响应。
三者具体分析过程与步骤有明显区别。
.1.1模态分析
1.模态分析应用
用模态分析可以确定一个结构的固有频率利振型,固有频率和振型是承受动态载荷结构设计中的重要参数。
如果要进行模态叠加法谐响应分析或瞬态动力学分析,固有频率和振型也是必要的。
可以对有预应力的结构进行模态分析,例如旋转的涡轮叶片。
另一个有用的分析功能是循环对称结构模态分析,该功能允许通过仅对循环对称结构的一部分进行建模,而分析产生整个结构的振型。
ANSYS产品家族的模态分析是线性分析,任何非线性特性,如塑性和接触(间隙)单元,即使定义也将被忽略。
可选的模态提取方法有6种,即BlockLanczos(默认)、Subspace、PowerDynamics、Reduced、Unsymmetric、Damped及QRDamped,后两种方法允许结构中包含阻尼。
2.模态分析的步骤
模态分析过程由4个主要步骤组成,即建模、加载和求解、扩展模态,以及查看结果和后处理。
(1)建模。
指定项目名和分析标题,然后用前处理器PREP7定义单元类型、单元实常数、材料性质及几何模型。
必须指定杨氏模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量),材料性质可以是线性或非线性、各向同性或正交各向异性,以及恒定或与温度有关的,非线性特性将被忽略。
(2)加载及求解。
在这个步骤中要定义分析类型和分析选项,施加载荷,指定加载阶段选项,并进行固频率的有限元求解。
在得到初始解后,应对模态进行扩展以供查看。
ANSYS提供的用于模态分析的选项如下。
NewAnalysis[ANTYPE]:
选择新的分析类型。
AnalysisType:
Modal[ANTYPE]:
指定分析类型为模态分析。
ModeExtractionMethod[MODOPT]:
可选模态提取方法如下。
BlockLanczosmethod(默认):
分块的兰索斯法,它适用于大型对称特征值求解问题,比子空间法具有更快的收敛速度;Subspacemethod:
子空间法,适用于大型对称特征值问题;PowerDynamicsmethod:
适用于非常大的模型(100000个自由度以上)及求解结构的前几阶模态,以了解结构如何响应的情形。
该方法采用集中质量阵(LUMPM,ON);Reduced(Householder)method:
使用减缩的系统矩阵求解,速度快。
但由于减缩质量矩阵识近似矩阵,所以相应精度较低;Unsymmetricmethod:
用于系统矩阵为非对称矩阵的问题,例如流体—结构相同作用;Dampedmethod:
用于阻尼不可忽略的问题;QRDampedmethod:
采用减缩的阻尼阵计算复杂阻尼问题,所以比Dampedmethod方法有更快的计算速度和更好的计算效率。
NumberofModestoExtract[MODOPT]:
除Reduced方法外的所有模态提取方法都必须设置该选项。
NumberofModestoExpand[MXPAND]:
仅在采用Reduced、Unsymmetric和Damped方法时要求设置该选项。
但如果需要得到单元的求解结果,则不论采用何种模态提取方法需要得到单元的求解结果,则不论采用何种模态提取方法都需选择Calculateelemresults复选框。
MassMatrixFormulation[LUMPML]:
使用该选项可以选定采用默认的质量矩阵形成方式(和单元类型有关)或集中质量阵近似方式,建议在大多数情况下应采用默认形成方式。
但对有些包含薄膜结构的问题,如细长梁或非常薄的壳,采用集中质量矩阵近似经常产生较好的结果。
另外,采用集中质量阵求解时间短,需要内存少。
PrestressEffectsCalculation[PSTRES]:
选用该选项可以计算有预应力结构的模态。
默认的分析过程不包括预应力,即结构是处于无应力状态的。
完成模态分析选项(ModalAnalysisOption)对话框中的选择后,单击
按钮。
一个相应于指定的模态提取方法的对话框将会出项。
对话框中给出如下选择域的组合。
FREQB,FREQE:
指定模态提取的频率范围,大多数情况无需设置。
PRMODE:
要输出的减缩模态数,只对Reduced方法有效。
Nrmkey:
关于振型归一化的设置,可选择相对于质量矩阵[M]或单位矩阵[I]进行归一化处理。
RIGID:
设置提取对已知有刚体运动结构进行子空间迭代分析时的零频模态,只对Subspace和PowerDynamics法有效。
SUBOPT:
指定多种子空间迭代选项,只对Subspace和PowerDynamics方法有效。
CEkey:
指定处理约束方程的方法,只对BlockLanczos方法有效。
(3)定义自由度。
使用Reduced模态提取法时要求定义自由度:
GUI:
MainMenu>Solution>MasterDOFs>-userSelected-Define。
命令:
M。
(4)在模型上加载荷。
在典型的模态分析中惟一有效的“载荷”是零位移约束,如果在某个DOF处指定了一个非零位移约束,则以零位移约束替代该DOF处的设置。
可以施加除位移约束之外的其他载荷,但它们将被忽略。
在未加约束的方向上,程序将解算刚体运动(零频)及高频(非零频)自由体模态。
载荷可以加在实体模型(点,线和面)上或加在有限元模型(点和单元)上。
(5)指定载荷步选项。
模态分析中可用的载荷步选项见表5-1。
阻尼只在用Damped模态提取法时有效,在其他模态提取法中将被忽略。
如果包含阻尼,且采用Damped模态提取法,则计算特征值时复数解。
表5-1模态分析中可用的载荷步选项
选项
命令
GUI路径
Alpha(质量)阻尼
ALPHAD
MainMenu>Solution>LoadStepOpts>Time/Frequenc>Damping
Beta(刚度)阻尼
BETAD
MainMenu>Solution>LoadStepOpts>Time/Frequenc>Damping
恒定阻尼比
DMPRAT
MainMenu>Solution>LoadStepOpts>Time/Frequenc>Damping
材料阻尼比
MP,DAMP
MainMenu>Solution>Other>ChangeMatProps>Polynomial
单元阻尼比
R
MainMenu>Solution>LoadStepOpts>Other>RealConstants>Add/Edit/Delete
PrintedOutput
OUTPR
MainMenu>Solution>LoadStepOpts>OutputCtrls>SoluPrintout
(6)开始求解计算。
GUI:
MainMenu>Solution>-Solve-CurrentLS。
命令:
SOLVE。
求解器的输出内容主要为写到输出文件及振型文件中的固有频率,也可以包含减缩的振型和参与因子表,这取决于设置的分析选项的输出控制。
由于振型现在尚未写到数据库或结果文件中,因此还不能对结果进行后处理。
如果采用Subspace模态提取法,则输出内容中可能包括警告:
STURMnumber=nshouldbem。
其中n和m为整数,表示某阶模态被漏掉或第m阶和第n阶模态的频率相同,而要求输出的只有第m阶模态。
如果采用Dmaped模态提取方法,求得的特征值和特征向量将是复数解。
特征值的虚部代表固有频率,实部为系统稳定性的量度。
(7)退出SOLUTION。
GUI:
MainMenu>Finish。
命令:
FINISH。
3.扩展模态
从严格意义上来说,扩展意味着将减缩解扩展到完整的DOF集上;而缩减解常用主DOF表达。
在模态分析中扩展指将振型写入结果文件,即扩展模态适用于Reduced模态提取方法得到的减缩振型和使用其他模态提取方法得到的完整振型。
因此如果需要在后处理器中查看振型,必须先将振型写入结果文件。
模态扩展要求振型、、及文件(如果采用Reduced方法)必须存在且数据库中必须包含和结算模态时所用模型相同的分析模型。
扩展模态的操作步骤如下。
(8)进入ANSYS求解器,可采用如下命令。
GUI:
MainMenu>Solution。
命令:
/SOLU。
在扩展处理前必须退出求解并重新进入(/SOLU)。
(9)激活扩展处理及相关选项(如表5-2所示)。
表5-2扩展处理选项
选项
命令
GUI路径
ExpansionPassOn/Off
EXPASS
MainMenu>Solution>AnalysisType>ExpansionPass
No.ofModestoExpand
MXPAND
MainMenu>Solution>LoadStepOpts>ExpansionPass>SingleExpand>ExpandModes
Freq.RangeforExpansion
MXPAND
MainMenu>Solution>LoadStepOpts>ExpansionPass>SingleExpand>ExpandModes
StressCalc.On/Off
MXPAND
MainMenu>Solution>LoadStepOpts>ExpansionPass>SingleExpand>ExpandModes
ExpansionPassOn/Off[EXPASS]:
选择ON(打开)。
NumberofModestoExpand[MXPAND,NMODE]:
指定要扩展的模态数。
记住,只有经过扩展的模态可在后处理中查看。
默认为不进行模态扩展。
FrequencyRangeforExpansion[MXPAND,,FREQB,FREQE]:
这是另一种控制要扩展模态数的方法。
如果指定一个频率范围,那么只有该频率范围内的模态会被扩展。
StressCalculationsOn/Off[MXPAND,,,,Elcalc]:
是否计算应力,默认为不计算。
模态分析中的应力并不代表结构中的实际应力,而只是给出一个各阶模态之间相对应力分布的概念。
(10)指定载荷步选项,模态扩展处理中惟一有效的选项是输出控制。
GUI:
MainMenu>Solution>LoadStep>OutputCtrls>DB/ResultsFile。
命令:
OUTRES。
(11)开始扩展处理,扩展处理的输出包括已扩展的振型,而且还可以要求包含各阶模态相对应的应力分布。
GUI:
MainMenu>Solution>CurrentLS。
命令:
SOLVE。
(12)如须扩展另外的模态(如不同频率范围的模态)重复步骤
(2)~(4),每次扩展处理的结果文件中保存为单步的载荷步。
(13)退山SOLUTION,可以在后处理器中查看结果。
GUI:
MainMenu>Finish。
命令:
FINISH。
4.查看结果和后处理
模态分析的结果(即扩展模态处理的结果)写入结构分析文件中,其中包括固有频率、己扩展的振型和相对应力和力分布(如果要求输出)可以在普通后处理器(/POST1)中查看模态分析结果。
查看结果数据包括读入合适子步的结果数据。
每阶模态在结果文件中保存为一个单独的子步。
如扩展了6阶模态,结果文件中将有6个子步组成的一个载荷步。
GUI:
MainMenu>GeneralPostproc>ReadResults>ByLoadStep>Substep。
命令:
SET和SBSTEP。
GUI:
MainMenu>GeneralPostproc>PlotResults>DeformedShape。
命令:
PLDISP。
.1.2谐响应分析
5.谐响应分析应用
谐响应分析是用于确定线性结构在承受随时间按正弦(简谐)规律变化的载荷时的稳态响应的一种技术。
分析的目的是计算结构在几种频率下的响应并得到一些响应值(通常是位移)对频率的曲线,从这些曲线上可找到“峰值”响应并进一步查看峰值频率对应的应力。
这种分析技术只计算结构的稳态受迫振动,发生在激励开始时的瞬态振动不在谐响应分析中考虑。
作为一种线性分析,该分析忽略任何即使已定义的非线性特性,如塑性和接触(间隙)单元。
但可以包含非对称矩阵,如分析在流体-结构相互作用问题。
谐响应分析也可用于分析有预应力的结构,如小提琴的弦(假定简谐应力比预加的拉伸应力小得多)。
6.3种求解方法
谐响应分析可以采用如下3种方法。
(14)Full方法(完全)。
该方法采用完整的系统矩阵计算谐响应(没有矩阵减缩),矩阵可以是对称或非对称的,其优点如下。
容易使用,因为不必关心如何选择主自由度和振型。
使用完整矩阵,因此不涉及质量矩阵的近似。
允许有非对称矩阵,这种矩阵在声学或轴承问题中很典型。
用单一处理过程计算出所有的位移和应力。
允许施加各种类刑的载荷,如节点力、外加的(非零)约束利单元载荷(压力和温度)。
允许采用实体模型上所加的载荷。
该方法的缺点是预应力选项不可用,并且采用Frontal方程求解器时通常比其他方法运行时间长。
但是采用JCG求解器或JCCG求解器时,该方法的效率很高。
(15)Reduced方法。
该方法通常采用主自由度和减缩矩阵来压缩问题的规模,计算主自由度处的位移后,解可以被扩展到初始的完整DOF集上,其优点如下。
在采用Frontal求解器时比Full方法更快。
可以考虑预应力效果。
该方法的缺点如下:
初始解只计算出主自由度的位移。
要得到完整的位移,应力和力的解则需执行被称为扩展处理的进一步处理,扩展处理在某些分析应用中是可选操作。
不能施加单元载荷(压力和温度等)。
所有载荷必须施加在则户定义的自由度上,限制了采用实体模型上所加的载荷。
(16)ModeSuperposition方法(模态叠加)。
该方法通过对模态分析得到的振型(特征向量)乘上因子并求和计算出结构的响应,其优点如下:
对于许多问题,比Reduced或Full方法更快。
在模态分析中施加的载荷可以通过LVSCALE命令用于谐响应分析中。
可以使解按结构的固有频率聚集,可产生更平滑且更精确的响应曲线图。
可以包含预应力效果。
允许考虑振型阻尼(阻尼系数为频率的函数)。
该方法的缺点如下:
不能施加非零位移。
在模态分析中使用PowerDynamics方法时,初始条件中不能有预加的载荷。
谐响应的3种方法有如下共同局限性。
所有载荷必须随时间按正弦规律变化。
所有载荷必须有相同的频率。
不允许有非线性特性。
不计算瞬态效应。
7.步骤
使用Full方法进行谐响应分析的过程的主要步骤为建模、加载并求解,以及查看结果及后处理。
8.建模
在该步骤中需指定文件名和分析标题,然后用PREP7来定义单元类型、单元实常数、材料特性及几何模型,需记住的要点如下:
(17)只有线性行为是有效的,如果有非线性单元,则按线性单元处理。
(18)必须指定杨氏模量EX(或某种形式的刚度)和密度DENS(或某种形式的质量)。
材料特性可为线性、各向同性或各向异性,以及恒定的或和温度相关的,忽略非线性材料特性。
9.加载并求解
在该步骤中定义分析类型和选项、加载、指定载荷步选项并开始有限元求解。
需要注意的是,峰值响应分析发生在力的频率和结构的固有频率相等时。
在得到谐响应分析解之前,应首先执行模态分析,以确定结构的固有频率。
(19)进入ANSYS求解器。
GUI:
MainMenu>Solution。
命令:
/SOLU。
(20)定义分析类型和分析选项,ANSYS提供的用于谐响应分析的选项见表5-3。
表5-3用于谐响应分析的选项
选项
命令
GUI路径
NewAnalysis
ANTYPE
MainMenu>Solution>AnalysisType>NewAnalysis
AnalysisType:
HarmonicResponse
ANTYPE
MainMenu>Solution>AnalysisType>NewAnalysis>Harmonic
SolutionMethod
HROPT
MainMenu>Solution>AnalysisType>AnalysisOptions
SolutionListingFormat
HROUT
MainMenu>Solution>AnalysisType>AnalysisOptions
MassMatrixFormulation
LUMPM
MainMenu>Solution>AnalysisType>AnalysisOptions
EquationSolver
EQSLV
MainMenu>Solution>AnalysisType>AnalysisOptions
NewAnalysis[ANTYPE]选择新分析,在谐响应分析中Restart不可用。
如果需要施加另外的简谐载荷,可以另进行一次新分析。
AnalysisType:
HarmomcResponse[ANTYPE]:
选择分析类型为HarmomcResponse(谐响应分析)。
SolutionMethod[HROPT]选择Full、Reduced或ModeSuperposition求解方法之一。
SolutionListingFormat[HROUT]:
确定在输出文件中谐响应分析的位移解如何列出,可选方式有realandimaginary(实部和虚部)(默认)和amplitudesandphaseangles(幅值和相位角)。
MassMatrixFormulation[LUMPM]:
指定采用默认的质量阵形成方式(取决于单元类型)或使用集中质量阵近似。
.
EquationSolver[EQSLV]:
可选求解器有Frontal(默认),SparseDirect(SPARSE)、JacobiConjugateGradient(JCG),以及IncompleteCholeskyConjugateGradient(ICCG)。
对大多数结构模型,建议采用Frontal或SPARSE求解器。
(21)在模型上加载。
根据定义,谐响应分析假定所施加的所有载荷随时间按简谐(正弦)规律变化。
指定一个完整的简谐载荷需输入3个数据,即Amplitude(振幅),phaseangle(相位角)和forcingfrequencyrange(强制频率范围)。
(22)指定载荷步选项,谐响应分析可用的选项见表5-4。
表5-4谐响应分析可用的选项
选项
命令
GUI途径
普通选项(GeneralOptions)
NumberofHarmonicSolution
NSUBST
MainMenu>Solution>LoadStepOpts>Time/Frequenc>FreqandSubsteps
SteppeorRampedLoads
KBC
MainMenu>Solution>LoadStepOpts>Time/Frequenc>Time>TimeSteporFreqandSubsteps
动力学选项(DynamicsOptions)
ForcingFrequencyRange
HARFRQ
MainMenu>Solution>LoadStepOpts>Time/Frequenc>FreqandSubsteps
Damping
ALPHAD,BETAD,DMPRAT
MainMenu>Solution>LoadStepOpts>Time/Frequenc>Damping
输出控制选项(OutputControlOptions)
PrintedOutput
OUTPR
MainMenu>Solution>LoadStepOpts>OutputCtrls>SoluPrintout
DatabaseandResultsFileOutput
OUTRES
MainMenu>Solution>LoadStepOpts>OutputCtrls>DB/ResultsFile
ExtrapolationofResults
ERESX
MainMenu>Solution>LoadStepOpts>OutputCtrls>IntegrationPt
普通选项如下。
NumberofHarmonicSolutions[NSUBST]:
请求计算任何数目的谐响应解,解(或子步)将均布于指定的频率范围内[HARFQR]。
例如,如果在30Hz~40Hz范围内要求出10个解,则计算在频率31Hz~40Hz处的响应,而不计算其他频率处。
SteppedorRampedLoads[KBC]:
载荷以Stepped或Ramped方式变化,默认为Ramped,即载荷的幅值随各子步逐渐增长。
如果用命令[KBC,1)设置了Stepped载荷,则在频率范围内的所有子步载荷将保持恒定的幅值。
动力学选项如下。
ForcingFrequencyRange[HARFRQ]:
在谐响应分析中必须指定强制频率范围(以周/单位时间为单位),然后指定在此频率范围内要计算处的解数。
Damping:
必须指定某种形式的阻尼,如Alpha(质量)阻尼[ALPHAD]、Beta(刚度)阻尼[BETAD]或恒定阻尼比[DMPRAT]否则在共振处的响应将无限大。
(23)开始求解。
GUI:
MainMenu>Solution>Solve>CurrentLS。
命令:
SOLVE。
(24)如果有另外的载荷和频率范围(即另外的载荷步),重复步骤(3)~(5)。
如果要做时间历程后处理(POST26),则一个载荷步和另一个载荷步的频率范同间不能存在重叠。
(25)退出SOLUTION。
GUI:
关闭Sol
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 动力学 分析