实体单元后处理中的求内力.docVIP

实体单元后处理中的求内力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,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 显示变形前和变形后的结构?????????? ???! 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，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-9pdef,sigmyz,s,yzplpath,sigmyzpcalc,intg,qfor,sigmyz,yg,1*get,qfor1,path,,last,qfor!跨中弯矩mfor1=差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,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.对于较复杂的截面如何进行呢？例如箱形截面，可否用面积直接积分呢？并且path是断的是否也可以呢？??请继续探讨，也不知对否，这里仅为抛砖引玉而已。 不好意思，差不多都忘记了。1.*get,mfor1中的mfor1是自己任意定义的变量