有限元分析大作业报告要点.docx
《有限元分析大作业报告要点.docx》由会员分享，可在线阅读，更多相关《有限元分析大作业报告要点.docx（34页珍藏版）》请在冰点文库上搜索。
有限元分析大作业报告要点
船海1004黄山U201012278
有限元分析大作业报告
试题1：
一、问题描述及数学建模
图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较：
（1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；
（2）分别采用不同数量的三节点常应变单元计算；
（3）当选常应变三角单元时，分别采用不同划分方案计算。
该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图所示。
二、采用相同单元数目的三节点常应变单元和六节点三角形单元计算
1、有限元建模
（1）设置计算类型：
两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences为Structural
（2）选择单元类型：
三节点常应变单元选择的类型是SolidQuad4node182；六节点三角形单元选择的类型是SolidQuad8node183。
因研究的问题为平面应变问题，故对Elementbehavior（K3）设置为planestrain。
（3）定义材料参数：
弹性模量E=2.1e11，泊松比σ=0.3
（4）建几何模型：
生成特征点；生成坝体截面
（5）网格化分：
划分网格时，拾取lineAB和lineBC，设定inputNDIV为15；拾取lineAC，设定inputNDIV为20，选择网格划分方式为Tri+Mapped，最后得到600个单元。
1
船海1004黄山U201012278
（6）模型施加约束：
约束采用的是对底面BC全约束。
大坝所受载荷形式为Pressure，作用在AB面上，分析时施加在L上，方向水平向右，载荷大小沿L由小到大均匀分布。
以ABABB为坐标原点，BA方向为纵轴y，则沿着y方向的受力大小可表示为：
P?
gh?
gyY}*{?
）?
?
98000?
9800（10?
2、计算结果及结果分析
（1）三节点常应变单元
三节点常应变单元的位移分布图
三节点常应变单元的应力分布图
2
船海1004黄山U201012278
（2）六节点三角形单元
六节点三角形单元的变形分布图
六节点三角形单元的应力分布图计算数据
单元类型
最小位移（mm）
最大位移（mm）
最小应力（Pa）
最大应力（Pa）
三节点
0
0.0284
5460.7
392364
六节点
0
0.0292
0.001385
607043
（4）结果分析
3
船海1004黄山U201012278
①最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况；
B②结果显示三节点和六节点单元分析出来的最大应力值相差较大，原因可能是点产生了虚假应力，造成了最大应力值的不准确性。
根据结果显示，最小三节点和六节点单元分析出来的最小应力值相差极为悬③
殊，结合理论分析，实际上A点不承受载荷，最小应力接近于零，显然六节点三角形单元分析在这一点上更准确。
④六节点的应力范围较大，所以可判断在单元数目相同的前提下，节点数目越多，分析精度就越大；但是节点数目的增多会不可避免地带来计算工作量增加和计算效率降低的问题。
三、分别采用不同数量的三节点常应变单元计算
1、有限元建模（单元数目分别为150和1350）
2、计算结果及结果分析
（1）单元数目为150的常应变三节点单元
4
船海1004黄山U201012278
（2）单元数目为1350的常应变三节点单元
（3）计算数据表
单元数
最大位移（mm）
最小应力（Pa）
最大应力（Pa）
91
0.0270
10923.5
301924
1350
0.0288
3640.16
452618
果分析结4）（单元数目的增加，最大位移变化不大，应力变化范围逐步增大；网格划分越密，分析
的结果准确度将会提高；单元数目的增加和节点数目的增加都会造成计算量的增加和计算速度的下降的问题。
四、当选常应变三角单元时，分别采用不同划分方案计算1、方案一
5
船海1004黄山U201012278
方案二2、
6
船海1004黄山U201012278
3、计算数据表
最大位移（mm）
最小应力（Pa）
最大应力（Pa）
方案一
0.0128
76772
147567
方案二
0.0107
50772.9
156173
4、数据分析
由以上分析结果可知，由于方案一和二都只有四个单元，所以在计算应力和位移的时结果的准确度较低。
分析应力图可知，方案二得出的最大应力不在坝底和水的交界处，不符合实际情况，而方案一的最大应力所在位置符合实际情况，所以总体来说，方案一的分析结果优于方案二。
7
船海1004黄山U201012278
试题3：
一、问题描述及数学建模
图示为一带圆孔的单位厚度（1M）的正方形平板，在x方向作用均布压力0.25Mpa，试用三节点常应变单元和六节点三角形单元对平板进行有限元分析，并对以下几种计算方案的计算结果进行比较：
（1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；
（2）分别采用不同数量的三节点常应变单元计算；
?
?
MPa75.?
?
0?
0.25MPa孔边应力的精确解为：
y在轴上，孔边应力的精确解为：
轴上，在，xxy
y
0．25MPa
25MPa
0．x
6mm
48m
由图可知，本题所研究问题为平面应力问题，又因此平板结构关于图示中X、Y轴对称，可以利用此对称性，取截面的四分之一进行分析计算。
二、采用相同单元数目的三节点常应变单元和六节点三角形单元计算
1、三节点常应变单元
8
船海1004黄山U201012278
2、六节点三角形单元
3、计算结果及分析
DMX
SMN
SMX
0.225三节点0.309e-040.390e-050.309e-04-0.708
0.309e-04
0.391e-05
六节点0.309e-04
-0.780.258
理论值0.25
-0.75
9
船海1004黄山U201012278
由上表可看出，在单元数目相同的情况，六节点常应变三角形单元的分析精度要高于三节点常应变三角形单元。
所以，当单元形状和大小相同时，高阶单元的计算精度要高于低阶单元。
三、采用不同数量的三节点常应变单元计算
1、第一次加密
2、第二次加密
10
船海1004黄山U201012278
3、计算结果及分析
DMX
SMN
SMX
28
100?
单元数
板厚t=15mm
9m
22250.2250.309e-40.391e-50.309e-4-0.701不加密
5020-0.7260.393e-50.309e-40.234第一次加密0.309e-4
91400.309e-4-0.7410.247第二次加密0.309e-40.395e-5
0.25-0.75理论值
由上表可知，虽然常应变三角形单元的计算结果没有高阶单元的计算结果精确，但是随着单元数目的增多，计算结果逐渐的接近高阶单元的计算结果。
11
船海1004黄山U201012278
加筋板建模：
一、问题描述及数学建模
加筋板的几何模型如图所示。
横向加强筋
L150?
100?
10
纵向加强筋200?
206m
?
0.1Mpa的作用，求变形和应力。
四边简支的板，受到均布压力和要求：
使用shell63beam188单元。
有限元建模二、几何建模1、
用工作平面把平板按照加强筋的位置分割成如下图所示的几何模型。
12
船海1004黄山U201012278
2、属性定义：
单元类型——板：
shell63加强筋：
beam188
材料属性——杨氏模量E=2.1e+11，泊松比0.3
实常数——板厚0.02m
14?
250L100?
100?
12）——纵向加强筋；横向加强筋梁截面（Section?
26?
803、梁的方向点：
每一根横向加强筋和纵向加强筋都要定义一个方向点
4、改变线的方向：
改变线的方向的目的是改变梁的方向
5、有限元模型
不显示梁截面的有限元模型
不显示梁截面的位移云图
13
船海1004黄山U201012278
不显示梁截面的应力云图
显示梁截面的有限元模型
18016?
10?
200?
L250（截面形状）纵向加强筋截面形状（）横向加强筋ab?
120?
2414
船海1004黄山U201012278
显示梁截面的载荷模型
显示梁截面的应力云图显示梁截面的位移云图
结果分析三、
最大位移
最大应力
无加筋板
2.301m
1080Mpa
有加筋板
0.004233m
46Mpa
0.0472Mpa
最小应8.09Mpa
从结果中可以看出有加筋板的最大位移，最大应力都相对没有加筋板的小很多，分析：
承受载荷的能力显著提高，即加筋板可以有效的增强板材的结构强度，使变形程度大大减小，所以，工程中带加筋板的板材性能优越，应用广泛。
15
升级会员 