1 問(wèn)題描述

如下圖所示鋼制圓盤(pán)轉(zhuǎn)子結(jié)構(gòu)，材料的彈性模量為200GPa，密度為7800kg/m3，泊松系數(shù)為0.3，長(zhǎng)度分別為L(zhǎng)1=0.2m，L2=0.3m，L3=0.5m，L4=0.3m，軸徑d=0.05m，圓盤(pán)半徑分別為R1=0.12m，R2=0.2m，R3=0.2m，厚度分別為t1=0.05m，t2=0.05m，t3=0.06m。軸的兩端為軸承支撐，軸承的動(dòng)力特性隨著轉(zhuǎn)速變化見(jiàn)下表。求該轉(zhuǎn)子結(jié)構(gòu)渦動(dòng)頻率、振型、臨界轉(zhuǎn)速。

假定軸承剛度為轉(zhuǎn)速（rad/s）的非線性函數(shù)，但阻尼不隨著轉(zhuǎn)速而改變（變阻尼的方法與變剛度一樣）。

2 結(jié)果對(duì)比

采用BEAM188單元模擬轉(zhuǎn)子結(jié)構(gòu)，COMBI214模擬軸承。分別用ANSYS APDL和ANSYS Workbench模擬該例子作對(duì)比，可以看出APDL和Workbench的結(jié)果是一致的。

APDL結(jié)果

Workbench結(jié)果

3 APDL模擬過(guò)程

當(dāng)軸承的動(dòng)力特性隨轉(zhuǎn)速變化時(shí)，必須采用COMBI214單元來(lái)模擬軸承，COMBI214單元的實(shí)常數(shù)采用表數(shù)組定義，注意定義表數(shù)組的變量名項(xiàng)（VAR1）需為“OMEGS”，以表明數(shù)組隨著轉(zhuǎn)速而改變。具體命令流如下：

*dim,k11,table,7,1,1,omegs

k11(1,0)=0,200,400,600,800,1000,1200 !rad/s

k11(1,1)=1.6e6,2.0e6,4.8e6,8.8e6,1.3e7,1.8e7,2.3e7 !N/m

*dim,k22,table,7,1,1,omegs

k22(1,0)=0,200,400,600,800,1000,1200 !rad/s

k22(1,1)=2e6,3e6,6e6,1.1e7,1.7e7,2.3e7,3e7 !N/m

*dim,k12,table,7,1,1,omegs

k12(1,0)=0,200,400,600,800,1000,1200 !rad/s

k12(1,1)=-2e4,-2.4e4,-6.8e4,-1.2e5,-1.9e5,-2.6e5,-3.6e5 !N/m

*dim,k21,table,7,1,1,omegs

k21(1,0)=0,200,400,600,800,1000,1200 !rad/s

k21(1,1)=0.6e4,4e4,1e5,1.7e5,2.5e5,3.4e5,4.3e5 !N/m

r,2,%k11%,%k22%,%k12%,%k21%,2e3,5e3

rmore,-4e3,3e3

4 Workbench模擬過(guò)程

---

轉(zhuǎn)子動(dòng)力學(xué)系列（十）：不平衡激勵(lì)下的啟動(dòng)過(guò)程瞬態(tài)轉(zhuǎn)子動(dòng)力學(xué)分析

轉(zhuǎn)子動(dòng)力學(xué)系列（九）：基于ANSYS Workbench的多軸轉(zhuǎn)子臨界轉(zhuǎn)速

轉(zhuǎn)子動(dòng)力學(xué)系列（八）：軸對(duì)稱實(shí)體單元Solid272/Solid273的應(yīng)用

轉(zhuǎn)子動(dòng)力學(xué)系列（七）：帶支承結(jié)構(gòu)的復(fù)雜轉(zhuǎn)子分析

轉(zhuǎn)子動(dòng)力學(xué)系列（六）：考慮預(yù)應(yīng)力的轉(zhuǎn)子動(dòng)力學(xué)分析

轉(zhuǎn)子動(dòng)力學(xué)系列（五）：隨轉(zhuǎn)速變剛度和變阻尼的模擬

轉(zhuǎn)子動(dòng)力學(xué)系列（四）：不同軸承單元對(duì)比（COMBIN14和COMBI214）

轉(zhuǎn)子動(dòng)力學(xué)系列（三）：不同建模單元對(duì)比（BEAM188與SOLID186）

轉(zhuǎn)子動(dòng)力學(xué)系列（二）：不平衡響應(yīng)分析

轉(zhuǎn)子動(dòng)力學(xué)系列（一）：臨界轉(zhuǎn)速與坎貝爾圖