薄壁箱梁剪力滞效应解读.doc

作业3 一、题 采用有限元方法对教材P31页算例进行计算，具体分两个工况进行： （1）跨中截面腹板位置作用一对对称集中竖向荷载，荷载大小为P/2=225.5KN； （2）跨中截面腹板位置作用一对反对称集中竖向荷载，荷载大小为P/2=225.5KN。 分别计算跨中截面、1/4跨位置截面上的正应力和剪应力分布，并绘制相应的正应力和剪应力分布曲线。 ANSYS软件命令流法 三、ANSYS命令流 （1）工况 finish /clear /title,the analysis of simply supported box-girder !********前处理模块******** /prep7 !建立几何模型 k,1,0,0,0 k,2,0,0,2.4 k,3,0,0,7.1 k,4,0,0,9.5 k,5,0,-2.12,2.4 k,6,0,-2.12,7.1 kgen,2,all,,,40 a,1,2,8,7 a,2,3,9,8 a,3,4,10,9 a,6,5,11,12 a,5,2,8,11 a,3,6,12,9 !定义单元属性 et,1,shell63 r,1,0.22 r,2,0.34 r,3,0.30 mp,ex,1,3.40e10 !弹性模量 mp,gyz,1,1.445e10 !剪切模量 !赋予相应的单元属性和材料特性 aatt,1,1,1,,1 aatt,1,1,1,,2 aatt,1,1,1,,3 aatt,1,2,1,,4 aatt,1,3,1,,5 aatt,1,3,1,,6 !网格划分 mshape,0,2d !采用四边形网格 mshkey,1 !采用映射网格 esize,0. amesh,all finish !**********求解模块*********** /solu antype,static !在跨中腹板位置施加集中荷载 allsel f,node(20,0,2.4),fy,-225500 f,node(20,0,7.1),fy,-225500 !边界条件 allsel dk,5,ux,,,,,uy,uz,roty dk,6,ux,,,,,uy,roty dk,11,uy,,,,,uz,roty dk,12,uy,,,,,roty sbctran !把实体单元模型的荷载和边界条件，转化到有限元几何模型中 solve finish !*********后处理模块************ /post1 !梁的变形 allsel pldisp,2 !查看跨中截面正应力 allsel nsel,s,loc,x,19.79,20.01 esln,s plnsol,s,x !路径方法得到跨中截面正应力分布曲线和数据 path,zengyingli_a,2,,50 ppath,1,,20,0,0 ppath,2,,20,0,9.5 pdef,sx,s,x plpath,sx path,zengyingli_b,2,,50 ppath,1,,20,0,0 ppath,2,,20,0,9.5 pdef,sx,s,x prpath,sx !查看跨中截面剪应力 allsel nsel,s,loc,x,19.79,20.01 esln,s plnsol,s,xy !路径方法得到跨中截面剪应力分布曲线和数据 path,jianyingli_a,2,,50 ppath,1,,20,0,0 ppath,2,,20,0,9.5 pdef,sxy,s,xy plpath,sxy path,jianyingli_b,2,,50 ppath,1,,20,0,0 ppath,2,,20,0,9.5 pdef,sxy,s,xy prpath,sxy !查看1/4跨截面正应力 allsel nsel,s,loc,x,9.79,10.01 esln,s plnsol,s,x !路径方法得到1/4跨正应力分布曲线和数据 path,zengyingli_a,2,,50 ppath,1,,10,0,0 ppath,2,,10,0,9.5 pdef,sx,s,x plpath,sx path,zengyingli_b,2,,50 ppath,1,,10,0,0 ppath,2,,10,0,9.5 pdef,sx,s,x prpath,sx !查看1/4跨截面剪应力 allsel nsel,s,loc,x,9.79,10.01 esln,s plnsol,s,xy !路径方法得到1/4跨截面剪应力分布曲线和数据 path,jianyingli_a,2,,50 ppath,1,,10,0,0 ppath,2,,10,0,9.5 pdef,sxy,s,xy plpath,sxy