华科大有限元分析题及大作业题答案――船海专业_图文.docx

姓名:学号:班级: 有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变 单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较 1分别采用相同单元数冃的三节点常应变单元和六节点三角形单元计算 ； 2分别采用不同数量的三节点常应变单元计算； 3当选常应变三角单元时，分别采用不同划分方案计算。  、几何建模与分析 图1-2力学模型 由于大坝长度 横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作 无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用于大坝 上的载荷平行丁横截面且沿纵向方向均匀分布，两端面不受力,满足平面应变问题的 载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方 向如图1-2所示,建立几何模型，进行求解。 假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11 泊松比b =0.3 三、第1问的有限元建模 本题将分别采用相同单元数冃的三节点常应变单元和六节点三角形单元计算 1设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取 Preferences为 Structural 2选择单元类型：三节点常应变单元选择的类型是 PLANE42(Quad 4node42该单 元属于是四节点单元类型， 在网格划分时可以对节点数冃控制使其蜕化为三节点单元 ；六节点三角形单元 选择的类型是PLANE183(Quad 8node183该单元属于是八节点单元类型，在网格划 分时可以对节点数冃控制使其蜕化为六节点单元。因研究的问题为平面应变问题 ， 故对 Element behavior(K3 设置为 plane strain。 3定义材料参数 4生成几何模 生成特征点 生成坝体截面 5网格化分:划分网格时，拾取所有线段设定input NDIV为10,选择网格划分方 式为Tri+Mapped，最后得到200个单元。 6模型施加约束： 约束采用的是对底面BC全约束。 大坝所受载荷形式为Pressure作用在AB面上,分析时施加在L AB上，方向水平 向右，载荷大小沿L AB由小到大均匀分布(见图1-2。以B为坐标原点,BA方向为纵 轴y,则沿着y方向的受力大小可表示为： p (1 gh p- =p g 10 (* } 98000 98000 (Y y 其中p为水的密度，取g为9.8m/s2,可知P max为98000N,P min为0。施加载荷 时只需对L AB插入预先设置的载荷函数(1即可。 网格划分及约束受载情况如图1-3(a和1-4(a所示。 7分析计算 8结果显示 四、计算结果及结果分析 4.1计算结果 (1三节点常应变单元(4 node 42 图1-3(a常应变三节点单元的网格划分及约束受载图 AN APR 19 2D13 Norn soLirricw Norn soLirricw AN JLPR 2013 图1-3(b常应变三节点单元的位移分布图 AN JLPR 2013 (2六节点三角形单元 图1-4(a六节点三角形单元网格划分及约束受载图 图1-4(b六节点三角形单元的变形分布图 根据以上位移和应力图，可以得出常应变三节点单元和六节点三角形单元的最 小最大位移应力如表1-1所示 单元类型 院小位移(mm) 最大位移(mm〉 最小成力(Pm) 最大应力(Pa) 常应变三节点单元 0 0.02$ 11468 34S6S6 六节点三角形单元 0 0.0315 0.11988 538035 4.2结果分析 由以上各图和数据表可知，采用三节点和六节点的三角形单元分析计算： (1最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近,即坝底和水 的交界处，且整体应力 和位移变化分布趋势相似，符合实际情况； (2结果显示三节点和六节点单元分析出来的最大应力值相差较大，原因可能是 B点产生了虚假应力，造 成了最大应力值的不准确性。 (3根据结果显示，最小三节点和六节点单元分析出来的最小应力值相差极为悬 殊,结合理论分析，实际 上A点不承受载荷，最小应力接近于零，显然六节点三角形单元分析在这一点上 更准确。 (4六节点的应力范围较大，所以可判断在单元数冃相同的前提下，节点数冃越多， 分析精度就越大；但 是节点数冃的增多会不可避免地带来计算工作量增加和计算效率降低的问题。 五、第2问的有限元建模及计算结果 此次分析选择的单元类型为常应变三节点单元。选用三种不同单元数冃情况进 行比较分析。具体做法如下：有限元建模步骤与第1小题类似，只是在划分网格时，依 次设置NDIV值为5,10,50所获得的单元数冃依次为23(图1-9(a、80(图1-10(a、 1850(图1-11(a;分别计算并得到位移变化图如