【讨论】初学乍练，恳请各位高手给会诊一下问题所在

Insomnia

刚刚开始学习ABAQUS,试着写一个INP文件，总是出错，恳请各位高手给会诊一下问题所在。
问题描述：
一个梁，长40in，两端弹簧支撑，k=10e4 lb/in，离左端5in处有一集中质量，mass=0.02lb-s^2/in。梁截面宽度为2in*2in，I=1.33in^4, 材料E=30*E6, POSSION=0.29, DENSITY=.0007 lb-s^2/in^4
  
预求该体系的前6阶振动频率。错误的INP文件如下：
  
*HEADING
  Determine the normal modes and natural frequencies of a spirng-beam-mass system
  Units - lb, -s, in
*NODE
  1,  0.00, 0.00
  9, 40.00, 0.00
10, 0.00, -1.00
11, 40.00, -1.00
*NGEN
  1, 9
*ELEMENT, TYPE = B21
1, 1, 2
*ELGEN, ELSET=BEAM
  1,8
*BEAM SECTION, SECTION=RECTANGULAR, ELSET=BEAM, MATERIAL=USER_DEFINED
  2., 2.
*ELEMENT, TYPE = MASS, ELSET=CMASS
9, 2
*MASS, ELSET=CMASS
0.02,
* ELEMENT, TYPE=SPRINGA, ELSET=SPRINGS
10, 1, 10
11, 9, 11
*SPRING, ELSET = SPRINGS
,10000
*MATERIAL, NAME=USER_DEFINED
*ELASTIC
  30.E6, 0.29
*DENSITY
  0.0007
* STEP, name=Frequency, PERTURBATION
*BOUNDARY
  10,ENCASTRE
  11,ENCASTRE
*FREQUENCY
  6,
*END STEP

ghostfp

我有一个桥梁振动模态的分析例子，和你的问题很像
改天给你传上来

Insomnia

不会吧？:O) 以斑竹的法力，应该可以直接看出问题所在吧。

seansheu

恕刪...
* ELEMENT, TYPE=SPRINGA, ELSET=SPRINGS
10, 1, 10
11, 9, 11

keyword應該是*ELEMENT，不能有空格
  

*SPRING, ELSET = SPRINGS
,10000

這裡有誤，elset的名字前面不能有空格
另外，springA單元定義第一行應為空白，第二行第一個值才是Spring stiffness
  

...
* STEP, name=Frequency, PERTURBATION

同理，*STEP不能有空格
  

*BOUNDARY
  10,ENCASTRE
  11,ENCASTRE
*FREQUENCY
  6,
*END STEP

*Boundary如要在STEP裡宣告，要擺到*Frequency後面
或者，擺到Step外的Model Data裡
  
經修改後計算得前六項Natural frequency為:
0.0000,  9.67450E-05,  55.187,  101.08,  271.29,  678.22     
  
完整input file
==========================================================================
  
*HEADING  
Determine the normal modes and natural frequencies of a spirng-beam-mass system  
Units - lb, -s, in  
*NODE  
1, 0.00, 0.00  
9, 40.00, 0.00  
10, 0.00, -1.00  
11, 40.00, -1.00  
*NGEN  
1, 9  
*ELEMENT, TYPE = B21  
1, 1, 2  
*ELGEN, ELSET=BEAM  
1,8  
*BEAM SECTION, SECTION=RECTANGULAR, ELSET=BEAM, MATERIAL=USER_DEFINED  
2., 2.  
*ELEMENT, TYPE = MASS, ELSET=CMASS  
9, 2  
*MASS, ELSET=CMASS  
0.02,  
*ELEMENT, TYPE =SPRINGA, ELSET=SPRINGS  
10, 1, 10  
11, 9, 11  
*SPRING, ELSET=SPRINGS  
  
10000  
*MATERIAL, NAME=USER_DEFINED  
*ELASTIC  
30.E6, 0.29  
*DENSITY  
0.0007  
*STEP, name=Frequency, PERTURBATION  
*FREQUENCY  
6,  
*BOUNDARY  
10,ENCASTRE  
11,ENCASTRE  
*END STEP  
  
===============================================================================

Insomnia

非常感谢seansheu 大侠的指点。受益非浅。BOW//
  
现在的计算结果中有两个接近零的频率，说明系统存在刚体位移。所以在你的文件的基础上，在*BOUNDARY中在添加以下命令，
  1,1
  9,1
  2,3,5
  
则频率计算值为： 100.64， 268.77， 458.81， 779.26， 1276.7，1879.4
事实上，本例取自：
Cifuentes A.Q., Using MSC/NASTRAN: Statics and Dynamics, Springer-Verlag, 1989
该书中用NASTRAN计算的结果为：
100.80，269.64，461.46，790.95，1313.17，1965.98

seansheu

你的問題順利解決，可喜可賀...
感謝你提出對照數據，這樣看來，ABAQUS還是值得信賴的，呵...

ghostfp

你太抬举我了！
呵呵，我可没有那么大本事！
另，最近忙于实验，版上的问题我很难抽出时间解答！
多多仰仗各位大侠了！
  

Insomnia wrote:
不会吧？:O) 以斑竹的法力，应该可以直接看出问题所在吧。

ghostfp

seansheu大侠才是真正的高手！
以后要多多向你请教了
：）
  

seansheu wrote:
[quote]
恕刪...
* ELEMENT, TYPE=SPRINGA, ELSET=SPRINGS
10, 1, 10
11, 9, 11

keyword應該是*ELEMENT，不能有空格
  
[[/quote]

willams_lee

seansheu 确实是高手，以后有空就到这里逛，学东西。