ABAQUS简支梁分析梁单元和实体单元.docx
- 文档编号:12244975
- 上传时间:2023-04-17
- 格式:DOCX
- 页数:17
- 大小:1,011.52KB
ABAQUS简支梁分析梁单元和实体单元.docx
《ABAQUS简支梁分析梁单元和实体单元.docx》由会员分享,可在线阅读,更多相关《ABAQUS简支梁分析梁单元和实体单元.docx(17页珍藏版)》请在冰豆网上搜索。
ABAQUS简支梁分析梁单元和实体单元
基于ABAQUS简支梁受力和弯矩的相关分析
(梁单元和实体单元)
对于简支梁,基于ABAQUS2016,首先用梁单元分析了梁受力作用下的应力,变形,剪力和力矩;对同一模型,并用实体单元进行了相应的分析。
另外,还分析了梁结构受力和弯矩作用下的剪力及力矩分析。
对于CAE仿真分析具体细节操作并没有给出详细的操作,不过在后面上传了对应的cae,odb,inp文件。
不过要注意的是本文采用的是ABAQUS2016进行计算,低版本可能打不开,可以自己提交inp文件自己计算即可。
可以到小木虫搜索:
“基于ABAQUS简支梁受力和弯矩的相关分析”进行相应文件下载。
对于一简支梁,其结构简图如下所示,梁的一段受固支,一段受简支,在梁的两端受集中载荷,梁的大直径D=180mm,小直径d=150mm,a=200mm,b=300mm,l=1600mm,F=300000N。
现通过梁单元和实体单元分析简支梁的受力情况,变形情况,以及分析其剪力和弯矩等。
材料采用45#钢,弹性模量E=2.1e6MPa,泊松比v=0.28。
图1简支梁结构简图
1.梁单元分析
ABAQUS2016中对应的文件为beam-shaft.cae,beam-shaft.odb,beam-shaft.inp。
在建立梁part的时候,采用三维线性实体,按照图1所示尺寸建立,然后在台阶及支撑梁处进行分割,结果如图2所示。
图2建立part并分割
接下来为梁结构分配材料,创建材料,定义弹性模量和泊松比,创建梁截面形状,如图3,非别定义两个圆,圆的直接分别为180和150mm。
然后创建两个截面,截面选择梁截面,再选择图2中的所有梁,定义梁的方向矢量为(0,0,-1)(点击图3中的n2,n1,t那个图标即可创建梁的方向矢量),最后把创建好的梁赋给梁结构。
图3创建梁截面形状
接下来装配实体,再创建分析步,在创建分析步的时候,点击主菜单栏的Output,编辑EditFieldOutputRequest,在SF前面打钩,这样就可以在结果后处理中输出截面剪力和力矩,如图4所示。
在Load加载中,在固支处剪力边界条件,约束x,y,z,及绕x和y轴的转动,如图5所示,同理,在固支另一处约束y,z,及绕x和y轴的转动。
在梁的两端添加集中力,集中力的大小为300000N。
最后对实体部件进行分网,采用B32梁单元,网格尺寸为10。
完成以上工作后,创建作业并提交分析。
(由于操作比较简单,故没有详细列下所有操作步骤。
)
图4Step中SF输出编辑
图5边界条件约束
图6为等效应力云图,可知最大应力为181.1MPa,最大位置出现在梁台阶处(梁直径变化处)。
根据材料力学,最大弯矩应力产生在C截面,同时根据材料力学知道AB段处的最大应力,其应力为
(1)
(2)
从图6和图7可以知道,梁的最大应力以及AB段的应力都与理论解一致。
图8为梁的等效应力图,可见最大位移出现在梁的两端,最大有1.639mm。
沿着梁的轴线建立路径,然后绘制出梁的变形,图9和图10分别给出了截面剪力和力矩沿路径的变化情况。
值得注意的是,图9中剪力图与材料力学的剪力图有区别,其并不是按照设正法画的剪力图,不过其数值的绝对值与材料力学上的一致。
图10的弯矩图也材料力分析一致,图11为等效位移沿路径的变化情况。
图6等效应力
图7中间段等效应力图
图8等效位移图
图9剪力SF2沿路径情况
图10弯矩SM1沿路径情况
图11等效位移U沿路径情况
2.实体单元分析
ABAQUS2016中对应的文件为beam-solid.cae,beam-solid.odb,beam-solid.inp。
按照图1建立相应的实体单元,然后在支撑处切割实体,再建立材料属性,装配,设置步,在创建分析步的时候,点击主菜单栏的Output,编辑EditFieldOutputRequest,在SF前面打钩,同梁设置操作一样。
图12实体单元结构
图13给出了力载荷及边界条件加载情况,在梁左右两端加载力载荷,可以选择对应的面,然后加载界面切应力(F/S=300000/(3.1414*75^2=16.9765MPa,S为圆的面积),也可以在加载面的中心建立参考点,然后把加载面与参考点相耦合,然后直接在参考点加载300000N的力即可。
在支撑处选择对应的面,然后加载对应的边界条件约束即可。
最后完成相应的分网并提交分析。
图13载荷及边界条件
图14等效应力图
图15等效位移图
图16输出剪力弯矩操作
图17输出剪力弯矩操作
图18剪力数据数据输出
图14给出了等效应力图,其数值结果与梁给出的结果一致,与理论计算也一致,但是图15给出的等效位移图与梁单元计算的等效位移图不同,因此在具体问题分析的时候,我们应该判断具体使用什么单元进行分析。
图16到图18给出了如何把沿轴向方向各个截面的剪力和弯矩输出到一个txt文件的操作。
在图16中打开主菜单Tools中的ViewCutManager,然后在轴线平面前面打钩,结果如16所示,然后单击图16中的Options,再单击Slicing进入到图17,设置实体切割的数目,然后单击OK即可。
再在主菜单中单击Report的FreeBodyCut,然后按照图18中设置,然后单击OK,这样就可以在ABAQUS的当前文件夹找到moment.txt文件,里面记录了各个截面的力矩和弯矩,把里面的数据复制到excel中就可以绘制出弯矩和剪力图。
3.梁同时受集中力和弯矩分析
ABAQUS2016中对应的文件为beam-force-moment.cae,beam-force-moment.odb,beam-force-moment.inp。
对于该分析,还给出用Workbench给出的对应分析,其文件为beam-force-moment.wbpj,Workbench版本为15.0。
图19为梁在中间受集中载荷和弯矩的结构简图,梁的长度l=1000mm,梁的直径d=100mm,Me=9e7N·mm,F=300000N,E=210000MPa,v=0.28。
现采用梁单元进行剪力和弯矩等相关分析。
图19梁受力和弯矩结构简图
其建模等分析过程与前面第一个实例的梁单元分析一样,在梁的中点部分进行切割。
为了对比分析只有集中力,弯矩以及集中力和弯矩同时作用在梁上的几种情形,如图20,对应Load处右键单击,然后选择Suppress,可以抑制载荷作用,选择Resume,可以激活抑制的载荷。
通过载荷抑制或者激活的设置,然后提交分析,就可以得到对应载荷情况下的结果。
图20载荷加载情况
(1)梁中部只受集中力情形
下面四个图给出了只受集中力载荷作用下的等效应力,等效位移,剪力和弯矩图,剪力和弯矩图与材料力学分析一致。
图21等效应力图
图22等效位移图
图23轴线方向剪力图
图24轴线方向弯矩图
(2)梁中部只受弯矩作用
下面四个图给出了在梁的中部只加载弯矩作用下的等效应力,等效位移,力矩和弯矩图。
从剪力和弯矩图可以知道,整个梁受到90000N(Me*l=9e7N·mm*1000mm=90000N)的剪力作用,这与材料力学分析不一致,这应该引起关注。
在材料力学分析中,只受到弯矩作用时,应该没有剪力作用,不过按照这思路也可以做出有弯矩作用下的剪力和弯矩图。
在弯矩图中可以知道,弯矩中部处出现了突变,有均匀剪力作用处,其弯矩呈现线性变化。
图25等效应力图
图26等效位移图
图27轴线方向剪力图
图29轴线方向弯矩图
(3)梁中部同时受集中力和弯矩作用
下面四图给出了梁在中部同时受到集中力和弯矩作用下的等效应力,等效位移,剪力和弯矩图。
从剪力图知道,当给梁加载弯矩后,改变了梁的剪力和弯矩,不过在梁的中部,梁的弯矩图出现了突变,这与材料力学的分析变化趋势相一致。
图30等效应力图
图31等效位移图
图32轴线方向剪力图
图33轴线方向弯矩图
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ABAQUS 简支梁 分析 单元 实体
![提示](https://static.bdocx.com/images/bang_tan.gif)