练习实例-结构分析(中间孔).docVIP

... ... 平面问题：板中圆孔的应力集中 实例[1]：如图所示为承受双向拉伸的板件，其 中心位置有一个小圆孔，尺寸 (mm)如图所示。 5N/mm 2， 弹性模量 E＝2× 10 泊松比 v=0.3 拉伸载荷 :q=20N/mm 平板的厚度： t=20mm 解题思路分析： 属于平面应力问题 中心带孔，应使用 8 节点四边形单 元或三角形单元 注意单位：尺寸 mm，力 N，故应力 N/mm2 最大变形约为 0.001mm（忽略孔的影响） ,最大应力在孔的顶部和底部，大 小约为 3.9N/mm2，即 3.9MPa。依次检验有限元的分析结果。 5 可使用 PLANE82单元（平面 4 节点）；或单元 PLANE42（平面 4 节点） KEYOPT,1,3，,3 !设置为带厚度的平面问题 7 R,1,20, ! 设定实常数 No.1，厚度 20 一. 相关设置 设置 jobname 为 bracket、Title 为”a Example for Bracke，t 且”过滤参数为 Structural; 二.建立网格模型 创建几何模型：在 XY平面内建立一个矩形和圆，并用布尔 subtract 得到几何模型； 定义单元类型：选择 Structural Solid选项下的 Quad 8node 82 单元，确定返回最后自动得 到单元类型为 PLANE82。它是 8 节点的四边形单元，是平面 4 节点单元 PLANE42的高阶形 式，更适合有曲线边界的模型。对于本问题，我们需要有厚度的平面应力单元，只需单击 单元类型表中的 options 按钮，弹出 PLANE82选项设置窗口，如图所示，在 K3 对应的方框 中选择 ”Plane strs w/thk，“使得可以设置板的厚度。 定义实常数： PreprocessorReal ConstantsAdd/Edit/Delete。弹出如下对话框 ,输入板的厚 度 20 定义材料特性： EX＝ 200000，PRXY＝0.3 定义网格尺寸：网格边长＝ 25，划自由网格。如下页所示 PreprocessorMeshingSize CntrlsManual SizeAll Areas PreprocessorMeshingMeshAreasFree 或用 MeshTool 保存工作 网格模型图约束、载荷模型图 三. 加载和求解 定义分析类型： SolutionAnalysis TypeNew Analysis，设为 Static，即结构静态分析； 施加约束： SolutionDefine LoadsApplyStructuralDisplacementon Line，s用鼠标点选模 型最左侧边，并全部约束 (All DOF) 施加载荷：板右侧边缘上有一个背离平板的 20N/mm 的均布线载荷，则均布压力＝线载 荷除以板厚 20mm＝ 1N/mm2 。对对模型右侧边施加－1 的均布表面压力。 SolutionDefine LoadsApplyStructuralPressureOn Lines 求解： SolutionSolveCurrent LS 四.检查计算结果，观察收敛情况 ,决定是否修改网格模型 节点最大应力检查 由于最大应力点的应力值有解析解，因此可以检查该点的应力值。首先显示节点编号， 找出孔部对应的节点： Utility MenuPlotNodes（如果没有显示节点编号，则 Utility MenuPlotCtrlsNumbering ，?打开节点编号），记下与圆顶部对应的节点编号。如图所示 列出应力值： General PostprocList ResultsNodal SolutionStre，ss检查所要考察节点 的 SEQV 值（等效应力值），查结果可知其大小为 2.67N/mm 2，与手工计算的结果 3.9Mpa 有较大的差别，因此需要在孔的周围采用更小的网格尺寸才能获得更为精确的解。 修改网格模型：网格模型的修改方式很多，这里选择将孔周围的单元进行网格细化。运 用 MeshTool 的 Refine 完成，细化级别＝ 1（稍作细化） 。 结果如下页图所示。 孔周围要细化的单元 细化后的网格 2 重新计算：这时检查孔顶部最大等效应力（编号可能会改变）＝ 3.39N/mm 四. 后处理 绘制变形图： General PostprocPlot ResultsDeformed ShapeDef+undeformed 从图中可以看出孔的变形情况、整体变形情况，并且从图中左上角说明得知，最大位移＝ mm 绘制等效应力云图： General PostprocPlot ResultsContour plotNodal Solution，在弹出 的窗口中选 von Mise