资源描述:
《有限元分析及应用大作业》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、有限元分析及应用作业报告II目录有限元分析及应用作业报告I目录II试题11一、问题描述1二、几何建模与分析2三、第1问的有限元建模及计算结果2四、第2问的有限元建模及计算结果7五、第3问的有限元建模及计算结果12六、总结和建议15试题516一、问题的描述16二、几何建模与分析16三、有限元建模及计算结果分析16四、总结和建议28试题629一、问题的描述29二、几何建模与分析29三、有限元建模及计算结果分析30五、总结和建议37II试题1一、问题描述图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元
2、对坝体进行有限元分析,并对以下几种计算方案进行比较:1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;2)分别采用不同数量的三节点常应变单元计算;3)当选常应变三角单元时,分别采用不同划分方案计算。图1-1模型示意图及划分方案35二、几何建模与分析图1-2力学模型由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝
3、所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3三、第1问的有限元建模本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。1)设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural2)选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类
4、型是PLANE183(Quad8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题,故对Elementbehavior(K3)设置为planestrain。3)定义材料参数:按以上假设大坝材料为钢,设定:ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic→inputEX:2.1e11,PRXY:0.3→OK4)生成几何模型:a.生成特征
5、点:ANSYSMainMenu:Preprocessor→Modeling→Create→Keypoints→InActiveCS→依次输入三个点的坐标:35input:1(0,0),2(6,0),3(0,10)→OKb.生成坝体截面:ANSYSMainMenu:Preprocessor→Modeling→Create→Areas→Arbitrary→ThroughKPS→依次连接三个特征点,1(0,0),2(6,0),3(0,10)→OK5)网格化分:划分网格时,拾取lineAB和lineBC进行SizeConrotls,设定input
6、NDIV为15;拾取lineAC,设定inputNDIV为20,选择网格划分方式为Tri+Mapped,最后得到600个单元。6)模型施加约束:约束采用的是对底面BC全约束。大坝所受载荷形式为Pressure,作用在AB面上,分析时施加在LAB上,方向水平向右,载荷大小沿LAB由小到大均匀分布(见图1-2)。以B为坐标原点,BA方向为纵轴y,则沿着y方向的受力大小可表示为:(1)其中ρ为水的密度,取g为9.8m/s2,可知Pmax为98000N,Pmin为0。施加载荷时只需对LAB插入预先设置的载荷函数(1)即可。网格划分及约束受载情况如
7、图1-3(a)和1-4(a)所示。7)分析计算ANSYSMainMenu:Solution→Solve→CurrentLS→OK(toclosethesolveCurrentLoadStepwindow)→OK8)结果显示ANSYSMainMenu:GeneralPostproc→PlotResults→DeformedShape…→selectDef+Undeformed→OK(backtoPlotResultswindow)→ContourPlot→NodalSolu…→select:DOFsolution,UX,UY,Def+Und
8、eformed,Stress,SX,SY,SZ,Def+Undeformed→OK四、计算结果及结果分析4.1计算结果(1)三节点常应变单元35图1-3(a)常应变三节点单元的网格划分及约束受