有限元仿真分析动力学-explicit总结
- 格式:docx
- 大小:12.20 KB
- 文档页数:4
有限元动⼒学分析知识点复习⽬录⼀、模型输⼊、建模A 输⼊⼏何模型1、两种⽅法:No defeaturing 和 defeaturing(Merge合并选项、Solid实体选项、Small选项)2、产品接⼝。
输⼊IGES ⽂件的⽅法虽然很好,但是双重转换过程CAD > IGES > ANSYS 在很多情况下并不能实现100%的转换.ANSYS 的产品接⼝直接读⼊“原始”的CAD ⽂件,解决了上⾯提到的问题.3、输⼊有限元模型。
除了实体⼏何模型外, ANSYS 也可输⼊由某些软件包⽣成的有限元单元模型数据(节点和单元)。
B 实体建模1、定义实体建模:建⽴实体模型的过程。
(两种途径)1)⾃上⽽下建模:⾸先建⽴体(或⾯),对这些体或⾯按⼀定规则组合得到最终需要的形状.开始建⽴的体或⾯称为图元.⼯作平⾯⽤来定位并帮助⽣成图元.对原始体组合形成最终形状的过程称为布尔运算1第1页总体直⾓坐标系 [csys,0] 总体柱坐标系[csys,1]总体球坐标系[csys,2] ⼯作平⾯ [csys,4]2)⾃下⽽上建模:按照从点到线,从线到⾯,从⾯到体的顺序建⽴模型。
B ⽹格划分1、⽹格划分三步骤:定义单元属性、指定⽹格的控制参数、⽣成⽹格2、单元属性(单元类型 (TYPE)、实常数 (REAL)、材料特性 (MAT))3、单元类型单元类型是⼀个重要选项,它决定如下单元特性:⾃由度(DOF)设置、单元形状、维数、假设的位移形函数。
1)线单元(梁单元、杆单元、弹簧单元)2)壳⽤来模拟平⾯或曲⾯。
3)⼆维实体⽤于模拟实体截⾯4)三维实体2第2页⽤于⼏何属性,材料属性,荷载或分析要求考虑细节,⽽⽆法采⽤更简单的单元进⾏建模的结构。
也⽤于从三维CAD系统转化⽽来的⼏何模型,⽽这些⼏何模型转化成⼆维模型或壳体会花费⼤量的时间和精⼒4、单元阶次与形函数单元阶次是指单元形函数的多项式阶次。
什么是形函数–形函数是指给出单元内结果形态的数值函数。
有限元动力学分析知识点复习目录一、模型输入、建模A 输入几何模型1、两种方法:No defeaturing 和 defeaturing(Merge合并选项、Solid实体选项、Small选项)2、产品接口。
输入IGES 文件的方法虽然很好,但是双重转换过程CAD > IGES > ANSYS 在很多情况下并不能实现100%的转换.ANSYS 的产品接口直接读入“原始”的CAD 文件,解决了上面提到的问题.3、输入有限元模型。
除了实体几何模型外, ANSYS 也可输入由某些软件包生成的有限元单元模型数据(节点和单元)。
B 实体建模1、定义实体建模:建立实体模型的过程。
(两种途径)1)自上而下建模:首先建立体(或面),对这些体或面按一定规则组合得到最终需要的形状.✓开始建立的体或面称为图元.✓工作平面用来定位并帮助生成图元.✓对原始体组合形成最终形状的过程称为布尔运算✓总体直角坐标系 [csys,0] 总体柱坐标系[csys,1]总体球坐标系[csys,2] 工作平面 [csys,4]2)自下而上建模:按照从点到线,从线到面,从面到体的顺序建立模型。
B 网格划分1、网格划分三步骤:定义单元属性、指定网格的控制参数、生成网格2、单元属性(单元类型 (TYPE)、实常数 (REAL)、材料特性(MAT))3、单元类型单元类型是一个重要选项,它决定如下单元特性:自由度(DOF)设置、单元形状、维数、假设的位移形函数。
1)线单元(梁单元、杆单元、弹簧单元)2)壳用来模拟平面或曲面。
3)二维实体用于模拟实体截面4)三维实体✓用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
✓也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数?–形函数是指给出单元内结果形态的数值函数。
ANSYS显示动力学分析实例ANSYS有限元仿真3月10日1410在仿真过程中遇到瞬态大变形,材料破坏失效等情况下可以借助ANSYS 的显示动力学分析来解决。
ANSYS显示动力学模块包括三种:Explicit Dynamics、ANSYS AUTODYN、ANSYS L S-DYNA。
本期通过一个实例来简单介绍下这三个模块的具体操作。
实例问题描述:一个金属圆柱体快速穿透金属板。
求解穿透过程中的最大应力和穿透方向的变形。
通过用不同仿真模块计算并比较仿真结果。
分析流程图如下。
其中A、B、C分别对应上面提到的三个模块。
这三个模块建立了数据共享,可避免重复的前处理操作,便于提高仿真效率。
分析树如下:1.Explicit Dynamics材料添加和几何建模略过...加载情况:固定约束金属板八条边、金属圆柱体运动速度300m/s。
注意分析设置Analysis setting 中的最大循环次数Maximum number of cycle和结束时间End Time应设置合理,不宜过大。
过大容易导致计算时间过长。
等效应力和变形求解结果如下图:最大等效应力为4.9e8,Z轴方向的最大变形为20.52m。
2.AUTODYNANSYS AUTODYN软件它有不同于Explicit Dynamics的交互式图形界面。
如下图所示主界面。
在AUTODYN软件中不需要再做其他前处理了!因为已经和Explicit Dynamics建立数据共享,只需要你轻轻点击RUN即可!这就是流程式分析的优点,大大的减少了工作量。
下面是后处理:求取应力数据。
按照图中步骤1.选择绘制云图类型contour 2.调出绘图变量contour variable 对话框 3.点击对于变量 4.勾选。
求取变形云图同理。
仿真结果截图如下:最大等效应力为3.8e8,Z轴方向的最大变形为19.88m。
用两个模块分析比较来看,等效应力相差稍微较大,而变形相差很小。
机械系统的动力学特性分析与仿真在现代工程领域中,机械系统的动力学特性分析与仿真是一项重要的技术。
它可以帮助工程师们更好地了解和预测机械系统的运动行为,为设计和优化机械系统提供可靠的依据。
本文将从理论与实践两个方面介绍机械系统的动力学特性分析与仿真。
一、动力学特性分析机械系统的动力学特性包括质量、惯性、刚度、阻尼等。
这些特性能直接影响机械系统的运动响应和稳定性。
在动力学特性分析中,常用的方法有质量矩阵法、阻尼矩阵法和刚度矩阵法等。
质量矩阵法利用质量矩阵描述机械系统各个部分的质量分布情况,并通过矩阵运算得到系统的动力学方程。
通过分析质量矩阵可以得知机械系统的质量分布情况,为设计优化提供依据。
阻尼矩阵法则通过对系统进行阻尼特性分析,得到系统的阻尼矩阵。
阻尼矩阵可以反映机械系统的阻尼分布和阻尼能力,对减少系统振动与噪音具有重要作用。
刚度矩阵法通过分析机械系统的刚度分布情况,得到系统的刚度矩阵。
刚度矩阵能反映机械系统的刚度分布和变形特性,为系统的优化设计提供依据。
二、动力学仿真动力学仿真是通过计算机建立机械系统的数学模型,并利用数值计算方法求解动力学方程,从而模拟机械系统的运动行为。
动力学仿真可以有效地预测机械系统的响应和稳定性,为系统的设计和优化提供指导。
在动力学仿真中,常用的方法有多体系统仿真和有限元分析。
多体系统仿真是通过建立各个部件之间的动力学联系,构建机械系统的数学模型。
通过求解模型的动力学方程,可以得到系统的运动轨迹和响应。
多体系统仿真在车辆动力学、机械振动与噪声等领域得到广泛应用。
有限元分析将机械系统分割成有限个单元,每个单元具有特定的材料和几何性质。
通过求解单元之间的力平衡方程,可以得到机械系统的运动行为。
有限元分析在结构强度、疲劳分析等方面具有重要应用。
三、实例分析以汽车悬挂系统为例,介绍动力学特性分析与仿真的应用。
汽车悬挂系统通过减震器和弹簧等组件,为车身提供舒适的行驶环境。
在悬挂系统的设计过程中,需要对系统的动力学特性进行分析与仿真。
复习目录一、模型输入、建模A 输入几何模型1、两种方法:No defeaturing 和 defeaturing(Merge合并选项、Solid实体选项、Small选项)2、产品接口。
输入IGES 文件的方法虽然很好,但是双重转换过程CAD > IGES > ANSYS 在很多情况下并不能实现100%的转换.ANSYS 的产品接口直接读入“原始”的CAD 文件,解决了上面提到的问题.3、输入有限元模型。
除了实体几何模型外, ANSYS 也可输入由某些软件包生成的有限元单元模型数据(节点和单元)。
B 实体建模1、定义实体建模:建立实体模型的过程。
(两种途径)1)自上而下建模:首先建立体(或面),对这些体或面按一定规则组合得到最终需要的形状.开始建立的体或面称为图元.工作平面用来定位并帮助生成图元.对原始体组合形成最终形状的过程称为布尔运算1第1页总体直角坐标系 [csys,0] 总体柱坐标系[csys,1]总体球坐标系[csys,2] 工作平面 [csys,4]2)自下而上建模:按照从点到线,从线到面,从面到体的顺序建立模型。
B 网格划分1、网格划分三步骤:定义单元属性、指定网格的控制参数、生成网格2、单元属性(单元类型 (TYPE)、实常数 (REAL)、材料特性 (MAT))3、单元类型单元类型是一个重要选项,它决定如下单元特性:自由度(DOF)设置、单元形状、维数、假设的位移形函数。
1)线单元(梁单元、杆单元、弹簧单元)2)壳用来模拟平面或曲面。
3)二维实体用于模拟实体截面4)三维实体2第2页用于几何属性,材料属性,荷载或分析要求考虑细节,而无法采用更简单的单元进行建模的结构。
也用于从三维CAD系统转化而来的几何模型,而这些几何模型转化成二维模型或壳体会花费大量的时间和精力4、单元阶次与形函数•单元阶次是指单元形函数的多项式阶次。
•什么是形函数–形函数是指给出单元内结果形态的数值函数。
机械设计中的仿真和模拟技术在机械设计领域中,仿真和模拟技术是非常重要的工具和方法。
借助这些技术,设计师可以在计算机中建立虚拟的模型,以模拟和预测机械系统的性能、行为和性质。
本文将介绍机械设计中常用的仿真和模拟技术,并分析其在实际应用中的优势和挑战。
一、仿真技术1. 数值仿真数值仿真是机械设计中常用的一种仿真技术。
通过建立数学模型和运用数值计算方法,可以模拟机械系统的运动、变形、热力学特性等。
数值仿真广泛应用于机械结构强度分析、流体力学仿真、热传导分析等方面。
借助数值仿真,设计师可以快速了解机械系统的行为,并优化设计方案,减少试验成本。
2. 有限元分析有限元分析是机械设计中一种常见的仿真技术。
通过将复杂的连续体分割成有限个单元,并建立节点之间的关系,可以在计算机中计算出结构的应力、变形、振动等情况。
有限元分析广泛应用于结构强度分析、材料力学性能评估、振动和声学分析等领域。
它可以快速评估设计的可行性,并指导设计的优化和改进。
3. 多体动力学仿真多体动力学仿真是研究机械系统运动学和动力学行为的技术。
通过建立机械系统各个零部件之间的连接关系和力学特性,可以模拟机械系统的运动规律、力学特性和能量传递情况。
多体动力学仿真广泛应用于机械系统的运动学分析、轨迹规划、运动控制等方面。
它可以帮助设计师更好地理解机械系统的工作原理,提高设计的准确性和可靠性。
二、模拟技术1. 动力学模拟动力学模拟是机械设计中的一种重要模拟技术。
借助物理和数学模型,可以模拟机械系统在不同工况下的运动和力学特性。
通过动力学模拟,设计师可以预测机械系统在不同负载下的响应,分析系统的稳定性、振动特性等。
动力学模拟广泛应用于机械系统的动态性能评估、操纵性分析、碰撞仿真和可靠性评估等方面。
2. 流体力学模拟流体力学模拟是一种模拟和预测流体流动和传热行为的技术。
借助流体力学模拟,设计师可以研究液体和气体在不同流动条件下的行为、压力分布和热传导情况。
动力学-abaqus/explict 总结
动力学分为:线性动力学和非线性动力学。
Standard 适合模拟与模型的振动频率相比响应 周期较长的问题:explicit :适合于模拟高速动力学问题。
线性动力学在abaqus/standard 中求解,是基丁•模态的分析方法。
应用有:
模态动力学:在时域内计算结构的线性动力学响应;可以使用直接积分
稳态动力学:计算由谐波激励引起的动态响应,可以使用直接积分。
响应谱分析:计算运动过程中的峰值响应;
随即响应分析:计算随即连续激励的响应,如地震波。
非线性动力学:需要对运动方程进行直接积分;abaqus/standard 中使用newmark 积分方 法,是隐式非线性直接积分法(无条件稳定,可以使用任意的时间增量,并11解仍然是有 界的)。
Abaqus/explicit 使用二阶精度的中心差分法(该方法是条件稳定的,只有在时间 增量小『•一定的临界值时才能给出有界的解)。
下面对explicit 使用过程中的一些细节作 简要的总结。
1. Abaqus/explicit :提供两种方案定义接触:
1.1 General contact:通用接触。
一般在模型中存在多个部件或复杂的拓扑结构情况下使用, 该功能强大,不需像在abaqus/standard 一样定义相互作用的接触对,在abaqus/explicit 里 会口动搜索相互作用的接触。
Examples
The following input specifies that the contact domain is based on self-contact of an all- inclusive, automatically generated surface but that contact (including self-contact in any overlap regions) should be ignored between the all-inclusive, automatically generated surface and surface CONTACT
CONTACT INCLUSIONS, ALL EXTERIOR 或ALL ELEMENT BASED
*CONTACT PROPERTY ASSIGNMENT
/,prop_l *alum__surf, steel__surf,
prop__2 *alum__sur f,
alum^sur f, prop_3 ★CONTACT
EXCLUSIONS
,surface^2 Either of the following methods can be used to exclude self-contact for surface_l from the contact domain:
(以全局的方式重新制定属
性) (局部修改) (局部修改)
★CONTACT EXCLUSIONS
surface 1,
or
^CONTACT1 EXCLUSIONS
surface^l, surface^l
1.2.接触问题中调整初始节点位置
Abaqus/explicit不允许接触表面的初始过盈。
在分析之前,接触面上的节点将被自动调整,以删除任何的初始过盈(不产生应变和应力〉。
在随后的分析步中,这样的调整将引起应变。
大量的调整将引起初始单元形状的严重扭曲。
如果发生这样的错误信息,最好运行datacheck 分析,并利用abaqus/viewer 查找问题。
即运行:abaqusjob=? Datacheck
2. 在大多数的abaqus/explicit分析中,用户只需要供总时间步长(该时间步长是有实际意义的)。
时间增量的大小由程序口动确定,这样可以满足稳定时间增量的限制。
对于特殊情况,可以人为的控制时间增量的大小。
一般t = Al/c,其中山为单元最小尺寸,c为波速=7E7P,E为弹性模量,P为密度。
显示动力学过程需要使用单元的集中质量矩阵。
因为使用abaqus/explicit时求解效率一般是非常重要的因素,所以可用的单元只有一阶减缩积分单元。
例外:修正的三角形和四面体单元(CPS6M Z CPE6M,C3D10M),二阶梁单元
(B22,B32),全积分膜单元(M3D4),全积分壳单元(S3,S4T),全积分一阶六面体单元(C3D&
C3D8I,C3D8T)
3. abaqus/explicit分析结果中经常会出现沙漏现象,沙漏现象是网格出现扭曲但不产生应变能导致结果不准确。
在某些荷载条件下,线性减缩积分单元可能会出现所谓沙漏(hourglassing)的非物
理变形。
考虑一个一阶减缩积分单元,模拟一小片承受纯弯荷載的材料的情况,
当单元变形时,单元中虚线的长度没有改变,并11他们之间的夹角也未改变。
因此,在单元单个积分点上的所有应力和应变分量都为零。
因为这种方式的单元扭曲并不产生应变能, 所以这种变形模式是一个零能量模式。
由于单元在此模式下没有刚度,因而不能阻止这种的变形。
在粗网格中,这种零能量模式可以通过网格扩展,从而产生不正确的结果。
在ABAQS/Explicit中,仅一阶减缩积分四边形和六面体单元具有沙漏模式。
沙漏有时可能扩展到整个网格。
-ABAQS/Explicit包含比较严密的控制以防止由沙漏现象导致大多数实际分析出现问题。
然而,控制沙漏的工作是通过校正作用力的大小实现的,有可能要通过儿个増量步才能控制沙漏。
在一些比较严重的情况下,在沙漏控制能够纠正问题以前,沙漏现象可能
己经扩展到了整个网格。
沙漏现象可以通过观察网格的变形情况來进行判断。
一个更为量化的途径是研究伪应变能,它是控制沙漏变形所耗散的主要能量。
如果伪应变能过高,说明过多的应变能可能被用來控制沙漏变形。
在ABAQS/Explicit中,变量ALLAE是伪应变能的能量耗散总和。
ALLAE包含粘性和弹性两项;然而,由于粘性项通常占主要地位,因此大部分转化伪伪应变能的能量是不可能恢复的。
当ALLAE的变化曲线斜率非常陡时,耗散的伪应变能也相应很多:反之,若ALLAE 的变化曲线趋于平缓,伪应变能的耗散则很少。
评判标准:ALLAE/ALLIE<0. 5%
4. Abaqus/explicit中的准静态分析,如金属成型过程;
准静态分析的目的是:在惯性力的影响较小的前提下,尽量缩短计算的时间周期,可考虑的因素有:
载荷速率:在模拟过程中,人为的增加准静态成型过程的速度是必要的,它可以让求解过程更经济。
一般推荐的载荷速率为材料中波速的1%,材料波速c=VE7po推荐方法:以不同的速率多次模拟,检查结果(变形形状,应力,应变,能量),分析不同载荷速率对结果的影响。
在分析步内,以光滑缓坡地方式,把把碰撞速度从零增加到所施加的碰撞速度可以得到更加精确的解。
(施加力以同样的方式)需定义SMOOTH STEP幅值曲线。
能量平衡:能量检査未Abaqus/explicit金属成型过程的结果是否反应了准静态提供了评估方法:在主要的成型过程中,变形材料的动能不可以超过内能的一小部分,ALLKE/ALLIE —般为1%~5%,使用光滑幅值曲线将改进早期响应。
质量缩放:人为的增加成型速度可以提高解的经济性。
同时,材料应变率以同样的速度增加。
(如果材料对于应变率是不敏感的,这是不相关的:但如果模型中考虑应变率敏感的, 将导致错误的结果。
)
人为的增加工具速度之后,质量缩放对惯性效应具有同样的影响,过多的质量缩放将导致非真实的解。
如果质量缩放用于完全的动态条件下,总质量的变化应该尽量小(小于1%)。
显示动力学过程稳定极限的估计公式为:At = Al/c,其中山为单元最小尺寸,c为波速=7E7P,E 为弹性模量,P为密度。
如果以f?的方式人为的增加材料密度:膨胀波速以f的方式减小,稳定时间增量以f的方式增加。
利用*FIXED MASS SCALING选项,可以在分析步开始时施加质量缩放。
句法:*FIXED MASS SCALING, ELSE=?, FACTOR=f2
网格自适应:在许多非线性模拟中,结构或加工过程中的材料承受非常大的变形,这些变形将扭曲有限元网格,过度扭曲的网格将不能得到准确地解,或由于数值原因导致分析提前终止。
在这些模拟中,必须使用网格自适应工具周期性的减小网格的扭曲。
在Abaqus/explicit中,口适应网格功能通过任意的拉格朗日■欧拉(ALE)方法实现。