支架的线性静力分析实例 刘海宝.docx
- 文档编号:11982162
- 上传时间:2023-04-16
- 格式:DOCX
- 页数:28
- 大小:856.46KB
支架的线性静力分析实例 刘海宝.docx
《支架的线性静力分析实例 刘海宝.docx》由会员分享,可在线阅读,更多相关《支架的线性静力分析实例 刘海宝.docx(28页珍藏版)》请在冰豆网上搜索。
支架的线性静力分析实例刘海宝
支架的线性静力分析实力
本算例的实体是支架,一段牢固的焊接在大型结构上,支架的圆孔中穿过一个相对较软的杆件,材料的弹性模量E=210000MPa,泊松比u=0.3。
本算例的支架有以下两种工作情况:
1)杆件的一端受到沿Y轴负方向上的集中力F=2kN,其大小随时间变化。
2)除上述载荷外,支架的自由端还在局部区域上受到匀布的剪力Ps=36MPa。
要求确定这两种工况下支架绕度随时间变化的情况,以及内圆角处的最大主应力。
根据分析结构来改进设计,以减小应力集中。
建模时应注意以下要点:
1)此问题研究的是结构的静态响应,所以分析步类型应为Static,General(使用ABAQUS/Standard作为处理器)。
2)由于关心的是应力集中部位的应力状态,所以在模型中使用C3D20R单元(20节点六面体二次减缩积分实体单元)。
3)基于结构和载荷的对称性,可以只取模型的1/2进行分析。
4)由于圆孔处的应力应变不是所关心的重点,可以简化简化杆件和圆孔之间的关系,不精细建模,而只在杆件一端的受力点和圆孔内表面之间建立分布耦和约束。
1创建二维平面图
选择part功能模块后,单击
,随后会进入草图绘制环境,根据尺寸要求绘制二维平面图,完成的二维平面图如图1-1所示。
图1-1支架的二维平面图
1)创建部件
然后再试图区中点击鼠标中键两次来退出绘图环境,在弹出的editbaseextrusion对话框中将depth(拉伸长度)设为15,然后点OK。
然后点击窗口顶部工具栏中的
,不易透视图方式显示部件,看到的三维部件如图1-2所示。
图1-2拉伸后达到的部件
2)倒圆角,切割圆孔
点击左侧工作区中的
,点击部件中需要倒角的边,在试图区中点击鼠标中键来确认,在窗口底部提示区中的radius后输入5,按回车键来确认,得到的部件如图1-3所示。
图1-3倒圆角后的部件
3)切割圆孔
选择左侧工具区中的
,点击基准面,并选择相应的边以后进入二维绘图环境。
绘制切割轮廓即半圆,定义好尺寸后单击鼠标中键两次后退出绘图环境得到如图1-4所示的模型。
图1-4切割圆孔后的部件
2.创建材料和界面属性
1)创建材料在窗口左上角的module列表中选择property功能模块,点击左侧工具区中的
,弹出editmaterial对话框,mechanical-elasticity-elastic,在数据表中设置young`smodulus为210000,poisson`sratio为0.3,点击ok,如图1-5所示。
图1-5创建材料
2)创建界面属性点击左侧工具区中
,点击continue,在弹出的editsection对话框中,保持默认参数不变(material:
material-1;planestress/strainthickness:
1)。
如图1-6所示。
图1-6创建截面属性
3)给部件赋予截面属性点击左侧工具中的
,再点击部件,在试图区中点击鼠标中键,然后点击ok。
部件变为绿色,表明已经被赋予了界面属性,如图1-7所示。
图1-7赋予部件截面属性
3.定义装配件
在module列表中选择assembly功能模块。
点击左侧工具区中的
,接受默认参数instancetype:
dependent(meshonpart),即类型为非独立实体,点击ok。
4.划分网格
1)进入mesh功能模块在module列表中选择mesh功能模块,在窗口顶部的环境栏中把object选项设为part:
bracket,即对部件bracket划分网格,而不是对整个装配件划分网格。
如图1-8所示,当前部件显示为橙色,表明无法使用默认的网格划分技术来生成网格,需要把部件划分为几个简单区域再在每个区域上分别生成结构化网格或扫略网格。
图1-8为分割的部件
2)分割部件点击左侧工具区中的
,然后点击窗口底部提示区中的point&normal(通过指定分割面上的一点和法线来定义分割面)选择相应线段和中点,在视图区域点击鼠标中键来确定,得到的部件如图1-9所示,圆孔附近区域显示为黄色,表示此区域已经可生成扫略网格。
图1-9第一次分割后的部件
3)继续分割部件点击左侧工具区中的
上面按住鼠标左键不放,在展开的按钮栏中点击(partitioncell:
extrude/sweepedges),然后点继续继续分割的橙色区域,在视图区中点击鼠标中间来确认。
点击将要拉伸的线段,再次在试图区域点击鼠标中键来确认。
点击窗口底部提示区中的extruealongdirection,再点击拉伸方向线段,然后在视图区域点击鼠标中键两次来确认,得到的部件如图1-10所示,这是整个部件显示为黄色,表示所有区域都可以生成扫略网格。
图1-10完成第二次分割后的部件
4)设置全局种子点击左侧工具区中的
,在弹出的globalseeds对话框中,在approximateglobalsize后面输入2,点击ok,视图区中显示出全部种子。
5)布置边上的种子在左侧工具区中的
上按住鼠标左键不放,在展开的按钮栏中点击seededge:
bynumber,然后点击相应圆弧,在视图区中点击鼠标中键来确认。
在窗口底部的提示区中输入圆弧上的单元数12,按回车键来确认,如图1-11所示。
图1-11生成种子
6)设置网格参数点击左侧工具区中的
,按住鼠标左键,在视图区中画一个矩形框来选中部件的所有区域,然后在视图区中点击鼠标中键。
在弹出的meshcontrols对话框中,默认参数为
Elementshape(单元形状):
hex(六面体)
Techniques(网格划分技术):
sweep(扫略网格划分技术)
Algorithm(算法):
medialaxis(中性轴算法),选中minimizethemeshtransition(最小化网格过渡)
点击ok来确认上述参数。
如图1-12所示。
图1-12设置网格参数
7)设置单元类型点击左侧工具区中的
,画一个矩形框来选中部件的所有区域,然后在视图区中点击鼠标中键。
在弹出的elementtype对话框中将geometricorder(几何阶次)设为quadratic(二次单元),保持对reducedintegration(减缩积分)的选择,即当前单元类型为c3d20r(20节点六面体二次减缩积分单元),点击ok。
如图1-13所示。
图1-13设置单元类型
8)划分网格点击左侧工具区中
,然后在视图区中点击鼠标中键,得到如图1-14所示的网格。
图1-14使用medialaxis(中性轴算法)生成的网格
图1-15使用advancingfront(进阶算法)算法生成的网格
可以看出,尽管在所关心的应力集中部位布置了较密的种子,但此处的单元形状太狭长,这会对计算精度有很大影响,若使用advancingfront算法,得到的网格(图1-15)优于medialaxis算法,但在应力集中区域单元形状仍不规则。
更好的方法是通过分割,在应力集中区域生成结构化网格,具体方法如下:
9)在应力集中区域分割面点击左侧工具区中的
,然后点击拐角L型区域,在视图区中点击鼠标中键来确认。
这时会弹出一个信息框,提示将删除已有的网格,点击ok即可。
点击相应线段后自动进入绘图环境。
在应力集中区域绘制九十度扇形区域(扇形外圆弧半径为6),然后在试图区域中连续点击鼠标中键来退出绘图环境得到如图1-16所示部件。
图1-16分割面后的部件
10)在应力集中部位分割部件点击左侧工具区中的(partitioncell:
extrude/sweepedges),然后点击扇形分割面所在区域,在视图区中点击鼠标中键来确认。
按住shift键,依次点击扇形区域各边,再次在试图区域中点击鼠标中键来确认。
点击窗口底部提示区中的extrudealongdirection,然后点击部件中任意一条沿3轴方向的线段,点击鼠标中键两次来确认。
这样通过拉伸扇形区域,完成了在应力集中部位的分割。
分割后的扇形区域显示为绿色,表示在此区域可以使用结构化网格划分技术,如图1-17所示。
图1-17完成应力集中部位分割后的部件
11)在应力集中区域布置种子点击左侧工具区中的
,按住shift键,依次选中图1-16中扇形区域的两条圆弧,在视图区中点击鼠标中键,在窗口底部的提示区中设置单元格数为12,点击提示区右侧的constraints按钮,选中donotallowthenumberofelementstochange(种子受完全约束,不允许单元数目改变),点击ok,然后再次在视图区域中点击鼠标中键。
用同样的方法,选中图1-16中扇形区域的水平线和竖直线,设置单元数为3,并点击提示区右侧的constraints按钮,选中donotallowthenumberofelementstochange(种子受完全约束,不允许单元数目改变),布置种子后的部件如图1-18所示。
图1-18在应力集中区域布置受完全约束的种子
12)设置网格划分方式点击左侧工具区中的
,点击与扇形区域相邻的黄色区域,在试图区中点击鼠标中键,弹出meshcontrol对话框,将algorithm设为advancingfront,然后点击ok。
如图1-19所示。
图1-19设置网格划分方式
13)划分网格点击左侧工具区中的
,在视图区中点击鼠标中键,得到的网格如图1-20所示,可以看出在应力集中区域生成了形状规则的细画网格,其单元形状优于图1-14和1-15所示的网格。
图1-20分割应力集中区域后,使用advancingfront算法生成的网格
14)检验网格质量在主菜单中选择mesh-verify,画一个矩形框来选中所有单元,在弹出的verifymesh对话框中,将type设为analysischecks,然后点击highlight。
在模型中没有单元显示为黄色或红色,说明网格划分没有问题,窗口底部信息区中显示了所选择区域的单元总数,如图1-21所示。
图1-21检验网格质量
5.设置分析步
下面定义两个分析步,分别用来定义施加在杆件一端的集中载荷和施加在支架自由端上的剪力。
1)创建分析步在module列表中选择set功能模块,点击左侧工具区中的
。
在name后输入step-load-1,其余参数保持默认值,点击continue。
在弹出的editstep对话框中,点击incrementation标签,设置incrementsize:
initial(初始增量步大小)为0.2,然后点击ok。
2)创建第二个分析步再次点击左侧工具区中的
。
在name后输入step-load-2,其余参数保持默认值,点击continue。
在弹出的editstep对话框中,保持各参数的默认值不变(初始增量步大小为1),点击ok。
3)设置场变量输出结果在主菜单中选择output-fieldoutputrequests-manager,在弹出的fieldoutputrequestsmanager对话框中可以看到,abaquus/cae已经自动创建了一个名为F-output-1的场变量输出控制,它在分析步step-load-1中开始起作用,并自动延续到step-load-2中。
点击edit,在弹出的editfieldoutputrequests对话框中,点击stresses前面的黑色三角,在下拉选项中取消对VS,stressintheelastic-viscousnetwork和PS,stressintheelastic-plasticnetwork的选择,然后取消对strain,forces/reactions和contact的选择。
这样,分析过程中将指数出两种变量:
应力结果S和位移结果U。
点击OK,再点击dismiss。
6.定义耦合约束
模型中将不对杆件建模,而只是在杆件一端的受力点处定义一个参考点,然后在此点和圆孔内表面之间建立分布耦合约束,从而模拟杆件和圆孔之间的连接关系。
1)定义参考点在module列表中选择interaction功能模块,在主菜单中选择tools-referencepoint,在窗口底部的提示区中输入此点的坐标:
(30,20,15),然后按回车键确认。
试图去中显示出名为RP-1的参考点。
2)为参考点创建集合在主菜单中选择tools-set-manager,点击create,弹出createset对话框。
在name后输入set-point,点击continue,点击参考点,在视图区域中点击鼠标中键来确认。
3)定义受约束的面在主菜单中选择tools-surface-manager,点击create,弹出createsurface对话框。
在name后输入surf-hole,点击continue,再点击圆孔的内表面。
然后在视图区域点击鼠标中键来确认。
4)定义耦合约束点击左侧工具区中的
,弹出createconstraint对话框。
在name后输入constraint-hole,选择约束类型为coupling,点击continue。
点击窗口底部提示区右侧的sets按钮,在弹出的regionselection对话框中雪中set-point作为耦合约束的控制点,点击continue。
点击窗口底部提示区中间的surface按钮,在再次弹出的regionselection对话框中选中surf-hole来作为被约束面,点击continue。
这时弹出如图1-22所示的editconstraint对话框,设置couplingtype(耦合类型)为distributing,点击ok。
视图区中的模型将如图1-23所示。
图1-22设置耦合约束
图1-23定义耦合约束后的模型
5)保存模型点击窗口顶部工具栏中的
来保存所建模型。
7.定义载荷
穿过圆孔的杆件在一端受到沿Y轴负方向上的集中力,力的大小随时间变化,下面首先定义此集中力。
1)定义载荷随时间变化的幅值在module列表中选择load功能模块,点击主菜单tools-amplitude(幅值)-create,接受默认的幅值名称amp-1和类型tabular(表格),点击continue。
在弹出的editamplitude对话框中,输入如图1-24所示的分析步时间和幅值,然后点击ok。
图1-24定义载荷随时间变化的幅值
2)定义集中力在主菜单中选择load-manager,在弹出的loadmanager对话框中点击creat,弹出creatload对话框。
在name后输入loadpoint,将step设为step-load-1,载荷类型为默认值,然后点击continue。
在弹出的regionselecation对话框中,选中set-point,点击continue。
在弹出的editload对话框中,在CF2后输入-1000,将amplitude设为amp-1,然后点击ok。
在loadmanager对话框中可以看到,名为load-point的载荷在分析步step-load-1中开始起作用,并延续到分析步step-load-2中。
下面定义在支架自由端在局部区域上受到的匀布剪力Ps。
3)分割支架的端面剪力只作用在支架端面的局部区域上,因此要首先使用左侧工具区中的
来把端面分割成两部分。
如果在load功能模块中直接点击左侧工具区中的
,会出现错误信息:
无法对独立实体进行编辑。
这时应首先切换到mesh模块,在窗口顶部的环境中把object选项设为part:
bracket,然后在进行以下分割操作:
点击左侧工具区中的
,然后点击图1-23中的分割端面,在试图区中点击鼠标中键来确认。
这时会弹出一个信息窗口,提示将删除已有网格,点击ok即可。
点击图1-23中的相应线段,abaqus会自动计入绘图环境,在矩形区域的正中绘制竖直线,然后在视图区域中连续点击鼠标中键来退出绘图环境。
点击左侧工具区中的
,在视图区中点击鼠标中键,,重新划分网格。
4)定义受剪力的面返回load功能模块,在主菜单中选择tools-surface-create,在弹出的creatsurface对话框中,在name后输入surf-load,然后点击continue,点击划分好的小矩形,在视图区中点击鼠标中键来确认。
5)定义载荷面在loadmanager对话框中再次点击creat,在name后输入load-surface,设置step为step-load-2,载荷类型为surface-traction,然后点击continue。
在弹出的regionselection对话框中,选中surf-load,点击continue。
在弹出的editload对话框中,点击vectorbeforeprojection后面的edit按钮,在窗口底部提示区中可以看到向量的起始点坐标为(0.0,0.0,0.0),按回车键确认,然后再提示区中输入向量的终点坐标(0.0,-10.0,0.0),再次按回车键确认。
再重新出现的editload对话框中,设置magnitude为36,点击ok,得到如图1-25所示部件。
图1-25施加载荷后的部件
在loadmanager对话框中可以看到,名为load-surface的载荷在分析步step-load-2中开始起作用。
点击dismiss来退出loadmanager对话框。
8.定义边界条件
1)为施加固支边界条件的区域创建集合在load功能模块中,选择主菜单tools-set-manager,点击creat,在name后输入set-fix,点击continue,按住shift键,依次点击图1-26中的两个区域,在视图区中点击鼠标中键来确认。
图1-26施加固支边界条件的区域
2)为施加对称边界条件的区域创建集合在setmanager对话框中再次点击creat,在name后输入set-sym,点击continue。
按住shift键,依次点击图1-27所示的各个区域和参考点RP-1,在视图区中点击鼠标中键来确认。
在setmanager对话框中,点击dismiss。
图1-27施加对称边界条件的区域
3)定义固支边界条件在主菜单中选择BC-manager,在弹出的boundaryconditionmanager对话框中点击creat,弹出creatboundarycondition对话框。
在name后面输入BC-fix,将step设为initial,其余各参数保持默认值,点击continue。
在弹出的regionselection对话框中,选中set-fix,点击continue。
在弹出的editboundarycondition对话框中,选中ENCASTRE(U1=U2=U3=UR1=UR2=UR3=0),然后点击ok。
4)定义对称边界条件在boundaryconditionmanager对话框中再次点击creat,在name后面输入BC-sym,点击continue。
在editboundarycondition对话框中,选中ZSYMM(U3=UR1=UR2=0),然后点击ok。
得到的部件如图1-28所示。
图1-28施加边界条件后的部件
5)保存模型。
9.提交分析作业
1)创建分析作业在moudle列表中选择job功能模块,点击左侧工具区中的
,弹出jobmanager对话框。
点击create,在name后输入bracket-c3d20r,点击continue。
在弹出的editjob对话框中,点击memory标签页,abaqus的默认使用内存是256M,可以根据自己计算机的内存情况,适当增大此参数值,以加快分析速度。
2)提交分析作业这时已返回jobmanager对话框,点击submit,然后点击monitor来监控分析作业的运行状态,如图1-29所示。
图1-29分析作业运行状态
当jobmanager对话框中的status变为completed时,表示对模型的分析已经成功完成,点击results,进入visualization功能模块。
10.后处理
如上所述,在job功能模块完成分析计算后直接进入visualization功能模块进行后处理。
1)显示变形图
点击左侧工具区中的
,显示出变形后的网格模型。
如图1-30所示。
图1-30部件的变形图
2)显示云纹图
点击左侧工具区中的
,显示出最后一个分析步结束时的mises应力云纹图如图1-31。
图1-31部件的云纹图
3)逐格显示各个时间增量步
点击窗口底部提示区中的
或
,可以逐格各个时间增量步下的云纹图。
点击
或
可以直接跳至当前分析步的开始或结束时刻。
4)改变场变量
点击主菜单results-fieldoutput,在fieldoutput对话框中选择输出变量为invariant:
max.principle,点击ok。
视图区中显示出最大主应力的云纹图。
如图1-32所示。
图1-32最大主应力云纹图
5)查询节点上分析结果
点击窗口顶部工具栏中的
,在弹出的query对话框中选择probevalues,点击ok。
在弹出的probevalues对话框中,将probe设为nodes,选择s,max.principle,然后将鼠标移至应力集中处,节点的最大主应力就会在probevalues对话框中显示出来。
6)把查询结果写入文件
点击需要写入文件的节点,其查询结果就会在probevalues对话框的底部显示出来。
点击writetofile,弹出reportprobevalues对话框,在name后输入希望保存的文件名,然后点击ok,此文件会被保存在默认的工作目录下。
7)通过切面视图来观察模型内部的分析结果
在主菜单中选择tools-viewcut-manager,在viewcutmanager对话框中可以看到,abaqus/cae已经建立了三个基于全局坐标系的切面视图:
x-plane,y-plane,z-plane,点击其中的某个切面视图前面的小方框,视图区中会相应的显示出切面的效果。
拖动viewcutmanager对话框底部的滚动条,可以看到切面的位置会随之移动,点击此对话框中的creat,可以创建新的切面视图。
8)显示矢量图
点击主菜单results-fieldoutput,在fieldoutput对话框中选择输出变量为name:
U,点击ok。
点击左侧工具区中的
,视图区中显示出位移的矢量图,见图1-33。
图1-33显示位移的矢量图
9)用X-Y曲线图来显示位移随时间的变化
点击左侧工具区中的
来恢复云纹图显示。
在主菜单中选择tools-XYdata-manager,点击create,在createXYdata对话框中选择ODBfieldoutput,然后点击continue。
在XYdatafromODBfieldoutput对话框中将默认position:
integrationpoint改为position:
uniquenodal,点击U:
spatialdisplacement旁边的三角形,然后选中magnitude。
点击此对话框中的elements/no
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 支架的线性静力分析实例 刘海宝 支架 线性 静力 分析 实例 刘海