ANSYS学习记录
基于教程
https://www.bilibili.com/video/BV1u4411B7Fo?p=3
核心:学软件是找到自己要的功能 不是每个功能都关注
学习过程问题记录:
- 打开下部的message查看警告。关于报错:pivot转轴;matrix矩阵,iterative积分,unconverged不收敛。
- An internal solution magnitude limit was exceeded. Please check your Environment for inappropriate load values or insufficient supports. Please see the Troubleshooting section of the Help System for more information.某个方向上的载荷超过强度。
- ansys界面混乱:view-reset。
- ansys不支持撤销命令,所以随时保存。
- At least one body has been found to have only 1 element in at least 2 directions along with reduced integration. This situation can lead to invalid results. consider changing to full integration element control or meshing with more elements.网格密度不够。可以使用goto-Bodies With One Element Through the Thickness定位。
- 结果中的等效应变和等效压力和等效变形是不同的,可以通过左上角unit量纲来判断。
- 参考资料[:](https://wenku.baidu.com/view/34709ff7182e453610661ed9ad51f01dc381577d.html?rec_flag=default&word=ansys troubleshooting section&fr=pc_oldview_relate-1001_1-4-wk_rec_doc2-1001_1-57270722192e45361066f505-34709ff7182e453610661ed9ad51f01dc381577d)这个整理的很好,但是38块钱,我没买哈哈哈。[https://wenku.baidu.com/view/34709ff7182e453610661ed9ad51f01dc381577d.html?rec_flag=default&word=ansys%20troubleshooting%20section&fr=pc_oldview_relate-1001_1-4-wk_rec_doc2-1001_1-57270722192e45361066f505-34709ff7182e453610661ed9ad51f01dc381577d](https://wenku.baidu.com/view/34709ff7182e453610661ed9ad51f01dc381577d.html?rec_flag=default&word=ansys troubleshooting section&fr=pc_oldview_relate-1001_1-4-wk_rec_doc2-1001_1-57270722192e45361066f505-34709ff7182e453610661ed9ad51f01dc381577d)
Day 1:
关于软件选择:我们主要使用Workbench部分 在sw等软件中完成建模而后直接分析 ANSYS自带的建模使用感受一般
关于版本选择:从12.0后界面基本不变了,所以可以使用临近的3-5个版本均可 例如:我现在最新的版本是19.2 那么推荐使用的是19/18/17三个系列 16/15两个系列也可以 再低版本不推荐
关于知识储备:材料力学基础(应力剪力等)+较好的英文基础(主要是掌握机械设计的专业术语 要不看不懂操作选项QAQ)
大体制作流程:三维建模→模型转化→条件附加→网格划分→求解及后处理
界面分配:主界面+DM设计界面+材料附加界面+分析界面(主要!!)
文件类型:sw导出的step/x_t/igs均可
中英操作对照:
https://wenku.baidu.com/view/3f2e4e4ee518964bcf847cca.html
我们能用到的基本都在这了
关于WORKBENCH入门基础:
- 完整流程:
- 界面说明:
- 分析流程:
大体来说即选择分析类型->工程数据(材料等,其中材料默认为结构钢)->几何模型(导入时类型选择所有,sw可以使用stl和-xt)->分析模型(网格,约束,载荷(大小和方向))->解算和结果
实例1:
分析实例:带孔矩形板
问题描述:带孔矩形板一侧面均布110000N载荷,另一侧固定,利用WORKBENCH求解矩形板的应力、应变和位移,并对分析结果进行解释。
材料:Structure Steel
屈服强度:250MPa
注意:generate mesh→force(方向默认法线)→fixed support(apply)→1:1显示
通过改变比例尺上的屈服强度来达到不同的颜色显示效果(后处理)
扩展知识:失稳起皱-在弯曲或挤压过程中部分区域过厚堆积引起皱褶
实例2:最大应力值位置及数值的统计 应力产生规律
分析实例:L型支架
问题描述:L形支架上端面固定,同时在下端面施加900N弯曲载荷,分别使用10mm、3mm和1mm网格求解模型,对比应力结果。
材料:Structure Steel
注意:
mesh→sizing→element size(工具栏units切换单位)
force选线是平行于线, 选面为面的法线。工具栏选maximum查看最大值点
改变网格大小时要新建工程算例不要在原分析上直接改,复制后修改网格大小即可
工程数量较多时要进行Archive打包操作 因为ANSYS默认只保存例程不保存数据
扩展知识:
网格无关性检查-一直减小网格直到计算结果收敛在允许的范围内即认为结果此时与网格大小无关
应力奇异-按照网格无关性检查不一定能找到结果进入收敛半径的网格大小,即应力趋于无穷,此时为应力奇异,常出现在不连续处如不加圆角的直角边,详情见
http://blog.sina.com.cn/s/blog_13b525e650102yejg.html
Day 2
实例2:圆角L型支架的应力结果
分析实例:L型支架
材 料:Structure Steel
问题描述:修改第二课的L形支架模型,如图位置增加半径2.54mm的圆角,上端面固定,同时在下端面施加900N弯曲载荷,分别使用10mm、3mm、1mm和0.5mm网格求解模型,对比应力结果。
注意:geometry→dm(design model)添加圆角→每一步都generate
数据连接而非复制的好处是能够只修改源例程而实现全部数据替换
分析结果中的geometry从all body改成其他点线面以实现部分结果预览 应用于已经对仿真结果有初步预估 (如对应力最大值点初步判断)可减少处理时间
菜单最上面有一栏元素过滤器 可选择线/面/体
在达到幅值的5%左右波动就可以认为满足了网格无关性判定 此时红色危险区域应覆盖两层基本单元
1mm结果 不可
0.5mm结果 可
实例4 应力奇异及位移结果对比
应力集中:
https://baike.baidu.com/item/%E5%BA%94%E5%8A%9B%E9%9B%86%E4%B8%AD/1953796?fr=aladdin
是一种很正常的现象
应力奇异是应力集中的极端现象 是一种数学算法问题 反推可得有限元分析是不适用于尖角/折弯处的红色≠应力奇异/应力集中 只要产生应力梯度就是应力集中
99%的物体都有应力奇异点 但是99%都在无危险处
当出现应力奇异点 我们不是考虑更改网格把他算准 而是考虑更改结构来避免 如给L型角件加圆角
更改网格大小不会影响位移的分析(误差在1%左右)因此优先对比位移结果
Day 3
课程5:初步了解单元基本知识
改变局部网格密度:mesh→sizing→geometry→选择面→generate mesh
此时大概率出现网格变为长方体 解决方法:mesh→sizing→选择包围该面的所有边界线→generate mesh
讲网格从六面体改为四面体:mesh→method→automatic改为tetrahedrons
使网格更加平滑:sizing→transition→slow
在同一网格精度下, 四面体堆积比六面体堆积更密集, 因而计算结果更接近最大值
因此目前的普遍观点是(仍在争论)四面体计算精度更高, 可以保留更多的物理细节,六面体划分更简单 其原理就是对物理特征进行了自动简化——这种简化是有益的,使其计算速度更高(囿于目前的电脑性能),所以这二者互补,这里列举一些老师的话帮助大家理解:
1.六面体网格划分有助于学习者提升对网格划分的理解;
2.随着模型的复杂程度越来越高,四面体网格在划分效率上优势会非常明显
课程6:装配体分析预备(重点!!)
分析实例:支架
问题描述:支架两端通过地脚螺栓固定在地面上,零件凹孔面承受5000N的法向推力,求解支架的应力和位移。
材料:Structure Steel
当出现应力奇异等其他难以区分的分析结果时,可以考虑分析位移结果,选择directional可选择x/y/z其中一个方向,添加位移结果后,可选择加入0位移线ISOlines,来验证是否有不合理位移,此时要将坐标轴随便选一点归零
Day 4
课程7:装配体分析基本思路
导入装配体后在connection中自动显示了接触部分 选择body view来预览接触物体 explode爆炸视图后可以自行添加connection
接触类型选择:
Bonded(绑定):默认接触形式,不允许面或线间有相对滑动或分离,可以将此区域看做被连接在一起。
Frictionless(无摩擦):这种接触形式代表单边接触,即如果出现分离则法向压力为零,同时假设摩擦系数为0。
Rough(粗糙):这种接触方式和无摩擦类似。但表现为完全的摩擦接触,即没有相对滑动,法向可分离,切向不滑动。
Frictional(有摩擦):这种情况下,在发生相对滑动前,两接触面可以通过接触区域传递一定数量的剪应力。法向可分离,切向摩擦滑动。
No Separation(不分离):这种接触形式和绑定类似,只适用于面。不允许接触区域的面分离,但是沿着接触面可以有小范围滑动,即法向不分离,切向可以有小位移。
做出结果后 要多次改变网格大小 没有结果崩溃即证明所做分析可靠。
课程8:刚体平移和弱弹簧
fixed support 三维全约束 displacement可按x/y/z分类约束 默认为free输入0约束
show undeformed model 可以显示平移前的位置
电脑计算存在这样的问题:给一块板左右分别加100N的力,我们知道这是平衡的,但是当两边网格不一致时,100N会被读取为两个力如100.2和99.8这其中就产生了微小力,这是一种误差,不可避免。
利用弱弹簧功能可以抵消微小力 analysis settings→Weak Springs→on
检测微小力的大小:solution→probe→force reaction→boundary condition→weak springs
Day 5:
课程 9:接触基础和接触设置
计算错误最常见的情况就是接触错误和整体模型不平均
在分析结果之前一定要看message里面的警告和错误
设置Frictional接触类型(摩擦系数0.2)可以解决非常规平移的问题
Lesson 10:分析结果解读
分析实例:L型支架
材料:Structure Steel
问题描述:如下图所示,支架红色端面施加箭头方向900N载荷,横杆两端面固定,假设横杆和支架为零配合状态,考察结构的应力和位移
分析结果与实验结果的比较:
1.最大变形:实验结果有指定测量位置和方向而分析结果只有全变形
2.测量结果:应力值不是直接获取,测试时必须给定弹性模量,才有对比性
Day 6
Lesson 11
圣维南原理和模型简化
圣维南原理:
https://baike.baidu.com/item/%E5%9C%A3%E7%BB%B4%E5%8D%97%E5%8E%9F%E7%90%86/2590043?fr=aladdin
我个人的理解是:有两点A和B, 当A和B距离足够远,B处只要合力和合力矩不变,那么B处有多少个力/方向如何(即力的复杂程度)都不会对A处的应力情况产生影响。
1.关键是找到这个“足够远”的距离是多少
2.在工程实例中广泛满足,但是没有理论验证
3.应力分析符合,位移分析不符合
将方杆
对L型支架的力通过圣维南转换成固定约束发现位移量有很大变化,这是因为图一显示的位移是方杆和L型支架共同的位移量
四大强度理论:
1、最大拉应力
理论
2、最大伸长线应变理论 →→脆性断裂
3、最大切应力理论
4、形状改变比能理论 →→塑性屈服(我们要用的)
强度条件:
$$
Von Mises应力=[(σ1-σ2)^2+(σ2-σ3)^2+(σ3-σ1)^2]≤[σ]
$$
式中σ1,σ2,σ3,为三大主应力,σ为屈服应力
安全系数:Tools→Stress Tool→Safety Factor 计算:屈服值/应力值
安全系数应大于1 我们要求应该在1-2左右
Lesson 12 有限元分析学习方法介绍
Day 7
Lesson 13 对称问题预备:水杯
分析实例:水杯
材料:ABS
问题描述:杯子放置在桌子上,内部装满水,求杯子底部径向扩张量。
材料修改:Material→Assighment→New Material或者返回Project界面进入Engineering Data
新建材料:命名→Density(密度)→Isotropic Elasticity(各向同性下的弹力属性)→Young’s Modulus(杨氏模量)→Poisson’s Ratio(泊松比)→返回Model界面File里Refresh
水压设置:Hydrostatic Pressure→Fluid Density(液体密度)→Magnitude(重力加速度)→Direction选择一个面,他的法线为所选方向→Location定义一个零压强面
需要杯子和桌子两个物体 ,限定为有摩擦接触,千万不可用全约束把杯底固定死
这里就要用到之前讲的局部网格细分,四面体网格划分,缓慢过渡,软弹簧等操作
上图分别为有摩擦接触/弱弹簧/全约束
Lesson 14 对称
对于水杯这种中心对称/轴对称物体我们可以只分析他的一部分(以1/4为例)
切除:在DM界面切掉3/4之后Tools→Symmetry(对称,即告诉系统你这里还有镜像的另一半只是不显示不计算)→Plane→Model Type→Partial Model
建立新plane:Create→New Plane→Form Face→Base Face选择切割面(plane必须与xy平面重合)
plane与face的区别:plane值得广义无穷大面或坐标轴平面 face指物体上的一个小端面
mesh分割:New Section Plane新建一个任意切割面
Day 8
Lesson 15 轴对称问题应用
对称问题:如果模型的几何形状及边界条件都对称于同一平面,则结构内各点的位移、应力及应变都对称于此平面,这类问题称之为对称问题。如之前所做的矩形板/L型支架/水杯
对于水杯这种绕轴旋转而成的物体,他的每一个切面都相同,所以可以简化到一个面,称为轴对称模型
因为从三维转到了二维所以要新建工程分析,不要复制之前的算例
新建后点geometry右侧属性栏中把analysis从3D改成2D(view→properties勾选上即可显示属性栏)
DM先把水杯切成1/4 然后thin/surface选中水杯截面和桌子截面两个面 FD1,thickness改为0
所选面必须在xy平面上且y为旋转轴
转换到xy平面:create→body transformation
在mesh前选择geometry→2D behavior→ axisymmetric切换到镜像问题模式
分析出2D结果,转换为3D显示:
①主界面→tools→options→appearance→beta options勾选上
②model→symmetry→type改为2D axisymmetric
③num repeat 平面复制的份数-1 Δθ每一份的角度 初始为37份10度
Lesson 16对称问题 平面问题 自由度
静力学要求每个物体的六个空间自由度全部被约束(如接触/摩擦/对称)
带孔矩形板就可变为 对称→平面→边界条件→求解
Day 9
Lesson18 简化模型正确性
分析实例:桌面放置物体
问题描述:考察桌子正中间放置一个质量为45.4kg的长方体,考察桌子的变形情况。
材料:结构钢
分别用两种方法判断:
①想象成均布载荷用force代替
②物体自重:inertial→standard earth gravity→+z direction
通过probe→force reaction 对四个桌脚的支反力检测可以确保两个方法其他条件均相同
求解时提示数据错误可以solution中清空数据重新求解
应变结果相差超过25% 对比应力结果:
对比发现桌子受压变形 ,物体给的并不是均布载荷 ,实际上均布载荷是一种理想化受力,只有很少情况出现
tools→contact tool→插入pressure→求解出实际应力分布 这是对简化的重要验证
结果:物体压力应等效为四个点作用而非均布载荷
Day 10
Lesson 19 力学计算对比
在DM界面把两个solid放在一起 form new part里可以把两个物体视作一个来分析
这里复习材料力学中梁的有关词汇:超静定问题,转角和挠度,弯矩等
Lesson 20
线性静态分析所涉及的材料参数:
必要参数:弹性模量、泊松比(金属在0.23-0.33之间)
一个特殊的材料参数:屈服强度
可能需要参数:密度(惯性力)、线膨胀系数(热边界条件)
添加了重力就必须添加密度/ 添加重力与不添加重力结果差10%就一定要添加重力
tools→safety factor安全系数才需要填写屈服强度
添加屈服强度:strength→tensile yield strength
三个材料参数对线性材料小变形计算结果的影响
弹性模量:应力结果不影响,变形结果影响;
屈服强度:应力结果变形结果均不影响;
泊松比:应力结果影响,变形结果影响,但是变化范围较小。
对比实验结果和仿真结果时不要添加重力 因为在实验进行调零时已经完成重力形变 所以重力已经被忽略