瞬态动力学分析汇总文档格式.docx

瞬态动力学分析汇总文档格式.docx

《瞬态动力学分析汇总文档格式.docx》由会员分享，可在线阅读，更多相关《瞬态动力学分析汇总文档格式.docx（21页珍藏版）》请在冰点文库上搜索。

NewAnalysisType

ANTYPE

MainMenu>

Solution>

AnalysisType>

NewAnalysis

Transient（瞬态）

SolutionMethod

TRNOPT

MainMenu>

AnalysisOptions

Full/Reduced/Modesuperpos’n

MassMatrixFormulation

LUMPM

“薄膜”结构使用如：

细长梁、薄壳

LargeDeformationEffects

NLGEOM

大变形效应如：

细长杆、金属成型问题

StressStiffeningEffects

SSTIF

小变形分析应力显著改变结构刚度；

大变形问题中为了加强收敛

Newton-RaphsonOption

NROPT

求解期间切线矩阵被刷新的频度，仅在存在非线性时使用

EquationSolver

EQSLV

选用适合求解器进行求解

ANSYS在进行瞬态动力学分析中可以采用三种方法，即Full（完全）法、Reduced（缩减）法和ModeSuperposition（模态叠加）法。

ANSYS提供了各种分析类型和分析选项，使用不同方法ANSYS软件会自动配置相应选择项目，常用的分析类型和分析选项如表16-1所示。

在瞬态分析中，时间总是计算的跟踪参数，在整个时间历程中，同样载荷也是时间的函数，有两种变化方式：

Ramped：

如图16-1（a）所示，载荷按照线性渐变方式变化。

Stepped：

如图16-1（b）所示，载荷按照解体突变方式变化。

表16-2常用的分析类型和分析选项

Full（完全）法

采用完整的系统矩阵计算瞬态响应。

功能最强大，允许包括非线性的类型。

Reduce（缩减）法

采用主自由度及缩减矩阵计算瞬态响应。

需定义字主自由度，计算速度快。

ModeSuperposition

（模态叠加）法

从模态分析中得到中模态乘参与因子并求和，计算瞬态相应。

速度快。

图16-1载荷增加方式渐变与突变

依据载荷变化方式可以将整个时间历程划分成多个载荷步（LoadStep），每个载荷步代表载荷发生一次突变或一次渐变阶段。

在每个载荷步时间内，载荷增量又可以划分多个子步（Substep），在子步载荷增量的条件下程序进行迭代计算即Iteriation，经过多个子步的求解实现一个载荷步的求解，进而求出多个载荷步的求解实现整个载荷时间历程的求解。

利用ANSYS进行瞬态动力学分析时可以在实体模型或有限元模型上施加下列载荷：

约束（Displacement）、集中力（Force）、力矩（Moment）、面载荷（Pressure）、体载荷（Temperature、Fluence）、惯性力（Gravity，Spinning，ect.）。

在ANSYS中，进行多载荷步加载的基本方法常用有三种：

（1）连续多载荷步加载法。

（2）定义载荷步文件批加载法。

（3）定义表载荷加载法。

第2节瞬态动力学分析实例

案例1——自由度弹簧质量系统瞬态分析

表16-3随时间变化载荷

TIME（s）

2

4

6.8

9.5

FORCE（N）

-500

-150

图16-2弹簧质量系统/载荷图

问题

如图16-2所示，单自由度的弹簧质量系统，试对质点M在变力FORCE作用下的瞬态分析，并绘出位移瞬态响应曲线。

条件

弹簧刚度50000N/m，长度0.2m,质量大小为150KG，质点受力如表16-3所示，忽略阻尼。

解题过程

以弹簧上部端点为坐标原点，建立直角坐标系。

制定分析方案。

分析类型为瞬态动力学分析；

模型类型为线、质点模型，由于结构简单可以直接创建节点和单元，弹簧部分选用Combin14单元，质量块部分简化为质点选用MASS21单元，边界条件为上端施加固定全约束，据图16-2中Force—Time图的特点采用外力以Ramped线性渐变方式加载，连续多载荷步加载方法。

瞬态分析的求解方法采用Reduce（缩减）法。

