apdl吧 关注:146贴子:300
  • 1回复贴,共1

一个关于几何非线性的增量计算的相关问题。求大神指点

只看楼主收藏回复

!主程序文件内容如下:
/clear
*DIM,A,,5,10
DF=-628.31853
!(以下)给出初始位置的节点坐标
X1=1 $ Y1=0.5*TAN(0.12217) $ X2=1.5 $ Y2=0
!(以下)设置节点的初始位移为零
UX1=0 $ UY1=0 $ UX2=0 $ UY2=0
!(以下)设置节点的初始增量位移为零
DUX1=0 $ DUY1=0 $ DUX2=0 $ DUY2=0
F=0
*DO,I,1,10
!(以下)计算节点的当前坐标
X1=X1+DUX1 $ Y1=Y1+DUY1 $ X2=X2+DUX2 $ Y2=Y2+DUY2
HONC
F=F+DF
!(以下)计算累计位移
UX1=UX1+DUX1 $ UY1=UY1+DUY1 $ UX2=UX2+DUX2 $ UY2=UY2+DUY2
A(1,I)=F,UX1,UY1,UX2,UY2
*ENDDO
*CFOPEN,OUT,dat
*DO,I,1,10,1
*VWRITE,A(1,I),A(2,I),A(3,I),A(4,I),A(5,I)
(5F12.5)
*ENDDO
*CFCLOSE
FINISH
!*******************************
!HONC.mac内容如下:
/PREP7
ANTYPE,STATIC
ET,1,SHELL61
R,1,2
MP,EX,1,3E7
MP,NUXY,1,0
N,1,X1,Y1
N,2,X2,Y2
E,1,2
D,2,UY,,,,,UZ
FINISH
/SOLU
F,1,FY,DF
SOLVE
FINISH
/POST1
*GET,DUX1,NODE,1,U,X
*GET,DUY1,NODE,1,U,Y
*GET,DUX2,NODE,2,U,X
*GET,DUY2,NODE,2,U,Y
FINISH
我的计算结果截图如下:

参考书中的计算结果截图如下:

问题:为什么我的计算结果与原书中的计算结果差别这么大,我仔细对照了命令流,没有发现错误,希望大神能给予指点?


1楼2017-07-07 11:24回复
    我已经找到问题的答案,这本书中把单位弄错了,所以F12.5这样的显示方式对于计算结果来说太短了,特别是小数部分,当把小数部分扩大到10,结果就有数据了。这个计算结果并不是很正确,因为作者把单位搞得很乱,但是命令流还是可以借鉴的


    2楼2017-07-07 11:59
    回复