资源描述:
《实体单元后处理中的求内力》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实体单元后处理中的求内力1.采用实体单元计算(有时不得不采用实体单元)有其优点,但因实体单元无法直接得到结构的内力(M,N,Q),所以对于混凝土结构想进行配筋计算就带来了难度,这是本题目提出的缘由。2.考察了本论坛和其它几个论坛,似乎大家都知道用后处理中的“积分”,编制一些APDL,但又没有给出比较详细的说明和方法,这是本题目讨论的必要性。 为方便讨论,这里做了个简支梁受均布荷载的分析。finish/clearb0=200h0=300l0=3000ec=3.3e5p0=0.2/prep7et,1,solid95mp,ex,1,ecmp,prxy,1,0.167blc4,,,b0
2、,h0,l0/view,1,1,1,1/ang,1vplotlsel,s,loc,x,1,b0-1lsel,r,loc,y,0dl,all,,uyksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxksel,r,loc,z,0dk,all,uzasel,s,loc,y,h0sfa,all,1,pres,p0allsel,allesize,50vsweep,allfinish/solusolvefinish/post1pldisp,2!显示变形的结构 ! kund:0仅显示变形后的结构 ! 1显示变形前和变形后的结构
3、 !2显示变形结构和未变形结构的边缘!------可以在下面添加语句完成(要求使用积分命令)!A.跨中弯矩!B.1/4跨剪力经过摸索,讨论如下:1.积分命令流如下:path,kz,2 ! 定义路径属性,PATH,NAME,nPts,nSets,nDivppath,1,,b0/2,0,l0/2! PPATH,POINT,NODE,X,Y,Z,CS(定义路径上的点)ppath,2,,b0/2,h0,l0/2pdef,sigmz,s,z! PDEF,Lab,Item,Comp,Avglab(将数据影射到路径上) plpath,sigmz! PLPATH,Lab1
4、,Lab2,Lab3,Lab4,Lab5,Lab6! (绘制路径上一数据项关于另一数据项的变化曲线图) !跨中轴力nfor1=-2.55e-5pcalc,intg,nfor,sigmz,yg,b0! PCALC,Oper,LabR,Lab1,Lab2,FACT1,FACT2,CONST ! (路径上数据的数学运算)INTG(积分)、(自然对数)、*get,nfor1,path,,last,nfor! *GET,Par,NODE,N,Item1,IT1NUM,Item2,IT2NUM! (得到数值并将其存储为标量参数或者数组中元素) !跨中剪力qfor1=-4.039e-9pde
5、f,sigmyz,s,yzplpath,sigmyzpcalc,intg,qfor,sigmyz,yg,1*get,qfor1,path,,last,qfor!跨中弯矩mfor1=-45251688误差0.6%pcalc,mult,md1,sigmz,s!相乘计算pcalc,intg,mfor,md1,s,b0!积分计算*get,mfor1,path,,last,mfor!1/4截面剪力qz(分条计算)=-30083.11误差0.3%*dim,qq,,20*do,i,1,20path,pname,2ppath,1,,(i-1)*10,0,l0/4ppath,2,,i*10,h0
6、,l0/4pdef,sigmyz,s,yzpcalc,intg,qfor,sigmyz,yg,10*get,qq(i),path,,last,qfor*enddoqz=0*do,i,1,20qz=qz+qq(i)*enddo*stat2.关于pcalc,intg,varname,lab1,lab2,fact1的探讨varname----积分结果变量lab1,lab2--变量,可取xg,yg,zg(总体坐标),s(当前路径)及定义的变量(如上述中的md1);fact1------系数即varname=fact1×∫lab1×d(lab2),积分路径为当前路径S的长度 3.对于较
7、复杂的截面如何进行呢?例如箱形截面,可否用面积直接积分呢?并且path是断的是否也可以呢? 请继续探讨,也不知对否,这里仅为抛砖引玉而已。不好意思,差不多都忘记了。1.*get,mfor1中的mfor1是自己任意定义的变量名;2.yg是整体坐标的标识符,可参见哪个讨论;轴力和剪力可沿Y轴积分得到;而M积分时,使用的是路径,所以用s符号表示。 该例子很简单,应该可以读懂,只要你仔细分析每一句即可明白。1.跨中弯矩计算中的几个语句解释pcalc,mult,md1,sigmz,s!相乘计算,其