1．ANSYS分析开始准备工作

（1）清空数据库并开始一个新的分析

选取UtilityMenu>

File>

Clear&

StartNew，弹出ClearsdatabaseandStartNew对话框，单击OK按钮，弹出Verify对话框，单击OK按钮完成清空数据库。

（2）指定新的工作文件名

指定工作文件名。

ChangeJobname，弹出ChangeJobname对话框，在EnterNewJobname项输入工作文件名，本例中输入的工作文件名为“Transientexample1”，单击OK按钮完成工作文件名的定义。

（3）指定新的标题

指定分析标题。

ChangeTitle，弹出ChangeTitle对话框，在EnterNewTitle项输入标题名，本例中输入“exercise1”为标题名，然后单击OK按钮完成分析标题的定义。

（4）重新刷新图形窗口

Plot>

Replot，定义的信息显示在图形窗口中。

（5）定义结构分析

运行主菜单MainMenu>

Preferences，出现偏好设置对话框，选中赋值分析模块为Structure结构分析模块，单击OK按钮完成分析模块的定义。

2．定义单元及材料

（1）新建单元类型

Preprocessor>

ElementType>

Add/Edit/Delete命令，弹出ElementTypes对话框，单击Add按钮新建单元类型，弹出LibraryofElementTypes对话框，先选择单元大类为StructuralMass，接着选择3DMass21（MASS21单元），单击Apply按钮，选择单元大类为Combination，接着选择Spring-damper14（COMBIN14单元）单击OK按钮，完成单元类型TYPE1和单元类型TYPE2的定义，如图16-3所示。

图16-3定义单元类型

选择TYPE1，单击Option按钮进入单元设置选项，弹出MASS21Elementtypeoptions对话框，在Rotaryinertiaoptions（K3）栏中更改选项为2-Dw/orotiner（2D无转动惯量），在单击OK关闭窗口，如图16-4所示，再选择TYPE2，单击Option按钮进入单元设置选项如图16-5所示，在DOFselectionfor2D+3D（K3）栏中更改选项为2-Dlongitudinal（定义弹簧单元自由度），单击OK按钮。

3．定义实常数

（1）新建实例常量

RealConstantsAdd/Edit/Delete命令，弹出实常数定义对话框，如图16-6所示。

16-4设置单元配置项

图16-5设置单元配置项

（2）输入实常数

单击Add按钮进入选择单元类型对话框如图16-7所示，选择TYPE1单击OK按钮。

进入实例常量输入对话框，如图16-8，输入质点质量150，单击OK按钮，出现如16-9左图所示。

如上步骤继续单击Add按钮，如图16-9、图16-10所示，完成实常数Set2的输入，在Springconstant（K）项，输入50000，在单击OK按钮。

（3）关闭实常数对话框

回到实例常量对话框，此时显示出新建编号为Set1和Set2的实例常量，单击Close按钮完成输入，如图16-11所示。

图16-6定义实常数对话框图16-7选择定义实常数的单元类型对话框

图

16-8实例常量Set2输入对话框

图16-9选择定义实常数的单元类型对话框

图16-10实例常量Set2输入对话框图16-11定义实常数对话框

4．建立有限元模型（由于本案例模型较为简单，可以直接创建节点和单元形成有限元模型。

）

（1）创建节点

Modeling>

Create>

Nodes>

InActiveCS命令，在对话框中分别在节点号NODE栏输入1，节点坐标X，Y，Z栏输入0，0，0，单击Apply按钮完成一次输入，如图16-12所示。

在对话框中分别在节点号NODE栏输入2，节点坐标X，Y，Z栏输入0，-0.2，0，单击OK按钮完成节点输入，单击OK按钮，如图16-13所示。

16-12当前坐标系下创建节点对话框

图16-13当前坐标系下创建节点对话框

（2）创建单元

Elements>

ElemAttributes命令，在[TYPE]单元类型序号栏中更改选项为1MASS21，在[REAL]实常数序号栏中更改选项为1，单击OK按钮，如图16-14所示。

AutoNumbered>

ThruNodes命令，如图16-15所示，输入节点号2，单击OK按钮，创建质点单元。

