第二章 谐响应分析文档格式.docx

1、节点力、外加的（非零）位移、单元载荷（压力和温度）。允许在实体模型上定义载荷。完全法的 缺点 是：预应力选项不可用。当采用Frontal方程求解器时这种方法通常比其它方法都开销大。但在采用JCG求解器或ICCG求解器时，完全法的效率很高。2.3.2缩减法缩减法通过采用主自由度和缩减矩阵来压缩问题的规模。主自由度处的位移被计算出来后，解可以被扩展到初始的完整DOF集上（参见“模态分析”中的“矩阵缩减技术”部分关于缩减技术的细节）。这种方法的 优点 是：在采用Frontal求解器时比完全法更快且开销小；可以考虑预应力效果。缩减法的 缺点 是：初始解只计算出主自由度处的位移。要得到完整的位移、应力和

2、力的解则需执行扩展过程（扩展过程在某些分析应用中是可选操作）；不能施加单元载荷（压力、温度等等）所有载荷必须施加在用户定义的主自由度上（限制了采用实体模型上所加载荷）。2.3.3模态叠加法模态叠加法通过对模态分析得到的振型（特征向量）乘上因子并求和来计算出结构的响应。它的 优点 是：对于许多问题，此法比Reduced或完全法更快且开销小；模态分析中施加的载荷可以通过LVSCALE命令用于谐响应分析中；可以使解按结构的固有频率聚集，便可得到更平滑、更精确的响应曲线图；可以包含预应力效果；允许考虑振型阻尼（阻尼系数为频率的函数）。模态叠加法的缺点是：不能施加非零位移；在模态分析中使用PowerDy

3、namics法时，初始条件中不能有预加的载荷。2.3.4三种方法共同的局限性谐响应分析的三种方法存在共同的 局限性 ：所有载荷必须随时间按正弦规律变化；所有载荷必须有相同的频率；不允许有非线性特性；不计算瞬态效应。可以通过进行瞬态动力学分析来克服这些限制，这时应将简谐载荷表示为有时间历程的载荷函数。“瞬态动力学分析”中描述了瞬态动力学分析的过程。2.4完全法谐响应分析2.4.1完全法谐响应分析过程下面首先将描述的是如何用完全法来进行谐响应分析，然后列出用缩减法和模态叠加法时有差别的步骤。完全法谐响应分析过程由三个主要步骤组成：1.建模；2.加载并求解；3.观察结果。2.4.2建模建模阶段需要指

4、定文件名和标题，然后进入前处理器PREP7定义单元类型、单元实常数、材料特性以及几何模型。该过程与其它分析基本相同，但必须注意下面两个要点：在谐响应分析中，只有线性行为是有效的。如果有非线性单元，它们将按线性单元处理。例如，如果分析中包含接触单元，则它们的刚度取初始状态值并在计算过程中不再发生变化。必须指定杨氏模量EX（或某种形式的刚度）和密度DENS（或某种形式的质量）。材料特性可以是线性的、各向同性的或各向异性的、恒温的或和温度相关的。非线性材料特将被忽略。2.4.3加载并求解该过程将指定分析类型及其相关选项、定义模型载荷以及指定载荷步选项，然后开始有限元求解，下面详细介绍每个步骤。注意：

5、峰值响应发生在力的频率和结构的固有频率相等时。在得到谐响应分析解之前，应该首先做一下模态分析（如“模态分析”中所述）以确定结构的固有频率。2.4.3.1进入ANSYS求解器命令：/SOLUGUI路径：MainMenuSolution2.4.3.2定义分析类型和分析选项ANSYS提供下表所示的用于谐响应分析的选项：表1分析类型和分析选项选项命令GUI 途径NewAnalysisANTYPEMain MenuSolution-Analysis Type-New AnalysisAnalysisType:HarmonicResponseAnalysis Type-New AnalysisHarmon

6、ic ResponseSolutionMethodHROPTAnalysisOptionsSolutionListingFormatHROUTMassMatrixFormulationLUMPMEquationSolverEQSLV下面将对各个选项进行详细解释。 选项：NewAnalysisANTYPE 选NewAnalysis（新分析）。在谐响应分析中Restart不可用；如果需要施加其他简谐载荷，可以另进行一次新分析。HarmonicResponseANTYPE选分析类型为HarmonicResponse（谐响应分析）。SolutionMethodHROPT选择下列求解方法中的一种：完全法

