Ansys經(jīng)典有很多值得學(xué)習(xí)的案例，但是經(jīng)典有個問題就是，很多案例的單元，現(xiàn)在已經(jīng)沒有了。

比如在《Ansys機械工程應(yīng)用50例》一書中，第13.2，有預(yù)應(yīng)力的模態(tài)分析-弦的橫向振動。

這個題好，有理論解，結(jié)構(gòu)簡單，很適合做例題。

按照文章進行練習(xí)：

在Element Type, Add/Edit/Delete,打開一看，LINK1在哪里？

Link下只有180和11，link1單元已經(jīng)作古，隨風(fēng)消逝了。

一般來說有些單元雖然作廢了，但是在命令流里還是可以繼續(xù)存在的。

在命令行輸入

/clear

/filename,examplee9

/prep7

Et,1,link1

然后編輯一下實常數(shù)，ansys已經(jīng)不支持此單元的實常數(shù)了。

那么只能用Link180單元來替換link1單元了。

Link1單元是平面2D單元，link180單元是3D單元，用180單元替換能夠得到一樣的結(jié)果嗎？

選擇了link180單元，然后選擇實常數(shù)，程序警告：

原來LINK180不支持實常數(shù)。

其實很正常，很多原來的2D單元，如梁，pipe單元都有實常數(shù)，隨著版本更新?lián)Q成了高階單元后，實常數(shù)都沒有了，Beam，pipe，link的截面特性都在Sections里定義。

選擇link后，add，可以填寫link area的數(shù)據(jù)。

其中Added Mass是指桿件額外增加的重量。如果材料的密度（或者等效密度）輸入正確的話，此處可以不再輸入。

按照書中的方法，設(shè)置左端點Ux，Uy=0，右端點Uy=0

算完后查詢計算結(jié)果，結(jié)果如下：

為什么1階為0呢，23階的數(shù)值一樣呢？

合理的猜想是：因為沒有約束Z軸的位移。link1單元是2D單元，默認(rèn)是約束Z軸方向的。而link180屬于3D單元，Z軸方向是沒有約束的，這樣算出來會有y和z方向的自振頻率。

那么如果將Z軸方向約束住，是否就能夠和link1一樣的效果呢？

約束Z軸方向后，再看結(jié)果：

從結(jié)果看，與書中的結(jié)果完全一致。

最后看一下?lián)軇拥那傧?，是否能夠撩動你的心弦?/p>

總結(jié)：

ansys中很多單元已經(jīng)被改進的新單元，在GUI操作用無法使用，但是一般來說不影響在命令流中的使用。

當(dāng)需要在GUI中使用改進的新單元來替換老單元時應(yīng)當(dāng)注意：

1. 新單元和老單元的區(qū)別。新單元是包含更多的自由度，那么模擬老單元需要將多余自由度約束。新單元往往將好幾種老單元的特點合并，那么需要區(qū)分老單元對應(yīng)的參數(shù)應(yīng)該如何在新單元中體現(xiàn)。

2. 如BEAM，PIPE，LINK等單元的慣性矩，面積等實常數(shù)，現(xiàn)在經(jīng)常在Sections對應(yīng)的單元下輸入實際的尺寸。

3. 盡可能的使用新單元，如果對新單元不了解可以先做幾個簡單例子試驗一下。

!改進命令流如下

/clear

/filename,link180

/prep7

Et,1,link180

KEYOPT,1,12,1

Mp,ex,1,2e11

Mp,prxy,1,0.3

Mp,dens,1,7800

SECTYPE,1,LINK, ,1

SECDATA,1e-6,

K,1,,,,

K,2,1,,

L,1,2

Lesize,1,,,50

Lmesh,1

Finish

/solu

!設(shè)置Z方向位移為0

DL,1, ,UZ,

Dk,1,ux

Dk,1,uy

Dk,2,uz

Dk,2,uy

Dk,2,uz

Fk,2,fx,2000

Pstres,on

Solve

Save

Finish

/solu

Antype,modal

Modopt,lanb,10

Mxpand,10

Dk,2,ux

Dk,2,uy

Dk,2,uz

Pstres,on

Solve

Finish

/post1

Set,list

Set,first

Pldi

Anmode,10,0.5,,0

Set,next

Pldi

Anmode,10,0.5,,0

finish