实体单元后处理求内力

lwd899

这是我做的一个实体单元后处理求内力的实际工程，发上来共同学习，共同探讨！绝对原创啊,希望版主加分！
!单位：kN,m
finish
/clear,start
/filname,'zhugong'
/prep7
et,1,solid45
mp,ex,1,3.20e7                        !C30
mp,prxy,1,0.2
mp,dens,1,25
/view,1,1,1,1
WPSTYLE,,,,,,,,1
wprot,0,90,0
BLC5,0,0,14,9.6,1
*get,maxkp,kp,,num,max
wpoff,0,0,1
CSYS,4
k,maxkp+1,6.8,4.6
k,maxkp+2,-6.8,4.6
k,maxkp+3,-6.8,-4.6
k,maxkp+4,6.8,-4.6
k,maxkp+5,4.35,3.10,1.50
k,maxkp+6,-4.35,3.10,1.50
k,maxkp+7,-4.35,-3.10,1.50
k,maxkp+8,4.35,-3.10,1.50
v,maxkp+1,maxkp+2,maxkp+3,maxkp+4,maxkp+5,maxkp+6,maxkp+7,maxkp+8
wpoff,0,0,1.5
*get,maxkp,kp,,num,max
k,maxkp+1,4.25,3.0
k,maxkp+2,2.25,3.0
k,maxkp+3,2.25,1.5
k,maxkp+4,4.25,1.5
         
k,maxkp+5,4.25,3.0,13.0
k,maxkp+6,2.25,3.0,13.0
k,maxkp+7,2.25,1.5,13.0
k,maxkp+8,4.25,1.5,13.0
v,maxkp+1,maxkp+2,maxkp+3,maxkp+4,maxkp+5,maxkp+6,maxkp+7,maxkp+8
*get,maxvolu,volu,,num,max
vgen,2,maxvolu,,,-6.5,
vsel,s,loc,z,0,20
vgen,2,all,,,,-4.5
wpoff,0,0,13.0
BLC5,0,0,10,8,2.5
allsel
vadd,all
!
vsbw,all
wpoff,0,0,-13
vsbw,all
wpoff,0,0,-1.5
vsbw,all
wpoff,0,0,17
CSYS,0
asel,s,loc,y,0
/pnum,area,1
*get,areanum,area,,num,max
hptcreate,area,areanum,,coord,-5.4,0,3.65   !HPTCREATE, TYPE, ENTITY, NHP, LABEL, VAL1, VAL2, VAL3(nhp=57)
hptcreate,area,areanum,,coord,-5.4,0,-3.65
hptcreate,area,areanum,,coord,5.4,0,3.65
hptcreate,area,areanum,,coord,5.4,0,-3.65
hptcreate,area,areanum,,coord,0,0,0
               
hptcreate,area,areanum,,coord,-2.25,0,1.5
hptcreate,area,areanum,,coord,-2.25,0,0
hptcreate,area,areanum,,coord,-2.25,0,-1.5
               
