版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、1.画出梁的剪力图和弯矩图某工字梁材料弹性模量2.1105MPa,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)
2、梁截面尺寸 图1 梁的有限元力学模型1.1 APDL程序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, , 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,2k,
3、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,1 esll,ssfbeam,all,1,pres,20e3 FK,4,FY,-20e3FK,2,MZ,160e3alls!solve!/post1!/CONT,1,8,AUTO!ETABLE,Fzi,SMISC, 5 ETABLE,fzj,SMISC, 18!-ETAB
4、LE,Mzi,SMISC, 2 ETABLE,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 ! 通过Mzmi
5、n得到最小值!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
6、,eq,S_MAX,thenNe_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,SUM ! 最大位移!-NSORT,ROT,Z,0,0,All *GET,Rot_Num,SORT,0,IMAX ! 最大转角节点编号 *GET,RotzMax ,NODE,rot_Num,ROT,Z ! 最大转角!*! 获取总重量density=7850 ! kg/m3!ETABLE,e
7、vol,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,
8、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
9、.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验
10、算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获取梁横截面上的Von Mises、变形分布云图图
11、5-7 Von Mises图图5-8 变形分布云图2、带孔眼钢板结构的应力及其变形有一带孔眼钢板结构尺寸及其力学模型见图2所示,其厚度为25mm,材料弹性模量为2.1E5MPa,泊松比为0.3,材料屈服应力650MPa。1) P=30MPa,求孔眼内壁sx的值(用路径曲线分析),钢板内sx、 sy、txy以及Von Mises的应力场分布,形成Word文件分析报告。2) 钢板开始屈服时P极限载荷为多少?钢板內最大Von Mises应力发生在什么位置?孔眼内壁附近的位移变化以及各种应力变化等值线云图,形成分析报告。要求1. 用APDL编程序,程序内标注拼音字母姓名和学号。2. 结果要求Word文
12、件分析报告。3. 报告中曲线用Excell生成。 APDL程序语言Finish/clear/filnam,homework2a=150 !mmb=150r=13h=25Ex=2.1e5 !MPaux=0.3p=30 !MPa/PREP7ET,1,PLANE42MPTEMP, MPTEMP,1,0 MPDATA,EX,1,ex MPDATA,PRXY,1,ux RECTNG,0,a,0,b, CYL4, , ,r ASBA,1, 2 LESIZE,10, , ,20,0.2, , , ,1 LESIZE,5, , ,10, , , , ,1 LESIZE,9, , ,20,0.2, , , ,1
13、 LESIZE,2, , ,20, , , , ,1 LESIZE,3, , ,20, , , , ,1 MSHAPE,0,2D MSHKEY,0AMESH, 3 !/eof /SOL DL,10, ,SYMM DL,9, ,SYMM, SFL,2,PRES,p,SOLVE FINISH 图2.1 1/4钢板 X方向所受应力x图2.2 P=30MPa时孔眼附近sx的路径曲线当P=30MPa时,钢板的受力图:图2.3 X方向所受应力图2.4 钢板Y方向所受应力图2.5 钢板所受剪切应力xy图2.6 Von Mises应力云图欲使钢板屈服,最大应力必须达到材料的屈服极限,因分析是在弹性范围内,应力
14、的变化服从线形分布,预测P为225.5Mpa。P=225.5Mpa带入程序,进入ANSYS运行分析得钢板的Von Mises应力分布云图:图2.7 Von Mises应力云图(P=225.5Mpa)图2.8 钢板开始屈服时位移变化(XY方向)图2.9 钢板开始屈服时X方向所受应力图2.10 钢板开始屈服时Y方向所受应力图2.11 钢板开始屈服时所受剪切应力图2.12 钢板开始屈服时Von Mises应力云图3、管道层流和紊流CFD分析二维管道层流和紊流CFD分析,其结构尺寸见图3所示(单位mm)。流体介质为空气,空气密度1.225kg/m3,空气粘度1.7894e-05 kg/m-s,出口压力
15、为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生成。4、两块板间受压的O-形橡胶密封圈的接触分析两块无摩擦板间受压的O-形橡胶密封圈的接触问题如图4所示,橡胶材
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 《柴油机燃料供给》课件
- 《畜禽免疫学》课件
- 中医诊断学课件-中医诊断学绪论
- 2024年高考历史总复习考前冲刺攻略 第4讲 高考应试能力的培养
- 单位管理制度合并汇编职员管理十篇
- 3号桥台(扩基U型桥台)计算书
- 《标书制作商务部分》课件
- 2016年高考语文试卷(新课标Ⅱ卷)(解析卷)
- 《寄小读者》导读课件
- 文明的历史足迹
- 东方电影学习通超星期末考试答案章节答案2024年
- 人教版四年级上册数学数学复习资料
- SB/T 10439-2007酱腌菜
- T∕ZZB 2665-2022 免洗手消毒凝胶
- 化粪池计算表格Excel(自动版)
- 2022年人美版美术六年级上册教案全一册
- 超外差调幅收音机课设报告——内蒙古工业大学
- 3.2熔化和凝固-人教版八年级上册课件(21张PPT)pptx
- 2017衢州新城吾悦广场开业安保方案
- 名师工作室考核评价表.doc
- 公司宣传品管理办法1
评论
0/150
提交评论