机械零件有限元分析--2-概述2应用步骤举例.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,七、有限元分析的基本步骤,下面以,ANSYS,系统为例，用几个简单的例子介绍应用有限元法的步骤。,例,1,、,GUI,方式,问题描述：,对材料为软钢（弹性模量,:210000MPa,，泊松比为,0.3,）的悬臂梁施加,100,牛顿的集中载荷，用有限元法考察梁的,Y,向位移情况和梁内应力分布情况。具体情况如图所示。,悬臂梁示意图,90mm,10mm,10mm,5mm,载荷,100N,断面,步骤,1,、建模,经分析，此问题可用平面模型解决，因此可以在,ANSYS,中，用,90mm5mm,的面积来表示此悬臂梁。如下图所示。,步骤,2,、确定单元类型,8,节点矩形单元,指定单元属性,平面应力,+,厚度,设置单元厚度常数,厚度为,10mm,设置材料属性,步骤,3,、划分有限元网格,指定单元尺寸：,2mm,自动划分网格结果,自动划分网格上的节点,自动生成的节点编号（梁左端放大）,步骤,4,、定义约束条件,梁左面的节点,X,、,Y,自由度全,约束模拟悬臂梁情况,定义载荷,步骤,5,、运算求解,步骤,6,、后置处理,Y,方向位移等值图,X,方向应力等值图,各节点的应力值,例,2,、命令与命令流方式,21,问题描述,y,H,=10 mm,P,=100 N,Q,=1,MPa,x,L,=60 mm,W,=6 mm,Material,E,=200,Gpa,n,=0.3,22,步骤,1,、,建立,实体,模型,/PREP7,BLOCK,0,6,0,-5,5,-3,3,/PREP7,K,1,0,-5,-3,K,2,0,5,-3,K,3,0,5,3,K,4,0,-5,3,K,5,60,-5,-3,K,6,60,5,-3,K,7,60,5,3,K,8,60,-5,3,V,1,2,3,4,5,6,7,8,/VIEW,1,2,3,VPLOT,or,23,ANSYS,命令格式,各项目间以逗号隔开，逗号前后空格无关紧要,不分大小,写,命令必,须,在一行,内,完成,Command,Arg1,Arg2,.,24,步骤,2,、,建立,单元属性表,ET,1,SOLID45,MP,EX,1,200000,MP,NUXY,1,0.3,单元属性包括：,element type,material properties,real constants,element coordinate system.,25,步骤,3,：单元,切割,（网格化）,TYPE,1,MAT,1,ESIZE,3,VMESH,ALL,EPLOT,FINISH,26,步骤,4,：负载,描述,/SOLU,NSEL,S,LOC,X,0,D,ALL,ALL,0,NSEL,S,LOC,Y,5,SF,ALL,PRES,1,NSEL,ALL,N1=NODE(60,-,5,-3),N2=NODE(60,-5,3),F,N1,FY,-50,F,N2,FY,-50,EPLOT,27,步骤,5,：,解題,结构分析的数值解会存在,Jobname.RST,文档中,，,比如,FILE.RST,热分析的数值解会存在,Jobname.RTH,文档中,，如,FILE.RTH,SOLVE,FINISH,28,步骤,6,：后处理,(1/3),/POST1,/VIEW,0,0,1,PLDISP,2,ANDSCL,29,后处理,(2/3),PLNSOL,S,X,ANCNTR,30,后处理,(3/3),/VIEW,0,1,0,PLNSOL,S,X,31,ANSYS,宏程序（命令流）,FINISH,/CLEAR,!Units:mm,N,L=60 !Length,H=10 !Height,B=6 !Width,E=200000 !Youngs modulus,NU=0.3 !Poissons ratio,SIZE=3 !Element size,Q=1 !Distributed load,P=100 !Point load,/PREP7,K,1,0,-H/2,-B/2,K,2,0,H/2,-B/2,K,3,0,H/2,B/2,K,4,0,-H/2,B/2,K,5,L,-H/2,-B/2,K,6,L,H/2,-B/2,K,7,L,H/2,B/2,K,8,L,-H/2,B/2,V,1,2,3,4,5,6,7,8,/VIEW,1,2,3,VPLOT,ET,1,SOLID45,MP,EX,1,E,MP,NUXY,1,NU,TYPE,1,MAT,1,ESIZE,SIZE,VMESH,ALL,FINISH,/SOLU,!Fixed end,NSEL,S,LOC,X,0,D,ALL,ALL,0,!Distributed load,NSEL,S,LOC,Y,H/2,SF,ALL,PRES,Q,NSEL,ALL,!Loads at free end,N1=NODE(L,-H/2,-B/2),N2=NODE(L,-H/2,B/2),F,N1,FY,-P/2,F,N2,FY,-P/2,SOLVE,FINISH,/POST1,/VIEW,0,0,1,PLDISP,2,PLNSOL,S,X,32,About ANSYS Command Format,一个命令不可以分成两行,两个,命令可以放在一行,，,中,间,用“,$”,分,开,所有在“,!”,后,面的文字,视为注释,适当地插入,/EOF,及*,ASK,命令,例,3,、,APDL,方式,内容：复合材料纤维强化规律,ANSYS,有限元仿真研究,任务：,通过有,ANSYS,限元计算机仿真系统，建立不同类型的纤维强化树脂的模型，并进行受力后的模型变形及应力计算，考查不同的纤维形状和分布对复合材料的机械性能的影响。,具体方法如下图所示，在,ANSYS,中建立纤维强化树脂的断面模型（二维模型，简化计算），指定模型中的材料属性，进行有限元网格的划分，指定边界约束条件，施加载荷，进行计算，记录计算结果，改变模型数据再进行计算，直到多组数据计算完毕。整理计算数据画出刚度曲线。,基本情况,模拟试验试件厚度,1.6mm,，长度,22mm,，宽度,10mm,，；试件由两种材料纤维和树脂复合而成，其中纤维的弹性模量为,26400MPa,、纤维粗,0.05mm-0.1mm,；树脂的弹性模量为,500MPa,，纤维强化的复合材料的性能将与纤维在复合材料中的密度，纤维长度，纤维的粗细，纤维的排列方向有关系。,1,、确定有限元模型,为了加快运算速度，减小建模难度，根据模型简化法则，考虑使用二维模型，即用试件断面的纤维分布情况代替整个试件的纤维分布，产生的误差用增加试件的数量取算术平均值的办法解决。另外纤维的形状简化为矩形，粗细完全一致，方向完全一致。,2,、确定建模方式分析方式,因试件数量众多，且每个试件内的纤维数量很大且前后位置随机分布，因此用手工方式建模显然是不可取的，故考虑采用,ANSYS,的宏程序（,APDL-ANSYS Parametric Design Language,）自动建模，并自动指定材料属性，自动划分有限元网格，自动指定约束关系和自动加载，自动计算，自动将后置处理的结果写入文件。然后再通过对多个试件的结果文件的处理得到复合材料不同密度和不同纤维长度下的刚度曲线。,3,、具体建模方法,1,）、随机生成纤维断面,用,APDL,编程实现，纤维长度随机，位置随机，密度一定，下图为密度,0.62,，纤维平均长度,5mm,，纤维粗,0.2mm,，纤维间水平间距均值,0.2mm,。,2,）、生成基体材料范围,生成基体材料范围（试件轮廓）。,3,）、整理纤维范围,通过“与”命令，在指定的范围内（试件轮廓范围内）整理纤维。,4,）、生成试件几何模型并指定材料,再画试件基体范围矩形，用面积重合命令“,aovlap,”,最终生成试件，这样纤维和基体共用相同的边而形成一个整体。分别指定纤维和基体的材料属性。,5,）、划分有限元网格,局部放大,6,）、施加约束并加载,三点加力试验下的约束加载情况。,7,）、解算，并生成结果文件,将解算结果中感兴趣的部分生成结果文件（图形文件或文本文件）,Y,方向的位移等值图,X,方向的拉应力等值图,文本文件,8,）、处理结果文件,按三点加力刚度公式算出每个试件的弹性模量，,E=FL,3,/,（,4bh,3,U,Y,）,式中：,F,载荷,L,支点间距离,b,横截面宽度,h,横截面高度,U,Y,加载点,Y,向位移,注意,计算用,U,Y,的确定要考虑点载荷加载与实际情况的误差：,U,Y,=,U,Y,取点,-,（,U,Y,取点,-,U,Y,min,）,/2,U,Y,取点,处理所有试件数据，得到弹性模量与纤维密度关系曲线（可以使用,EXCEL,电子表格工具）。,改变纤维长度，可得弹性模量,纤维长度曲线,