Ansys Beam單元詳解——BEAM44(一)

2013-08-14  by:有限元優(yōu)化設(shè)計  來源:仿真在線

Ansys Beam單元詳解——BEAM44

 

BEAM44

3-D 彈性 變截面(Tapered)非對稱 梁

BEAM44 Element Description

BEAM44 是一種具有承受拉、壓、扭轉(zhuǎn)和彎曲能力的單軸梁。單元每個節(jié)點有6個自由度:x、y、z方向的平移和x、y、z軸向的轉(zhuǎn)動。這個單元允許具有不對稱的端面結(jié)構(gòu),并且允許端面節(jié)點偏離截面形心位置。如果你并不需要這些特性那么可以選用均質(zhì)對稱的BEAM4單元。同時也可以選用同類型2D單元(BEAM54)。對于非線性材料,使用BEAM188和BEAM189來代替BEAM44。

對于剪切變形的作用是可選的。例外一個可選是可顯示在單元坐標(biāo)系下,力作用單元上。同時包括應(yīng)力強化和大變形的特性??梢圆殚啞禔NSYS理論手冊》得到關(guān)于本單元更詳細的說明。

BEAM44可以使用SECTYPE, SECDATA, SECOFFSET, SECWRITE, and SECREAD命令來建立任何形狀的橫截面。但是只有當(dāng)沒有定義實常數(shù)時,上面的定義截面的命令才有效。

Ansys Beam單元詳解——BEAM44(一)ansys workbanch圖片1

 

BEAM44 輸入數(shù)據(jù)