ElemAttributes命令，在[TYPE]单元类型序号栏中更改选项为2COMBIN14，在[REAL]实常数序号栏中更改选项为2，单击OK按钮，如图16-16所示。

ThruNodes命令，如图16-17所示，输入节点号1，2（中间用“，”间隔）或用鼠标按顺序点选1、2节点，单击OK按钮，创建弹簧单元。

图16-14单元属性对话框图16-15以节点创建单元对话框

图16-16单元属性对话图16-17以节点创建单元对话框

运行功能菜单UtilityMenu>

PlotCtrls>

Numbering命令，弹出PlotNumberingControl菜单在NODENodenumbers栏中鼠标点击选项为On，在Elem/attribnumbering栏中更改选项为Elementnumbers单击OK按钮，如图16-18所示。

5．定义分析类型和分析选项并加载

（1）定义分析类型及选项

AnalysisType>

NewAnalysis命令，弹出NewAnalysis对话框，选择Transient瞬态分析选项，单击OK按钮，弹出TransientAnalysis对话框，选择Reduced选项，单击OK按钮，如图16-19、图16-20所示。

图16-18显示编号对话框

图16-19定义分析类型对话框

图16-20定义瞬态分析方法对话框

NewAnalysis>

AnalysisOptions命令，弹出ReducedTransientAnalysis对话框，在[TRNOPT]Dampingeffects栏中更改选项为Ignore（忽略阻尼），单击OK按钮，如图16-21所示。

图16-21缩减法瞬态分析选项对话框

（2）定义主自由度

MasterDOFs>

UserSelected>

Define命令，弹出DefineMasterDOFs对话框，输入节点号2，单击OK按钮，如图16-22所示，又弹出DefineMasterDOFs对话框，如图16-23所示，在Lab1栏中更改选项为UY，单击OK按钮完成质点单元之自由度的定义即UY为其主自由度。

图16-22定义主自由度图16-23定义主自由度对话框

（3）定义载荷步选项

LoadStepOpts>

Time/Frequenc>

Time—TimeStep命令，弹出TimeandTimeStepOptions对话框，在[DELTIM]Timestepsize栏中输入0.02，在[KBC]stepSteppedorrampedb.c.栏鼠标点选Ramped，在[DELTIM]Minimumtimestepsize栏中输入0.002，在[DELTIM]Maximumtimestepsize栏中输入0.2，单击OK按钮，如图16-24所示。

16-24定义载荷步对话框

（4）定义输出选项

OutputCtrls>

SoluPrintout命令，弹出SolutionPrintoutControls对话框，在ItemItemforprintoutcontroled栏中更改选项为Allitems，在FREQPrintfrequency栏鼠标点选Everysubstep，单击OK按钮，如图16-25所示。

图16-25定义打印输出控制对话框

DB/ResultsFile命令，弹出ControlsforDatabaseandResultsFileWriting对话框，在ItemItemtobecontroled栏中更改选项为Allitems，在FREQFilewritefrequency栏鼠标点选Everysubstep，单击OK按钮。

其它选项默认如图16-26所示。

图16-26定义打印输出控制对话框

（5）施加约束

DefineLoads>

Apply>

Structural>

Displacement>

OnNodes命令，出现拾取菜单，如图16.27拾取节点1，单击OK按钮确定，出现约束定义对话框如图16.28，选择AllDOF约束所有自由度，在DisplacementValue选项输入0，再单击OK按钮完成约束定义。

图16-27拾取对话框图16-28施加节点约束对话框

单击ANSYSToolbar中的SAVE__DB进行存盘。

（6）对第一载荷步施加载荷

Force/Moment>

OnNodes命令，出现拾取菜单如图16-29所示，选择如图16-30中节点2，单击OK按钮确定，出现载荷定义对话框如图16-31所示，在Lab栏中更改选项为FY，在VALUE栏输入0，单击OK按钮确定。

Time—TimeStep命令，出现拾取菜单如图16-32所示，在[TIME]栏输入0，其它默认，单击OK按钮确定。

图16-29定义节点载荷路径图16-31定义节点载荷方向和值对话框

图16-30选择节点对话框图16-32定义载荷步对话框

（7）对第一载荷步进行求解

