有限元分析及应用大课后复习.docxVIP

,. 有限元分析及应用作业报告 ,. 目 录 有限元分析及应用作业报告 I 目 录 II 试题 1 1 一、问题描述 1 二、几何建模与分析 2 三、第 1 问的有限元建模及计算结果 2 四、第 2 问的有限元建模及计算结果 7 五、第 3 问的有限元建模及计算结果 13 六、总结和建议 16 试题 5 17 一、问题的描述 17 二、几何建模与分析 18 三、有限元建模及计算结果分析 18 四、总结和建议 26 试题 6 27 一、问题的描述 27 二、几何建模与分析 27 三、有限元建模及计算结果分析 27 五、总结和建议 35 ,. 试题 1 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应 变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比 较： 分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 分别采用不同数量的三节点常应变单元计算； 当选常应变三角单元时，分别采用不同划分方案计算。 M 0 1 6M 不同划分方案试例 图 1-1 模型示意图及划分方案 ,. 二、几何建模与分析 图 1-2 力学模型 由于大坝长度 横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看 作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用 于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应 变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分 布情况及方向如图 1-2 所示，建立几何模型，进行求解。 假设大坝的材料为钢，则其材料参数：弹性模量 E=2.1e11 ，泊松比σ =0.3 三、第 1 问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。 ）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取 Preferences  为  Structural 2 ） 选择单元类型：三节点常应变单元选择的类型是  PLANE42（ Quad 4node42  ），该 单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三 节点单元；六节点三角形单元选择的类型是  PLANE183（  Quad 8node183  ），该单 元属于是八节点单元类型，在网格划分时可以对节点数目控制使其蜕化为  六节 点单元。因研究的问题为平面应变问题，故对  Element behavior  （ K3 ）设置为 plane strain  。 3 ）定义材料参数：按以上假设大坝材料为钢，设定： Preprocessor → Material Props → Material Models → Elastic → Isotropic → input EX:2.1e11, PRXY:0.3  ANSYS Main Menu: → Structural → OK  → Linear ）生成几何模型： a. 生成特征点：  ANSYS Main Menu: Preprocessor  → Modeling  → Create  → Keypoints  → In Active CS  →依次输入三个点的坐标： ,. input:1(0,0),2(6,0),3(0,10) b. 生成坝体截面： → Areas → Arbitrary  → OK ANSYS Main Menu: Preprocessor → Through KPS →依次连接三个特征点，  → Modeling  → Create 1(0,0),2(6,0),3(0,10)  → OK 5 ）网格化分：划分网格时，拾取  lineAB  和  lineBC  进行  Size Conrotls  ，设定 input NDIV  为 15；拾取  lineAC  ，设定  input NDIV  为 20，选择网格划分方式为 Tri+Mapped  ，最后得到  600 个单元。 6) 模型施加约束： 约束采用的是对底面 BC 全约束。 大坝所受载荷形式为 Pressure ，作用在 AB 面上，分析时施加在 L AB 上，方向水 平向右，载荷大小沿 LAB 由小到大均匀分布（见图 1-2 ）。以 B 为坐标原点， BA 方 向为纵轴 y，则沿着 y 方向的受力大小可表示为： P gh  g(10  y)  98000  98000 * { Y} （ 1） 其中ρ为水的密度，取 g 为 9.8m/s 需对 L AB 插入预先设置的载荷函数（网格划分及约束受载情况如图  2 ，可知 Pmax 为 98000N ， Pmin 1）即可。