关于充气膜找形问题

题目：气承式膜结构平面形状为一边长25.4m的正方形，与x轴平行的两边固定，而与y轴平行的两边界可以沿z方向自由变形。膜材厚度t=0.38mm，膜的弹性模量为1.055*10^9N/m^2，泊松比0.4.内压为51pa；设定利用内压将薄膜充成一初应力为1.174*10^7N/m^2的充气结构。

我的程序如下：
/prep7
et,1,shell41
r,1,0.38e-3
mp,ex,1,1.055
mp,prxy,1,0.4
mp,alpx,1,6.5232
blc5,0,0,25.4,25.4
lesize,1,,,10
lesize,3,,,10
lesize,2,,,10
lesize,4,,,10
mshape,1,2d
mshkey,1
amesh,1
finish
/solu
antype,static
nlgeom,on
nropt,full
eqslv,front
dl,1,,all
dl,2,,ux
dl,4,,ux
dl,3,,all
tunif,-1
solve
finish
/solu
kbc,1
nsubst,10
sfa,all,,pres,51e-6
solve

原题给出结果如下：

66号节点显示的为最高点位移，而我算出来的节点位移和题目不一致，请问一下，我的程序哪里有问题，谢谢各位！！