ansys动力学瞬态分析详解(20220301011126).pdf
《ansys动力学瞬态分析详解(20220301011126).pdf》由会员分享,可在线阅读,更多相关《ansys动力学瞬态分析详解(20220301011126).pdf(50页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、3.1 瞬态动力学分析的定义瞬态动力学分析(亦称时间历程分析)是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。载荷和时间的相关性使得惯性力和阻尼作用比较重要。如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。瞬态动力学的基本运动方程是:其中:M=质量矩阵C=阻尼矩阵K=刚度矩阵=节点加速度向量=节点速度向量u=节点位移向量在任意给定的时间,这些方程可看作是一系列考虑了惯性力(M)和阻尼力(C)的静力学平衡方程。ANSYS 程序使用 Newmark 时间积分方法在
2、离散的时间点上求解这些方程。两个连续时间点间的时间增量称为积分时间步长(integration time step)。3.2 学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂,因为按“工程”时间计算,瞬态动力学分析通常要占用更多的计算机资源和更多的人力。可以先做一些预备工作以理解问题的物理意义,从而节省大量资源。例如,可以做以下预备工作:1.首先分析一个较简单模型。创建梁、质量体和弹簧组成的模型,以最小的代价深入的理解动力学认识,简单模型更有利于全面了解所有的动力学响应所需要的。2.如果分析包括非线性特性,建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。在某些场合,动力学分析
3、中是没必要包括非线性特性的。3.掌握结构动力学特性。通过做模态分析计算结构的固有频率和振型,了解这些模态被激活时结构的响应状态。同时,固有频率对计算正确的积分时间步长十分有用。4.对于非线性问题,考虑将模型的线性部分子结构化以降低分析代价。中将讲述子结构。3.3 三种求解方法瞬态动力学分析可采用三种方法:完全(Full)法、缩减(Reduced)法及模态叠加法。ANSYS/Professional 产品中只允许用模态叠加法。在研究如何实现这些方法之前,让我们先探讨一下各种方法的优点和缺点。3.3.1完全法完全法采用完整的系统矩阵计算瞬态响应(没有矩阵缩减)。它是三种方法中功能最强的,允许包括各
4、类非线性特性(塑性、大变形、大应变等)。注如果并不想包括任何非线性,应当考虑使用另外两种方法中的一种。这是因为完全法是三种方法中开销最大的一种。完全法的优点 是:容易使用,不必关心选择主自由度或振型。允许各种类型的非线性特性。采用完整矩阵,不涉及质量矩阵近似。在一次分析就能得到所有的位移和应力。允许施加所有类型的载荷:节点力、外加的(非零)位移(不建议采用)和单元载荷(压力和温度),还允许通过TABLE数组参数指定表边界条件。允许在实体模型上施加的载荷。完全法的主要缺点是它比其它方法开销大。3.3.2模态叠加法模态叠加法通过对模态分析得到的振型(特征值)乘上因子并求和来计算结构的响应。此法是A
5、NSYS/Professional 程序中唯一可用的瞬态动力学分析法。模态叠加法的优点 是:对于许多问题,它比缩减法或完全法更快开销更小;只要模态分析不采用PowerDynamics方法,通过 LVSCALE 命令将模态分析中施加的单元载荷引入到瞬态分析中;允许考虑模态阻尼(阻尼比作为振型号的函数)。模态叠加法的缺点 是:整个瞬态分析过程中时间步长必须保持恒定,不允许采用自动时间步长;唯一允许的非线性是简单的点点接触(间隙条件);不能施加强制位移(非零)位移。3.3.3缩减法缩减法通过采用主自由度及缩减矩阵压缩问题规模。在主自由度处的位移被计算出来后,ANSYS 可将解扩展到原有的完整自由度集
6、上。(参见“模态分析”中的“矩阵缩减”部分对缩减过程的详细讨论。)缩减法的优点 是:比完全法快且开销小。缩减法的缺点 是:初始解只计算主自由度的位移,第二步进行扩展计算,得到完整空间上的位移、应力和力;不能施加单元载荷(压力,温度等),但允许施加加速度。所有载荷必须加在用户定义的主自由度上(限制在实体模型上施加载荷)。整个瞬态分析过程中时间步长必须保持恒定,不允许用自动时间步长。唯一允许的非线性是简单的点点接触(间隙条件)。3.4 完全法瞬态动力学分析首先,讲述完全法瞬态动力学分析过程,然后分别介绍模态叠加法和缩减法与完全法不相同的计算步骤。完全法瞬态动力分析(在ANSYS/Multiphsi
7、cs、ANSYS/Mechauioal及 ANSYS/Structural中可用)由以下步骤组成:1.建造模型2.建立初始条件3.设置求解控制4.设置其他求解选项5.施加载荷6.存储当前载荷步的载荷设置7.重复步骤 3-6 定义其他每个载荷步8.备份数据库9.开始瞬态分析10.退出求解器11.观察结果3.4.1建造模型在这一步中,首先要指定文件名和分析标题,然后用PREP7 定义单元类型,单元实常数,材料性质及几何模型。这些工作在大多数分析中是相似的。详细地说明了如何进行这些工作。对于完全法瞬态动力学分析,注意下面两点:可以用线性和非线性单元;必须指定杨氏模量EX(或某种形式的刚度)和密度DE
8、NS(或某种形式的质量)。材料特性可以是线性的或非线性的、各向同性的或各向异性的、恒定的或和温度有关的。划分合理的网格密度:网格密度应当密到足以确定感兴趣的最高阶振型;对应力或应变感兴趣的区域比只考察位移的区域的网格密度要细一些;如果要包含非线性特性,网格密度应当密到足以捕捉到非线性效应。例如,塑性分析要求在较大塑性变形梯度的区域有合理的积分点密度(即要求较密的网格);如果对波传播效果感兴趣(例如,一根棒的末端准确落地),网格密度应当密到足以解算出波动效应。基本准则是沿波的传播方向每一波长至少有20个单元。3.4.2建立初始条件在执行完全法瞬态动力学分析之前,用户需要正确理解建立初始条件和正确
9、使用载荷步。瞬态动力学分析顾名思义包含时间函数的载荷。为了定义这样的载荷,用户需要将载荷时间关系曲线划分成合适的载荷步。载荷时间曲线上的每个“拐角”对应一个载荷步,如图3.1 所示。图 3.1 载荷时间关系曲线第一个载荷步通常被用来建立初始条件,然后为第二和后继瞬态载荷步施加载荷并设置载步选项。对于每个载荷步,都要指定载荷值和时间值,同时指定其它的载荷步选项,如采用阶梯加载还是斜坡加载方式施加载荷以及是否使用自动时间步长等。然后,将每个载荷步写入载荷步文件,最后一次性求解所有载荷步。施加瞬态载荷的第一步是建立初始条件(即零时刻时的情况)。瞬态动力学分析要求给定两种初始条件(因为要求解的方程是两
10、阶的):初始位移()和初始速度()。如果没有进行特意设置,和都被假定为 0。初始加速度()一般假定为 0,但可以通过在一个小的时间间隔内施加合适的加速度载荷来指定非零的初始加速度。下面的段落描述了如何施加不同组合形式的初始条件。3.4.2.1零初始位移和零初始速度这是缺省的初始条件,即如果=0,则不需要指定任何条件。在第一个载荷步中可以加上对应于载荷/时间关系曲线的第一个拐角处的载荷。3.4.2.2非零初始位移及/或非零初始速度可以用 IC 命令设置这些初始条件。命令:IC GUI:Main MenuSolution-Loads-ApplyInitial Condit nDefine注意:不要
11、定义矛盾的初始条件。例如,在某单一自由度处定义了初始速度,则在所有其它自由度处的初始速度将为0.0,潜在地会产生冲突的初始条件。在大多数情形下要在模型的每个未约束自由度处定义初始条件。如果这些条件对各自由度是不同的,那么就可以较容易地明确指定初始条件,如下所述。关于 TIMINT 和 IC 命令的说明参见。3.4.2.3零初始位移和非零初始速度非零速度是通过对结构中需指定速度的部分加上小时间间隔上的小位移来实现的。比如如果=0.25,可以通过在时间间隔0.004 内加上 0.001 的位移来实现,命令流如下:.TIMINT,OFF!Time integration effects off D,
12、ALL,UY,.001!Small UY displ.(assuming Y-direction velocity)TIME,.004!Initial velocity=0.001/0.004=0.25 LSWRITE!Write load data to load step file(Jobname.S01)DDEL,ALL,UY!Remove imposed displacements TIMINT,ON!Time integration effects on.3.4.2.4非零初始位移和非零初始速度和上面的情形相似,不过施加的位移是真实数值而非“小”数值。比如,若=1.0 且=2.5,则
13、应当在时间间隔0.4 内施加一个值为 1.0 的位移:.TIMINT,OFF!Time integration effects off D,ALL,UY,1.0!Initial displacement=1.0 TIME,.4!Initial velocity=1.0/0.4=2.5 LSWRITE!Write load data to load step file(Jobname.S01)DDELE,ALL,UY!Remove imposed displacements TIMINT,ON!Time integration effects on.3.4.2.5非零初始位移和零初始速度需要用两
14、个子步 NSUBST,2来实现,所加位移在两个子步间是阶跃变化的KBC,1。如果位移不是阶跃变化的(或只用一个子步),所加位移将随时间变化,从而产生非零初速度。下面的例子演示了如何施加初始条件=1.0,=0.0:.TIMINT,OFF!Time integration effects off for static solution D,ALL,UY,1.0!Initial displacement=1.0 TIME,.001!Small time interval NSUBST,2!Two substeps KBC,1!Stepped loads LSWRITE!Write load data
15、 to load step file(Jobname.S01)!transient solution TIMINT,ON!Time-integration effects on for transient solution TIME,.!Realistic time interval DDELE,ALL,UY!Remove displacement constraints KBC,0!Ramped loads(if appropriate)!Continue with normal transient solution procedures.3.4.2.6非零初始加速度可以近似地通过在小的时间
16、间隔内指定要加的加速度ACEL实现。例如,施加初始加速度为 9.81 的命令如下:.ACEL,9.81!Initial Y-direction acceleration TIME,.001!Small time interval NSUBST,2!Two substeps KBC,1!Stepped loads LSWRITE!Write load data to load step file(Jobname.S01)!transient solution TIME,.!Realistic time interval DDELE,.!Remove displacement constraint
17、s(if appropriate)KBC,0!Ramped loads(if appropriate)!Continue with normal transient solution procedures.参见中关于命令 ACEL、TIME、NSUBST、KBC、LSWRITE、DDELE 和 KBC 的论述。3.4.3设置求解控制设置求解控制涉及定义分析类型、分析选项以及载荷步设置。执行完全法瞬态动力学分析,可以使用最新型的求解界面(称为求解控制对话框)进行这些选项的设置。求解控制对话框提供大多数结构完全法瞬态动力分析所需要的缺省设置,即用户只需要设置少量的必要选项。完全法瞬态动力分析建议采
18、用求解控制对话框,本章将详细进行介绍。如果完全瞬态动力分析需要初始条件,必须在分析的第一个载荷步进行,然后反复利用求解控制对话框为后续荷步设置载载荷步选项(即重复求解的3-6步)。如果不喜欢使用求解控制对话框(Main MenuSolution-Analysis Type-Soln Control),仍然可以沿用标准ANSYS 求解命令及其对应的菜单路径(Main MenuSolutionUnabridged Menuoption)。求解控制对话框一般形式参见 ANSYS 基本分析指南的针对确定的结构分析类型选用特定的求解控制。3.4.3.1使用求解控制对话框选择菜单路径 Main MenuS
19、olution-Analysis Type-Soln Control,就弹出求解控制对话框。下面将详细讲述求解控制对话框各页片夹中的选项。想要知道设置各选项的细节,选择感兴趣的页片夹,然后单击Help 按钮。本章还会讲述相关非线性结构分析的一些细节问题。3.4.3.2使用页片夹求解控制对话框包含5 各页片夹,各页片夹中分组设置控制选项,并将大多数基本控制选项设置在第一个页片夹中,其他页片夹提供更高级的控制选项。通过各页片夹,轻松达到控制求解过程。打开求解控制对话框,基本页片夹总是处于激活状态,只包含ANSYS 分析所需要设置的最少选项。如果基本页片夹已经满足控制要求,其他高级选项只有缺省状态不
20、符合求解控制才需要进一步进行调整。一旦单击任何页片夹中的OK按钮,所有求解控制对话框中选项设置都定义到ANSYS 数据库中,同时关闭求解控制对话框。可以是用基本页片夹设置下表中的选项。打开求解控制对话框,选择Basic页片夹,进行设置。基本页片夹选项选项更多信息参见:指定分析类型 ANTYPE,NLGEOM ANSYS基本分析指南中的定义分析类型和分析设置ANSYS结构分析指南中的非线性结构分析ANSYS基本分析指南中的重启动分析控制时间设置,包括:载荷步终点时间 TIME,自动时间步长 AUTOTS 以及载荷步内的子步数 NSUBST 或 DELTIM ANSYS基本分析指南中的跟踪中的时间
21、作用ANSYS基本分析指南中的一般选项指定写入数据库的求解(结果)数据 OUTRES ANSYS基本分析指南中的输出控制在瞬态动力学中,这些选项的特殊考虑有:1)当设置 ANTYPE 和 NLGEOM 时,如果执行一个新分析希望忽略大位移效应,如大变形、大转角和大应变,就选择小位移瞬态。如果希望考虑大变形(如弯曲的长细杆件)或大应变(如金属成型),就选择大位移瞬态。如果希望重启动一个失败的非线性分析,或者前面完成一个静态预应力分析或完全法瞬态动力分析,而后希望继续下面的时间历程计算,就可以选择重启动当前分析。2)当设置AUTOTS时,记住该载荷步选项(瞬态动力学分析中也称为时间步长优化)基于结
22、构的响应增大或减小积分时间步长。对于多数问题,建议打开自动时间步长与积分时间步长的上下限。通过 DELTIM和 NSUBST 指定积分步长上下限,有助于限制时间步长的波动范围;更多信息参见Automatic Time Stepping。缺省值为不打开自动时间步长。3)NSUBST和DELTIM是载荷步选项,用于指定瞬态分析积分时间步长。积分时间步长是运动方程时间积分中的时间增量。时间积分增量可以直接或间接指定(即通过子步数目)。时间步长的大小决定求解的精度:它的值越小,精度就越高。使用时应当考虑多种因素,以便计算出一个好的积分时间步长,详情参见积分时间步长章节。4)当设置 OUTRES 时,记
23、住下面注意事项:在完全法瞬态动力分析,缺省时只有最后子步(时间点)写入结果文件(Jobname.RST)为了将所有子步写入,需要设置所有子步的写入频率。同时,缺省时只有 1000 个结果序列能够写入结果文件。如果超过这个数目(基于用户的 OUTRES 定义),程序将认为出错终止。使用命令/CONFIG,NRES可以增大限制数(参见 ANSYS 基本分析指南中的内存和配置章节)。3.4.3.3使用瞬态页片夹利用瞬态页片夹设置其中的瞬态动力选项。有关设置这些选项的具体信息,打开求解控制对话框,选择瞬态页片夹,然后单击Help 按钮。瞬态页片夹选项选项更多信息参见:指定是否考虑时间积分效应 TIMI
24、NT ANSYS 结构分析指南中“执行非线性瞬态分析”指定采用渐变加载还是突变加载方式 KBC ANSYS 基本分析指南中“突变渐变载荷”ANSYS 基本分析指南中“突变或渐变载荷”指定质量和刚度阻尼 ALPHAD,BETAD ANSYS 结构分析指南中“阻尼”定义积分参数 TINTP ANSYS,Inc.Theory Reference 在完全法瞬态动力学中,这些选项的特殊考虑有:1)TIMINT 是动力载荷步选项,用于指定是否打开时间积分效应 TIMINT 。对于需要考虑惯性和阻尼效益的分析,必须打开时间积分效应(否则当作静力进行求解),所以缺省值为打开时间积分效应。进行完静力分析之后接着
25、进行瞬态分析时,该选项十分有用;也就是说,前面的载荷步必须关闭时间积分效应。2)ALPHAD (alpha 或 mass,damping)和 BETAD (beta 或 stiffness,damping)是动力载荷步选项,用于指定阻尼。大多数结构中都存在某种形式的阻尼,必须在分析中考虑进来。3)TINTP 是动力载荷步选项,用于指定瞬态积分参数。瞬态积分参数控制Newmark 时间积分技术,缺省值为采用恒定的平均值加速度积分算法。3.4.3.4使用求解选项页片夹求解选项页片夹选项用于完全法瞬态分析的具体设置完全与结构分析指南静力分析中一致。详情参见结构分析中使用求解选项页片夹。3.4.3.5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ansys 动力学 瞬态 分析 详解 20220301011126
限制150内