workbench瞬态动力分析.ppt

Click to edit Master text styles,,Second level,,Third level,,Fourth level,,Fifth level,,Training Manual,DYNAMICS 11.0,Training Manual,Click to edit Master text styles,,Second level,,Third level,,Fourth level,,Fifth level,,第四章瞬态动力分析,瞬态动力分析总论,定义：,,确定结构在任意随时间变化载荷作用下系统瞬态响应特性的技术输入数据：,,最一般形式是载荷为时间的任意函数；,,输出数据：,,随时间变化的位移和其它的导出量，如：应力和应变运动方程,基本运动方程,,,这是动力学最通常的方程形式，载荷可以是任意随时间变化的,.,,按照求解方法，,ANSYS,允许在瞬态动力分析中包括各种类型的非线性,——,,大变形、接触、塑性等等,.,,求解方法,求解运动方程,直接积分法,模态叠加法,隐式积分,显式积分,完整矩阵法,缩减矩阵法,完整矩阵法,缩减矩阵法,求解方法,两种求解运动学方程方法,:,,模态叠加法,,直接积分法,,运动方程可以直接对时间按步积分。

在每个时间点（,time = 0, Dt , 2Dt, 3Dt,….),，需求解一组联立的静态平衡方程（,F=ma）；,,需假定位移、速度和加速度是如何随时间而变化的，,(,积分方案选择,),,有多种不同的积分方案，如中心差分法，平均加速度法,, Houbolt, Wilson,Q,, Newmark,等,.,,求解方法,时间积分方案,–,,两种积分方案,Newmark,和,HHT.,缺省为,Newmark,,,,,,不同的,a,,和,d,,造成积分方案的变化,(,隐式,/,显式,/,平均加速度,).,,Newmark,是隐式积分方案,.,,ANSYS/LS-DYNA,利用显式积分方案,.,求解方法,时间积分方案,- HHT,方法,:,,Newmark,方法是求解,t,n+1,时刻的运动方程,HHT,方法,–,求解中间时间点的运动方程然后外推到,t,n+1.,,(Note:,缺省,HHT,方法,a,m,= 0,),求解方法,时间积分方案,-,时间积分参数,, γ, a,,d,, af, am,,通过求解控制选项输入,,,TRNOPT, FULL ,,, ,, NMK|HHT !,缺省,Newmark,,[TINTP,GAMMA,ALPHA,DELTA,THETA ,,, ,,, ALPHAF,ALPHAM],,,指定,GAMMA,或,ALPHAF/ALPHAM,0,<,,a,f,<,0.5,a,m,,<,,a,f,,求解方法,时间积分方案,–,为了稳定性与精度要求,,,下列关系需满足,.,(HHT,方法退化成,Newmark,当,a,f,与,a,m,＝,0,时）,HHT,法可以通过简单指定,GAMMA,值或指定,ALPHAF,与,ALPHAM,可以得到其他的方法,Hilber, Hughes and Taylor,,(HHT),Wood, Bossack and Zienkiewicz,Chung and Hulbert,缩减/完整结构矩阵,求解时既可用缩减结构矩阵，也可用完整结构矩阵；,,缩减矩阵：,,用于快速求解；,,不允许非线性因素存在,,根据主自由度写出[,K]、[C],和[,M],等矩阵，主自由度是完全自由度的子集；,,缩减的,[K],是精确的，但缩减的,[C],和 [,M],是近似的。

完整矩阵：,,不进行自由度缩减，采用完整的[,K]、[C],和[,M],矩阵；,,下面的讨论都是基于此种方法积分时间步长,积分时间步长（亦称为,ITS,或,,D,t,）是时间积分法中的一个重要概念,,ITS =,两个时刻点间的时间增量,D,t ；,,积分时间步长决定求解的精确度，因而其数值应仔细选取对于缩减矩阵法与模态叠加法瞬态分析,ANSYS,只允许,ITS,常值,.,,完全法瞬态分析，,ANSYS,可以自动调整时间步大小在用户指定的范围内,积分时间步长,ITS,小到足够获取下列动力学现象：,,响应频率,,载荷突变,,接触频率,,波传播效应,,响应频率,响应频率,,不同类型载荷激发系统不同的响应频率；,,ITS,小到足够获取所关心的最高响应频率（最低响应周期）；,,每个循环中有20个时刻点应是足够的，即：,,,D,t = 1/20f,,式中,，f,是所关心的最高响应频率响应周期,载荷突变,载荷突变,,ITS,小到足够获取载荷突变现象,Load,t,Load,t,接触频率,接触频率,,当两个物体发生接触，间隙或接触表面通常用刚度（间隙刚度）来描述；,,ITS,小到足够获取间隙,“,弹簧,”,频率；,,建议每个循环三十个点，才足以获取两物体间的动量传递。

更小的,ITS,会造成能量损失，并且冲击可能不是完全弹性的波传播,波传播,,由冲击引起在细长结构中更为显著（如下落时以一端着地的细棒）,,需要很小的,ITS ，,并且在波传播方向需要精细的网格,,显式积分法（在,ANSYS-LS/DYNA,采用）可能对此更为适用,非线性响应,非线性响应,,全瞬态分析可包括任何非线性类型,.,,更小的,ITS,通常有助于平衡迭代收敛,.,,塑性、蠕变及摩擦等非线性本质上是非保守的，需要精确地遵循载荷加载历程,.,小的,ITS,通常有助于精确跟踪载荷历程,.,,小的,ITS,可跟踪接触状态的变化,.,积分时间步长,如何选择,ITS?,,推荐打开自动时间步长选项,(AUTOTS),,并设置初始时间步长,D,t,initial,和最小时间步长,D,t,min,,、最大时间步长,D,t,max,. ANSYS,会利用自动时间步长功能来自动决定最佳时间步长,D,t.,,例如,:,如果,AUTOTS,是打开的， 并且,D,t,initial,= 1 sec,,D,t,min,= 0.01 sec, and,,D,t,max,= 10 sec;,那,ANSYS,起始采用,ITS= 1 sec,，并依据结构的响应允许其在,0.01,和,10,之间变动,.,AUTOTS,对于全瞬态分析缺省是打开的,.,对于缩减法和模态叠加法，是不可用的,.,,AUTOTS,会减小,ITS (,直到,D,t,min,),在下列情况,:,,在响应频率处，小于,20,个点,,求解发散,,求解需要大量的平衡迭代（收敛很慢）,,,塑性应变在一个时间步内累积超过,15%,,蠕变率超过,0.1,,如果接触状态要发生变化,(,决大多数接触单元可由,KEYOPT(7),控制,),积分时间步长,分析过程,讨论完全法瞬态分析过程,.,,五个主要步骤,:,,建立模型,,选择分析类型和选项,,指定边界条件和初始条件,,施加载荷历程并求解,,查看结果,,模型：,,所有的非线性因素可允许注意要求密度！,分析选项,进入求解阶段，并选择瞬态分析,.,,选择完全法,,求解选项,,阻尼,,求解方法,,完整矩阵方法为缺省方法。

允许下列非线性选项：,,大变形,,应力硬化,,Newton-Raphson,解法,,集中质量矩阵,,主要用于细长梁和薄壁壳或波的传播,,方程求解器,,由程序自行选择,分析选项,求解选项,,选择大位移瞬态分析或小变形瞬态分析,,.,,当不确定时，就选择大变形瞬态分析,,自动时间步长,(discussed next),指定载荷步结束时间,指定初始、最大、最小时间步长,D,t.,,输出控制,controls (discussed next),分析选项,自动时间步长,,在瞬态分析过程中，可自动计算正确的时间步长,.,,推荐激活该选项同时指定最大与最小积分步长,.,,如果有非线性因素，选择,,“,Program Chosen,”,选项,,注意: 在,ANSYS,中，总体求解器控制开关,[SOLCONTROL],的缺省状态为开， 建议保留这一状态，更为重要的是，不要在载荷步之间打开或关闭此开关,,,,,分析选项,输出控制,,用来控制写到结果文件的内容,.,,使用命令,OUTRES,,或选择,Solution > Sol’n Control.. > Basic,,通常的选项用来将每个子步的结果写到结果文件中去,.,,可光滑绘制结果与时间的关系曲线,.,,可能造成结果文件庞大,.,,分析选项,瞬态效应,on/off,,用来设置初始条件,,阶跃或渐进载荷,,指定阻尼,,使用缺省积分参数值,,,分析选项,阻尼,,,α,和,b,阻尼均可用；,,在大多数情况下，忽略,α,阻尼（粘性阻尼），仅指定,b,阻尼（由滞后造成的阻尼）：,,b,= 2,,/,w,,,式中,x,,为阻尼比，,w,,为主要响应频率,（rad/sec）。

典型命令:,,ALPHAD，…,,BETAD，…,分析选项,求解器选择,,缺省,ANSYS,选择稀疏求解器,,对于大自由度问题,(>100000 dofs),使用,PCG,法,初始条件,初始条件,,时间,t = 0,时的条件：,u,0 ，,v,0，,a,0,,它们的缺省值为,， u,0,= v,0,= a,0,= 0,,可能要求非零初始条件的实例：,,飞机着陆,(v,0,,0),,高尔夫球棒击球,(v,0,,0),,物体跌落试验,(a,0,,0),施加初始条件的两种方法,以静载荷步开始,,当只需在模型的一部分上施加初始条件时，例如，用强加的位移将悬臂梁的自由端从平衡位置,“,拨,”,开时，这种方法是有用的；,,用于需要施加非零初始加速度时,使用,IC,命令,,Solution > Apply > Initial Condit,’,n > Define +,,当需在整个物体上施加非零初始位移或速度时,IC,,命令法是有用的,零初始位移和零初始速度,,是缺省的初始条件，即如果,u,0,= v,0,= 0,,，,则不需要指定任何条件在第一个载荷步中可以加上对应于载荷—时间关系曲线的第一个拐角处的载荷。

非零初始位移及/或非零初始速度─可以用,IC,命令设置这些初始条件命令：,IC,,,GUI：MainMenu>Solution>-Loads-Apply>Initial Condit’n> Define,,零初始位移和非零初始速度,,非零速度是通过对结构中需指定速度的部分加上小时间间隔上的小位移来实现的比如如果,,v,0,=0.25，可以通过在时间间隔0.004内加上0.001的位移来实现，命令流如下：,,...,,TIMINT,OFF ! Time integration effects off,,D,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,,...,,非零初始位移和非零初始速度,和上面的情形相似，不过施加的位移是真实数值而非“小”数值。

比如，若,,u,0,,= 1.0且,v,0,= 2.5，则应当在时间间隔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,,...,非零初始位移和零初始速度,,需要用两个子步,[,NSUBST,2],来实现，所加位移在两个子步间是阶跃变化的[,KBC,1]。

如果位移不是阶跃变化的（或只用一个子步），所加位移将随时间变化，从而产生非零初速度下面的例子演示了如何施加初始条件,u,0,= 1.0，,v,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 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,,...,,非零初始加速度,,可以近似地通过在小的时间间隔内指定要加的加速度,[,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 constraints (if appropriate),,KBC,0 ! Ramped loads (if appropriate),,! Continue with normal transient solution procedures,,...,施加时间-历程载荷方法,施加时间-历程载荷,,时间- 历程载荷是随时间变化的载荷,,这类载荷有两种施加方法：,,函数工具,,表输入法,,多载荷步施加法,Load,t,Load,t,Load,t,函数法,函数工具,,允许施加复杂的边界条件,.,可通过函数编辑器,Solution > Define Loads > Apply > Functions > Define/Edit,,建议,:,,如果边界条件可直接用表格输入，不推荐采用函数法,.,,表输入法,表输入法,,允许定义载荷随时间变化的表（用数组参数）并采用此表作为载荷；,,尤其是在同时有几种不同的载荷，而每种载荷又都有它自己的时间历程时很方便；,,例如，要施加下图所示的力随时间变化曲线：,,1.,选择,,Solution > Apply > Force/Moment > On Nodes,，,然后拾取所需节点,0.5,Force,t,22.5,10,1.0,1.5,表输入法,2. 选择力方向和,,“,新表,New table,”,，,然后确定（,OK）；,,3.,输入表名和行数（时间点的数量），然后确定（,OK）；,,4.,填入时间和载荷值，然后,File > Apply/Quit；,表输入法,5. 规定终止时间和积分时间步长,,Solution > Time/Frequenc > Time - Time Step…,,不必指定载荷的分步或线性条件，这已包含在载荷曲线中,,6.,激活自动时间步，规定输出控制，然后求解（稍后讨论）,。

典型命令:,,TIME，… !,终点时间,,DELTIM，0.002，0.001，0.1 !,起始,最小和最大,ITS,,AUTOTS，ON,,OUTRES，…,,SOLVE,多载荷步法,多载荷步法,,允许将载荷－时间历程采用多个载荷步,.,,不需要数组参数,.,只需简单地施加每段载荷并求解或写成载荷步文件,(LSWRITE).,多载荷步法,例如,,,要施加如图力时间曲线,:,,1.,,需要三个载荷步,:,一个是上升渐进载荷,,,一个是下降渐进载荷,,,另一个是阶跃载荷,.,Force,t,22.5,10,0.5,1.0,1.5,2.,定义载荷步,1:,,在期望的节点施加力值,22.5,单位,.,,指定结束时间,(0.5),,积分时间步长和渐进载荷,.,,激活自动时间步长，指定输出控制，或求解或写成载荷步 文件,.,,,多载荷步法,3,.,定义载荷步,2:,,改变力值为,10.0.,,指定结束时间,(1.0).,不需指定积分时间步长或渐进载荷,.,,求解或写入载荷步文件,.,,,4.,定义载荷步,3:,,删除载荷或设置其值为,,指定结束时间,(1.5),和阶跃载荷,.,,求解或写入载荷步文件,求解,,利用,SOLVE,命令,(,或,LSSOLVE,,如果是写载荷步文件,).,,在每个时间步，,ANSYS,基于载荷时间曲线计算载荷值,.,,多载荷步法,多载荷步法,多载荷步过程文件：,,,,TIME,... ! Time at the end of 1st transient load step,,Loads ... ! Load values at above time,,KBC,... ! Stepped or ramped loads,,LSWRITE ! Write load data to load step file,,,TIME,... ! Time at the end of 2nd transient load step,,Loads ... ! Load values at above time,,KBC,... ! Stepped or ramped loads,,LSWRITE ! Write load data to load step file,,,TIME,... ! Time at the end of 3rd transient load step,,Loads ... ! Load values at above time,,KBC,... ! Stepped or ramped loads,,LSWRITE ! Write load data to load step file,,…,观察结果,由三步构成：,,绘制结构中某些特殊点的结果-时间曲线,,确定临界时间点,,察看在这些临界时间点处整个结构上的结果,,采用,POST26，,时间-历程后处理器,,采用,POST1，,通用后处理器,POST26：,绘制结果-时间曲线,绘制结果-时间曲线：,,首先定义,POST26,的变量,,节点或单元数据列表,,用一个,, 2,的编号来识别,,变量,1,含有各时间点，并且是预见定义了的,POST26：,绘制结果-时间曲线,定义变量,（,接上页）,,拾取那些变形最大的节点，然后选择自由度的方向,,更新变量定义表,POST26：,绘制结果-时间曲线,变量一旦已定义，就可以对它们绘图或列表,典型命令:,,/POST26,,NSOL，…,,PLVAR，...,POST26：,确定临界点,确定临界点,,采用各种极值列表（,List Extremes）,菜单,,记下发生最小和最大值时的时间点,典型命令,:,,EXTREM，…,,FINISH,,,POST1：,临界时间点时的总体结果,察看在各临界时间点处整个结构上的结果,,进入,POST1，“By Time/Freq...”,读出结果，并输入近似时间值,,绘制变形后的形状和应力等值线图,POST1：,临界时间点时的总体结果,POST1：,临界时间点时的总体结果,。