[轉(zhuǎn)載]基于ANSYS經(jīng)典界面的光盤的轉(zhuǎn)子動(dòng)力學(xué)分析

 

我們經(jīng)常通過(guò)光驅(qū)來(lái)讀取光盤數(shù)據(jù)。光盤是一個(gè)很薄的盤,當(dāng)讀取數(shù)據(jù)時(shí),它的轉(zhuǎn)速很高。本篇博文試圖對(duì)光盤的轉(zhuǎn)子動(dòng)力學(xué)問(wèn)題進(jìn)行分析。

該光盤的結(jié)構(gòu)如下,它的內(nèi)圈被固定,而外圈自由。其轉(zhuǎn)速是120HZ.

已知:

幾何參數(shù): 厚度:0.8mm;內(nèi)徑:16.5mm,外徑:47.5mm

物理參數(shù):彈性模量 7.e10MPa;泊松比:0.3;密度:2800kg/m3.

【問(wèn)題分析】

1. 用SHELL181單元建模。這是一個(gè)殼單元。并通過(guò)實(shí)常數(shù)指明厚度。

2. 分析時(shí)取3個(gè)轉(zhuǎn)速:0.1, 該光盤轉(zhuǎn)速的一半,以及光盤的實(shí)際轉(zhuǎn)速來(lái)進(jìn)行三次模態(tài)分析。從而可以繪制拷貝圖。

3. 本篇博文用命令流方式說(shuō)明步驟。讀者只需要把相應(yīng)的命令流拷貝到ANSYS經(jīng)典界面的命令窗口,執(zhí)行即可。

【求解過(guò)程】

首先打開ANSYS的經(jīng)典界面

一 建模

1.定義參數(shù)

在命令窗口中輸入下列命令

pi = acos(-1)xa = 47.5e-3xb = 16.5e-3zh = 0.8e-3spin = 120*2*pi第1行定義了圓周率的大小第2行是光盤的外徑第3行是光盤的內(nèi)徑第4行是在計(jì)算光盤的轉(zhuǎn)速,用弧度每秒為單位。

2.定義單元類型和實(shí)常數(shù)

繼續(xù)輸入

/prep7et,1,181r,1,zh第1行指明進(jìn)入前處理器第2行定義光盤的單元類型是SHELL181第3行說(shuō)明該單元的厚度是0.8e-3

3.定義材料模型

繼續(xù)輸入

mp,ex,,7.2e+10mp,nuxy,,.3mp,dens,,2800.這三行分別說(shuō)明新建材料的彈性模量,泊松比和密度。

4.創(chuàng)建幾何模型

繼續(xù)輸入

cyl4,,,xb,0,xa,360這里繪制出一個(gè)圓環(huán),這就是光盤的幾何模型。此時(shí)主窗口中顯示如下圖

繼續(xù)輸入

esize,0.0025amesh,all第1行指明了單元的尺寸是2.5mm第2行劃分面網(wǎng)格。此時(shí)主窗口中顯示如下圖 可見,網(wǎng)格并不怎么樣。雖然有更好的處理方式,但是本文的主要目的是說(shuō)明轉(zhuǎn)子動(dòng)力學(xué)的計(jì)算步驟,就忽略這些細(xì)節(jié)問(wèn)題了。

6.設(shè)定邊界條件

繼續(xù)輸入

lsel,,,,5,8dl,all,1,allallselfini第1行說(shuō)明要選擇內(nèi)圈的四根圓弧第2行說(shuō)明對(duì)內(nèi)圓弧上所有節(jié)點(diǎn)全固定第3行說(shuō)明選擇所有對(duì)象第4行退出前處理器此時(shí)主窗口中顯示如下圖 可見,內(nèi)圈已經(jīng)被固定,而外圈自由。這正確的模擬了光盤裝在光驅(qū)內(nèi)部的情況。

二 仿真

1.設(shè)置分析類型及選項(xiàng)

繼續(xù)輸入

/soluantype,modalmodopt,qrdamp,30,,,onmxpand,30qrdopt,on第1行說(shuō)明進(jìn)入到求解器第2行說(shuō)明進(jìn)行模態(tài)分析第3行設(shè)置模態(tài)分析的選項(xiàng):使用QRDAMP算法,提取前30階模態(tài)。第4行指明要展開這30階模態(tài)第5行說(shuō)明,后一個(gè)載荷步要重用前一個(gè)載荷步的Block Lanczos的特征值解,以加快計(jì)算速度。2.開啟科里奧利效果繼續(xù)輸入coriolis,on,,,on說(shuō)明要打開科里奧利效果,而且是在固定坐標(biāo)系中求解。

3.第1次模態(tài)分析

繼續(xù)輸入

omega,,,0.1solve第1行指明該光盤繞軸轉(zhuǎn)動(dòng)的角速度是0.1弧度每秒第2行是求解的命令

4.第2次模態(tài)分析

繼續(xù)輸入

omega,,,spin/2solve第1行指明該光盤繞軸轉(zhuǎn)動(dòng)的角速度是最終轉(zhuǎn)速的二分之一第2行是求解的命令

5.第3次模態(tài)分析

繼續(xù)輸入

omega,,,spinsolvefinish第1行指明該光盤繞軸轉(zhuǎn)動(dòng)的速度是題目中給定的速度第2行是求解的命令第3行退出求解器

三 后處理

1. 輸出坎布爾圖

繼續(xù)輸入

/post1/yrange,500,1500plcamp第1行說(shuō)明進(jìn)入到通用后處理器第2行指明下面要畫的圖形的縱坐標(biāo)的范圍是500-1500.第3行指明請(qǐng)繪制坎貝爾圖結(jié)果如下

繼續(xù)輸入

prcampfinish第1行指明請(qǐng)?jiān)谳敵龃翱谥写蛴】藏悹枅D第2行是退出通用后處理器則在輸出窗口發(fā)現(xiàn)了下面的數(shù)據(jù) 它列出了對(duì)于3個(gè)給定轉(zhuǎn)速(0.1,376.991.753.982)時(shí)的前10階固有頻率。