自己做的ansys实例-建筑楼板分析

欢迎各位大师指正，抛砖引玉，讨论讨论深层次的问题。
!板梁内力计算
finish
/clear
/prep7
et,1,shell63
Mp,ex,1,3e10
Mp,prxy,1,0.3
r,1,0.15
et,2,beam189
sectype,1,beam,rect$secdata,0.2,0.5,4,8
sectype,2,beam,rect$secdata,0.20,0.6,4,8
sectype,3,beam,rect$secdata,0.20,0.65,4,8
sectype,4,beam,rect$secdata,0.1,0.1,2,2

k,1$k,2,,6.6$k,3,5.7,6.6$k,4,5.7,3
k,5,8.4,3$k,6,5.7,1.5$k,7,8.4,1.5$k,8,5.7,0$k,9,8.4,0
k,10,0,3
a,1,10,4,8$a,10,2,3,4$a,4,5,7,6$a,8,6,7,9
wprota,,,90$wpoffs,,,1.85$asbw,all
wpoffs,,,2$asbw,all
wprota,,90$wpoffs,,,-1.2$asbw,all
wpoffs,,,-1.2$asbw,all
aglue,all
wpcsys,-1$/eshape,1$/pnum,line,1
lsel,s,loc,x,0$latt,1,,2,,,,3$cm,cmsel1,line
lsel,s,loc,x,5.7$latt,1,,2,,,,1$cm,cmsel2,line
lsel,s,loc,x,8.4$latt,1,,2,,,,1$cm,cmsel3,line
lsel,s,loc,y,0$latt,1,,2,,,,1$cm,cmsel4,line
lsel,s,loc,y,6.6$latt,1,,2,,,,2$cm,cmsel5,line
lsel,s,line,,8,10,2$latt,1,,2,,,,1$cm,cmsel6,line
allsel,all
cmsel,u,cmsel1,line$cmsel,u,cmsel2,line$cmsel,u,cmsel3,line
cmsel,u,cmsel4,line$cmsel,u,cmsel5,line$cmsel,u,cmsel6,line
latt,1,,2,,,,4$allsel,all
lesize,all,0.05$lmesh,all$amesh,all

dk,1,all,0$dk,2,all,0$dk,3,all,0$dk,4,all,0
dk,5,all,0$dk,7,all,0$dk,8,all,0$dk,9,all,0
dl,7,,all,0$dl,8,,all,0
allsel,all
sfa,all,,pres,-10e3$sftran
/psf,pres,norm,2,0,on$/pnum,kp,1

nn1=node(1.85,1.2,0)$nn2=node(3.85,1.2,0)
nn3=node(1.85,2.2,0)$nn4=node(3.85,2.2,0)

f,nn1,fz,10000$f,nn2,fz,10000$f,nn3,fz,10000$f,nn4,fz,10000
allsel,all
/solu
solve