西南石油大学CAE大作业.docx

1.画出梁的剪力图和弯矩图 某工字梁材料弹性模量2.1×105MPa，v=0.3，密度7850kg/m3，其结构尺寸及其受力见图1所示。用有限元法进行力学分析。要求： 1.用APDL编程序，程序内标注拼音字母姓名和学号。 2.优化工字梁的横截面尺寸，保证梁的最大变形小于1mm，梁的重量控制在5t左右。 3.验算Ansys获取梁的重量与精确计算的重量的准确性。 4.获取A、C点反力？画出剪力图、弯矩图？ 5.获取最大转角、最大弯矩位置及其大小？ 6.梁上最大应力、最大位移分析及其位置？ 7.获取梁横截面上的Von Mises、变形分布云图。 8.要求完成Word文件分析报告。 (a) 梁结构尺寸及其载荷分布 (b) 梁截面尺寸 图1 梁的有限元力学模型 1.1 APDL程序 fini /clear /filnam, Beam188_exe /PLOPTS, DEFA !数值标示放图形窗口右边，不是放底部 /prep7 ! et,1,beam188,,,2 ! mp,ex,1,2.1e11 mp,nuxy,1,0.3 ! SECTYPE, 1, BEAM, I, , 2 SECOFFSET, CENT SECDATA,0.759,0.759,0.915,0.05,0.05,0.04 ! 原始数据 !SECDATA,0.759/2,0.759/2,0.915/1.2,0.05/2,0.05/2,0.05 ! 优化数据 ! k,1, k,2,2 k,3,10 k,4,12 k,10,0,1 ! L,1,2 L,2,3 L,3,4 ! esize,0.2 ! latt,1,,1,,10,,1 ! lmesh,all ! /eshape,1 /solu DK,1, , , ,0,UX,UY,UZ,ROTX,ROTY, , DK,3, , , ,0,UY,UZ,ROTX,ROTY, , , !-------------------- lsel,s,,,2,3,1 esll,s sfbeam,all,1,pres,20e3 FK,4,FY,-20e3 FK,2,MZ,160e3 alls ! solve ! /post1 ! /CONT,1,8,AUTO ! ETABLE,Fzi,SMISC, 5 ETABLE,fzj,SMISC, 18 !------------------------- ETABLE,Mzi,SMISC, 2 ETABLE,mzj,SMISC, 15 ! SMULT,Fzi,FZI, ,-0.001,1, ! 剪力单位，kN SMULT,Fzj,FZJ, ,-0.001,1, !PLLS,FZI,FZJ,1,0 !------------------------------- SMULT,Mzi,MZI, ,-0.001,1, ! 弯矩单位，kN.m SMULT,Mzj,MZJ, ,-0.001,1, !-------------------------- PLLS,MZI,MZJ,1,0 ! ESORT,ETAB,Mzi,1,0 !按照单元表Mzi的大小排序 *GET,Mzimax,SORT,0,MAX ! 通过Mzmax得到最大值 *GET,Mzimin,SORT,0,Min ! 通过Mzmin得到最小值 ! ESORT,ETAB,Mzj,1,0 !按照单元表Mzi的大小排序 *GET,Mzjmax,SORT,0,MAX ! 通过Mzmax得到最大值 *GET,Mzjmin,SORT,0,Min ! 通过Mzmin得到最小值 !--------------------------------- *GET,S_MAX,SECR,ALL,S,EQV,MAX ! 整个梁上的最大Von_Mises应力Beam188 /eof *get,nnmax,node,,count *get,nemax,elem,,num,max ! 最大单元数 *get,nemin,elem,,num,min ! *do,i,1,nemax,1 *get,vs,secr,i,s,eqv,max *if,vs,eq,S_MAX,then Ne_Smax=i !整个梁上的最大应力Beam188对应的单元编号 *endif *enddo !----------- ALLSEL NSORT,U,SUM,0,0,ALL *GET,Umax_Num,SORT,0,IMAX ! 最大位移节点编号 *GET,Umax_U,NODE,Umax_Num,U,S