有限元分析题及大作业题答案.docx
- 文档编号:28465817
- 上传时间:2023-07-14
- 格式:DOCX
- 页数:25
- 大小:712.40KB
有限元分析题及大作业题答案.docx
《有限元分析题及大作业题答案.docx》由会员分享,可在线阅读,更多相关《有限元分析题及大作业题答案.docx(25页珍藏版)》请在冰豆网上搜索。
有限元分析题及大作业题答案
有限元分析及应用作业报告
有限元分析及应用作业报告
一、问题描述
图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较:
1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;
2)分别采用不同数量的三节点常应变单元计算;
3)当选常应变三角单元时,分别采用不同划分方案计算。
4)
不同划分方案试例
、几何建模与分析
图1-2力学模型
满足平
两端
由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,面不受力,满足平面应变问题的载荷条件。
因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况
E=2.1e11,泊松比d=0.3
及方向如图1-2所示,建立几何模型,进行求解。
假设大坝的材料为钢,则其材料参数:
弹性模量
三、第1问的有限元建模
本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。
1)设置计算类型:
两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural
2)选择单元类型:
三节点常应变单元选择的类型是PLANE42(Quad4node42),该单元属于是四节点单元类型,
在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE18(Quad
8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。
因
研究的问题为平面应变问题,故对Elementbehavior(K3)设置为planestrain。
3)定义材料参数
4)生成几何模
a.生成特征点
b.生成坝体截面
5)网格化分:
划分网格时,拾取lineAB和lineBC进行SizeConrotls,设定inputNDIV为15;拾取lineAC,
设定inputNDIV为20,选择网格划分方式为Tri+Mapped,最后得到600个单元。
6)模型施加约束:
约束采用的是对底面BC全约束。
大坝所受载荷形式为Pressure,作用在AB面上,分析时施加在Lab上,方向水平向右,载荷大小沿Lab由小
到大均匀分布(见图1-2)。
以B为坐标原点,BA方向为纵轴y,则沿着y方向的受力大小可表示为:
P=ph二二g(10-y)=98000-98000*{Y}
(1)
其中P为水的密度,取g为9.8m/s2,可知Pmax为98000N,Pmin为0。
施加载荷时只需对Lab插入预先设置的载荷函数
(1)即可。
网格划分及约束受载情况如图1-3(a)和1-4(a)所示。
7)分析计算
8)结果显示
四、计算结果及结果分析
4.1计算结果
(1)三节点常应变单元
图1-3(a)常应变三节点单元的网格划分及约束受载图
1
毗呛肛aOLUT-ICiW
are=J
iwwx=■盘的E> .^GE-O1? ..15SE-C4 .231E-04 .2S4E-04 图1-3(b)常应变三节点单元的位移分布图 (2) 图1-3(c)常应变三节点单元的应力分布图 六节点三角形单元 ANSYS DEC2010 图1-4(a)六节点三角形单元网格划分及约束受载图 DEC262010 iotasi3« SUE=1 T,IWE.=J ASiiS=O [*&<=u2^2K-Q4 WH律=,292E-04 ..M5E-05..13CE-(H..155E-H 3^5E-C5-&7+E-C5..lt5E-01,.^37E-C^ 图1-4(b)六节点三角形单元的变形分布图 WC 8'CE.P=1 3U&=1 SEQV IM兰二.迂9齟一□理 SMW=bO(I1JB5 SWC=€070-11 -001385134^326^79T«Q4696&3^594 *744^20234S3^724172145607Q4^ 图1-4(c)六节点三角形单元的应力分布图 1-1所示。 根据以上位移和应力图,可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表 表1-1计算数据表 单元类型 最小位移(mm) 最大位移(mm) 最小应力(Pa) 最大应力(Pa) 常应变三节点单兀 0 0.0284 5461 392364 六节点三角形单元 0 0.0292 0.001385 607043 4.2结果分析 由以上各图和数据表可知,采用三节点和六节点的三角形单元分析计算: (1)最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力 和位移变化分布趋势相似,符合实际情况; (2)结果显示三节点和六节点单元分析出来的最大应力值相差较大,原因可能是B点产生了虚假应力,造成了最大应力值的不准确性。 (3)根据结果显示,最小三节点和六节点单元分析出来的最小应力值相差极为悬殊,结合理论分析,实际上A点不承受载荷,最小应力接近于零,显然六节点三角形单元分析在这一点上更准确。 (4)六节点的应力范围较大,所以可判断在单元数目相同的前提下,节点数目越多,分析精度就越大;但是节点数目的增多会不可避免地带来计算工作量增加和计算效率降低的问题。 五、第2问的有限元建模及计算结果 此次分析选择的单元类型为常应变三节点单元。 选用三种不同单元数目情况进行比较分析。 具体做法如下: 有限元建模步骤与第1小题类似,只是在划分网格时,对Lbc的AERA: SizeControls不同: 依次设置单元边长 度Elementedgelength为1.6、0.4、0.1,所获得的单元数目依次为96(图1-9(a))、1536(图1-10(a))、20880 (图1-11(a));分别计算并得到位移变化图如图1-9(b)、1-10(b)、1-11(c)所示;分别计算并得到应力变化云 图如图1-9(c)、1-10(c)、1-11(c)所示。 (1)单元数目为96的常应变三节点单元 ETiEJMfFHT煜 DEC262010 25i5Qs34 图1-9(a)单元数目为96的网格划分及约束受载图 1 总gALHC'LUTIGB sre.p^i 密帕=1 吃IMG,二工 USUU^i RSKS=O 3K==2ii2E-a^ SMCC=a2€^E-a4i ANSYS DEC26201025153: 33 „116£-D4..3J3E-(J4 胡孚1上一巧872£-0*5J45E-O电..JD4E-G4丄曲 8-re.p=jSUB=1 二J3EGVEWX=B.SMJ=2 =2 12SS6 DEC262010 120022 100^31159114 24S28S 275480 42747 图1-9(c)单元数目为96的应力分布图 (2)单元数目为1536的常应变三节点单元 1 EJLEMEHTfS DEC262010 ia: 4i: 4Q 图1-10(a)单元数目为1536的网格划分及约束受载图 即PEA]SHE=1TUWE.^1□soitf H3¥3=-n -.28? E--a4 泅工三・2S? &-0: 4 3.31E-05 ..&G4L-C5 ..225L-04 图1-10(b)单元数目为1536的位移分布图 ANSYS DEC262010 13: 01: 40 2730 S794Q 414013 49^624 (3)单元数目为20880的常应变三节点单元 ANSYS DEC262010 “詡0: 曲 图1-11(a)单元数目为20880的网格划分及约束受载图 HC-MLiaC-L>UT.JIC-N ^TT.P=1 邑HE=1 抽址=■・肥 BjS64BE-D5..IW-fi*_1«£-04胡5無-04 „.324^-05-9I73E-C5..2ME-C4 g're.p=i DEC262010it: ai: 12 銘5”203 804133 TUHE=1 3E^i^工毘呻 WX=,2S2&-04弓阳灿=925.-203訥垃=71fe4? .2 15^40il(3355477S70515935 239441839*4153SS747771641^2 图1-11(c)单元数目为20880的应力分布图 由以上不同单元数目的位移应力分布图可以看出,大坝截面所受位移和应力的变化趋势是相同的,最大应力都发生在坝底和水的交界点附近,最小应力发生在大坝顶端;最大变形位移也是发生在坝顶。 不同单元数目下计算的数据如表1-2所示。 表1-2不同单元数目下计算数据表 序号 单元数 最大位移(mm) 最小应力(Pa) 最大应力(Pa) 1 96 0.0262 13656 275480 2 1536 0.0289 2730 499624 3 20880 0.0292 925 716492 (4)结果分析 由以上分析结果可知: (1)随着单元数目的增加,最大位移变化不大,应力变化范围逐步增大; (2)随着单元数目的增加,即网格划分越密,分析的结果准确度将会提高;但是单元数目的增加和节点数目的增加都会造成计算量的增加和计算速度的下降的问题。 (3)对于本次计算结果,仍可能存在虚假应力,应力的准确值无法准确得出,只是网格划分越密,计算结果越精确。 所以减少虚假应力影响的措施之一就是增加单元的数目,提高网格划分的密度。 五、第3问的有限元建模及计算结果 由图1-1所示的划分方案可知,需采用手动划分网格: 首先创建6个节点,然后采用不同的方式连接节点创 建单元,从而分别得到两种不同的网格划分方式,见下图1-12所示。 对底边的三个节点施加全约束;载荷建立方程式并创建table;其他的处理方式与第1小题相同。 图1-12方案一和二的划分方案图 1-13(a)、1-13(b)、 有限元模型建立完成后进行求解,则可得到方案一和方案二的的位移图和应力图,如图 1-14(a)、1-14(b)所示。 ANSYS 3TrSP=J8UB=1 US1W R3¥3=01*利=f.107E-Q4鱼砒=.107R-Q4 1EF ;119£«C5 -.35IE-05 -5KE"C5 B3aE-05 DEC262010 39: 24 ..1C7E-04 图1-13(a)方案一网格划分方式下的位移图 3C-LUTICIU 3DB=1 T11ME.=1 dwt=r307e--0--i却=507735JHW=J5«173 ^0773741959751712104014^52 £24^4詰10^3291227SI1S6172 图1-13(b)方案一网格划分方式下的应力图 图1-14(a)方案二网格划分方式下的位移图 ANSYS DEC262010 17: 17: 47 »C>DALiSOIAJTIC'fel S'rEF=l3DB=1TEWE.-l 3EQV恨妬) UC€--J2SE.-fl42HN=7*772 SWX=3475€7 1QS23712^591W0L 1003711151031318351417567 图1-14(b)方案二网格划分方式下的应力图 由以上两种方案的位移和应力图可得出的最大位移和最小最大应力如表1-3所示: 表1-3方案一和方案二计算数据表 最大位移(mm) 最小应力(Pa) 最大应力(Pa) 万案一 0.0107 50773 156173 万案二 0.0128 76772 147567 由以上分析结果可知,由于方案一和二都只有四个单元,所以在计算应力和位移的时结果的准确度较低。 分析应力图可知,方案二得出的最大应力不在坝底和水的交界处,不符合实际情况,而方案一的最大应力所在位 置符合实际情况,所以总体来说,方案一的分析结果优于方案二。 原因是方案一具有整体几何保形性的单元数目多于方案二的数目。 六、总结和建议 通过以上分析情况可以看出,如果要使分析结果较为精确,单元的类型选择要恰当。 由第 (1)小问计算结 果可知,不同的单元类型会造成结果的不同,节点较多可以保证计算精度较高;由第 (2)小问的计算结果可知, 划分网格时,单元数目也不能太少,单元数目的增加也可以提高计算的精度;但是对于实际工程而言,采用较多节点的单元反而会增加计算的工作量,影响工作效率和经济性。 因此在保证网格划分大小适当和均匀的前提下,使应力集中处划的密集些,这样也能得到较为精确的结果。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元分析 作业题 答案