直墙拱形有限元分析

1、有限元法上机实习作业(一)弹性力学平面问题有限元分析一、设计目的掌握用有限元方法分析弹性力学平面问题的步骤和技巧，熟练使用编制的程序解决工程 实际问题，并分析计算结果的正确性，分析应力分布和位移分布规律。二、设计任务某地下岩石洞室断面为直墙拱形，如图所示。岩石的弹性参数e = 1.0*105mpa, u = 0.25,在远处作用有水平均布荷载q”和竖向均布荷载qy,不计岩石重力。取5 (r+h)计算 范围，圆心以上取5r,圆心以下取5h按平面应变问题取二分之一洞室计算岩石的位移、应 变和应力，右边界和上边界为荷载，左边界和下边界为约束，并对计算结果进行必要的讨论。三、实习成果1设计条件图sii

2、.3直墙拱形断面洞室示意图序号学号r (m)h (m)a、(mpa)ck, (mpa)803.45.43.86.22.单元网格划分：输入tgen1.dat的数据见tgen1.dat,执行tgenl后的输 出结果见tgen3.out,绘制网格程序所需文件tgen3.dat见tgen3.dat,执行绘 图程序后的结果见网格图。3输入数据文件：输入simple.dat的数据见simple.dato4执行 simpl.exe 后的结果见 simple.outo5 分析讨论：(1)正确性验证：a.水平外边界上单元水平应力与水平分布荷载相近(误差不大于0.5mpa)o单元号水平应力分布荷载误差值49-3.

3、82226-3.8-0.0222650-3.81026-3.8-0.0102659-3.8013-3.8-0.001360-3.80785-3.8-0.078569-3.782093.80.0179170-3.8159-3.8-0.015979-3.765933.80.0340780-3.77337-3.80.02663893.73643-3.80.06366903.76839-3.80.0316199-3.77184380.02816100-3.835423.8-0.035421093.80903-3.8-0.009031103.85761-3.8-0.05761b竖直外边界上单元竖向应力与

4、竖向分布荷载相近（误差不大于0.5mpa）o单元号竖向应力分布荷载误差值119-6.20574-6.2-0.00574120-6.27578-6.2-0.07578129-6.15768620.04232130-6.2569762-0.05697139-6.13201620.06799140-6.1423-6.20.0577c静力平衡：计算模型竖向约束反力之和等于作用于模型上表面竖向荷载之和（误差不大于0.005mn）onodeycomp613.002521226.056683013.25127竖向约束反力之和105.404139竖向荷载之和105.4误差0.004139d.静力平衡：计算模型

5、水平约束反力之和等于作用于模型右表面水平荷载之和（误差不大于0.005mn）。nodex-comp7.19354 16.55773 18.15817 19.9283722.6834312.17408857.703248611.885878712.354098813.659328915.80728909.12363约束反力之和 167.20393水平荷载之和 167.20.00393误差（2）分析应力、位移分布乩绘出临近水平轴各单元水平应力与竖向应力对单元横向坐标的关系曲线。水平轴上方形心坐标/m200201015pxpy形心坐标/m20一io l0 5 0 5- 1 1- edlal'fpz水平轴下方px-»-py见图表1b.绘出临近竖直轴各单元水平应力与竖向应力对单元横向坐标的关系曲线。心以下px-»-py应力/mpa见图表2c绘出临近水平轴各节点水平位移与节点横向坐标的关系曲线。水平位移y水平位移ya淺活 >系横向坐标x见图表3d.绘出临近水平轴各节点水平位移与节点横向坐标的关系曲线。见图表4竖向位移30-200.0qe+0q-2.00e04 4-4.00e-04-6.00e-04竖向位移y-1.00e-03竖向位移y5-90节点竖向位移竖向坐标xo.ooe+oo -5.00e-8p