完整word版有限元分析大作业报告要点.docxVIP

船海 船海1004 黄山 U201012278 有限元分析大作业报告 试题2： 一、问题描述及数学建模 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变 单元和六节点三角形单元对坝体进行有限元分析，并对以下儿种计?算方案进行比 较： （1） 分别采用相同单元数LI的三节点常应变单元和六节点三角形单元讣算； （2） 分别釆用不同数量的三节点常应变单元计算； （3） 当选常应变三角单元时，分别采用不同划分方案计算。 该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图所示。  二、采用相同单元数目的三节点常应变单元和六节点三角形单元计 2、有限元建模 （1） 设置讣算类型：两者因儿何条件和载荷条件均满足平面应变问题，故均取 Preferences 为 Structural （2） 选择单元类型：三节点常应变单元选择的类型是Solid Quad 4 nodel82； 六节点三角形单元选择的类型是Solid Quad 8 nodel83o因研究的问题为平面 应变问题，故对 Element behavior （K3）设置为 plane strain。 （3） 定义材料参数:弹性模量E=2. lell,泊松比o=0.3 （4） 建儿何模型：生成特征点；生成坝体截面 （5） 网格化分：划分网格时，拾取lineAB和lineBC,设定input NDIV为15；拾取lineAC,设定input NDIV为20,选择网格划分方式为Tri+Mapped,最后 得到600个单元。 得到600个单元。 U201012278 模型施加约束：约束采用的是对底面BC全约束。大坝所受载荷形式为 Pressure,作用在AB面上，分析时施加在L上，方向水平向右，载荷大小沿L 由小到大均匀分布。以sB为坐标原点，BA方向为纵轴y,则沿着y方向的受力 大小可表示为： P?gh?gyY} * {?) ??98000?9800 (10? 2、计算结果及结果分析 AN(1)三节点常应变单元 AN NODAL SOLUTIOH 5IEE=1 SUB =1 TIKE-1U5UM (AVG) RS?S=0EMX -.234E-C4SMX =.234E-C4IPR 11 2013 5IEE=1 SUB =1 TIKE-1 U5UM (AVG) RS?S=0 EMX -.234E-C4 SMX =.234E-C4 IPR 11 2013 21:47:35 三节点常应变单元的位移分布图 NODAL 50LJII0N 5IEP=1 SUB =1 TIME=1 SEQV (AVG) EMX =.284E-04 SMK =5460.7 SMX =392364 APR 14 2013 21:48:23 三节点常应变单元的应力分布图 黄山船海1004 黄山 U201012278 （2）六节点三角形单元 NODAL 5OLUTIO1T STEP SOT =1 IIME=1 USCM (AVG| RSYS=0 =.292E-C4 SMX ?.292E-C4 MX AN APR 14 2313 21:59:15 0 .6492-05 .130E-04 .19SE-04 .260E-04 .32SE-0S .974H-0S .162E-04 .227H-04 .292E-04 六节点三角形单元的变形分布图 NODAL 50LJII0K 57EP=1 SUB =1 TIMZ=1 3EQV (AV3) DHX -.252Z-04 SKN =.001384 =607013 MN AN APR 14 2013 21：59：48 .002384 二 34399 269797 4X696 67449.3 202343 337246 472145 607043 （ （4）结果分析 单元类型 最小位移（mm） 最大位移（mm） 最小应力（Pa） 最大应力（Pa） 三节点 0 0. 0284 5460. 7 392364 六节点 0 0. 0292 0.001385 607043 六节点三角形单元的应力分布图 （3）计算数据表 3 黄山船海1004 黄山 U201012278 ①最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和 水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况； B②结果显示三节点和六节点单元分析岀来的最大应力值相差较大，原因可能是 点产生了虚假应力，造成了最大应力值的不准确性。 根据结果显示，最小三节点和六节点单元分析出来的最小应力值相差极为悬③ 殊，结合理论分析，实际上A点不承受载荷，最小应力接近于零，显然六节点三 角形单元分析在这一点上更准确。 ④ 六节点的应力范圉较大，所以可判断在单元数H相同的前提下，节点数H越 多，分析精度