《ansys动力学瞬态分析详解》由会员分享，可在线阅读，更多相关《ansys动力学瞬态分析详解（88页珍藏版）》请在金锄头文库上搜索。
1、瞬态动力学分析的定义瞬态动力学分析（亦称时间历程分析）是用于确定承受任意的随时间变化载 荷结构的动力学响应的一种方法。可以用瞬态动力学分析确定结构在稳态载荷、 瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。 载荷和时间的相关性使得惯性力和阻尼作用比较重要。如果惯性力和阻尼作用不 重要，就可以用静力学分析代替瞬态分析。瞬态动力学的基本运动方程是：网冋+琲+附心即）其中：M=质量矩阵C二阻尼矩阵K二刚度矩阵茂=节点加速度向量茂=节点速度向量u二节点位移向量在任意给定的时间匚 这些方程可看作是一系列考虑了惯性力CM）和 阻尼力（C m）的静力学平衡方程。ANSYS程序使用Ne
2、wmark时间积分方法 在离散的时间点上求解这些方程。两个连续时间点间的时间增量称为积分时间步 长 （integration time step）。学习瞬态动力学的预备工作瞬态动力学分析比静力学分析更复杂，因为按工程时间计算，瞬态动力学 分析通常要占用更多的计算机资源和更多的人力。可以先做一些预备工作以理解 问题的物理意义，从而节省大量资源。例如，可以做以下预备工作：1. 首先分析一个较简单模型。创建梁、质量体和弹费组成的模型，以最小的 代价深入的理解动力学认识，简单模型更有利于全面了解所有的动力学响应所需 要的。2. 如果分析包括非线性特性，建议首先利用静力学分析掌握非线性特性对结 构响应的
3、影响规律。在某些场合，动力学分析中是没必要包括非线性特性的。3. 掌握结构动力学特性。通过做模态分析计算结构的固有频率和振型，了解 这些模态被激活时结构的响应状态。同时，固有频率对计算正确的积分时间步长 十分有用。4. 对于非线性问题，考虑将模型的线性部分子结构化以降低分析代价。vv高 级技术分指南 中将讲述子结构。三种求解方法瞬态动力学分析可采用三种方法：完全（Full）法、缩减（Reduced）法及模 态叠加法。ANSYS/Professional产品中只允许用模态叠加法。在研究如何实现这 些方法之前，让我们先探讨一下各种方法的优点和缺点。33.1完全法完全法采用完整的系统矩阵计算瞬态响应
4、（没有矩阵缩减）。它是三种方法 中功能最强的，允许包括各类非线性特性（塑性、大变形、大应变等）。油如果并不想包括任何非线性,应当考虑使用另外两种方法中的一种。这 是因为完全法是三种方法中开销最大的一种。完全法的优点是：容易使用，不必关心选择主自由度或振型。允许各种类型的非线性特性。采用完整矩阵，不涉及质量矩阵近似。在一次分析就能得到所有的位移和应力。允许施加所有类型的载荷：节点力、外加的（非零）位移（不建议釆用） 和单元载荷（压力和温度），还允许通过TABLE数组参数指定表边界条件。允许在实体模型上施加的载荷。完全法的主要缺点是它比其它方法开销大。3.3.2态叠加法模态叠加法通过对模态分析得到
5、的振型（特征值）乘上因子并求和来计算结 构的响应。此法是ANSYS/Professional程序中唯一可用的瞬态动力学分析法。模态叠加法的优点是：对于许多问题，它比缩减法或完全法更快开销更小；只要模态分析不采用PowerDynamics方法，通过LVSCALE命令将模态分析 中施加的单元载荷引入到瞬态分析中；允许考虑模态阻尼（阻尼比作为振型号的函数）。模态叠加法的缺点是：整个瞬态分析过程中时间步长必须保持恒定，不允许釆用自动时间步长；唯一允许的非线性是简单的点点接触（间隙条件）；不能施加强制位移（非零）位移。333缩减法缩减法通过采用主自由度及缩减矩阵压缩问题规模。在主自由度处的位移被 计算出
6、來后，ANSYS可将解扩展到原有的完整自由度集上。（参见“模态分析中 的“矩阵缩减部分对缩减过程的详细讨论。）缩减法的优点是：比完全法快且开销小。缩减法的缺点是：初始解只计算主自由度的位移，第二步进行扩展计算，得到完整空间上的 位移、应力和力；不能施加单元载荷（压力，温度等），但允许施加加速度。所有载荷必须加在用户定义的主自由度上邙艮制在实体模型上施加载荷）。整个瞬态分析过程中时间步长必须保持恒定，不允许用自动时间步长。唯一允许的非线性是简单的点一点接触（间隙条件）。完全法瞬态动力学分析首先，讲述完全法瞬态动力学分析过程，然后分别介绍模态叠加法和缩减法 与完全法不相同的计算步骤。完全法瞬态动力
7、分析（在ANSYS/Multiphsics、 ANSYS/Mechauioal 及 ANSYS/Structural 中可用）由以下步骤组成：1. 建造模型2. 建立初始条件3. 设置求解控制4. 设置其他求解选项5. 施加载荷6. 存储当前载荷步的载荷设置7. 重复步骤3-6定义其他每个载荷步8. 备份数据库9. 开始瞬态分析10. 退出求解器11. 观察结果3.4.1建造模型在这一步中，首先要指定文件名和分析标题，然后用PREP7定义单元类型， 单元实常数，材料性质及几何模型。这些工作在大多数分析中是相似的。ANSYS 建模与网格指南 详细地说明了如何进行这些工作。对于完全法瞬态动力学分析
8、，注意下面两点：可以用线性和非线性单元；必须指定杨氏模量EX （或某种形式的刚度）和密度DENS （或某种形式的 质量）。材料特性可以是线性的或非线性的、各向同性的或各向异性的、恒定的 或和温度有关的。划分合理的网格密度：网格密度应当密到足以确定感兴趣的最高阶振型；对应力或应变感兴趣的区域比只考察位移的区域的网格密度要细一些；如果要包含非线性特性，网格密度应当密到足以捕捉到非线性效应。例如, 塑性分析要求在较大塑性变形梯度的区域有合理的积分点密度（即要求较密的网 格）；如果对波传播效果感兴趣（例如，一根棒的末端准确落地），网格密度应 当密到足以解算出波动效应。基本准则是沿波的传播方向每一波长至
9、少有20个 单元。342建立初始条件在执行完全法瞬态动力学分析之前，用户需要正确理解建立初始条件和正确 使用载荷步。瞬态动力学分析顾名思义包含时间函数的载荷。为了定义这样的载荷，用户 需要将载荷一时间关系曲线划分成合适的载荷步。载荷一时间曲线上的每个“拐 角对应一个载荷步，如图所示。图载荷一时间关系曲线第一个载荷步通常被用來建立初始条件，然后为第二和后继瞬态载荷步施加载荷并设置载步选项。对于每个载荷步，都要指定载荷值和时间值，同时指定其它的载荷步选项，如采用阶梯加载还是斜坡加载方式施加载荷以及是否使用自动时间步长等。然后，将每个载荷步写入载荷步文件，最后一次性求解所有载荷步。施加瞬态载荷的第一
10、步是建立初始条件（即零时刻时的情况）。瞬态动力学 分析要求给定两种初始条件（因为要求解的方程是两阶的）：初始位移（） 和初始速度（肌）。如果没有进行特意设置，吗和肌都被假定为0。初始加 速度（肌）一般假定为0,但可以通过在一个小的时间间隔内施加合适的加速 度载荷來指定非零的初始加速度。下面的段落描述了如何施加不同组合形式的初始条件。3.4.2.1零初始位移和零初始速度这是缺省的初始条件，即如果= 0,则不需要指定任何条件。在第一个载荷步中可以加上对应于载荷/时间关系曲线的第一个拐角处的载荷。3.4.2.2非零初始位移及/或非零初始速度可以用IC命令设置这些初始条件。命令：ICGUI： Main
11、 MenuSolution-Loads-Applylnitial ConditznDefine注意:不要定义矛盾的初始条件。例如,在某单一自由度处定义了初始速度, 则在所有其它自由度处的初始速度将为,潜在地会产生冲突的初始条件。在大 多数情形下要在模型的每个未约束自由度处定义初始条件。如果这些条件对各 自由度是不同的,那么就可以较容易地明确指定初始条件,如下所述。关于TIMINT和IC命令的说明参见ANSYS命令参考手册。3.4.23零初始位移和非零初始速度非零速度是通过对结构中需指定速度的部分加上小时间间隔上的小位移來实 现的。比如如果肌 =,可以通过在时间间隔内加上的位移来实现，命令流如
12、下： TIMINTQFF! Time integration effects offD,ALL,UY,001! Small UY displ. (assuming Y-directio n velocity)TIME,.004! Initial velocity =LSWRITE! Write load data to load step fileDDELALLUY! Remove imposed displacementsTIMINTQN! Time integration effects on 3424非零初始位移和非零初始速度和上面的情形相似，不过施加的位移是真实数值而非小数值。比如，若
13、血二且珈 =,则应当在时间间隔内施加一个值为的位移： TIMINTQFF! Time integration effects offD,ALL,UYInitial displacement =TIME,.4! Initial velocity =LSWRITE! Write load data to load step fileDDELE,ALL,UY! Remove imposed displacementsTIMINTQN! Time integration effects on3425非零初始位移和零初始速度需要用两个子步NSUBST,2來实现，所加位移在两个子步间是阶跃变化的 KBCz
14、lo如果位移不是阶跃变化的(或只用一个子步)，所加位移将随时间变化,从而产生非零初速度。下面的例子演示了如何施加初始条件=,肌= : TIMINTQFF! Time integration effects off for static solutionD,ALL,UYJ Initial displacement =TIME,.001! Small time intervalNSUBST2! Two substepsKBC,1! Stepped loadsLSWRITE! Write load data to load step fileItransient soluti ontTIMINTON
15、! Time-integrati on effects on for transient solutionTIME,.! Realistic time intervalDDELE,ALL,UY! Remove displacement constraintsKBCO! Ramped loads (if appropriate)Continue with nor mal transient solution procedures 3.4.2.6非零初始加速度可以近似地通过在小的时间间隔内指定要加的加速度ACEL实现。例如，施 加初始加速度为的命令如下：ACEL,J Initial Y-direction accelerationTIME,.001!Small time intervalNSUBS2!Two substepsKBCl
