西南石油大学CAE大作业

1.画出梁的剪力图和弯矩图某工字梁材料弹性模量2.1×105MPa，v=0.3，密度7850kg/m3，其结构尺寸及其受力见图1所示。用有限元法进行力学分析。要求：1.用APDL编程序，程序内标注拼音字母姓名和学号。2.优化工字梁的横截面尺寸，保证梁的最大变形小于1mm，梁的重量控制在5t左右。3.验算Ansys获取梁的重量与精确计算的重量的准确性。4.获取A、C点反力？画出剪力图、弯矩图？5.获取最大转角、最大弯矩位置及其大小？6.梁上最大应力、最大位移分析及其位置？7.获取梁横截面上的VonMises、变形分布云图。8.要求完成Word文件分析报告。(a)梁结构尺寸及其载荷分布(b)梁截面尺寸图1梁的有限元力学模型单位：mmMB=160kNmP=20kNq=20kN/mABCD8m2m2m7597595040915501.1APDL程序fini/clear/filnam,Beam188_exe/PLOPTS,DEFA!数值标示放图形窗口右边，不是放底部/prep7!et,1,beam188,,,2!mp,ex,1,2.1e11mp,nuxy,1,0.3!SECTYPE,1,BEAM,I,,2SECOFFSET,CENTSECDATA,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,2k,3,10k,4,12k,10,0,1!L,1,2L,2,3L,3,4!esize,0.2!latt,1,,1,,10,,1!lmesh,all!/eshape,1/soluDK,1,,,,0,UX,UY,UZ,ROTX,ROTY,,DK,3,,,,0,UY,UZ,ROTX,ROTY,,,!--------------------lsel,s,,,2,3,1esll,ssfbeam,all,1,pres,20e3FK,4,FY,-20e3FK,2,MZ,160e3alls!solve!/post1!/CONT,1,8,AUTO!ETABLE,Fzi,SMISC,5ETABLE,fzj,SMISC,18!-------------------------ETABLE,Mzi,SMISC,2ETABLE,mzj,SMISC,15!SMULT,Fzi,FZI,,-0.001,1,!剪力单位，kNSMULT,Fzj,FZJ,,-0.001,1,!PLLS,FZI,FZJ,1,0!-------------------------------SMULT,Mzi,MZI,,-0.001,1,!弯矩单位，kN.mSMULT,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,thenNe_Smax=i!整个梁上的最大应力Beam188对应的单元编号*endif*enddo!-----------ALLSELNSORT,U,SUM,0,0,ALL*GET,Umax_Num,SORT,0,IMAX!最大位移节点编号*GET,Umax_U,NODE,Umax_Num,U,SUM!最大位移!------------------------------------NSORT,ROT,Z,0,0,All*GET,Rot_Num,SORT,0,IMAX!最大转角节点编号*GET,RotzMax,NODE,rot_Num,ROT,Z!最大转角!!**!获取总重量density=7850!kg/m^3!ETABLE,evol,VOLU,SSUM*GET,vtot,SSUM,,ITEM,EVOLweight=Density*(vtot)/1000!梁重量(t)!--------------------------------------------------------N1=node(kx(1),ky(1),kz(1))!A节点N3=node(kx(3),ky(3),kz(3))!C节点*get,Fay,NODE,N1,RF,Fy!提取A节点y方向支座反力*get,Fcy,NODE,N3,RF,Fy!提取C节点y方向支座反力!----------------------------------Fay=fay/1000!关键点A的反力，kNFcy=fcy/1000!关键点C的反力，kN!-----------------------------Mz_max=Mzimax*if,Mzjmax,GE,Mzimax,thenMz_max=Mzjmax!最大弯矩，kN.m*endif!-------------Mz_min=Mzimin*if,Mzjmin,LE,Mzimin,thenMz_min=Mzjmin!最小弯矩，kN.m*endif!---------------S_MAX=nint(1000*S_MAX*(1e-6))/1000!最大应力MPaUmax_U=nint(1000*Umax_U*1000)/1000!最大位移，mmRotzMax=RotzMax*180/3.14!最大转角,度!-(删掉不用的变量)------Mzimax=Mzimin=Mzjmax=Mzjmin=I=vs=NEMIN=1.2结果分析图1.1剪力图图1.2弯矩图图1.3梁的Von-Mises应力图图1.4位移云图1.2.1优化工字梁的横截面尺寸，保证梁的最大变形小于1mm，梁的重量控制在5t左右优化前工字梁的横截面尺寸：W1=0.759，W2=0.759，W3=0.915，t1=0.05，t2=0.05，t3=0.04。此时：梁的最大变形为1.302mm，梁的重量为10.22t。优化后工字梁的横截面尺寸：W1=0.759/4，W2=0.759/4，W3=0.99/1，t1=0.05/2，t2=0.05/2，t3=0.05。此时：梁的最大变形为0.923mm，梁的重量为5.32t（注：后面计算所用梁的参数仍采用优化前的数据）。1.2.2验算Ansys获取梁的重量与精确计算的重量的准确性梁重Ansys计算(t)实际值(t)误差优化前10.220710.22070优化后5.325.320误差基本上为0，说明Ansys计算结果精度高。1.2.3获取A、C点反力表5-1反力节点ACFx00Fy721481.2.4剪力图、弯矩图图5-5剪力图图5-6弯矩图1.2.5最大转角、最大弯矩位置及其大小表5-2转角与弯矩名称位置大小最大转角节点571.48e-2度最大弯矩B点144KN.M1.2.6梁上最大应力、最大位移分析表5-3应力与位移名称大小最大应力0最大位移0.391mm1.2.7获取梁横截面上的VonMises、变形分布云图图5-7VonMises图图5-8变形分布云图2、带孔眼钢板结构的应力及其变形有一带孔眼钢板结构尺寸及其力学模型见图2所示，其厚度为25mm，材料弹性模量为2.1E5MPa，泊松比为0.3，材料屈服应力650MPa。1)P=30MPa，求孔眼内壁x的值(用路径曲线分析)，钢板内x、y、xy以及VonMises的应力场分布，形成Word文件分析报告。2)钢板开始屈服时P极限载荷为多少？钢板內最大VonMises应力发生在什么位置？孔眼内壁附近的位移变化以及各种应力变化等值线云图，形成分析报告。APDL程序语言Finish/clear/filnam,homework2a=150!mmb=150r=13h=25Ex=2.1e5!MPaux=0.3p=30!MPa/PREP7ET,1,PLANE42MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,exMPDATA,PRXY,1,,uxRECTNG,0,a,0,b,CYL4,,,rASBA,1,2LESIZE,10,,,20,0.2,,,,1LESIZE,5,,,10,,,,,1LESIZE,9,,,20,0.2,,,,1LESIZE,2,,,20,,,,,1LESIZE,3,,,20,,,,,1MSHAPE,0,2DMSHKEY,0AMESH,3!/eof/SOLDL,10,,SYMMDL,9,,SYMM,SFL,2,PRES,p,SOLVEFINISHyxR=13mm300mm300mmpp要求1.用APDL编程序，程序内标注拼音字母姓名和学号。2.结果要求Word文件分析报告。3.报告中曲线用Excell生成。图2.11/4钢板X方向所受应力σx图2.2P=30MPa时孔眼附近x的路径曲线当P=30MPa时，钢板的受力图：-100-80-60-40-200200510152025路径距离（mm）X方向应力（MPa）SX图2.3X方向所受应力图2.4钢板Y方向所受应力图2.5钢板所受剪切应力τxyσy(MPa)σy(MPa)图2.6VonMises应力云图欲使钢板屈服，最大应力必须达到材料的屈服极限，因分析是在弹性范围内，应力的变化服从线形分布，预测P为225.5Mpa。P=225.5Mpa带入程序，进入ANSYS运行分析得钢板的VonMises应力分布云图：图2.7VonMises应力云图(P=225.5Mpa)图2.8钢板开始屈服时位移变化（XY方向）图2.9钢板开始屈服时X方向所受应力图2.10钢板开始屈服时Y方向所受应力图2.11钢板开始屈服时所受剪切应力图2.12钢板开始屈服时VonMises应力云图3、管道层流和紊流CFD分析二维管道层流和紊流CFD分析，其结构尺寸见图3所示(单位mm)。流体介质为空气，空气密度1.225kg/m3，空气粘度1.7894e-05kg/m-s，出口压力为0。其中出口段长度L可以任意变化，问题：1)当L=200mm时，层流模拟，模拟雷诺数为90的层流问题，入口速度为0.025m/s。分析二维管道内的速度场分布。2)增加管道长度，即L=860mm时，紊流模拟，入口速度到1.27m/s，计算雷诺数大于4000时，用紊流模型求解，分析其速度场的变化和分布。要求1.直接在ICEM模块中建立实体模型，用结构化模块进行网格划分，边界层加密。2.在ICEM中建立CFD模型的边界条件，并存入Fleunt模型数据格式。3.将以上模型调入Ansys-Fluent模块中进行流体动力学分析。4.分析速度场和压力场变化规律。5.获取出口边界的速度变化曲线。6.结果要求Word文件分析报告，报告中曲线用Excell生成。100502565L外壁面内壁面进口过渡出口4、两块板间受压的O-形橡胶密封圈的接触分析两块无摩擦板间受压的O-形橡胶密封圈的接触问题如图4所示，橡胶材料参数：C10=0.01055!橡胶材料参数C10C01=0.01408!橡胶材料参数C01NU1=0.49967!橡胶泊松比DD=(1-2*NU1)/(C10+C01)要求：1.用APDL编程序编写计算模型。2.分析Disp=(0~8