原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实体单元在后处理中积分求内力
实体单元在后处理中积分求内力??[2004-12-10] ?说明：本程序转自中华钢结构论坛三雨兄写的一个例子，觉得有用转过来研究了一下。 混凝土简支梁均布受载时，求跨中截面弯矩和1/4截面处剪力，模型计算命令流如下： finish /clear b0=200 h0=300 l0=3000 ec=3.3e5 p0=0.2 /prep7 et,1,solid95 mp,ex,1,ec mp,prxy,1,0.167 blc4,,,b0,h0,l0 /view,1,1,1,1 /ang,1 vplot lsel,s,loc,x,1,b0-1 lsel,r,loc,y,0 dl,all,,uy ksel,s,loc,x,0 ksel,r,loc,y,0 dk,all,ux ksel,r,loc,z,0 dk,all,uz asel,s,loc,y,h0 sfa,all,1,pres,p0 allsel,all esize,50 vsweep,all finish /solu solve finish /post1!进入后处理模块 pldisp 积分求内力的方法如下： 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=差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
文档评论（0）