單元的幾何模型、節(jié)點位置和坐標(biāo)系如圖44.1所示。單元由參考坐標(biāo)(x', y', z')和偏移量來定位。這個參考坐標(biāo)由I、J和K節(jié)點或方向角來定義,如圖44.1所示。梁的主軸為單元坐標(biāo)系(x, y, z)中沿X方向,并經(jīng)過橫截面中心(C.G.)。

單元的X軸的方向是指從I節(jié)點(end 1)到J節(jié)點(end 2)。如果只給了兩個節(jié)點參數(shù),默認(θ = 0°)那單元Y軸的方向自動確定為平行于系統(tǒng)坐標(biāo)系下的X-Y平面。有關(guān)示例見上圖。當(dāng)單元坐標(biāo)的X軸平行于整體坐標(biāo)系下的Z軸(包括0.01%的偏差在內(nèi)),單元Y軸的方向是平行于總體坐標(biāo)系下的Y軸。用戶可以通過給定θ角或定義第三個節(jié)點的方法來控制單元的方向。如果前面的兩個參數(shù)同時給定時,則以給定第三點的控制為準(zhǔn)。第三點一經(jīng)給出就意味著定義了一個由I,J,K三點定義的平面且該平面包含了單元坐標(biāo)的X與Z軸。當(dāng)本單元用于大變形分析時,那么給定的第三節(jié)點(K)或旋轉(zhuǎn)角(θ)僅用來確定單元的初始狀態(tài)。(有關(guān)方向節(jié)點及單元劃分的詳細信息參見《實體單元分網(wǎng)》及《ansys建模與分網(wǎng)指南》。K節(jié)點的只能生成可以查看LMESH和LATT命令描述。

用實常數(shù)來描述梁的橫截面,需要定義面積、面積慣性矩、端點到形心的距離,形心的偏移,和剪切系數(shù)。IZ和IY兩個慣性矩是在單元主軸的側(cè)面。截面 1(IX1)的扭轉(zhuǎn)慣性矩如果沒有特別說明,那其值就默認為在截面1的極慣性矩(IZ1+IY1)。截面2 的轉(zhuǎn)動慣性矩(IX2,IY2和IZ2),如果空白,就默認為對應(yīng)截面1的值。單元的扭轉(zhuǎn)剛度隨著IX的值而減小。

偏移量常數(shù)(DX, DY, DZ)由截面中心位置相對節(jié)點位置來定義。沿著單元坐標(biāo)系正向的位置為正。所有截面2 處實常數(shù)(除了中心偏移實常數(shù)DX, DY, 和 DZ)如果為0,則都默認為其對應(yīng)截面 1 處的值。截面1節(jié)點的上層厚度(the "bottom" thicknesses),TKZT1 和 TKYT1,分別默認為截面1節(jié)點的下層厚度(the "bottom" thicknesses),TKZB1 and TKYB1。截面 2節(jié)點的上層厚度(the "bottom" thicknesses),TKZT2 和 TKYT2,也分別默認為截面 2節(jié)點的下層厚度(the "bottom" thicknesses),TKZB2 and TKYB2。上層厚度為頂點到截面中心的距離,下層厚度為底點到截面中心的距離。

剪切變形系數(shù)(SHEARZ和SHEARY)只有在考慮剪切變形時才使用。某個方向上剪切系數(shù)為0,一般用在忽略該方向上的剪切變形。查看剪切變形細節(jié)。

如果沒有實常數(shù)定義,使用SECTYPESECDATA命令來分別定義橫截面(細節(jié)查看《ANSYS結(jié)構(gòu)分析指南》中的梁分析和橫截面)。注意使用SECTYPESECDATA命令建立的截面可以在同一個模型中被BEAM44, BEAM188, 和 BEAM189中的任何組合使用。截面與單元用截面號(SECNUM)來關(guān)聯(lián),截面號是獨立的單元屬性。

KEYOPT(2)允許縮減質(zhì)量矩陣(reduced mass matrix)表達(刪除轉(zhuǎn)動方向的自由度)。這個選項有助于改善在分析細長桿在質(zhì)量荷載(mass loading)下計算得到彎曲應(yīng)力結(jié)果。

KEYOPT(7) 和 KEYOPT(8)允許在單元坐標(biāo)系上節(jié)點上的單元節(jié)點約束釋放。在做自由運動時,節(jié)點約束則不能釋放,同時將會有pivot警告和錯誤信息。同時,應(yīng)力剛度矩陣的平移自由度上不能進行節(jié)點約束釋放。作用在節(jié)點約束釋放方向上的荷載將被忽略。對于大變形,注意到單元的節(jié)點約束釋放是沿著單元的方向,但是在節(jié)點耦合處不釋放。在沒有節(jié)點約束釋放的模型中添加柔性(低彈性模量)梁單元有助于提高解的穩(wěn)定性。

剪切面積(ARES_ _)和扭轉(zhuǎn)應(yīng)力系數(shù)(TSF_) 在非零時,將起作用。剪切面積僅僅是為了計算剪切應(yīng)力,一般小于實際橫截面的面積。扭慣性矩乘上扭轉(zhuǎn)應(yīng)力系數(shù)可計算得到扭轉(zhuǎn)剪切應(yīng)力。扭轉(zhuǎn)剪切系數(shù)一般可在結(jié)構(gòu)手冊上查到。對于圓截面,TSF=直徑/(2*IX)。

對于有些梁的截面,剪切中心可以與形心不重合。非零的剪切中間偏移(DSC_ _)應(yīng)該按圖44.1“梁幾何模型”所示的輸入。剪切中心偏移以形心沿單元坐標(biāo)系正向為正。截面 2節(jié)點的偏移值如果為0,則默認為截面1節(jié)點的數(shù)值。如圖44.2“梁的應(yīng)力輸出”所示,如果從Y1到Z4的常數(shù)給定,梁每端可以輸出給定點4個點的附加應(yīng)力。

彈性基礎(chǔ)剛度(EFS_)是基礎(chǔ)產(chǎn)生單位法向變形所需的壓力。當(dāng)EFS_為0時,可以忽略這個性能。單元的初始應(yīng)變(ISTRN)通過Δ/L給定,這里的Δ是單元長度L(由節(jié)點I和J的坐標(biāo)所決定)與零應(yīng)變時的長度之差。參數(shù)ADDMAS要輸入的值是每單位長度的附加質(zhì)量。

“節(jié)點與單元荷載”一節(jié)對“單元荷載”有專門介紹??梢栽诒締卧谋砻媸┘用婧奢d,如上圖中帶圈數(shù)字所示,其中箭頭指向為面荷載作用正向。橫向均布壓力的單位為力每單位長度,端點作用的壓力應(yīng)以集中力的形式輸入。KEYOPT(10)用來控制線性變化的橫向壓力相對單元節(jié)點的偏移量??稍趩卧獛缀螆D形的八個角上設(shè)定溫度值,其被當(dāng)作體荷載處理。第一個角上的溫度T1的默認值為TUNIF,如其它角的溫度未給定時其默認值等于第一個角的溫度,如給定了T1和T2則T3的默認值為T2,T4的默認值為T1;T5到T8的值默認與T1到T4的值相對應(yīng)。對于其他輸入方法,沒有定義溫度,則都默認TUNIF。

KEYOPT(9)用來控制兩節(jié)點中間部分相關(guān)值的輸出情況,值是按平衡條件得出的。但在下列情況下這些值不能得到:

l 考慮應(yīng)力硬化時[SSTIF,ON];

l 一個以上的部件作用有角速度時[OMEGA];

l 通過命令CGOMGA, DOMEGA, or DCG/MG作用了角速度或加速度時。

概要的單元輸入在“BEAM44 輸入概要”中有介紹。單元輸入的一般描述在“單元輸入”中有介紹。

BEAM44 輸入概要

節(jié)點 I, J, K (K 方向節(jié)點可選)

自由度 UX, UY, UZ, ROTX, ROTY, ROTZ

實常數(shù)

AREA1, IZ1, IY1, TKZB1, TKYB1, IX1,

AREA2, IZ2, IY2, TKZB2, TKYB2, IX2,

DX1, DY1, DZ1, DX2, DY2, DZ2,

SHEARZ, SHEARY, TKZT1, TKYT1, TKZT2, TKYT2,

ARESZ1, ARESY1, ARESZ2, ARESY2, TSF1, TSF2,

DSCZ1, DSCY1, DSCZ2, DSCY2, EFSZ, EFSY,

Y1, Z1, Y2, Z2, Y3, Z3,

Y4, Z4, Y1, Z1, Y2, Z2,

Y3, Z3, Y4, Z4, THETA, ISTRN,

ADDMAS

查閱 Table 44.1: "BEAM44 Real Constants" 獲得實常數(shù)的描述

材料屬性 EX, ALPX (or CTEX or THSX), DENS, GXY, DAMP

面荷載

壓力 --

face 1 (I-J) (-Z normal direction)

face 2 (I-J) (-Y normal direction)

face 3 (I-J) (+X tangential direction)

face 4 (I) (+X axial direction)

face 5 (J) (-X axial direction)

(負值表示反方向面荷載)

體力

溫度 --

T1, T2, T3, T4, T5, T6, T7, T8

特性

應(yīng)力強化

大變形

單元生死

KEYOPT(2) 集總質(zhì)量矩陣:

0 -- 協(xié)調(diào)

1 -- 縮減

KEYOPT(6) 輸出力和力矩:

0 -- 沒有力的輸出

1 -- 輸出在單元坐標(biāo)系下的力和力矩

KEYOPT(7) I節(jié)點約束釋放:

1 -- 單元Z軸約束釋放

10 -- 單元Y軸約束釋放

100 -- 單元X軸約束釋放

1000 -- 單元Z方向約束釋放

10000 -- 單元Y方向約束釋放

100000 -- 單元X方向約束釋放

為了組合這些約束釋放, 可以輸入這些數(shù)字的和(例如輸入11表示Z軸和Y軸轉(zhuǎn)動).

KEYOPT(8) 與 KEYOPT(7)相同 只表示在J節(jié)點:

KEYOPT(9) 控制單元節(jié)點間計算點的輸出情況

N -- 輸出計算點數(shù)(N=0,1,3,5,7,9)

KEYOPT(10) 用于控制用SFBEAM指令輸入荷載偏移量的輸入方法:

0 -- 以長度為單位,直接輸入壓力相對于I,J節(jié)點的偏移量

1 -- 以偏移量與單元長度的比值為單位(0~1),用相對值形式輸入偏移

注意:

SHEARZ 與 IZ相關(guān). 如果 SHEARZ = 0.0, 表示在Y方向上沒有剪切變形.

SHEARY與 IY相關(guān). 如果 SHEARY = 0.0, 表示在Z方向上沒有剪切變形

表 44.1 BEAM44 實常數(shù)

序號

名稱

描述

1

AREA1

截面1的面積

2, 3

IZ1, IY1

截面1慣性矩

4, 5

TKZB1, TKYB1

截面1的Z方向和Y方向的下層厚度

6

IX1

截面1扭轉(zhuǎn)慣性矩

7

AREA2

截面2的面積

8, 9

IZ2, IY2

截面2慣性矩

10, 11

TKZB2, TKYB2

截面2的Z方向和Y方向的下層厚度

12

IX2

截面1扭轉(zhuǎn)慣性矩

13, 14, 15

DX1, DY1, DZ1

截面1的X, Y, and Z方向的形心偏移量

16, 17, 18

DX2, DY2, DZ2

截面2的X, Y, and Z方向的形心偏移量

19, 20

SHEARZ, SHEARY

Z方向和Y方向上的剪切變形系數(shù)

21,22

TKZT1, TKYT1

截面1的Z方向和Y方向的上層厚度

23, 24

TKZT2, TKYT2

截面2的Z方向和Y方向的上層厚度

25, 26

ARESZ1, ARESY1

截面1上Z方向和Y方向的剪切面積

27, 28

ARESZ2, ARESY2

截面2上Z方向和Y方向的剪切面積

29. 30

TSF1, TSF2

每個截面的扭轉(zhuǎn)應(yīng)力系數(shù)

31, 32

DSCZ1, DSCY1

截面1的剪切中心偏移量

33, 34

DSCZ2, DSCY2

截面2的剪切中心偏移量

35, 36

EFSZ, EFSY

基礎(chǔ)剛度

37, 38

Y1, Z1

截面1附加應(yīng)力輸出坐標(biāo)值系列1

39, 40

Y2, Z2

截面1附加應(yīng)力輸出坐標(biāo)值系列2

41, 42

Y3, Z3

截面1附加應(yīng)力輸出坐標(biāo)值系列3

43, 44

Y4, Z4

截面1附加應(yīng)力輸出坐標(biāo)值系列4

45, 46

Y1, Z1

截面2附加應(yīng)力輸出坐標(biāo)值系列1

47, 48

Y2, Z2

截面2附加應(yīng)力輸出坐標(biāo)值系列2

49, 50

Y3, Z3

截面2附加應(yīng)力輸出坐標(biāo)值系列3

51, 52

Y4, Z4 (at end J)

截面2附加應(yīng)力輸出坐標(biāo)值系列4

53

THETA

X軸角度

54

ISTRN

初始應(yīng)變

55

ADDMAS

附加質(zhì)量/單位長度

BEAM44 輸出數(shù)據(jù)

單元的求解結(jié)果總體為以下兩部分:

l 全部節(jié)點解包括節(jié)點位移解;

l單元解,細則見“BEAM4單元輸出數(shù)據(jù)一覽表”

圖44.2“BEAM44應(yīng)力輸出”已經(jīng)圖示了幾種情況。在每個截面計算結(jié)果包括軸向應(yīng)力和4個彎曲應(yīng)力。這是5個應(yīng)力值組合起來以矩形截面形式得到最大和最小應(yīng)力。如果從Y1到Z4的實數(shù)給定,在任何點的組合應(yīng)力也可以計算得到,如果圖44.2“BEAM44應(yīng)力輸出”所示。對于KEYOPT(6) = 1的單元,12個彎距和軸力(每邊6個)也可以輸出(在單元坐標(biāo)系上)。單元的X軸通過截面的重心。如果實常數(shù)25~30(ARES_ _, TSF_)給定,平均剪切應(yīng)力和扭轉(zhuǎn)應(yīng)力也可以輸出。如果他們都為0,剪切應(yīng)力和扭轉(zhuǎn)應(yīng)力的輸入將無效。如果設(shè)置了KEYOPT(9),還可以輸出2個端點之間的點的值。有關(guān)求解輸出的總體性介紹請見“Solution Output”,而有關(guān)結(jié)果內(nèi)容察看的一些方法見“ANSYS Basic Analysis Guide”

注意:

如果輸入/ESHAPE,1 ,3D顯示的BEAM44單元只僅限在ANSYS前處理器中,后處理器將無法顯示3D效果。

Ansys Beam單元詳解——BEAM44(一)ansys workbanch圖片2

單位輸出定義表使用下面符號:

在name 列的冒號表示該項目可以通過構(gòu)成名字的方法來獲得[ETABLE, ESOL]。第0 列表示該項有效的說明在文件Jobname.OUT 中。R 列表示該項的結(jié)果顯示在results 文件中。無論在0 還是R 列中,Y 表示該項一直是可用的。數(shù)值表示描述哪里該項是選擇性提供的腳注,-表示該項不提供

表 44.2 BEAM44 單元輸出定義

名稱

定義

O

R

EL

單元號

Y

Y

NODES

節(jié)點 - I, J

Y

Y

MAT

材料號

Y

Y

VOLU:

體積

-

Y

XC, YC, ZC

輸出結(jié)果坐標(biāo)

Y

5

TEMP

溫度 T1, T2, T3, T4, T5, T6, T7, T8

Y

Y

PRES

壓力P1 at nodes I, J; 偏移量1 at I, J; 壓力P2 at I, J; 偏移量2 at I, J; 壓力P3 at I, J; 偏移量3 at I, J; 壓力P4 at I; 壓力P5 at J

Y

Y

SDIR

軸向正應(yīng)力

1

1

SBYT

梁內(nèi)單元+Y面上的彎曲應(yīng)力

1

1

SBYB

梁內(nèi)單元-Y面上的彎曲應(yīng)力

1

1

SBZT

梁內(nèi)單元+Z面上的彎曲應(yīng)力

1

1

SBZB

梁內(nèi)單元-Z面上的彎曲應(yīng)力

1

1

SMAX

最大應(yīng)力(正應(yīng)力+彎曲應(yīng)力)

1

1

SMIN

最小應(yīng)力(正應(yīng)力-彎曲應(yīng)力)

1

1

EPELDIR

端部軸向彈性應(yīng)變

1

1

EPELBYT

梁內(nèi)單元+Y面上的彎曲彈性應(yīng)變

1

1

EPELBYB

梁內(nèi)單元-Y面上的彎曲彈性應(yīng)變

1

1

EPELBZT

梁內(nèi)單元+Z面上的彎曲彈性應(yīng)變

1

1

EPELBZB

梁內(nèi)單元-Z面上的彎曲彈性應(yīng)變

1

1

EPTHDIR

端點軸向熱應(yīng)變

1

1

EPTHBYT

梁內(nèi)單元+Y面上的彎曲熱應(yīng)變

1

1

EPTHBYB

梁內(nèi)單元-Y面上的彎曲熱應(yīng)變

1

1

EPTHBZT

梁內(nèi)單元+Z面上的彎曲熱應(yīng)變

1

1

EPTHBZB

梁內(nèi)單元-Z面上的彎曲熱應(yīng)變

1

1

EPINAXL

單元初始軸向應(yīng)變

1

1

S(XY, XZ, YZ)

平均剪切Y方向,平均剪切Z方向,扭轉(zhuǎn)應(yīng)力

2

2

S(AXL1, AXL2, AXL3, AXL4)

定義點的組合應(yīng)力

3

3

MFOR(X, Y, Z)

元坐標(biāo)系中X,Y,Z方向的力

4

Y

MMOM(X, Y, Z)

單元坐標(biāo)系中X,Y,Z方向的彎矩

4

Y

  1. 這些項目在I截面重復(fù), 中間計算點(查看 KEYOPT(9)), 和J截面
  2. 只在實常數(shù)25-30給定時輸出
  3. 只在實常數(shù)37-52給定時輸出
  4. 如果KEYOPT(6) = 1
  5. 在中心點使用 *GET 命令.

開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:Ansys Beam單元詳解——BEAM44(一) Ansys靜力分析 有限元熱分析 Ansys疲勞分析 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復(fù)電話,謝謝!
驗證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com