搜索
支架的线性静力分析实力本算例的实体是支架，一段牢固的焊接在大型结构上，支架的圆孔中穿过一个相对较软的杆件，材料的弹性模量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，弹