【有限元分析】有限元仿真分析与解析解的结果对比——以阶梯轴的静力分析为例
(1)对一个阶梯轴零件进行基于材料力学的理论计算,求解最大应力值;
(2)在WORKBENCH中对该阶梯轴零件进行有限元仿真,实行两种仿真方案,分别是1.梁模型建模+梁单元网格划分;2.实体模型建模+六面体单元网格划分,观察两种仿真结果并与理论计算结果的对比,对比结果发现解析解与仿真解相差很小。
(3)可以借此算例学习WB中的梁单元静力分析、三维实体静力分析、理解并施加若干种边界条件,举一反三即可了解此类轴系中轴零件的强度分析。
在进行阶梯轴零件设计的时候一般会对其进行强度校核,校核方式主要有理论计算和仿真分析两种。轴零件的强度校核计算方式已经标准化,查阅手册即可,仿真分析可使用有限元仿真软件,本文算例将在ANSYS WORKBENCH 进行。
本文的算例来自于《ANSYS Workbench 工程实例详解》,以校核阶梯轴强度问题为例,探讨使用解析解解法和有限元分析解的差异。
一、算例描述及其解析解
图1为阶梯轴的简图,现校核其受载后的静强度,已知直径d1=180mm{d_1=180mm}d1=180mm,d2=150mm{d_2=150mm}d2=150mm,a=300mm{a=300mm}a=300mm,b=200mm{b=200mm}b=200mm,L=1000mm{L=1000mm}L=1000mm,F=300kN{F=300kN}F=300kN,材料为45,弹性模量E=2.1e11Pa{E=2.1e11Pa}E=2.1e11Pa,泊松比v=0.28{v=0.28}v=0.28,屈服应力δs=355MPa{δ_s=355MPa}δs=355MPa。在AB段,轴只受弯矩MAB{M_{AB}}MAB,而外伸到加载处的这一段,既受弯矩又有剪力,属于横力弯曲。根据材料力学分析,最大正应力应该产生在C截面的圆边缘处,强度为:
δmax=δc=MC/WC=32Fb/πd23=181.083MPa{δ_{max}=δ_c={M_C}/{W_C}={32Fb}/{π{d_2}^3}=181.083MPa}δmax=δc=MC/WC=32Fb/πd23=181.083MPa
同理AB段的最大应力大小为:δAB=MAB/WAB=32Fa/πd13=157.19MPa{δ_{AB}={M_{AB}}/{W_{AB}}={32Fa}/{π{d_1}^3}=157.19MPa}δAB=MAB/WAB=32Fa/πd13=157.19MPa
图1 算例的理论解法
二、有限元仿真分析结果
为了简化仿真分析难度,考虑到目前ANSYS Workbench已经普及,且其流程化的操作方式也被越来越多的机械工程师所接受,故本文使用该仿真平台。
在有限元分析的操作过程中,流程可简化为**建模→网格划分→设置边界条件→求解→结果后处理。**就重要性来说,前处理过程包括建模,网格划分和设置边界条件都是非常关键的步骤。网格划分需要考虑网格的类型、形状和尺寸等因素,而在设置边界条件时需确保对模型施加的边界条件与实际加载工况一致,三者均需保证准确无误,否则会导致计算结果与实际情况大相径庭,误导未来的进一步设计。
梁单元静力学分析:
当结构长度对横截面的比率超过10:1,沿长度方向的应力为主要分析对象,且横截面始终保持不变时,在WB中默认为铁摩辛柯梁单元,即beam188和beam189,可计算弯曲、轴向、扭转、和横向剪切变形。beam188和beam189两者的区别是形函数种类不同。beam189的精度更高,计算消耗内存也多,所以在仿真时需要权衡计算精度和时间。本文均使用beam188。
建模:
很多人在最初学习WB的时候,已经掌握了一种或多种三维绘图软件(例如SW、PRO/E等),并且认为只需使用这些三维软件绘制三维图再转化相应格式并将其导入WB中即可。但是这种导入的模型,在进行前处理划分网格的时候只能使用三维实体单元划分网格,而不能使用梁单元、二维平面单元和三维壳单元。例如本文中的算例,若要使用梁单元划分网格,必须在WB中DM模块中建立line body模型,具体参见《ANSYS Workbench 工程实例详解》。
划分网格:
使用line body生成的模型,可以选择自动划分网格的,WB会将以梁单元将其划分完成,例如图2所示。
图2梁模型建模后划分的网格
使用DM建立的solid模型或者通过其它三维软件绘制并导入的模型,可以划分为六面体网格。例如图3所示。
图3 三维实体建模后划分的网格
设置边界条件:
梁模型设置边界条件:需考虑零件的受载和约束,受载即为该轴在两端外侧分别受到大小为300kN的力,方向沿-Y方向。约束需考虑图1中所示的A和B两点:A点是简支约束,在受载后只会产生绕Z轴的旋转,其它自由度均被锁死,故需要在A点施加约束Simply Support和Fixed Rotation(X和Y设置成Fixed,Z设置成Free);B点有X方向上有平移自由度,和绕Z轴的转动,故需在B点施加约束Fixed rotation(X和Y设置成Fixed,Z设置成Free)和Displacement(X方向设置为free,Y和Z设置为0)。如图4所示。注意下图中的字母标识为软件自动生成,与上述算例描述无关。
图4 梁模型施加边界条件
三维实体模型设置边界条件:可在两端面加载大小为300kN的力,设置约束时可使用远程边界条件-远程位移,来设置A和B两点的自由度,在添加远程边界条件时可基于remote point。如图5-图7所示,添加基于A和B所在面与轴相交的外圆线生成的remote point。(关于远程边界条件和远程点,感兴趣的同学可以参看ANSYS Workbench Help 文件)
图5 三维实体模型施加的边界条件
图6 A点remote point的设置
图7 B点remote point的设置
结果和后处理:
梁模型的后处理需导入 Beam Tool→Maximum Combined Stress,求解零件的在各个位置的最大应力,结果如下所示,如预期所料最大应力值出现在C截面,最大值如左侧的颜色条所示:
图8 梁模型应力分布
结论:最大应力值为181.46MPa,理论值为181.08MPa,误差为0.1%。
三维实体模型的von-mises应力云图如图9所示,但是求解的最大的应力值为258.44MPa,与理论计算值不符,这是因为由于有限元计算的特点,在该处会出现应力集中的现象。所以再分析C截面的应力值。
图9 实体模型应力云图
图10 C截面上的路径
图11 应力值随路径位置的变化
由于应力集中导致峰值应力相差很大,所以反映的总应力相差也很大。
结论:C截面的弯曲应力为173.15MPa,理论值为181MPa,误差小于5%。
【有限元分析】有限元仿真分析与解析解的结果对比——以阶梯轴的静力分析为例相关推荐
- ANSYS Workbench18.0有限元仿真分析视频教程
ANSYS Workbench18.0有限元仿真分析视频教程 链接:https://pan.baidu.com/s/1oTZNcf7_ZWMVgKiWM1l0Wg?pwd=it2v 提取码:it2v ...
- 有限元仿真分析技术中网格划分的类型与步骤
对于有限元分析来说,网格划分是其中最关键的一个步骤,网格划分的好坏直接影响到解算的精度和速度.网格划分有三个步骤:定义单元属性(包括实常数).在几何模型上定义网格属性.划分网格.今天给大家分享:有限元 ...
- Creo5.0结构有限元仿真分析视频教程
Creo5.0结构有限元仿真分析视频教程 链接:https://pan.baidu.com/s/1dgM3gnz4SlA5PRot94h-LA 提取码:pjw4
- 贝叶斯分析好坏_贝叶斯估计对比于经典估计的优势分析与其局限性
贝叶斯估计对比于经典估计的优势分析与其局限性 经典估计和贝叶斯估计 经典估计理论是通过一个随机抽样过程, 从总体中随机抽取一定数量的样本, 再结合总体分布或总体分布族提供的的信息, 推断出总体分布或总 ...
- ANSYS有限元仿真分析:大变形(Large Dflection)开关的意义
为什么我们想做ANSYS有限元仿真的时候,老是听别人说打开大变形开关,是否转念一想,这会啥玩意?那为什么要打开?打开关闭的好处是什么呢?背后的机理又是什么呢? 如上图所示,如果未打开大变形开关输出力矩 ...
- ANSYS有限元仿真分析:边界非线性 (接触Contact)
一.非线性的概念 二.边界非线性的接触 在接触行为中主要分成:黏死(Bonded).不分离(No Separation).无摩擦(Frictionless).无滑动(Rough).有摩擦(Fricti ...
- 【设计模式】装饰者模式 ( 概念 | 适用场景 | 优缺点 | 与继承对比 | 定义流程 | 运行机制 | 案例分析 )
文章目录 I . 装饰者模式概念 II . 装饰者模式适用场景 III . 装饰者模式优缺点 IV . 装饰者模式与继承对比 V . 装饰者模式相关设计模式 VI . 装饰者模式四个相关类 VII . ...
- 一个C#和C++执行效率对比的简单实例
这里用一个算法题进行比较. 原题是见http://acm.hdu.edu.cn/showproblem.php?pid=4090,登载在http://blog.csdn.net/woshi250hua ...
- 云南高性能云桌面搭建解决方案、云桌面与传统PC优势对比,云桌面适用场景分析
一.传统PC面临的问题 1.难管理:个体复杂度带来的管理困难 企业内部PC分布广.数量多.型号杂,难以标准化及统一管理 软件应用.业务系统安装部署繁琐.业务上线效率低 无法集中管控,维护跟着物理机跑等 ...
- 【论文解读】NAACL 2021 对比自监督学习的上下文和一般句子表征:以篇章关系分析为例
前言 本文是NAACL 2021的论文<Contextualized and Generalized Sentence Representations by Contrastive Self-S ...
最新文章
- 众多新点子推动比特币现金BCH普及
- JS字符串的下划线命名和驼峰命名转换
- CentOS7安装和配置samba
- Python多线程学习教程
- python 中numpy dot函数的使用方法
- 项目管理一般知识:单个项目的管理过程
- python入门之迭代器
- python的格式_python 代码格式
- jQuery在同级节点中查找
- 《Artifact》的得与失:成功的游戏工业品,却与主流背道而驰
- rpmbuild FILE
- 快速傅里叶变换应用之二 hdu 4609 3-idiots
- 积微论坛--用微生物组时序数据重现生物膜装配动态过程PPT对应的讲解和提问
- 苹果手机软件升级密码_微软:将向安卓和苹果iOS平台推出杀毒软件Defender
- 计算机科学导论教学大纲,lbrack;机器人学导论rsqb;课程教学大纲
- 数据库学习与应用之什么是数据库
- matlab群速度色散系数,第三章-群速度色散.ppt
- python处理word文档格式_python处理word文档(docxtpl 中文说明文档)
- Linux系统关闭virbr0
- ActiveMQ 停止 关闭 报错 ERROR: No or outdated process id in
热门文章
- 全站仪坐标计算机公式,全站仪使用方法及坐标计算,详细讲解!建议收藏!
- IPv4与IPv6的区别
- ArcEngine代码 地图制图标注之 TextElement标注,Annotation注记,MapTips显示
- 戴尔linux恢复镜像,如何从官网下载并使用系统恢复映像
- 颜色的原理,三基色原理以及HLS(色相、亮度、饱和度)原理
- 计算机页面格式和编排,论文开题格式
- linux应用项目(一)数码相框(2)数码相框之字符编码与字符的点阵显示
- 五阶段--使用 Kibana 操作 ES/ 搜索
- CGAL学习之路(三):CGAL读写点云
- html转化pdf文字无法显示,HTML 转 PDF 排雷 之 wkhtmltopdf 转换得到的PDF没有内容