abaqus中地动态分析报告方法.doc
《abaqus中地动态分析报告方法.doc》由会员分享,可在线阅读,更多相关《abaqus中地动态分析报告方法.doc(32页珍藏版)》请在课桌文档上搜索。
1、wordABAQUS线性动态分析如果你只对结构承受载荷后的长期响应感兴趣,静力分析static analysis是足够的。然而,如果加载时间很短例如在地震中或者如果载荷在性质上是动态的例如来自旋转机械的荷载,你就必须采用动态分析dynamic analysis。本章将讨论应用ABAQUS/Standard进展线性动态分析;关于应用ABAQUS/Explicit进展非线性动态分析的讨论,请参阅第9章“非线性显式动态分析。引言动态模拟是将惯性力包含在动力学平衡方程中:其中M结构的质量。结构的加速度。在结构中的内力。P所施加的外力。在上面公式中的表述是牛顿第二运动定律F=ma。在静态和动态分析之间最
2、主要的区别是在平衡方程中包含了惯性力M。在两类模拟之间的另一个区别在于内力的定义。在静态分析中,内力仅由结构的变形引起;而在动态分析中,内力包括源于运动例如阻尼和结构的变形的贡献。7.1.1固有频率和模态最简单的动态问题是在弹簧上的质量自由振动,如图7-1所示。图71 质量弹簧系统在弹簧中的内力给出为,所以它的动态运动方程为这个质量弹簧系统的固有频率natral frequency单位是弧度/秒rad/s给出为如果质量块被移动后再释放,它将以这个频率振动。假如以此频率施加一个动态外力,位移的幅度将剧烈增加,这种现象即所谓的共振。实际结构具有大量的固有频率。因此在设计结构时,非常重要的是防止使可
3、能的载荷频率过分接近于固有频率。通过考虑非加载结构在动平衡方程中令的动态响应可以确定固有频率。如此运动方程变为对于无阻尼系统,因此有这个方程的解具有形式为将此式代入运动方程,得到了特征值eigenvalue问题其中。该系统具有个特征值,其中是在有限元模型中的自由度数目。记是第个特征值;它的平方根是结构的第阶模态的固有频率natural frequency,而是相应的第阶特征向量eigenvector。特征向量也就是所谓的模态mode shape也称为振型,因为它是结构以第阶模态振动的变形形状。在ABAQUS/Standard中,应用频率的提取过程确定结构的振型和频率。这个过程应用起来十分容易,
4、你只要指出所需要的振型数目或所关心的最高频率即可。7.1.2振型叠加在线性问题中,可以应用结构的固有频率和振型来定性它在载荷作用下的动态响应。采用振型叠加modal superposition技术,通过结构的振型组合可以计算结构的变形,每一阶模态乘以一个标量因子。在模型中的位移矢量定义为其中是振型的标量因子。这一技术仅在模拟小变形、线弹性材料和无接触条件的情况下是有效的,换句话说,即线性问题。在结构的动力学问题中,结构的响应往往被相对较少的几阶振型控制,在计算这类系统的响应时,应用振型叠加成为特别有效的方法。考虑一个含有10,000个自由度的模型,对动态运动方程的直接积分将在每个时间点上同时需
5、要联立求解10,000个方程。如果通过100个振型来描述结构的响应,如此在每个时间增量步上只需求解100个方程。更重要的是,振型方程是解耦的,而原来的运动方程是耦合的。在计算振型和频率的过程中,开始时需要一点本钱,但是,在计算响应时将会节省大量的计算花费。如果在模拟中存在非线性,在分析中固有频率会发生明显的变化,因此振型叠加法将不再适用。在这种情况下,只能要求对动力平衡方程直接积分,它所花费的时间比振型分析昂贵得多。必须具备如下特点的问题才适合于进展线性瞬态动力分析: 系统应该是线性的:线性材料行为,无接触条件,以与没有非线性几何效应。 响应应该只受相对少数的频率支配。当在响应中频率的成分增加
6、时,诸如是打击和碰撞的问题,振型叠加技术的效率将会降低。 载荷的主要频率应该在所提取的频率X围之内,以确保对载荷的描述足够准确。 应用特征模态,应该准确地描述由于任何突然加载所产生的初始加速度。 系统的阻尼不能过大。阻尼如果允许一个无阻尼结构做自由振动,如此它的振幅会是一个常数。然而在实际中,能量被结构的运动耗散,振动的幅度减小直至振动停止。这种能量耗散被称为阻尼damping。通常假定阻尼为粘滞的,或者正比于速度。包含阻尼的动力平衡方程可以重新写为其中C是结构的阻尼矩阵是结构的速度。能量耗散来自于诸多因素,其中包括结构连接处的摩擦和局部材料的迟滞效应。阻尼是一种很方便的方法,它包含了重要的能
7、量吸收而又无需模拟具体的效果。在ABAQUS/Standard中,特征模态的计算是关于无阻尼系统的。然而,大多数工程问题都包含某种阻尼,尽管阻尼可能很小。对于每个模态,在有阻尼和无阻尼的固有频率之间的关系是其中是阻尼特征值,是临界阻尼比,是该振型的阻尼,是临界阻尼。对于的较小值,有阻尼系统的特征频率非常接近于无阻尼系统的相应值;当增大时,无阻尼系统的特征频率成为不太准确的;而当接近于1时,采用无阻尼系统的特征频率就成为无效的。如果结构是处于临界阻尼,在任何扰动后,结构不会有摆动而是尽可能迅速地恢复到它的初始静止构形。见图7-2图72 阻尼.1在ABAQUS/Standard中阻尼的定义对于瞬时
8、模态分析,在ABAQUS/Standard中可以定义一些不同类型的阻尼:直接模态阻尼direct modal damping,瑞利阻尼Rayleigh damping和复合模态阻尼posite modal damping。阻尼是针对模态动力学过程定义的,阻尼是分析步定义的一局部,每阶模态可以定义不同量值的阻尼。直接模态阻尼应用直接模态阻尼可以定义与每阶模态相关的临界阻尼比,其典型的取值X围是在临界阻尼的1%到10%之间。直接模态阻尼允许用户准确地定义系统的每阶模态的阻尼。Rayleigh阻尼在Rayleigh阻尼中,假设阻尼矩阵是质量和刚度矩阵的线性组合,其中和是由用户定义的常数。尽管阻尼是正
9、比于质量和刚度矩阵的假设没有严格的物理根底,实际上我们对于阻尼的分布知之甚少,也就不能保证其它更为复杂的模型是正确的。一般的,这个模型对于大阻尼系统不可靠;即超过临界阻尼的大约10%。相对于其它形式的阻尼,你可以准确地定义系统的每阶模态的Rayleigh阻尼。对于一个给定模态i,临界阻尼值为,而Rayleigh阻尼值和的关系为复合阻尼在复合阻尼中,对于每种材料定义一个临界阻尼比,这样就得到了对应于整体结构的复合阻尼值。当结构中有多种不同的材料时,这一选项是有用的。在本指南中将不对复合阻尼做进一步的讨论。.2 选择阻尼值在大多数线性动力学问题中,恰当地定义阻尼对于获得准确的结果是十分重要的。但是
10、,在某种意义上阻尼只是近似地模拟了结构吸收能量的特性,并非试图去模拟引起这种效果的物理机制。因此,在模拟中确定所需要的阻尼数据是很困难的。偶尔,你可以从动态试验中获得这些数据,但是,你不得不通过查阅参考资料或者经验获得这些数据。在这些情况下,你必须十分慎重地解释模拟结果,并通过参数分析研究来评估模拟对于阻尼值的敏感性。单元选择事实上,ABAQUS的所有单元均可用于动态分析,选取单元的一般原如此与静力分析一样。但是,在模拟冲击和爆炸载荷时,应该选用一阶单元,因为它们具有集中质量公式,这种公式模拟应力波的效果优于二阶单元采用的一致质量公式。动态问题的网格剖分当你正在设计应用于动态模拟的网格时,你需
11、要考虑在响应中将被激发的振型,并且使所采用的网格能够充分地反映出这些振型。这意味着能够满足静态模拟的网格,不一定能够计算由于加载激发的高频振型的动态响应。例如,考虑图7-3所示的板。一阶壳单元的网格对于板受均布载荷的静力分析是适合的,并也适合于一阶振型的预测。但是,该网格是明显地过于粗糙以至于不能够准确地模拟第六阶振型。图73板的粗网格图7-4显示了同样的板采用了一阶单元的精细网格的模拟。现在,第六阶振型的位移形状看起来明显变好,对于该阶振型所预测的频率更加准确。如果作用在板上的动态载荷会显著地激发该阶振型,如此必须采用精细的网格;采用粗网格将得不到准确的结果。图74 板的精细网格 例题:货物
12、吊车动态载荷这个例子采用在第节“例题:货物吊车中已分析过的同样的货物吊车,现在要求研究的问题是当10kN的载荷在0.2秒的时间中落到吊车挂钩上所引起的响应。在A, B, C和D点见图7-5处的连接仅能够承受的最大拉力为100kN。你必须决定这些连接的任何一个是否会断裂。图75 货物吊车加载的持续时间很短意味着惯性效应可能是很重要的,根本上要进展动态分析。这里没有提供关于结构的阻尼的任何信息。由于在桁架和交叉支撑之间采用的是螺栓连接,因此由摩擦效应引起的能量吸收可能是比拟显著的。因此,基于经验可以对每一阶振型选择5的临界阻尼。施加载荷的值与时间的关系,如图7-6所示。图76 载荷时间特性“Car
13、go crane dynamic loading提供了输入文件。当通过ABAQUS/CAE运行这个输入文件时,将创建关于该问题的完整的分析模型。根据下面给出的指导如果你遇到困难,或者如果你希望检查你的工作,如此可以运行这个输入文件。在附录A“Example Files中,给出了如何提取和运行输入文件的指导。如果你没有进入ABAQUS/CAE“Example:Cargo crane dynamic loading。.1 修改模型打开模型数据库文件,将Static模型复制成一个名为Dynamic的模型。除了下面描述的修改之外,动态分析的模型根本上与静力分析的模型一样。材料在动态分析中,必须给定每种
14、材料的密度,这样才能形成质量矩阵。在吊车中钢的密度为7800 kg/m3。在这个模型中,材料属性是作为截面特性定义的一局部给出的。所以需要在Property模块中编辑BracingSection和MainMemberSection截面定义来指定密度。在Edit Beam Section编辑梁截面对话框的Section material density截面材料密度域中,为每个截面输入密度值为7800。注意:如果材料数据的定义是独立于截面属性的,通过编辑材料定义可以将密度包括在内,即在Edit Material对话框中,选择General-Density。分析步应用于动态分析的分析步定义与静力分析
15、的分析步定义具有本质上的不同。因此,两个新的分析步将取代前面所建立的静力分析步。“非线性动态分析。ABAQUS/Standard提供了Lanczos和子空间迭代subspace iteration的特征值提取方法。对于具有很多自由度的系统,当要求大量的特征模态时,一般来说Lanczos方法的速度更快。当需要仅少数几个少于20特征模态时,如此应用子空间迭代法的速度可能更快。在这个分析中,我们使用Lanczos特征值求解器并求解前30个的特征值。除了指定所要提取模态的数目,也可以指定所感兴趣的最小和最大频率,因此,一旦ABAQUS/Standard已经提取了在这个指定X围内的所有特征值,就会完毕该
16、分析步。也可以指定一个变换点shift point,距离这个变换点最近的特征值将被提取。在默认情况下,不使用最小或最大的频率或变换点。如果没有约束结构的刚体模态,必须设置变换值为一个小的负值,以防止由于刚体运动产生的数值问题。采用频率提取分析步代替静态分析步:1从主菜单栏中,选择Step-Replace-Tip Load。在ReplaceStep替换分析步对话框中,从Linear perturbation线性摄动过程的列表中,选择Frequency频率。将删除不能转换的模型参数。在本例中删除了集中力,因为在频率提取分析中不能应用它们。但是,频率提取分析步继承了与静态分析步相关的边界条件和输出需
17、求。2在Edit Step编辑分析步对话框的Basic根底页中,输入分析步描述First 30 modes;承受Lanczos特征求解器选项;并要求前30阶特征值。3将分析步重新命名为Extract Frequencies。在结构动态分析中,响应通常地是与低阶模态有关。但是,应该提取足够的模态以便较好地表达结构的动态响应。检查是否已经提取了足够数量的特征值的一种方法是查看在每个自由度上的全部有效质量,它明确了在所提取模态的每个方向上激活了多少质量。在数据文件.dat的特征值输出中,给出了有效质量的列表。在理想情况下,对于每个振型在每个方向上,有效质量的总和应当至少占总质量的90“模态数目的影响
18、中将给出进一步的讨论。应用模型动态过程进展瞬时动态分析。瞬时响应是基于在第一个分析步中提取全部的模态;在全部的30阶模态中均采用了5%的临界阻尼。创建瞬时模型动态分析步:1 从主菜单栏中,选择Step-Create。从Linear perturbation过程列表中选择Modal dynamics,并命名分析步为Transient modal dynamics。在上面定义的频率提取分析步之后插入这个分析步。2 在Edit Step对话框的Basic页中,输入分析步的描述Simulation of Load Dropped on Crane,并指定分析步的时间为和时间增量time increme
19、nt为。在动态分析中,时间是一个真实的物理量。3 在Edit Step对话框的Damping页中,指定直接模态阻尼direct modal,并对第1阶至第30阶的模态输入临界阻尼比为。如果使用了模态阻尼,必须指定在基于模态的动态过程中使用的特征模态。ABAQUS/CAE默认自动地选择所有可能得到的特征模态。当然, 如果你希望改变默认的选择,也可以应用Keywords Editor关键词编辑器编辑*SELECT EIGENMODES块。在这个问题中,承受默认的选择。输出应用Field Output Request Manager场变量输出管理器,对于Extract Frequencies分析步,
20、修改场变量输出设置,因此,选择了Preselected defaults预选默认值。在默认情况下,ABAQUS/Standard将振型写入到输出数据库文件.odb,以便应用Visualization模块绘制振型图。对于每阶振型的节点位移都是经过单位化的,所以最大的位移为1单位。因此,这些结果和对应的应力和应变是没有物理意义的:它们仅能够用于相互的比拟。完成动态分析通常比静态分析需要更多更多的增量步。做为结果,来自动态分析的输出量可能是非常大,你应该控制输出要求以确保输出文件具有一个合理的量。在本例中,要求在每第5个增量步完毕时,向输出数据库文件中输出一次位移形状。在分析步中有100个增量步;因
21、此,有20组场变量输出。另外在每个增量步,将在模型加载端例如,Tip-a集合的位移和在固定端Attach集合的约束反力作为历史数据写入到输出数据库文件中,以便从这些数据中得到更好的解答。在动态分析中,我们也关心在模型中的能量分布以与能量采用的形式。在模型中表现出的动能是质量运动的结果;表现出的应变能是结构位移的结果;通过阻尼也耗散了能量。在默认情况下,对于模型动态过程,整个模型的能量将作为历史数据写入到.odb文件中。对瞬时模型动态分析步中的输出请求:1 从主菜单栏中,选择Output-Field Output Requests-Manager。在标记Transient modal dynam
22、ics的列中可能需要拉大这列表格才能看见完整的分析步名称,选择标有Created的单元格。2 编辑场变量输出要求,使得仅将每第5个增量步的节点的位移写入到.odb文件中。3 从主菜单栏中,选择Output-History Output Requests-Manager。在标记Transient modal dynamics的分析步中创建两个新的输出要求。在第一个中,输出集合Tip-a在每个增量步完毕时的位移;在第二个中,输出集合 Attach在每个增量步完毕时的约束反力。载荷和边界条件边界条件与在静力分析中的条件一样。由于在分析步替换过程中保存了这些条件,无需再定义新的边界条件。在吊车的端部施
23、加一个集中力,它的量级是与时间相关的,如图7-6所示。与时间相关的载荷可以应用幅值曲线amplitude curve进展定义。通过幅值曲线上的值乘以载荷的值10,000N,可以获得当时的任意点处施加荷载的实际值。指定与时间相关的载荷:1首先定义幅值。从在Load模块的主菜单栏中,选择Tools- Amplitude-Create,命名幅值为Bounce,并选择Tabular数据表类型。在Edit Amplitude编辑幅值对话框中,输入在表7-1中所示的数据。承受默认的Step time分析步时间的选择作为时间跨度,并输入作为光滑参数值。注意:点击鼠标键3,进入表格选项。表71 幅值曲线数据时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- abaqus 动态 分析 报告 方法
链接地址:https://www.desk33.com/p-6323.html