hptcreate,area,areanum,,coord,2.25,0,1.5
hptcreate,area,areanum,,coord,2.25,0,0
hptcreate,area,areanum,,coord,2.25,0,-1.5
allsel
NUMCMP,ALL
vatt,1,,1,0
MSHAPE,1,3d
MSHKEY,0
VMESH,all
nsel,s,loc,x,-5.4                          !定义组以方便后面加载
nsel,r,loc,y,0
nsel,r,loc,z,3.65
cm,cm1,node
nsel,s,loc,x,-5.4
nsel,r,loc,y,0
nsel,r,loc,z,-3.65
cm,cm2,node
nsel,s,loc,x,5.4
nsel,r,loc,y,0
nsel,r,loc,z,-3.65
cm,cm3,node
nsel,s,loc,x,5.4
nsel,r,loc,y,0
nsel,r,loc,z,3.65
cm,cm4,node
nsel,s,loc,x,0
nsel,r,loc,y,0
nsel,r,loc,z,0
cm,cm00,node
nsel,s,loc,x,-2.25
nsel,r,loc,y,0
nsel,r,loc,z,1.5
cm,cm01,node
nsel,s,loc,x,-2.25
nsel,r,loc,y,0
nsel,r,loc,z,0
cm,cm02,node
nsel,s,loc,x,-2.25
nsel,r,loc,y,0
nsel,r,loc,z,-1.5
cm,cm03,node
nsel,s,loc,x,2.25
nsel,r,loc,y,0
nsel,r,loc,z,1.5
cm,cm04,node
nsel,s,loc,x,2.25
nsel,r,loc,y,0
nsel,r,loc,z,0
cm,cm05,node
nsel,s,loc,x,2.25
nsel,r,loc,y,0
nsel,r,loc,z,-1.5
cm,cm06,node
asel,s,loc,y,-18
da,all,all                        !DA,P51X,ALL,
allsel
save
finish
/solu
antype,static
time,1                            !自重
allsel
acel,,9.8,
lswrite,01
time,2                            !双重主力
fdele,all
acel,,,
f,cm1,fy,-9430
f,cm2,fy,-9079
f,cm3,fy,-9079
f,cm4,fy,-9430
f,cm00,fx,594.0
f,cm00,fy,-3281.0
allsel
lswrite,02
time,3                             !单孔重载（主力）
fdele,all
acel,,,
f,cm1,fy,-9874
f,cm2,fy,-6560
f,cm3,fy,-6560
f,cm4,fy,-9874
f,cm00,fx,422
f,cm00,fy,-2653
allsel
lswrite,03
time,4                             !特载作用
fdele,all
acel,,,
f,cm1,fy,-6560
f,cm2,fy,-6560
f,cm3,fy,-6560
f,cm4,fy,-6560
f,cm00,fy,-5307.1
f,cm01,fy,-320
f,cm02,fy,-320
f,cm03,fy,-320
f,cm04,fy,-320
f,cm05,fy,-320
f,cm06,fy,-320
allsel
lswrite,04
lssolve,1,4
finish
/post1
lcdef,1,1                                !自重
lcdef,2,2                                !双重主力
lcdef,3,3                                !单孔重载
lcdef,4,4                                !特载   
lcfact,1,1.0
lcfact,2,1.0
lcfact,3,1.0
lcfact,4,1.0
lcase,1
lcoper,add,2
lcwrite,05
lcase,1
lcoper,add,3
lcwrite,06
lcase,1
lcoper,add,4
lcwrite,07
LCASE,5,                         !读取工况5
vsel,s,,,3
eslv,s                              
wpoff,2.25,3.00,-2.51
!asel,s,,,16
sucr,sufy1,cplane,3
sumap,s_x,s,x
sumap,s_y,s,y
sumap,s_z,s,z
sumap,s_xy,s,xy
sumap,s_yz,s,yz
sumap,s_xz,s,xz
supl,sufy1
supl,sufy1,s_
supl,sufy1,s_y
sueval,xforce,s_xy,intg
sueval,yforce,s_y,intg
sueval,zforce,s_yz,intg
sueval,a,da,sum
sueval,az,gcz,intg                         !gcz*da
az=az/a                                    !面重心至整体坐标轴x的距离
sueval,ax,gcx,intg                         !gcx*da
ax=ax/a                                    !面重心至整体坐标轴z的距离
sucalc,sygcz,s_y,mult,gcz
sueval,mz,sygcz,intg                       !绕整体坐标x轴的弯矩
sucalc,sygcx,s_y,mult,gcx
sueval,mx,sygcx,intg                       !绕整体坐标z轴的弯矩
mz0=mz-yforce*az
mx0=mx-yforce*ax
*stat

pfsr

感谢分享。如果能再加以说明，我给你加技术分。

civileol

如果不符合平截面假定时候，截面是一个曲面的时候 又该怎么样出内力呢

zimuyao

原帖由 civileol 于 2008-5-22 00:54 发表
如果不符合平截面假定时候，截面是一个曲面的时候 又该怎么样出内力呢

路径？

lwd899

我发的这个命令流就是我自己学习琢磨来的，用midas建模型验证过，二者得出的结果基本相同。最近由于太忙，好长时间没有来了，civileol兄的问题我也不太会，等到过段时间闲下来再研究研究，到时候我们共同探讨。

yz2003

非常感谢，对我帮助很大

jbsu

谢谢分享

bob410

为何加密？考下来是乱码嘛