NASTRAN_动力分析指南
nastran动力学培训
（续）
SOL 109 TIME 30 CEND TITLE= TRANSIENT RESPONSE WITH TIME DEPENDENT PRESSURE AND POINT LOADS SUBTITLE= USE THE DIRECT METHOD ECHO= PUNCH SPC= 1 SET 1= 11, 33, 55 DISPLACEMENT= 1 SUBCASE 1 DLOAD= 700 $ SELECT TEMPORAL COMPONENT OF TRANSIENT LOADING （必须） LOADSET= 100 $ SELECT SPACIAL DISTRIBUTION OF TRANSIENT LOADING（可选） TSTEP= 100 $ SELECT INTEGRATION TIME STEPS （必须） $
（续）
PARAM，LFREQ 给出保留模态频率的下限 PARAM，HFREQ 给出保留模态频率的上限 PARAM，LMODES 给出被保留的最低频率模态的数量。 载荷卡片：
时变载荷 TLOAD1（与表联用）
载荷类型（力、力矩或者强迫位 移、速度、加速度）
DAREA或SPCD入口
指定表的号
二、瞬态响应分析
分析目的：计算时变激励载荷作用下结构的动力行为。 载荷的形式可以是外力或强迫运动。 两种数值方法：直接法和模态法。直接法对全部耦合的 运动方程进行直接数值积分来求解；而模态法则是利用 结构的振型来对耦合的运动方程进行缩减和解耦，然后 再由单个模态响应的叠加得到问题的最终解答。 求解器：直接法 SOL 109；模态法 SOL 112 直接瞬态响应中的阻尼
时间步数 作用时间
三、频率响应分析
频率响应分析是计算在稳态振动激励作用下结构动力响应 的一种方法（比如偏心旋转部件在一组转动频率下的旋转 分析）。 在频率响应分析中，激励载荷是在频域中明确定义的，所 有外力在每一个指定的频率上都是已知的。而力的形式可 以是外力、也可以是强迫运动。 与瞬态分析一样，也有两种方法供选用：直接法和模态法。 对应的求解器为SOL108、SOL111。
nastran操作指南
“机械工程有限元分析基础”本科生课程有限元分析软件MSC.NASTRAN2005r2ed操作指南南京航空航天大学机电学院设计工程系陈剑张保强郭勤涛2007年11月有限元结构静力与动态分析详细步骤南京航空航天大学机电学院设计工程系陈剑张保强郭勤涛一、分析目的有限元分析(FEA)是对物理现象（几何及载荷工况）的模拟，是对真实情况的数值近似。
通过划分单元，求解有限个数值来近似模拟真实环境的无限个未知量。
借助有限元分析软件进行结构静力与结构动力分析可以节省大量的时间。
通过本分析可以熟悉有限元软件patran与nastran的使用。
二、分析内容1、使用nastran进行一个悬臂梁的静力分析和动力分析2、使用nastran进行直齿圆柱齿轮的静力分析三、使用软件简单介绍MSC.Patran作为一个优秀的前后之处理器，具有高度的集成能力和良好的适用性：自动有限元建模: MSC.Patran的新产品中不断增加了很多更灵活更方便的智能化工具, 同时提供了自动网格及工业界最先进的映射网格划分功能, 使用户快速完成他们想做的工作。
同时也提供手动和其它有限元建模方法，一满足不同的需求。
分析的集成：MSC.Patran提供了众多的软件接口，将世界上大部分著名的不同类型分析软件和技术集于一体，为用户提供一个公共的环境。
这样可以使用户不必担心不同软件之间的兼容问题，在其它软件中建立的模型，在MSC.Patran 中仍然可以正常使用，非常灵活。
用户也能够根据多种类型的仿真结果对产品的整体设计给出正确的判断, 进行相应的改进，这就大大的提高了工作效率。
用户可自主开发新的功能：用户可将MSC.Patran作为自己的前后置处理器, 并利用其强大的PCL（Patran Command Language ）语言和编程函数库把自行开发的应用程序和功能及针对特殊要求开发的内容直接嵌入MSC.Patran的框架系统, 或单独使用或与其它系统联合使用。
NASTRAN_动力分析指南
NASTRAN_动⼒分析指南第⼀章动⼒学分析⽅法及NX NASTRAN基本使⽤介绍1.1 有限元分析⽅法介绍计算机软硬件技术的迅猛发展，给⼯程分析、科学研究以⾄⼈类社会带来急剧的⾰命性变化，数值模拟即为这⼀技术⾰命在⼯程分析、设计和科学研究中的具体表现。
数值模拟技术通过汲取当今计算数学、⼒学、计算机图形学和计算机硬件发展的最新成果，根据不同⾏业的需求，不断扩充、更新和完善。
近三⼗年来，计算机计算能⼒的飞速提⾼和数值计算技术的长⾜进步，诞⽣了商业化的有限元数值分析软件，并发展成为⼀门专门的学科－计算机辅助⼯程CAE（Computer Aided Engineering）。
这些商品化的CAE软件具有越来越⼈性化的操作界⾯和易⽤性，使得这⼀⼯具的使⽤者由学校或研究所的专业⼈员逐步扩展到企业的产品设计⼈员或分析⼈员，CAE在各个⼯业领域的应⽤也得到不断普及并逐步向纵深发展，CAE⼯程仿真在⼯业设计中的作⽤变得⽇益重要。
许多⾏业中已经将CAE分析⽅法和计算要求设置在产品研发流程中，作为产品上市前必不可少的环节。
CAE仿真在产品开发、研制与设计及科学研究中已显⽰出明显的优越性：●CAE仿真可有效缩短新产品的开发研究周期；●虚拟样机的引⼊减少了实物样机的试验次数；●⼤幅度地降低产品研发成本；●在精确的分析结果指导下制造出⾼质量的产品；●能够快速的对设计变更作出反应；●能充分的和CAD模型相结合并对不同类型的问题进⾏分析；●能够精确的预测出产品的性能；●增加产品和⼯程的可靠性；●采⽤优化设计，降低材料的消耗或成本；●在产品制造或⼯程施⼯前预先发现潜在的问题；●模拟各种试验⽅案，减少试验时间和经费；●进⾏机械事故分析，查找事故原因；●等等当前流⾏的商业化CAE软件有很多种，国际上早20世纪在50年代末、60年代初就投⼊⼤量的⼈⼒和物⼒开发具有强⼤功能的有限元分析程序。
其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。
P023-用MSC.Nastran进行流固耦合系统的动力学分析
用MSC.Nastran进行流固耦合系统的动力学分析王安平刘兵山中国科学院光电研究院北京 100190摘要：本文用Nastran2005对一个流固耦合系统进行了模态分析，结合一个密闭的薄壁结构模型，给出了分析的一般过程和需要注意的问题，也给出了该薄壁结构的模态频率、空腔系统的声学模态频率，以及耦合系统中，结构和空腔的声学模态频率和振型的变化。
关键词：Nastran，流固耦合，声学Modal Analysis Using MSC.Nastran for CoupledFluid-Structure SystemWANG Anping, LIU BingshanAcademy of Opto-Electronics, Chinese Academy of Sciences, Beijing, 100190 ChinaAbstract：The paper introduced the modal analysis method for the coupled fluid-structure (CFS) system using MSC.Nastran2005. Combinedthe model of a sealed laminated structure, analytical approach andwatchful items are presented. And making use of the MSC software, thestructure modal analysis and the cavity acoustic modal analysis of the CFSsystem are simulated.Keywords: Nastran, Coupled fluid-structure, Acoustics0 前言流固耦合法广泛应用于声学和噪声控制领域，如发动机的噪声控制。
对空腔结构(比如汽车车室、宇宙飞船船舱)进行流固耦合分析，可以知道耦合作用对系统模态的影响，可为研究耦合系统的声学特性提供可靠的理论和试验依据。
NX_Nastran__超单元用户指南
因为处理单个超单元对计算机资源的需求比不使用超单元的整个模型要少，因此可以使用快速序列提交各超单元的处理 (或者甚至在不同的计算机上处理)，而不是通宵等待整个问题的一次性求解。
减小风险
处理一个不使用超单元的模型，是一种要么全有要么全无的做法。如果出现错误，必须在改正错误后重新处理整个模型。在使用超单元时，每个超单元只需处理一次，除非需要修改超单元以重新处理。如果在处理中发生错误，只有受影响的超单元和残余结构 (最后处理的超单元) 需要重新处理。无错误的超单元不需要重新处理，除非对该超单元进行修改。
安全性
许多公司进行私人或保密的项目。包括为了竞争而开发的新设计，到高度机密的国防项目。即使在安全程序下工作，有时也需要发送模型给合作者，以便他们进行部件的耦合分析。使用外部超单元允许用户发送减缩后的边界矩阵，其中不包含实际部件的几何信息，只有从边界上看到的质量、刚度、阻尼和载荷。在收到一组可以由 NX Nastran 读出的任意格式的减缩矩阵后，工程师可以用这些矩阵定义一个外部超单元，并将外来结构附加到自己的模型上。
另一类超单元是外部超单元，模型的一个零件由一组源于外部的矩阵所代表 (该矩阵也可以来自另一个 NX Nastran 分析的结果)。对这些矩阵，没有任何内部的几何信息可用，只有这些矩阵所附着的节点是已知的。图 1-3 显示了一个外部超单元的例子。图中，左边是有限元模型，右边虚线表示外部超单元。
在静力分析中，用于超单元处理的理论是精确的。在动力分析中，对质量、阻尼和载荷矩阵的减缩是近似的。在第 9 章 “动力分析超单元介绍” 中对这些近似作了解释，并说明了改善近似性的方法 (特别是部件模态综合方法)。
■介绍
■如何对情况控制进行内部分区和使用
■多个载荷情况的例子
MSCNASTRAN 颤振分析模块使用说明
1.MSC/NASTRAN 颤振分析模块使用说明1.1.颤振分析模块颤振分析模块考虑结构气动弹性问题的动力稳定性。
它可以分析亚音速或超音速流，提供五种不同的气动力理论，包括用于亚音速的Doublet Lattice理论、Strip 理论以及用于超音速的Machbox理论、Piston理论、ZONA理论等。
对于稳定性分析，系统提供三种不同的方法：二种美国方法(K法，KE法)和一种英国方法(PK 法)，输出结果包括阻尼、频率和每个颤振模态的振型。
本说明仅以亚音速Doublet Lattice理论为例。
1.2.建模的一般流程其中结构有限元建模技术较为普及，不予说明。
升力面建模和颤振分析文件以填卡较为实用，大致包括：1)建立气动坐标系；2)设定影响体；3)选择颤振解法；4)给出飞行环境；5)给出马赫数和减缩频率系列；6)设定求解参数，如参与耦合的频率范围或模态数；7)选择适当的气动理论，定义升力面几何及分网信息。
至此完成升力面建模，下一步定义结构结点与升力面单元的耦合，即选择适当的样条将升力面结点同结构结点联系起来。
其中升力面结点是在定义升力面后由系统自动生成的，定义样条时直接引用升力面单元号；所以我们需要做的是将参与耦合的结构结点定义为一个集合，以便在样条定义中引用。
1.3.数据文件组织形式颤振分析模型数据文件遵循固定格式：设定求解时间、标题等；设置求解采用的特征值解法和颤振解法；输入模型数据即结构刚度和质量数据，还有升力面模型数据。
结构模型和升力面模型可以分别是独立的数据文件，只在颤振分析文件中将其包括进来。
下面以一个简单的例子(HA145B)来实现上述过程，并对颤振分析常用的卡片做简略介绍。
1.3.1.升力面模型文件$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$定义气动坐标系， 其X轴正向为来流方向（即将被AERO卡片引用）。
从hypermesh到nastran——模态和瞬态动力学分析关键步骤设置
hypermesh——nastran——模态分析。
模态分析关键步骤：1. 创建一个load collector, card image选择EIGRL(LANCZOS方法)。
然后editV1 –V2为频率范围，ND为阶数及方程组解的个数。
两者随意选择一个。
2. 创建loadstep，type为normal modes, method选中刚才创建的load collector。
3. 在control cards的sol选择nomal modes，, 如果想生成op2文件，把post也选上值为-1.4. 导出成bdf文件，启动nastran进行分析。
瞬态动力学分析如果激励是力比较好作，如果是强迫位移，老版本的需要用大质量或大刚度法把位移转换成力的载荷。
nastran 2001版以后可以直接加位移，关键步骤如下：1. 定义随时间历程曲线，创建load collectors，card image为Tabled12. 创建瞬态相应的时间步长和时间，load collectors, card image为Tstep3. 创建一个load collectors，card image为DAREA(如果是强迫位移不能用DAREA）4. 创建一个load collectors，card image为Tload1, excited选择DAREA，TID选择TSTEP，注意TYPE的选择。
5. 创建一个subcase，类型选择直接瞬态分析，DLOAD和TSTEP选择刚才创建的两个相对应的load collectors6. 导出成bdf文件，提交nastran进行分析。
如果是强迫位移，还要多两个卡，就是SPCD, LSEQ详细步骤跟以上差不多，只要把各个卡片弄懂了就很容易了。
nastran动力学培训
四、强迫运动
用于分析带有地基加速度、位移和速度的输入的受约束结 构。 直接指定法 例：一端固支的矩形结构，在地基上受到沿Z方向频率为 250HZ的单位正弦脉冲加速度作用，使用直接方法，确定 该结构的瞬态响应。在地基上施加1000lb的大质量，使用 的结构阻尼系数：g=0.06，并将此阻尼转化为在250HZ下 的等效粘性阻尼。
动力学培训内容介绍
1.模态分析 2.瞬态响应分析 3.频率响应分析 4.强迫运动 5.随机响应分析
结构动力学分析
一、模态分析 求解器：103 质量矩阵形式：MSC认为耦合质量比集中质量更精确，在 动力分析里出于对计算速度的考虑，更倾向于使用集中质 量。 使用方法：用PARAM，COUPMASS，1选择耦合质量； 缺省为集中质量。 求解方法：推荐的Lanczos方法。 EIGRL卡片
（续）
激励的定义：与瞬态响应中TLOAD对应，在频率响应中 为RLOAD。其中RLOAD1是按照实部与虚部的形式来定 义频变载荷；RLOAD2按幅值和相位的形式来定义频变载 荷。 几点考虑：如果激励的最高频率比系统的最低谐振频率小 得多，那么使用静态分析就足够了；阻尼很小的结构在激 励频率接近于谐振频率的时候，会表现出很大的动力响应。 在这样的问题中，模型上一个小的改动（或仅换一台电脑 来计算）都可能产生响应的明显变化；如果希望对峰值响 应进行充分的预测，必须使用足够好的频率步长（Δ f）。 对每个半能带宽至少使用5个点。
（续）
SOL 109 TIME 30 CEND TITLE = TRANSIENT RESPONSE WITH BASE EXCITATION SUBTITLE = USING DIRECT TRANSIENT METHOD, NO REDUCTION ECHO = UNSORTED SPC = 200 SET 111 = 23, 33 DISPLACEMENT (SORT2) = 111 VELOCITY (SORT2) = 111 ACCELERATION (SORT2) = 111 SUBCASE 1 DLOAD = 500 TSTEP = 100 $
NASTRAN动力学分析教程
MSC.NASTRAN的分析功能
MSC.NASTRAN的分析功能作为世界CAE工业标准及最流行的大型通用结构有限元分析软件, MSC.NASTRAN的分析功能覆盖了绝大多数工程应用领域,并为用户提供了方便的模块化功能选项,MSC.NASTRAN的主要功能模块有:基本分析模块(含静力、模态、屈曲、热应力、流固耦合及数据库管理等)。
动力学分析模块、热传导模块、非线性分析模块、设计灵敏度分析及优化模块、超单元分析模块、气动弹性分析模块、DMAP用户开发工具模块及高级对称分析模块。
除模块化外, MSC.NASTRAN还按解题规模分成10,000节点到无限节点,用户引进时可根据自身的经费状况和功能需求灵活地选择不同的模块和不同的解题规模, 以最小的经济投入取得最大效益。
MSC.NASTRAN及MSC的相关产品拥有统一的数据库管理,一旦用户需要可方便地进行模块或解题规模扩充, 不必有任何其它的担心。
MSC.NASTRAN以每年一个小版本, 每两年一个大版本的速度更新, 用户可不断获得当今CAE发展的最新技术用于其产品设计。
目前MSC.NASTRAN的最新版本是1999年发布的V70.5版。
新版本中无论在设计优化、P单元、热传导、非线性还是在数值算法、性能、文档手册等方面均有大幅度的改进或突出的新增功能。
以下将就MSC.NASTRAN不同的分析方法、加载方式、数据类型或新增的一些功能做进一步的介绍:⒈静力分析静力分析是工程结构设计人员使用最为频繁的分析手段, 主要用来求解结构在与时间无关或时间作用效果可忽略的静力载荷(如集中/分布静力、温度载荷、强制位移、惯性力等)作用下的响应, 并得出所需的节点位移、节点力、约束(反)力、单元内力、单元应力和应变能等。
该分析同时还提供结构的重量和重心数据。
MSC.NASTRAN支持全范围的材料模式,包括: 均质各项同性材料,正交各项异性材料, 各项异性材料,随温度变化的材料。
方便的载荷与工况组合单元上的点、线和面载荷、,热载荷、强迫位移,各种载荷的加权组合，在前后处理程序MSC.PA TRAN中定义时可把载荷直接施加于几何体上。
nastran动力学分析
S8-19
TSTEPNL – 时间步和结果
1
2
3
4
TSTEPNL
ID
NDT
5
6
7
8
DT
NO
9
10
• DT, NDT – 时间增量和总时间步
• NO – 结果输出间隔增量步
– 如果为正，输出位置由DT控制，默认值为1，每步输出。 – 如果为负，在NO个实际计算增量步后输出(SOL 129风格)
NAS400, Section 8, August 2011 Copyright 2011 MSC.Software Corporation
K * ( M, B, K, t ) U
等效动力 刚度
=
P * ( t, Ů, Ü, M, B, P )
等效动力 载荷矢量
• 等效动力刚度和载荷矢量由积分策略决定。 • 例如，使用平均加速度策略，也称之为trapezoidal法则或Newmark策略( =
1/2, = ¼)
NAS400, Section 8, August 2011 Copyright 2011 MSC.Software Corporation
S8-22
TSTEPNL –结果输出(续)
• 用户可以使用使用OTIME进行附加时间控制，例如：
OTIME = 99
SET 99 = 0.025, 0.035
additional output
$
SUBCASE 10
STEP 1
LOAD = 1
NLPARM = 110
STEP 2
ANALYSIS = NLTRAN
S8-11
用户界面
• 载荷定义
TLOAD1 TLOAD2 DAREA LSEQ NOLIN1 NOLIN2 NOLIN3 NOLIN4
NASTRAN动态分析
3、 ω与ωn大小相当时，H(ω)的大小主要取决于阻尼比ξ，此时称为阻尼控制区
在很多情况下，通过模态分析，得到结构的自然频率和模态振型就足够了。在 这种情况下，结构的自然频率应该与外界对结构的激励频率错开，避免共振。
动态分析的基本类型有： 1、实特征值分析即模态分析（无阻尼自由振动） 2、线性频响分析 3、线性瞬态响应分析
Nastran还包括其他高级的动态响应分析功能，例如强迫运动（enforced motion), 随机响应分析（random response analysis），响应谱分析（response spectrum analysis）等，这些都是在前面三个类型的基础上进行的。
F sint
得到的位移响应为： X x•sint
速度响应为： V • x• cost
加速度响应为： A • • x•sint
动刚度为：
K
F x
1
x
加速度导纳为： J 2 • x() 2
K ()
通常使用加速度导纳曲线来表征结构的动态特性。
பைடு நூலகம்Z方向出现 峰值
粘滞阻尼： 结构阻尼：
其中： 其中：
结构阻尼的输入： 1、在MATi的GE域中输入，模型中与此材料相关的单元形成阻尼力； 2、PARAM,G,r，输入全局结构阻尼。 粘滞阻尼的输入： 1、通过阻尼单元CDAMPi、CVISC、CBUSH输入； 2、模态瞬态响应分析与模态频响分析的粘滞阻尼通过TABDMP1输入。
频响分析所求出的响应也是在频率域上的，比如位移的形式可以表示为：
可见表达响应的结果中需要包含两方面的信息：幅值和相位 一种表达响应的简洁的方式是使用复数表达，如下
nastran动力学分析中的矩阵组集
nastra n动力学分析中的矩阵组集1概述1）在瞬态响应分析、频率响应分析、复模态分析中，方法Nastran提供了直接方法和模态2）分析类型和方法的不同，动力学矩阵组集也不一样2阻尼矩阵 2.1阻尼概述1）阻尼反映结构内部能量的耗散2）阻尼产生的机理a）粘性效应（粘性阻尼器，振动减振器）b）外摩擦（结构连接处的相对滑动）c）内摩擦（材料粘性）d）结构非线性（Plasticity ）3）阻尼的模拟a）粘性阻尼力f v = eum u + bu + k u = pb）结构阻尼力= igilcu where i = 1g = structural dampino coefffclerit2.2结构阻尼与粘性阻尼假设结构简谐响应为u = ue itot对粘性阻尼mu + bu + ku Pd)m(—oj?Oe + t>Cia>ue"<o^) + kue"= pft)uj2mue i£nt + ibcnue iwt + kue i<nt = p(t)对结构阻尼mu + (1 + ig)ku = p( t)m(—+ (1 + ig)kuei°jt= p(t)+ igkDe*to^ + kue'<o^ —p( t)可以得到gk = bw -> b =—如果但因为b = 2mwc n得到b 尸gb c= 2f = critical damping ratio (percent critical damping) g = Q= structural damping factorQ = quality factor or magnification factor其中，结论1) 粘性阻尼与速度成比例2) 结构阻尼与位移成比例3) 临界阻尼比b/b cr=；4) 品质因子与能量耗散成反比Q = 3/2Q = 1/(20Q = 1/g2.3阻尼输入1)结构阻尼a) MATi卡片b) PARAM,G, factor (Default = 0.0)用结构阻尼系数乘整个系统刚度矩阵c) PARAM,W3, factor (Default = 0.0)将结构阻尼转化为等效粘性阻尼d) PARAM,W4, factor (Default = 0.0)将单元结构阻尼转化为等效粘性阻尼e) W3,W4的单位为rad/unit timef) 如果使用PARAM,G贝U PARAM,W3勺factor必须大于0,否则，瞬态响应分析中将忽略PARAM,G2)标量粘性阻尼CDAMP1两自由度间的数值阻尼器，需属性卡PDAMPCDAMP2两自由度间的数值阻尼器，不需属性卡PDAMPCDAMP3两自由度间的数值阻尼器，需属性卡PDAMPCDAMP4两自由度间的数值阻尼器，不需属性卡PDAMPCVISC :两节点间的单元阻尼器，需属性卡PVISCCBUSH广义弹簧和阻尼器单元，允许与频率相关3) 模态阻尼TABDMP1 Lists damping value versus frequencyIn Case Control selects TABDMP1 in Bulk Data SDAMP = n3直接法直接法中使用的动力学方程为'M dd p2+B dd p + K dd l{u d } = {p d }其中，p = a derivative operatoru d = the union of the 族 n a lysis set u a and extra points u c对频率响应和复特征值分析，动力学矩阵为4模态法模态法中的动力学方程为l IM hhP 2*B hhP + K hh |{u h> = %}模态坐标与物理坐标间变换为= %山弭 叫} =%叫}其中，性i °%旷 o -eeI u h }=对频率响应和复特征值分析，动力学矩阵为'K dd'【％】I M dd I对瞬态响应，动力学矩阵为d 2d[K+ I I dd2d 1d B d 2d4d[K ■.^1d d d d jo d K B M 丨rlr-d d d 1d1d1d K B Md d d 2d K B[I n ng%d 1d KIKhh"旳屮怙』佃【昭(讣和7K爲|)|愉1 l B hh' = l b il + l*dh'T(l B ddl +l B ddl)t*dhl[M hh l = [m i] + [|dh lT[l M2<J][ldh]如果，KDAMP = -1,贝S% = m Sibjj = O岭=(1 +ig(w j))k jj对瞬态响应，动力学矩阵为l K hhl =l k i' + ^dhl T t K dd^dhlEhl = I b i I +I*dh 'T(B dd + B dd + K dd I + K dd l)$dh I [%】=lm ii l + [4>dh]T[M2<J][0(Jh]。
Autodesk Nastran In-CAD Dynamics 动态分析手册说明书
1. Activate the Nastran environment 2. Double click on the Analysis 1 3. In the Model State tab change the
LOD to Two Brackets 4. Press OK
5. Note that the CAD materials are
2. Change the name to Mounts 3. Select the Associated Geometry box 4. Select the two brackets to complete
the last property definition 5. Click OK 6. Your tree should now look like this:
© 2017 Autodesk
3
Open the Part
1. Open the file Muffler&Brackets.iam (Get Started > Launch > Open)
8/11/2017
© 2017 Autodesk
Use Level of Details in Nastran In-Cad
Transient Response Analysis
User-defined Time Steps Control Solution Accuracy + Results Insight
© 2017 Autodesk
8/11/2017 2
8/11/2017
Example 4.1: Muffler Modal Analysis
box 6. In the Selected Quadrilaterals:
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一章动力学分析方法及NX NASTRAN基本使用介绍1.1 有限元分析方法介绍计算机软硬件技术的迅猛发展，给工程分析、科学研究以至人类社会带来急剧的革命性变化，数值模拟即为这一技术革命在工程分析、设计和科学研究中的具体表现。
数值模拟技术通过汲取当今计算数学、力学、计算机图形学和计算机硬件发展的最新成果，根据不同行业的需求，不断扩充、更新和完善。
近三十年来，计算机计算能力的飞速提高和数值计算技术的长足进步，诞生了商业化的有限元数值分析软件，并发展成为一门专门的学科－计算机辅助工程CAE（Computer Aided Engineering）。
这些商品化的CAE软件具有越来越人性化的操作界面和易用性，使得这一工具的使用者由学校或研究所的专业人员逐步扩展到企业的产品设计人员或分析人员，CAE在各个工业领域的应用也得到不断普及并逐步向纵深发展，CAE工程仿真在工业设计中的作用变得日益重要。
许多行业中已经将CAE分析方法和计算要求设置在产品研发流程中，作为产品上市前必不可少的环节。
CAE仿真在产品开发、研制与设计及科学研究中已显示出明显的优越性：●CAE仿真可有效缩短新产品的开发研究周期；●虚拟样机的引入减少了实物样机的试验次数；●大幅度地降低产品研发成本；●在精确的分析结果指导下制造出高质量的产品；●能够快速的对设计变更作出反应；●能充分的和CAD模型相结合并对不同类型的问题进行分析；●能够精确的预测出产品的性能；●增加产品和工程的可靠性；●采用优化设计，降低材料的消耗或成本；●在产品制造或工程施工前预先发现潜在的问题；●模拟各种试验方案，减少试验时间和经费；●进行机械事故分析，查找事故原因；●等等当前流行的商业化CAE软件有很多种，国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。
其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。
该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。
从那时到现在，世界各地的研究机构和大学也发展了一批专用或通用有限元分析软件，除了Nastran以外，主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABAQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。
虽然软件种类繁多，但是万变不离其宗，其核心求解方法都是有限单元法，也简称为有限元法（Finite Element Method）。
1.1.1 有限单元法的基本思路有限元法的基本思路可以归结为：将连续系统分割成有限个分区或单元，对每个单元提出一个近似解，再将所有单元按标准方法加以组合，从而形成原有系统的一个数值近似系统，也就是形成相应的数值模型。
下面用在自重作用下的等截面直杆来说明有限元法的思路。
等截面直杆在自重作用下的材料力学解答：图1-1 受自重作用的等截面直杆图1-2 离散后的直杆受自重作用的等截面直杆如图1-3所示，杆的长度为L ，截面积为A ，弹性模量为E ，单位长度的重量为q ，杆的内力为N 。
试求：杆的位移分布，杆的应变和应力。
)()(x L q x N -=EA dxx L q EA dx x N x dL )()()(-==⎰-==x x Lx EA q EA dx x N x u 02)2()()(（1- 1）)(x L EA q dx du x -==ε)(x L A qE x x -==εσ等截面直杆在自重作用下的有限元法解答： 1）连续系统离散化如图1-4所示，将直杆划分成n 个有限段，有限段之间通过公共点相连接。
在有限元法中，我们将两段之间的公共连接点称为节点，将每个有限段称为单元。
节点和单元组成的离散模型就称为对应于连续系统的‘有限元模型’。
有限元模型中的第i 个单元，其长度为Li ，包含第i ，i+1个节点。
2）用单元节点位移表示单元内部位移第i 个单元中的位移用所包含的节点位移来表示，)()(1i iii i x x L u u u x u --+=+（1- 2）其中iu 为第i 节点的位移，ix 为第i 节点的坐标。
第i 个单元的应变为i ε，应力为i σ，内力为iN ：iii i L u u dx du -==+1ε （1- 3）ii i i i L u u E E )(1-==+εσ （1- 4）ii i i i L u u EA A N )(1-==+σ （1- 5）3）把外载荷归集到节点上把第i 单元和第i+1单元重量的一半2)(1++i i L L q ，归集到第i+1节点上。
图1-3 集中单元重量4）建立节点的力平衡方程对于第i+1节点，由力的平衡方程可得：2)(11+++=-i i i i L L q N N（1- 6）令1+=i ii L L λ，并将（1- 8）代入得：221)11(2)1(iii i i i i L EA q u u u λλλ+=-++-++（1-7）根据约束条件，01=u 。
对于第n+1个节点，2nn qL N =EA qL u u n n n 221=+-+（1-8）建立所有节点的力平衡方程，可以得到由n+1个方程构成的方程组，可解出n+1个未知的节点位移。
1.1.2 有限元法的计算步骤有限元法的计算步骤归纳为以下三个基本步骤：网格划分，单元分析，整体分析。
1）网格划分有限元法的基本做法是用有限个单元体的集合来代替原有的连续体。
因此首先要对弹性体进行必要的简化，再将弹性体划分为有限个单元组成的离散体。
单元之间通过节点相连接。
由单元、节点、节点连线构成的集合称为网格。
通常把三维实体划分成4面体或6面体单元的实体网格，平面问题划分成三角形或四边形单元的面网格。
图1-4四面体四节点单元图1-5 六面体8节点单元图1-6 四边形4节点单元2）单元分析对于弹性力学问题，单元分析就是建立各个单元的节点位移和节点力之间的关系式。
由于将单元的节点位移作为基本变量，进行单元分析首先要为单元内部的位移确定一个近似表达式，然后计算单元的应变、应力，再建立单元中节点力与节点位移的关系式。
以平面问题的三角形3节点单元为例，见图1-7。
单元有三个节点I、J、M，每个节点有两个位移u、v和两个节点力U、V。
图1-7 三角形3节点单元单元的所有节点位移、节点力，可以表示为节点位移向量（vector）：节点位移{}⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=mmjjiievuvuvuδ节点力{}⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=mmjjiieVUVUVUF单元的节点位移和节点力之间的关系用张量（tensor）来表示，{}[]{}eee KFδ=（1-9）3）整体分析对由各个单元组成的整体进行分析，建立节点外载荷与节点位移的关系，以解出节点位移，这个过程称为整体分析。
同样以弹性力学的平面问题为例，如图1-8所示，在边界节点i上受到集中力iyixPP,作用。
节点i是三个单元的结合点，因此要把这三个单元在同一节点上的节点力汇集在一起建立平衡方程。
图1-8 整体分析i 节点的节点力：∑=++ee i i i i U U U U )()3()2()1(∑=++ee i i i i V V V V )()3()2()1(i 节点的平衡方程：⎪⎭⎪⎬⎫=∑∑=i y ee iei x e iP V P U )()(（1-10）1.1.3有限元法的进展与应用有限元法不仅能应用于结构分析，还能解决归结为场问题的工程问题，从二十世纪六十年代中期以来，有限元法得到了巨大的发展，为工程设计和优化提供了有力的工具。
当今国际上FEA 方法和软件发展呈现出以下一些趋势特征：● 从单纯的结构力学计算发展到求解许多物理场问题。
有限元分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值分析方法。
而且从理论上也已经证明，只要用于离散求解对象的单元足够小，所得的解就可足够逼近于精确值。
所以近年来有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算，最近又发展到求解几个交叉学科的问题。
例如比较常见的是将温度场和结构场之间进行耦合计算，确定由于温度场分布不均匀引起的结构应力和变形等。
● 由求解线性工程问题进展到分析非线性问题随着科学技术的发展，线性理论已经远远不能满足设计的要求。
例如建筑行业中的高层建筑和大跨度悬索桥的出现，就要求考虑结构的大位移和大应变等几何非线性问题；航天和动力工程的高温部件存在热变形和热应力，也要考虑材料的非线性问题；诸如塑料、橡胶和复合材料等各种新材料的出现，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性有限元算法才能解决。
众所周知，非线性的数值计算是很复杂的，它涉及到很多专门的数学问题和运算技巧，很难为一般工程技术人员所掌握。
为此近年来国外一些公司花费了大量的人力和投资开发求解非线性问题的分析功能，并广泛应用于工程实践。
●增强可视化的前后处理功能。
早期有限元分析软件的研究重点在于推导新的高效率求解方法和高精度的单元。
随着数值分析方法的逐步完善，尤其是计算机运算速度的飞速发展，整个计算系统用于求解运算的时间越来越少，而准备数值模型和处理计算结果的时间占整个分析工程的比例越来越高。
据统计，整个分析流程中，前处理占用的工作时间大致在80%，而加上后处理部分，占用的时间就要超过95%。
因此目前几乎所有的商业化有限元程序系统都有功能很强的前后处理模块与之相配合。
在强调"可视化"的今天，很多程序都建立了对用户非常友好的GUI（Graphics User Interface）,使用户能以可视图形方式直观快速地进行网格自动划分，生成有限元分析所需数据，并按要求将大量的计算结果整理成变形图、等值分布云图，便于极值搜索和所需数据的列表输出。
●与CAD软件的无缝集成。
当今有限元分析系统的另一个特点是与通用CAD软件的集成使用——即,在用CAD软件完成部件和零件的造型设计后，自动生成有限元网格并进行计算，如果分析的结果不符合设计要求则重新进行造型和计算，直到满意为止，从而极大地提高了设计水平和效率。
今天，工程师可以在集成的CAD和FEA软件环境中快捷地解决一个在以前无法应付的复杂工程分析问题。
所以当今所有的商业化有限元系统商都开发了和著名的CAD软件（例如Unigraphics、Pro/ENGINEER、SolidEdge、SolidWorks等）的接口。
1.2动力学分析概述1.2.1动力分析过程图1-9 动力分析过程1.2.2单自由度系统1)动力学方程其中，m为质量（惯性），b为阻尼（能量耗散），k为刚度（恢复力）n为非线性恢复力p为作用力u为位移u 为加速度u 为速度通常，作用力p、位移u、速度u 、加速度u 为时间函数，m、u, 的函数。