7、缩减法模态叠加法SolutionListingFormatHROUT此选项确定在输出文件Jobname.Out中谐响应分析的位移解如何列出。可选的方式有“realandimaginary（实部与虚部）”（缺省）形式和“amplitudesandphaseangles（幅值与相位角）”形式。MassMatrixFormulationLUMPM此选项用于指定是采用缺省的分布质量矩阵（取决于单元类型）还是集中质量矩阵。建议在大多数应用中采用缺省的分布质量矩阵。但对于某些包含“薄膜”结构的问题，如细长梁或非常薄的壳，集中质量近似矩阵经常能产生较好的结果。另外，集中质量近似矩阵可以减少运行时间并降低内存

8、要求。在设置完HarmonicAnalysisOption对话框的MassMatrixFormulation项后，单击OK则弹出一个名为HarmonicAnalysis的对话框，用于选择方程求解器。EquationSolverEQSLV可选的求解器有：Frontal求解器（缺省）、JacobiConjugateGradient（JCG）求解器及IncompleteCholeskyConjugateGradient（ICCG）求解器。建议对大多数结构模型用Frontal求解器。2.4.3.3在模型上施加载荷1. 谐响应分析的载荷描述方式根据定义，谐响应分析假定所施加的所有载荷随时间按简谐（正弦）

9、规律变化。指定一个完整的简谐载荷需要输入三条信息：Amplitude（幅值）、phaseangle（相位角）和forcingfrequencyrange（强制频率范围）（见图2）。图2实部/虚部分量和振幅/相位角间的关系Amplitude（幅值）指载荷的最大值，可用表2中所示的命令指定。phaseangle（相位角）指载荷滞后（或领先）于参考时间的量度。在复平面上（见图2），相位角是以实轴为起始的角度。当同时要定义多个相互间存在相位差的简谐载荷时，必须分别指定相位角。例如，图3显示的不平衡旋转天线将在它的四个支撑点处产生垂直的异步载荷。相位角不能直接指定，而是应该用加载命令的VALUE和VAL

10、UE2域来指定有相位角载荷的实部和虚部。压力、分布载荷和体载荷只能指定0相位角（即不能定义载荷的虚部）。图2显示了计算实部和虚部的公式。Forcingfrequencyrange（强制频率范围）指简谐载荷（以周/单位时间为单位）的频率范围。在后面描述载荷步选项命令HARFRQ时将提到它。图3非平衡旋转天线在它的四个支撑点处产生异步垂直载荷注：谐响应分析不能计算频率不同的多个强制载荷同时作用时的响应。这种情形的实例是两个具有不同转速的机器同时运转时的情形。但在POST1中可以对两种载荷状况进行叠加以得到总体响应。2. 载荷类型 表2概括了谐响应分析中可施加的载荷。除惯性载荷外，可以在实体模型（由

11、关键点，线，面组成）或有限元模型（由节点和单元组成）上定义载荷。关于实体模型载荷与有限元载荷比较参见的第二章。在分析过程中，可以施加、删除载荷或对载荷进行操作或列表。表2谐响应分析中可施加的载荷载荷形式范畴Displacement（UX,UY,UZ,ROTX,ROTY,ROTZ）约束D MainMenuSolutionLoads-ApplyStructuralDisplacementForce,Moment（FX,FY,FZ,MX,MY,MZ）力F Loads-ApplyStructural-Force/MomentPressure（PRES）面载荷SF Structural-Pressure

12、Temperature（TEMP）Fluence（FLUE）体载荷BF Structural-TemperatureGravity,Spinning,等惯性载荷Structural-Other3. 利用命令施加载荷 表3列出了在谐响应分析中所有可以用来进行载荷操作的命令。表3谐响应分析中的加载命令载荷类型实体/FE 模型图素施加删除列表操作施加设置Displace-ment实体模型关键点DKDKDELE DKLIST DTRAN 线DLDLDELE DLLIST 面DADADELEDALISTDTRANFE节点DDDELEDLISTDSCALEDSYMDCUMForceFKFKDELEFKLI

13、STFTRANFFDELEFLISTFSCALEFCUMPressureSFLSFLDELESFLLISTSFTRANSFGRADSFASFADELESFALISTSFSFDELESFLISTSFSCALESFCUM单元SFESFEDELESFELISTSFBEAMSFFUNTemperature，fluenceBFKBFKDELEBFKLISTBFKTRAN有限元模型BFBFDELEBFLISTBFSCALEBFCUMBFEBFEDELEBFELISTInertiaACELOMEGADOMEGACGLOCCGOMGADCGOM4. 利用GUI 施加载荷所有的载荷操作（除列表外；见下）都可通过

