Abaqus 中显示动力学分析步骤.docx
- 文档编号:23818177
- 上传时间:2023-05-21
- 格式:DOCX
- 页数:17
- 大小:33.31KB
Abaqus 中显示动力学分析步骤.docx
《Abaqus 中显示动力学分析步骤.docx》由会员分享,可在线阅读,更多相关《Abaqus 中显示动力学分析步骤.docx(17页珍藏版)》请在冰豆网上搜索。
Abaqus中显示动力学分析步骤
准静态分析——ABAQUS/Explicit
准静态过程(guasi—staticprocess)
在过程进行的每一瞬间,系统都接近于平衡状态,以致在任意选取的短时间dt内,状态参量在整个系统的各部分都有确定的值,整个过程能够看成是由一系列极接近平衡的状态所构成,这种过程称为准静态过程。
无限缓慢地压缩和无限缓慢地膨胀过程可近似看作为准静态过程、准静态过程是一种理想过程,实际上是办不到的、
准静态原为一个热力学概念,在这个地方引用主要是指模型在加载的过程中任意时刻所经历的中间状态都可近似地视为静力状态,因此当加载过程进行得无限缓慢时,在各个时刻模型所处的状态就可近似地看作是静态,该过程便是准静态过程。
准静态啮合过程仿真主要考虑的是弧齿锥齿轮副在加载时的接触状态,以及齿面和齿根的应力变化规律,其前提是不考虑齿轮副惯性的影响。
ABAQUS/Explicit准静态分析
显式求解方法是一种真正的动态求解过程,它的最初发展是为了模拟高速冲击问题,在这类问题的求解中惯性发挥了主导性作用。
当求解动力平衡的状态时,非平衡力以应力波的形式在相邻的单元之间传播、由于最小稳定时间增量一般地是特别小的值,因此大多少问题需要大量的时间增量步。
在求解准静态问题上,显式求解方法差不多证明是有价值的,另外ABAQUS/Explicit在求解某些类型的静态问题方面比ABAQUS/Standard更容易。
在求解复杂的接触问题时,显式过程相关于隐式过程的一个优势是更加容易。
此外,当模型特别大时,显式过程比隐式过程需要较少的系统资源。
将显式动态过程应用于准静态问题需要一些特别的考虑。
依照定义,由于一个静态求解是一个长时间的求解过程,因此在其固有的时间尺度上分析模拟常常在计算上是不切合实际的,它将需要大量的小的时间增量。
因此,为了获得较经济的解答,必须采取一些方式来加速问题的模拟。
然而带来的问题是随着问题的加速,静态平衡的状态卷入了动态平衡的状态,在这个地方惯性力成为更加起主导作用的力。
目标是在保持惯性力的影响不显著的前提下用最短的时间进行模拟。
准静态(Quasi—static)分析也能够在ABAQUS/Standard中进行。
当惯性力能够忽略时,在ABAQUS/Standard中的准静态应力分析用来模拟含时间相关材料响应(蠕变、膨胀、粘弹性和双层粘塑性)的线性或非线性问题。
关于在ABAQUS/Standard中准静态分析的更多信息,请参阅ABAQUS分析用户手册(ABAQUSAnalysisUser'sManual)的第6、2。
5节“Quasi-staticanalysis”、
1、显式动态问题类比
假设两个载满了乘客的电梯、在缓慢的情况下,门打开后您步入电梯。
为了腾出空间,邻近门口的人慢慢地推他周围的人,这些被推的人再去推他周围的人,如此接着下去。
这种扰动在电梯中传播,直到靠近墙边的人表示他们无法移动为止。
一系列的波在电梯中传播,直到每个人都到达了一个新的平衡位置。
假如您稍稍加快速度,您会比前面更用力地推动您周围的人,然而最终每个人都会停留在与缓慢的情况下相同的位置。
在快速情况下,门打开后您以特别高的速度冲入电梯,电梯里的人没有时间挪动位置来重新安排他们自己以便容纳您。
您将会直截了当地撞伤在门口的两个人,而其他人则没有受到影响。
关于准静态分析,实际的道理是同样的。
分析的速度经常能够提高许多而可不能严重地降低准静态求解的质量;缓慢情况下和有一些加速情况下的的最终结果几乎是一致的。
然而,假如分析的速度增加到一个点,使得惯性影响占主导地位时,解答就会趋向于局部化,而且结果与准静态的结果是有一定区别的。
2。
加载速率
一个物理过程所占用的实际时间称其为它的固有时间(naturetime)。
关于一个准静态过程在固有时间中进行分析,我们一般能得到准确的静态结果、毕竟,假如实际事件真实地发生在其固有时间尺度内,并在结束时其速度为零,那么动态分析应该能够得到如此的事实,即分析实际上差不多达到了稳态。
您能够提高加载速率使相同的物理事件在较短的时间内发生,只要解答保持与真实的静态解答几乎相同,而且动态效果保持是不明显的、
2。
1光滑幅值曲线
关于准确和高效的准静态分析,要求施加的载荷尽估计地光滑、突然、急促的运动会产生应力波,它将导致振荡或不准确的结果。
以估计最光滑的方式施加载荷要求加速度从一个增量步到下一个增量步只能改变一个小量、假如加速度是光滑的,随其变化的速度和位移也是光滑的。
ABAQUS有一条简单、固定的光滑步骤(smoothstep)幅值曲线,它自动地创建一条光滑的载荷幅值。
当您定义一个光滑步骤幅值曲线时,ABAQUS自动地用曲线连接每一组数据对,该曲线的一阶和二阶导数是光滑的,在每一组数据点上,它的斜率都为零。
由于这些一阶和二阶导数都是光滑的,您能够采纳位移加载,应用一条光滑步骤幅值曲线,只用初始的和最终的数据点,而且中间的运动将是光滑的。
使用这种载荷幅值允许您进行准静态分析而可不能产生由于加载速率不连续引起的波动。
2、2结构问题
在静态分析中,结构的最低模态通常控制着结构的响应。
假如已知最低模态的频率和相应的周期,您能够估计出得到适当的静态响应所需要的时间。
为了说明如何确定适当的加载速率,考虑在汽车门上的一根梁被一个刚性圆环从侧面侵入的变形,实际的实验是准静态的、采纳不同的加载速率,梁的响应变化特别大。
以一个极高的碰撞速度为400m/s,在梁中的变形是高度局部化的。
为了得到一个更好的准静态解答,考虑最低阶的模态。
最低阶模态的频率大约为250Hz,它对应于4ms的周期、应用在ABAQUS/Standard中的特征频率提取过程能够容易地计算自然频率、为了使梁在4ms内发生所希望的0、2m的变形,圆环的速度为50m/s。
尽管50m/s好像仍然像是一个高速碰撞速度,而惯性力相关于整个结构的刚度差不多成为次要的了,变形形状显示了特别好的准静态响应、尽管整个结构的响应显示了我们所希望的准静态结果,但通常理想的是将加载时间增加到最低阶模态的周期的10倍以确保解答是真正的准静态。
为了更进一步地改进结果,刚环的速度估计会逐渐增大,例如应用一条光滑步骤幅值曲线,从而减缓初始的冲击、
2、3金属成形问题
为了获得低成本的求解过程,人为地提高成型问题的速度是必要的,然而,我们能够把速度提高多少仍能够获得可接受的静态解答呢?
假如薄金属板毛坯的变形对应于其最低阶模态的变形形状,能够应用最低阶结构模态的时间周期来指导成型的速度。
然而在成型过程中,刚性的冲模和冲头能够以如此的方式约束冲压,使坯件的变形估计与结构的模态无关。
在这种情况下,一般性的建议是限制冲头的速度小于1%的薄金属板的波速、关于典型的成型过程,冲头速度是在1m/s的量级上,而钢的波速大约为5000m/s。
因此依照这个建议,一个50的因数为冲头提高速度的上限。
为了确定一个可接受的冲压速度,建议的方法包括以各种变化的冲压速度运行一系列的分析,这些速度在3m/s至50m/s的范围内。
由于求解的时间与冲压的速度成反比,运行分析是以冲压速度从最快到最慢的顺序进行。
检查分析的结果,并感受变形形状、应力和应变是如何随冲压速度而改变的。
冲压速度过高的一些表现是与实际不符的、局部化的拉伸与变薄,以及对起皱的抑止。
假如您从一个冲压速度开始,例如50m/s,并从某处减速,在某点上从一个冲压速度到下一个冲压速度的解答将成为相似的,这说明解答开始收敛于一个准静态的解答。
当惯性的影响成为不明显时,在模拟结果之间的区别也是不明显的。
随着人为地增加加载速率,以逐渐和平滑的方式施加载荷成为越来越重要的方式。
例如,最简单的冲压加载方式是在整个成型过程中施加一个定常的速度。
在分析开始时,如此加载会对薄金属板坯引起突然的冲击载荷,在坯件中传播应力波并估计产生不希望的结果。
当加载速率增加时,任何冲击载荷对结果的影响将更加明显。
应用光滑步骤幅值曲线,使冲压速度从零逐渐增加能够使这些不利的影响最小化。
2。
4回弹
回弹经常是成型分析的一个重要部分,因为回弹分析决定了卸载后部件的最终形状。
尽管ABAQUS/Explicit十分适合于成型模拟,对回弹分析却遇到某些特别的困难。
在ABAQUS/Explicit中进行回弹模拟最主要的问题是需要大量的时间来获得稳态的结果。
特别是必须特别小心地卸载,同时必须引入阻尼以使得求解的时间比较合理。
幸运的是,在ABAQUS/Explicit和ABAQUS/Standard之间的紧密联系允许一种更有效的方法、
由于回弹过程不涉及接触,而且一般只包括中度的非线性,因此ABAQUS/Standard能够求解回弹问题,同时比ABAQUS/Explicit求解得更快。
因此,关于回弹分析更偏爱的方法是将完整的成型模型从ABAQUS/Explicit输入(import)到ABAQUS/Standard中进行。
3、质量放大
质量放大(massscaling)能够在不需要人为提高加载速率的情况下降低运算的成本。
关于含有率相关材料或率相关阻尼(如减震器)的问题,质量放大是惟一能够节约求解时间的选择。
在这种模拟中,不要选择提高加载速度,因为材料的应变率会与加载速率同比例增加。
当模型的参数随应变率变化时,人为地提高加载速率会人为地改变了分析的过程。
人为地将材料密度增加因数倍,则波速就会降低因数f倍,从而稳定时间增量将提高因数f倍。
注意到当全局的稳定极限增加时,进行同样的分析所需要的增量步就会减少,而这正是质量放大的目的、然而,放大质量对惯性效果与人为地提高加载速率恰好具有相同的影响。
因此,过度地质量放大,正像过度地加载速率,估计导致错误的结果。
为了确定一个可接受的质量放大因数,所建议的方法类似于确定一个可接受的加载速率放大因数。
两种方法的唯一区别是与质量放大相关的加速因子是质量放大因数的平方根,而与加载速率放大相关的加速因子是与加载速率放大因数成正比。
例如,一个为100倍的质量放大因数恰好对应于10倍的加载速率因数。
通过使用固定的或可变的质量放大,能够有多种方法来实现质量放大编程、质量放大的定义也能够随着分析步而改变,允许有特别大的灵活性。
详细的内容请参阅ABAQUAS分析用户手册第7、15、1节“Massscaling”。
4。
能量平衡
评估模拟是否产生了正确的准静态响应,最具有普遍意义的方式是研究模型中的各种能量、下面是在ABAQUS/Explicit中的能量平衡方程:
Etotal=EI+EV+EKE+EFD+EW
式中,EI是内能(包括弹性和塑性应变能),EV是粘性耗散吸收的能量,EKE是动能,EFD是摩擦耗散吸收的能量,EW是外力所做的功,Etotal是在系统中的总能量。
假如模拟是准静态的,那么外力所做的功是几乎等于系统内部的能量、除非有粘弹性材料、离散的减震器、或者使用了材料阻尼,否则粘性耗散能量一般地是特别小的。
由于在模型中材料的速度特别小,因此在准静态过程中,我们差不多确定惯性力能够忽略不计、由这两个条件能够推论,动能也是特别小的。
作为一般性的规律,在大多数过程中,变形材料的动能将可不能超过它的内能的一个小的比例(典型的为5%到10%)。
当比较能量时,请注意ABAQUS/Explicit报告的是整体的能量平衡,它包括了任何含有质量的刚体的动能。
由于当评价结果时我们只对变形体感兴趣,当评价能量平衡时我们应在Etotal中扣除刚体的动能。
例如,假如您正在模拟一个采纳滚动刚体模具的传输问题,刚体的动能估计占据模型整个动能的特别大部分。
在这种情况下,您必须扣除与刚体运动有关的动能,然而才估计做出与内能有意义的比较。
5、例题:
ABAQUS/Explicit凹槽成型
修改由ABAQUS/Standard分析所创建的模型,如此才能在ABAQUS/Explicit中运行它。
这些修改包括在材料模型中增加密度,改变单元库,并改变分析步。
为了获得正确的准静态响应,在运行ABAQUS/Explicit分析前,您将应用在ABAQUS/Standard的频率提取过程来确定所需要的计算时间。
5、1前处理——应用ABAQUS/Explicit重新运算模型
关于一个准静态过程,假如我们明白了坯件的最低阶固有频率,即基(fundamental)频,我们就能够确定分析步时间的一个大致的下限。
一种获得这个信息的方法是在ABAQUS/Standard中运行频率分析。
在这个成型分析中,冲压对坯件产生的变形类似于它的最低阶模态。
因此,假如您想模拟整个结构而并非局部的变形,选择第一个成型时期的时间是大于或等于坯件最低阶模态的周期是十分重要的。
运行一个固有频率提取过程:
①。
将已存在的模型复制成为一个新的模型,命名为Frequency,并对Fequency模型进行如下全面的修改:
在频率提取分析中,您将用一个单独的频率提取分析步取代现在所有的分析步、此外,您将删除所有的刚性工具和接触相互作用;它们与确定毛坯的基频无关、
②、在Property模块中,为Steel材料模型增加一个7800的密度、
③、在Assembly模块中,删除冲模、冲头和夹具部件的实体。
关于频率分析并不需要这些刚体部件。
(提示:
您能够从工具箱中采纳Delete工具删除这些部件。
)
④、进入Step模块,用一个单独的频率提取分析步替代现存的所有分析步。
a。
在StepManager(分析步管理器)中,删除分析步RemoveRight Constraint、HolderForce、EstablishContactII和MoveRunch、
b、选择分析步Establish ContactI,并点击Replace、
c、在RepalceStep(替换分析步)对话框中,从LinearParturbation过程列表中选择Frequency,键入分析步描述为Frequencymodes;选择Lanczos特征值选项,并要求五个特征值。
重新命名分析步为ExtractFrequencies。
d、取消DOF Monitor(自由度监视器)选项、
(注意:
由于频率提取分析步是一个线性扰动过程,将忽略材料的非线性性质、在这个分析中,坯件的左端约束沿x-方向的位移和绕法线的转动;然而,没有约束沿y-方向的位移、因此,提取的第一阶模态将是刚体模态、关于在ABAQUS/Explicit中的准静态分析,第二阶模态的频率将确定合适的时间段。
)
⑤、在Interaction模块,删除所有的接触相互作用、
⑥。
进入Load模块,在BCManager(边界条件管理器)中检查在ExtractFrequencies分析步中的边界条件。
除了边界条件名称CenterBC以外,删除所有的边界条件。
将这个留下的采纳了对称边界条件的毛坯约束施加到左端。
⑦、在创建和提交作业前,假如有必要则重新剖分网格、
⑧。
进入Job模块,创建一个作业,命名为Forming—Frequency,采纳如下的作业描述:
Channelforming–- frequencyanalysis。
提交作业进行分析,并监控求解过程、
⑨、当分析完成时,进入Visualization模块,并打开由这个作业创建的输出数据库文件。
从主菜单栏中,选择Plot-—>DeformedShape;或者应用在工具箱中的工具。
绘制出一阶屈曲模态的模型变形形状。
进一步绘出毛坯的二阶模态,将未变形的模型形状叠加在模型变形图上、
频率分析表明坯件有一个140Hz的基频,对应的周期为0。
00714s、关于成型分析,我们现在明白最短的分析步时间为0。
00714s。
创建ABAQUS/Explicit成型分析
成型过程的目标是采纳0、03m的冲头位移准静态地成型一个凹槽。
在选择准静态分析的加载速率时,建议您在开始时用较快的加载速率,并依照需要减小加载速率,更快地收敛到一个准静态解答、然而,假如您希望在您的第一次分析尝试中就增加能够得到准静态结果的估计性,您应当考虑分析步时间是比相应的基频缓慢10到50倍的因数。
在这个分析中,关于成型分析步,您将从0、007s的时间开始。
这是基于在ABAQUS/Standard中进行的频率分析,它显示出毛坯具有140Hz的基频,对应于0。
00714s的时间周期。
这个时间周期对应于4、3m/s的常数冲头速度。
您将认真地检查动能和内能的结果,以检验结果中并没有包含显著的动态影响。
将Standard模型复制成一个新模型,命名为Explicit、假如必要,通过从位于工具栏下方的Model(模型)列表中选择Explicit模型作为当前的模型。
使所有接下来的模型改变成为Explicit模型。
在ABAQUS/Standard分析中,在冲头和坯件之间模拟一个初始的缝隙以便于接触计算、在ABAQUS/Explicit分析中则不需要采取这种预防措施。
因此,在Assembly模块中,沿U2方向平移冲头-0。
001m。
在警告对话框中出现的关于相对和绝对约束中,点击Yes、
在毛坯夹具上施加一个集中力,为了计算夹具的动态反应,必须在刚性体的参考点上赋予一个点质量、夹具的实际质量是不重要的;而重要的是它的质量必须与毛坯的质量(0、78kg)具有同一个数量级,以使在接触计算中的振荡最小化。
选择数值为0、1kg的点质量。
在Property模块中,创建一个点的截面定义,命名为Pointmass。
在EditSection对话框的InertialProperties域中,键入0。
1点质量的值、在参考点RigidRefHolder应用这个截面定义、此外,编辑Steel材料定义来包括7800 kg/m3的质量密度、
进入Step模块、您需要为ABAQUS/Explicit分析创建两个分析步。
在第一个分析步中施加夹具力;在第二个分析步中施加冲头压下力、除了命名为Establish ContactI的分析步之外,删除所有其他的分析步,并用一个单一的显式动态分析步替换这个分析步。
键入分析步描述为Apply holderforce,并指定0、0001s的分析步时间。
这个时间关于施加夹具载荷是适合的,因为它是足够长以幸免了动态效果,而且又足够短以防止了对整个作业运行时间的明显冲击。
将分析步重新命名为Holderforce。
创建第二个显式动态分析步,命名为Displacepunch,分析步的时间为0。
007s,键入Applypunchstroke作为分析步的描述、
为了帮助确定分析是如何接近于准静态假设,研究各种能量的历史是特别有用的、特别有用的是比较动能和内部应变能。
能量历史默认地写入了输出数据库文件、
在这个金属成型分析的第一次尝试中,关于施加的夹具力和冲头压力,您将应用具有默认的光滑参数的表格形式的幅值曲线、进入Load模块,为施加的夹具力创建一个名为Ramp1的表格形式的幅值曲线。
在表1中输入幅值数据。
为冲头压力定义第二个表格形式的幅值曲线,命名为Ramp2。
在表2中输入幅值数据。
在LoadManager(载荷管理器)中,在命名为Holderforce的分析步中创建一个集中力,命名为RefHolderForce,在施加的点上指定RefHolder和一个沿着CF2方向大小为—440000的力。
关于这个载荷,改变幅值定义为Ramp1。
在BoudaryConditionManager(边界条件管理器)中,删除命名为MidLeftBC和MidRightBC的边界条件。
编辑RefDieBC边界条件,如此在Holder force分析步中沿着U2方向的约束为零,不改变其他方向的约束。
关于RefHolderBC边界条件,解除沿着U2方向的约束,而其他方向的约束保持不变。
在Displace Punch分析步中,改变位移边界条件RefPunchBC,使沿着U2方向的位移为-0、03m。
关于这个边界条件,应用幅值曲线Ramp2。
监视自由度的值。
在这个模型中,您将在整个分析步中监视冲头的参考节点的竖向位移(自由度2)。
在ABAQUS/Standard成型分析中,由于差不多设置了DOFMonitor监视RefPunch的竖向位移,因此您无需做出任何改变。
创建网格和定义作业。
在网格Mesh模块中,将用于剖分坯件网格的单元族改变为Explicit,并指定增强沙漏控制,并剖分坯件网格。
因为差不多将工具模拟成了解析刚性表面,因此无需将它们剖分网格。
在Job模块中创建一个作业,命名为Forming-1,给予作业如下的描述:
Channelforming—-attempt1、
在运行成型分析前,您估计希望明白该分析将需要多少个增量步,进而了解该分析需要多少计算机时间、您能够通过运行数据检查(datacheck)分析来获得关于初始稳定时间增量的近似值。
在这个例题中,从一个增量步到下一个增量步的稳定时间增量可不能有太大的变化,因此明白了稳定时间增量,您能够确定完成成型时期的分析需要多少个增量步。
一旦分析开始,您就能够明白每一个增量步需要多少CPU时间,进而明白整个分析需要多少CPU时间。
将模型保存到模型数据库文件中,并提交作业进行分析、监视求解过程;改正任何检测到的模拟错误,并调查任何警告信息的原因。
完成整个分析估计需要运行10分钟或更长的时间、
一旦分析开始运行,在另一个视图窗中会显示出您选择来监视(冲头的竖向位移)的自由度值的X—Y曲线图。
从主菜单栏中,选择Viewport--〉JobMonitor:
Forming-1,在分析运行的整个时间中跟踪沿着2—方向冲头位移的发展进程、
评价结果的策略。
在查看我们最关怀的结果之前,诸如应力和变形形状,我们需要确定结果是否是准静态的。
一个好的方法是比较动能与内能的历史、在金属成型分析中,大部分的内能是由于塑性变形产生的。
在这个模型中,坯件是动能的主要因素(忽略夹具的运动,没有与冲头和模具相关的质量)。
为了确定是否差不多获得了一个可接受的准静态解答,坯件的动能应该小于其内能的几个百分点。
关于更高的精确度,特别地是对回弹应力感兴趣时,动能应该是更低的。
这个方法是特别有用的,因为它应用于所有类型的金属成型过程,而且不需要任何直观地理解在模型中的应力;许多成型过程估计是过于复杂,以至于不允许对结果有一个直观的判断。
尽管是衡量准静态分析的良好和重要的证明,仅凭动能与内能的比值还不足以确任解的质量。
您还必须对这两种能量进行独立地评估,以确定它们是否是合理的、当需要准确的回弹应力结果时,这一部分的评估是更增加了重要性,因为一个高度精确的回弹应力解答是高度地依赖于准确的塑性结果。
即使动能是特别小的量,假如它包含了高度的振荡,则模型也会经历显著的塑性。
一般说来,我们希望光滑加载以产生光滑的结果;假如加载是光滑的,然而能量的结果是振荡的,则结果估计是不合适的。
由于一个能量的比值无法显示这种行为,因此您也必须研究动能本身的历史以观察是否是光滑的依然振荡的。
假如动能不能显示出准静态的行为,在某些节点上观察速度的历史估计是有用的,以帮助理解在各个区域中模型的行为。
这种速度历史能够表明在模型的哪些区域是振荡的,并产生大量的动能。
评估结果。
进入Visulization模块,并打开由这个作业(Forming-1。
odb)创建的输出数据库。
绘制动能和内能。
创建能量历史的曲线:
①、从主菜单栏中,选择Plot-->HistoryOutput、显示出整个模型的伪应变能历史曲线
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Abaqus 中显示动力学分析步骤 显示 动力学 分析 步骤