关于实体单元后处理中的求内力讨论

xialibaren19

1.采用实体单元计算(有时不得不采用实体单元)有其优点，但因实体单元无法直接得到结构的内力(M,N,Q)，所以对于混凝土结构想进行配筋计算就带来了难度，这是本题目提出的缘由。
2.考察了本论坛和其它几个论坛，似乎大家都知道用后处理中的“积分”，编制一些APDL，但又没有给出比较详细的说明和方法，这是本题目讨论的必要性。
3.我也不明白究竟如何处理，望高手讨论一番。
经过摸索，讨论如下：
1.积分命令流如下：
path,kz,2
ppath,1,,b0/2,0,l0/2
ppath,2,,b0/2,h0,l0/2
pdef,sigmz,s,z
plpath,sigmz
!跨中轴力nfor1=-2.55e-5
pcalc,intg,nfor,sigmz,yg,b0
*get,nfor1,path,,last,nfor
!跨中剪力qfor1=-4.039e-9
pdef,sigmyz,s,yz
plpath,sigmyz
pcalc,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,20
path,pname,2
ppath,1,,(i-1)*10,0,l0/4
ppath,2,,i*10,h0,l0/4
pdef,sigmyz,s,yz
pcalc,intg,qfor,sigmyz,yg,10
*get,qq(i),path,,last,qfor
*enddo
qz=0
*do,i,1,20
qz=qz+qq(i)
*enddo
*stat
2.关于pcalc,intg,varname,lab1,lab2,fact1的探讨
varname----积分结果变量
lab1,lab2--变量，可取xg,yg,zg(总体坐标),s(当前路径)及定义的变量(如上述中的md1);
fact1------系数
即varname=fact1×∫lab1×d(lab2)，积分路径为当前路径S的长度
3.对于较复杂的截面如何进行呢？例如箱形截面，可否用面积直接积分呢？并且path是断的是否也可以呢？

以上是三月雨在中华钢结构中的帖子

  我试了一下，如果不用xyz坐标，而用node号，更为简便：
path ,node ,m!m为个数，对于一个截面=2，只需给出两个端点号码即可！
ppath,1,78
ppath,2,1274
pdef,sigmyz,s,1
!以上定义截面起始终点号码分别为78，1274
pcalc,mult,md1,sigmz,s !相乘计算
pcalc,intg,mfor,md1,s,b0 !积分计算
即可

我的问题是：
对于一个隧洞断面，往往截面与xy不同向，我们积分采用的是pdef,sigmyz,s,1 ，还是：pdef,sigmyz,s,x或者pdef,sigmyz,s,y?
还有，得到的弯矩是不是我们计算所需要的弯矩，有人说，这个弯矩应该是对于中性轴的，但我看了一下规范，好像是关于截面受拉侧钢筋的（近似等于端点）
这两个问题，静听讨论

freddieaid

xialibaren19兄:
没用过，帮你顶一下