下面的所有東西,希望樓主拷貝到記事本裏面查看,不要選擇自動換行!?
/PREP7!?進入前處理器
ET,1,BEAM4!?定義單元beam4,是個三維梁單元
KEYOPT,1,2,0?!?定義單元的關鍵選項,如果後面是0,代表默認的,可以先不用理解,
KEYOPT,1,6,0
KEYOPT,1,7,1?!?這裏定義第七個關鍵選項,定義為編號1,也就是計算陀螺阻尼矩陣方程,要求
!?IYY?等於?IZZ,也就是兩個轉動慣量要相等,這兩個量要在實常數中定義 !?也就是下面的命令RKEYOPT,1,9,0
KEYOPT,1,10,0?
*SET,p,acos(-1)?!?定義三個參數,分別是派(4.1315)、第壹個半徑R1,第二個半徑r2
*SET,R1,5
*SET,R2,60!?半徑的單位
!?定義單元的實常數,有兩個,因為有兩個半徑,分別就是下面的R,1 !?和r,2 !?r命令的定義中需要根據使用單元beam4來壹壹對應,不同的單元R命令 !?定義的意思是不壹樣的,具體每個意思,下圖看,壹定要壹壹對應 !?這裏就先是截面積,Z向轉動慣量,Y像轉動慣量,是壹樣的,從上面的 !?關鍵選項定義中可以看的出來,KEYOPT,1,7,1R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1,?,?
RMORE,?,p*R1**4/2,?,?,2175,?,?!?這個也是定義實常數,因為命令只能定義6個數,從第七個
!?就要使用這個命令,編號可以從下圖中對應下R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2,?,!?定義第二個實常數R,2
RMORE,?,p*R2**4/2,?,?,2175,?,?
MPTEMP,,,,,,,,!?這裏來定義材料屬性,目前樓主提供的是GUI操作以後的
!?,這個需要相應簡化,因為GUi操作中壹般包括了溫度的考慮,所以使用了幾個溫度的命令 !?目前看來只是有壹個溫度,所以不用考慮溫度,簡化後的,在後面的註釋中,註釋開始是!(感嘆號)MPTEMP,1,0?
MPDATA,EX,1,,2e5
MPDATA,PRXY,1,,.3?
MPTEMP,,,,,,,,?
MPTEMP,1,0?
MPDATA,DENS,1,,1e-10
!?簡化後的材料屬性
!?mp,ex,1,2e5?!?分別定義彈性模量,泊松比,密度,其中的1代表第壹個材料屬性
!?mp,PRXY,1,0.3
!?mp,DENS,1,1e-10
MPTEMP,,,,,,,,!?這裏就是定義第二個材料屬性了!同樣包括彈性模量,泊松比,密度
MPTEMP,1,0?
MPDATA,EX,2,,2E5
MPDATA,PRXY,2,,.3?
MPTEMP,,,,,,,,?
MPTEMP,1,0?
MPDATA,DENS,2,,8E-8?
K,?,,,,?!?定義兩個關鍵點,第壹個是0,第二個坐標是:(100,0,0)
K,?,100,,,?
TYPE,?1?!?選擇第壹個單元,也就是beam4
MAT,?1!?選擇第壹個材料屬性
REAL,?1?!?選擇第壹個實常數
ESYS,?0?!?選擇默認的單元坐標系統
LSTR,?1,?2!?用剛才生成的兩個關鍵點,建立壹條直線
LESIZE,ALL,?,?,200,?,1,?,?,1,!?把建立的這條直線劃分為200個
LMESH,?1?!?根據劃分的線的段數,網格這條線,從而生成有限元模型
D,1,UX!?對於第壹個節點(也就是原點位置的那個節點)進行約束,三個平移方向都約束,但是轉動方向不約束,
!?因為單元beam4是有六個自由度的D,1,UY
D,1,UZ
D,102,UY?!?對於第102個節點,這樣來算他的位置,總***長度是100,劃分為200分,那麽102是那個位置,樓主應該知道了吧
!?選擇了Y方向和Z方向的平移約束D,102,UZ?
FLST,2,1,2,ORDE,1!?後面的兩個命令結合起來就是選擇了第200個節點,但是沒有具體操作,
!?只是有了拾取選擇這個點的操作,這兩個命令壹般是GUI直接生成的,寫命令流的壹般不使用這兩個操作 !?所以這兩個命令其實什麽任務都沒有做目前FITEM,2,200
對上面的命令壹壹般是需要簡化的,簡化後的就是:
/PREP7?
ET,1,BEAM4?
KEYOPT,1,7,1
*SET,p,acos(-1)
*SET,R1,5
*SET,R2,60?
R,1,p*R1**2,p*R1**4/4,p*R1**4/4,2*R1,2*R1,?,?
RMORE,?,p*R1**4/2,?,?,2175,?,?
R,2,p*R2**2,p*R2**4/4,p*R2**4/4,2*R2,2*R2,?,?
RMORE,?,p*R2**4/2,?,?,2175,?,?
mp,ex,1,2e5?
mp,PRXY,1,0.3
mp,DENS,1,1e-10
mp,ex,2,2e5?
mp,PRXY,2,0.3
mp,DENS,2,8e-8
K,1,,,,?
K,2,100,,,?
MAT,?1
REAL,?1?
LSTR,?1,?2?
LESIZE,ALL,?,?,200
LMESH,?1
D,1,UX
D,1,UY
D,1,UZ
D,102,UY
D,102,UZ
FLST,2,1,2,ORDE,1?
FITEM,2,200
簡化後的命令樓主結合本人說的,可以理解下
這段話因為只是在建模,所以沒有執行什麽計算的,思路其實我已經說的很明確了:
首先定義單元,指定壹個關鍵選項包括了陀螺阻尼效果,然後再定義兩個實常數,然後定義了兩個材料屬性,然後建立壹個壹個線,使用第壹個實常數,和第壹個材料屬性對這個線進行了網格!其次再對這個線進行了約束的施加,原點位置三個平移方向全約束,對102個節點位置,進行了兩個方向的約束