14、一系列等效的下拉菜单实现。从求解器菜单中可选择菜单载荷操作（施加、删除等），然后选择载荷类型（位移、力等），然后再选择施加载荷的对象（关键点、线、节点等）。例如，要在一条线上施加位移载荷，可用如下GUI途径：GUI：-Loads-Apply-Structural-DisplacementOnLines5. 载荷列表要列出现有载荷，可用如下GUI途径：UtilityMenuListLoadsLoadtype2.4.3.4指定载荷步选项谐响应分析可用的选项如下表所示：表4载荷步选项普通选项（GeneralOptions）NumberofHarmonicSolutionsNSUBST-LoadSte

15、pOpts-Time/FrequencyFreq&SubstepsSteppedorRampedLoadsKBCTime/FrequencTime&TimeStep/Freq&Substeps动力学选项（DynamicOptions）ForcingFrequencyRangeHARFREQTime/FrequencyFreq&DampingALPHADBETADDMPARTDamping输出控制选项（OutputControlOptions）PrintedOutputOUTPRSolution-LoadStepOpts-OutputCtrlsSoluPrintoutDatabaseandRes

16、ultFileOutputOUTRESDB/ResultsFileExtrapolationofResultsERESXIntegrationPt1. 普通选项 NumberofHarmonicSolutionsNSUBST可用此选项请求计算任何数目的谐响应解。解（或子步）将均布在指定的频率范围内HARFRQ。例如，如果在3040HZ范围内要求出10个解，程序将计算出在频率为31，32，33，39和40Hz处的响应，而不计算频率范围低端处（30Hz）的响应。 SteppedorRampedLoadsKBC载荷可以是Stepped或Ramped方式变化的。缺省时方式为Ramped。即载荷的幅值随

17、各子步逐渐增长。而如果用命令KBC,1设置了Stepped载荷，则在频率范围内的所有子步中载荷将保持恒定的幅值。2. 动力学选项 ForcingFrequencyRangeHARFRQ在谐响应分析中必须指定强制频率范围（以周/单位时间为单位）。然后要指定在此频率范围内要计算出的解的数目（参见“普通选项”）。 Damping必须指定某种形式的阻尼，否则在共振频率处的响应将无限大。命令ALPHAD和BETAD指定的是和频率相关的阻尼系数，而DMPRAT指定的是对所有频率为恒定值的阻尼比。参见“瞬态动力学分析”中关于阻尼的更详细的描述。Alpha（质量）阻尼ALPHADBeta（刚度）阻尼BETAD

18、恒定阻尼比DMPRAT在直接积分谐响应分析（用完全法或缩减法）中如果没有指定阻尼，程序将缺省采用零阻尼。3. 输出控制选项Printed Output OUTPR此选项用于指定输出文件Jobname.OUT中要包含的结果数据。 DatabaseandResultsFileOutputOUTRES此选项用于控制结果文件Jobname.RST中包含的数据。 ExtrapolationofResultsERESX用于设置采用将结果复制到节点处方式而非缺省的外插方式得到单元积分点结果。2.4.3.5将数据库保存到一个命名的备份文件中。将数据库保存到一个命名的备份文件中，这样在重新进入ANSYS程序后用

19、RESUME命令便可恢复以前建的模型。SAVEFileSaveas2.4.3.6开始求解。SOLVE-Solve-CurrentLs2.4.3.7重复执行加载求解得到其他载荷步重复以上过程，计算其他载荷和频率范围（即另外的载荷步）的结果。如果希望进行时间历程后处理（在POST26中），载荷步之间的频率范围不能存在重叠。处理多步载荷还有一种方法：将载荷步保存到文件中，然后用一个宏进行一次性求解。该方法在DefineVariables可用FORCE命令选择合力，合力的静力分量、阻尼分量、惯性分量。2.绘制变量频率或其它变量的关系曲线。然后用PLCPLX指定是用幅值/相位角方式还是实部/虚部方式表示

20、解。PLVAR,PLCPLXGraphVariableSettingsGraph3.对变量值进行列表。如果只要求列出极值，可用EXTREM命令。PRVAR,EXTREM,PRCPLXMainMenuListVariables/ListExtremesList在POST26中还可以使用许多其它后处理功能，如在变量间进行数学运算（复数运算），将变量值传递给数组元素，将数组元素值传递给变量等。细节参见ANSYS基本分析过程指南第六章。通过观察整个模型中关建点处的时间历程结果，可以得到用于进一步POST1后处理的频率值。2.4.4.4使用POST11.读入所需谐响应分析结果。可以用SET命令来读入结果，但它将读入实部或者虚部，不能同时将二者都读入。结果的实际大小由实部和虚部的SRSS值（平方和取平方根）给出（见图2），在POST26中可得到模型中的指定点处的真实结果。2.显示结构的变形形状、应力、应变等的等值线或者矢量图PLVECT。要得到数据的列表表格，请用PRNSOL,