有限元分析与应用详细例题

此为本人有限元作业，供学习交流用，严禁商业用途1《有限元分析与应用》详细例题试题1：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较：1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；2）分别采用不同数量的三节点常应变单元计算；3）当选常应变三角单元时，分别采用不同划分方案计算。一．问题描述及数学建模无限长的刚性地基上的三角形大坝受齐顶的水压作用可看作一个平面问题，简化为平面三角形受力问题，把无限长的地基看着平面三角形的底边受固定支座约束的作用，受力面的受力简化为受均布载荷的作用。二．建模及计算过程1.分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算下面简述三节点常应变单元有限元建模过程（其他类型的建模过程类似）：1.1进入ANSYS【开始】→【程序】→ANSYS10.0→ANSYSProductLauncher→changetheworkingdirectory→JobName:shiti1→Run1.2设置计算类型ANSYSMainMenu:Preferences→selectStructural→OK1.3选择单元类型单元是三节点常应变单元，可以用4节点退化表示。ANSYSMainMenu:Preprocessor→ElementType→Add/Edit/Delete→Add→selectSolidQuad4node42→OK(backtoElementTypeswindow)→Options…→selectK3:PlaneStrain→OK→Close(theElementTypewindow)1.4定义材料参数此为本人有限元作业，供学习交流用，严禁商业用途2材料为钢，可查找钢的参数并在有限元中定义，其中弹性模量E=210Gpa，泊松比v=0.3。ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic→inputEX:2.1e11,PRXY:0.3→OK1.5生成几何模型生成特征点ANSYSMainMenu:Preprocessor→Modeling→Create→Keypoints→InActiveCS→依次输入四个点的坐标：input:1(0,0),2(3,0),3(6,0),4(3,5),5(0,10),6(0,5)→OK生成坝体截面ANSYSMainMenu:Preprocessor→Modeling→Create→Areas→Arbitrary→ThroughKPS→依次连接1,2,6;2,3,4;2,4,6;4,5,6这三个特征点→OK1.6网格划分ANSYSMainMenu:Preprocessor→Meshing→MeshTool→(SizeControls)Global:Set→inputNDIV:1→OK→(backtothemeshtoolwindow)Mesh:Areas,Shape:Tri,Free→Mesh→PickAll(inPickingMenu)→Close(theMeshToolwindow)1.7模型施加约束分别给下底边和竖直的纵边施加x和y方向的约束ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Displacement→Onlines→选择底边→OK→select:ALLDOF→OK给斜边施加x方向的分布载荷ANSYS命令菜单栏:Parameters→Functions→Define/Edit→1)在下方的下拉列表框内选择x,作为设置的变量；2)在Result窗口中出现{X},写入所施加的载荷函数：1000*{X}；3)FileSave(文件扩展名：func)→返回：Parameters→Functions→Readfromfile：将需要的.func文件打开，任给一个参数名，它表示随之将施加的载荷→OK→ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Pressure→OnLines→拾取斜边；OK→在下拉列表框中，选择：Existingtable(来自用户定义的变量)→OK→选择需要的载荷参数名→OK1.8分析计算ANSYSMainMenu:Solution→Solve→CurrentLS→OK(toclosethesolveCurrentLoadStepwindow)→OK1.9结果显示确定当前数据为最后时间步的数据ANSYSMainMenu:GeneralPostproc→ReadResult→LastSet查看在外力作用下的变形ANSYSMainMenu:GeneralPostproc→PlotResults→DeformedShape→selectDef+Undeformed→OK查看节点位移分布情况ContourPlot→NodalSolu…→select:DOFsolution→Displacementvctorsum→Def+Undeformed→OK此为本人有限元作业，供学习交流用，严禁商业用途3查看节点应力分布情况ContourPlot→NodalSolu…→select:Stress→XYshearstress→Def+Undeformed→OK查看节点应力分布情况ContourPlot→NodalSolu…→select:Strain→XYshearstress→Def+Undeformed→OK1.10退出系统ANSYSUtilityMenu:File→Exit…→SaveEverything→OK三．结果分析三节点常应变单元（6个节点，4个单元）几何模型图变形图，节点位移图，节点应力图，节点应变图此为本人有限元作业，供学习交流用，严禁商业用途4六节点常应变单元（6个节点，4个单元）几何模型图变形图，节点位移图，节点应力图，节点应变图此为本人有限元作业，供学习交流用，严禁商业用途5分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算结果比较由实际情况可推知坝体X向的变形和所受应力都为正，上面的结果与实际结果基本相符。分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算结果比较可知，三节点常应变单元的结果好些。2.分别采用不同数量的三节点常应变单元计算三节点常应变单元（13个节点，14个单元）几何模型图单元划分方案变形大小应力大小应变大小值的比较分析三节点三角形单元DMX:0.109E-05SMX:0.109E-05DMX:0.109E-05SMN:2778SMX:8749DMX:0.109E-05SMN:0.344E-07SMX:0.108E-061.最大变形值小；2.最大应力值小；3.最大应变值小。六节点三角形单元DMX:0.289E-05SMX:0.289E-05DMX:0.289E-05SMN:-976.202SMX:11598DMX:0.289E-05SMN:-0.121E-07SMX:0.144E-061.最大变形值大；2.最大应力值大；3.最大应变值小。此为本人有限元作业，供学习交流用，严禁商业用途6变形图，节点位移图，节点应力图，节点应变图分别采用不同数量的三节点常应变单元计算结果由实际情况可推知坝体X向的变形和所受应力都为正，上面的结果与实际结果基本相符。分别采用不同数量的三节点常应变单元计算结果比较可得，节点多的比节点少的精确。单元划分方案位移大小应力大小应变大小值的比较分析三节点三角形单元（6节点4单元）DMX:0.109E-05SMX:0.109E-05DMX:0.109E-05SMN:2778SMX:8749DMX:0.109E-05SMN:0.344E-07SMX:0.108E-061.最大变形值小；2.最大应力值大；3.最大应变值大。三节点三角形单元（13节点14单元）DMX:0.178E-05SMX:0.178E-05DMX:0.178E-05SMN:990.156SMX:8000DMX:0.178E-05SMN:0.123E-07SMX:0.990E-071.最大变形值大；2.最大应力值小；3.最大应变值小。此为本人有限元作业，供学习交流用，严禁商业用途73.当选常应变三角单元时，分别采用不同划分方案三节点常应变单元（6个节点，4个单元）（另一种方案）几何模型图变形图，节点位移图，节点应力图，节点应变图此为本人有限元作业，供学习交流用，严禁商业用途8当选常应变三角单元时，分别采用不同划分方案结果比较由实际情况可推知坝体X向的变形和所受应力都为正，上面的结果与实际结果基本相符。分别采用不同方案相同节点相同单元的三节点常应变单元计算结果比较可得，第二种划分的方案精确。单元划分方案位移大小应力大小应变大小值的比较分析三节点三角形单元（6节点4单元）DMX:0.109E-05SMX:0.109E-05DMX:0.109E-05SMN:2778SMX:8749DMX:0.109E-05SMN:0.344E-07SMX:0.108E-061.最大变形值小；2.最大应力值大；3.最大应变值大。三节点三角形单元（6节点4单元）DMX:0.130E-05SMX:0.130E-05DMX:0.130E-05SMN:2529SMX:8582DMX:0.130E-05SMN:0.313E-07SMX:0.106E-061.最大变形值大；2.最大应力值小；3.最大应变值小。此为本人有限元作业，供学习交流用，严禁商业用途9试题2：图示为带方孔（边长为80mm）的悬臂梁，其上受部分均布载荷（p=10Kn/m）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置（即方位）进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为1mm，材料为钢）一．问题描述及数学建模悬臂梁受集中载荷以及均布载荷的作用可看作一个平面问题，简化为平面四边形受力问题，悬臂梁固定在墙上的部分看作是受全约束的作用，悬臂梁受力面的受集中载荷以及均布载荷的作用。二．建模及计算过程有限元建模选用Solid单元的8节点82单元建模，材料为钢，可查找钢的参数并在有限元中定义，其中弹性模量E=210Gpa，泊松比v=0.3。悬臂梁的左侧受全约束作用，同时梁上受集中载荷以及均布载荷的作用。三．结果分析带斜方孔的悬臂梁（450节点130单元）几何模型图此为本人有限元作业，供学习交流用，严禁商业用途10变形图，节点位移图，节点应力图，节点应变图带正方孔的悬臂梁（441节点127单元）几何模型图此为本人有限元作业，供学习交流用，严禁商业用途11变形图，节点位移图，节点应力图，节点应变图带圆形孔的悬臂梁（423节点121单元）几何模型图此为本人有限元作业，供学习交流用，严禁商业用途12变形图，节点位移图，节点应力图，节点应变图将三种方案进行比较由实际情况可推知坝体X向的变形和所受应力都为正，上面的结果与实际结果基本相符。根据以上分别采用不同方案的计算结果比较可得，带圆孔的悬臂梁的变形应力应变都最小，可得出先用圆孔设计最合理。模型方案位移大小应力大小应变大小值的比较分析带方孔（斜置）的悬臂梁DMX:0.003182SMX:0.003182DMX:0.003182SMN:-127762SMX:16794DMX:0.003182SMN:-0.158E-05SMX:0.208E-061.最大变形值中；2.最大应力值中；3.最大应变值中。带方孔（正置）的悬臂梁DMX:0.003193SMX:0.003193DMX:0.003193SMN:-144314SMX:56804DMX:0.003193SMN:-0.179E-06SMX:0.703E-061.最大变形值大；2.最大应力值大；3.最大应变值大。带圆孔的悬臂梁DMX:0.003163SMX:0.003163DMX:0.003163SMN:-1