发一个双支点平衡舵的有限元命令流！

xiaofanfu

/filname,rudder
/title,rudder

/PREP7  
ET,1,BEAM3  
R,1,63.617e-4,322.062e-8,0.09
R,2,8320e-6,6646.499e-8,0.3  
R,3,0.0132,1401.985e-8,0.13  
R,4,9.50331e-3,718.688e-8,0.11
R,5,27007722  

N,1,,,,,,,  
N,3,0.05,,,,,,  
N,12,2.2,,,,,,  
N,18,2.85,,,,,,
N,30,4.85,,,,,,
N,31,0,-0.6,,,,,

FILL,1,3
FILL,3,12
FILL,12,18
FILL,18,30

ET,1,beam3
ET,2,COMBIN14   

MP,EX,1,2.07e11
MP,PRXY,1,0.3  
MP,DENS,1,7800

TYPE, 1   
MAT, 1
REAL,1   
ESYS,0
e,1,2
e,2,3   

TYPE,   1   
MAT,       1
REAL,       2   
ESYS,       0   
e,3,4
egen,9,1,3

TYPE,   1   
MAT,       1
REAL,       3   
ESYS,       0
e,12,13
egen,6,1,12

type,1
mat,1
real,4
esys,0
e,18,19
egen,12,1,18

TYPE,   2   
REAL,       5   
ESYS,       0
e,1,31  

/sol
antype,static

d,31,all
d,18,ux
d,18,uy
d,30,ux
d,30,uy

*do,i,3,8
sfbeam,i,1,pres,26507.5,26507.5
*enddo

solve
finish
  
/POST1  
SET,LAST
ETABLE,imoment,SMISC,6  
ETABLE,jmoment,SMISC,12
ETABLE,ishear,SMISC,2   
ETABLE,jshear,SMISC,8   

PLLS,IMOMENT,JMOMENT,1,0
PRRSOL,
SAVE

jorgenhojmann

楼上的
你用的是BEAM3的单元
只支持常实数
建议改成BEAM188或者是BEAM189
这样可以支持截面
使用起来更好
而且计算精度更好
不是么？

而且
如果在开始的部分采用变量的形式
更好
更方便
别人也比较容易读懂
不是么？

xiaofanfu

我用四种方法算过，包括beam188，别外用vb,matlab编程，结果都一样。

jorgenhojmann

是的
但是应该说是结果都很接近
而不是一样吧

关键是改成变量
和截面后
别人用起来更方便
不是么？
而且
增加的代码也不是很多
既然都做了好人
就索性做到底吧