CurrentLS命令，出现菜单中单击OK按钮确定，计算机开始进行求解，求解完成后出现“Solutionisdone”提示表示求解完成，单击Close按钮完成第一载荷步求解，如图16-33所示。

图16-33求解对话框

（8）对第二载荷步到第五载荷步施加载荷

重复以上（6）～（7）步骤，不同之处仅在图16-31VALUE一栏输入载荷数值第二载荷步到第五载荷步输入表16-3数值分别为-500、-150、-150、0。

在如图16-32所示，在[TIME]一栏输入载荷步时间值为第二载荷步到第五载荷步输入表16-3时间值分别为2、4、6.8、9.5。

分别对第二载荷步到第五载荷步进行求解（ansys缺省的荷载处理为替代方式，可用fcum,sfcum命令修改，可选择三种方式：

替代（repl）、累加（add）、忽略（ignore）。

当采用缺省时，对于同一自由度处的荷载，后面施加的荷载替代了前面施加的荷载（或覆盖），故本例前一步骤载荷被覆盖而不用另外删除，否则应删除上步载荷。

）第五载荷步求解完成后出现“Solutionisdone”提示表示求解完成，单击Close按钮完成所有载荷步的求解，如图16-33所示。

Finish命令，结束整个求解过程。

（9）定义扩展并求解

ExpansionPass命令，弹出ExpansionPass对话框，在EXPASS栏鼠标点选ON，单击OK按钮确定，如图16-34所示。

图16-34扩展求解对话框

ExpansionPass>

SingleExpand>

ByTime/Freq命令，弹出ExpandSingleSolutionbyTime/Frequency对话框，在TIMFRQ栏输入7.8（预扩展解时间值），单击OK按钮确定，如图16-35所示。

Finish命令，求解。

在ANSYSToolbar上单击SAVE保存。

图16-35定义扩展求解时间对话框

5．观察计算结果

TimeHistPostpro命令，进入时间历程后处理器，弹出如图16-36TimeHistroyVariables－\file.rdsp窗口（若为TimeHistroyVariables－\file.rst窗口可以运行如16-36图file>

打开TimeHistroyVariable－\file.rdsp文件），单击“＋”按钮，弹出如图16-36所示AddTime—HistroyVariables窗口，单击选中NodalSolution>

Y－Componentofdisplacement菜单，单击OK按钮，弹出如图16-37所示，输入节点号2（或用鼠标在窗口点选），单击OK按钮，返回TimeHistroyVariables－\file.rdsp窗口，界面显示出选择的节点2处的时间历程位移值如图16-38所示，单击Graph按钮，ANSYS图形窗口显示节点2处的时间历程位移图示，如图16-38所示。

图16-36定义时间历程变量对话框图16-37选取节点对话框

单击Listdata按钮，弹出如图16-39所示的FRVARCommand窗口显示出选择的节点2处的时间历程位移值列表。

图16.38节点2处的时间历程位移值图示

图16.39节点2处的时间历程位移值列表

6．保存并退出ANSYS

在ANSYSToolbar上单击QUIT，以选定方式保存。

第3节本章小结

1．瞬态分析中时间的意义

瞬态分析中，与模态、谐响应分析不同，时间是真正意义的时间概念，而且加载时仅输入载荷步的最后时间点的值，如本例图16-32所示，在[TIME]栏输入0，2，4，6.8，9.5，ANSYS会自动默认载荷步开始时间为上一步的结束时间。

2．瞬态分析中初始条件的处理

瞬态分析中往往初始条件不为零，如：

时间t=0时，位移u0，v0，a0不为零，处理如此问题可以分两个载荷步去解决。

载荷步1：

施加初始条件可以以静载荷步开始采用较小的时间间隔，如0.001秒，施加初始条件用MainMenu>

Displacement命令或运行菜单MainMenu>

Inertia>

AngularAccel>

Global命令或用MainMenu>

InitialCondit'

n>

Define命令施加，若物体初始条件静止可以施加强制约束保持物体静止，采用突变方式2个子步，并关闭瞬态效应并求解。

载荷步2：

打开瞬态效应，释放物体的强制约束，规定终止时间，连续进行瞬态分析。