瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析

2017-04-26  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

瞬態(tài)動力學分析(亦稱時間歷程分析)是用于確定承受任意的隨時間變化載荷結(jié)構(gòu)的動力學響應(yīng)的一種方法。可以用瞬態(tài)動力學分析確定結(jié)構(gòu)在穩(wěn)態(tài)載荷、瞬態(tài)載荷和簡諧載荷的隨意組合作用下的隨時間變化的位移、應(yīng)變、應(yīng)力及力。載荷和時間的相關(guān)性使得慣性力和阻尼作用比較重要。如果慣性力和阻尼作用不重要,就可以用靜力學分析代替瞬態(tài)分析。

瞬態(tài)動力分析用在以下的設(shè)計中:

• 承受各種沖擊載荷的結(jié)構(gòu),如:汽車中的門和緩沖器、建筑框架以及懸掛系統(tǒng)等。

• 承受各種隨時間變化載荷的結(jié)構(gòu),如:橋梁、地面移動裝置以及其它機器部件。

• 承受撞擊和顛簸的家庭和辦公設(shè)備,如:便攜式電話、筆記本電腦和真空吸塵器等。

下面將展現(xiàn)最簡單的單自由度系統(tǒng)的自由振動的瞬態(tài)動力學ANSYS分析,重點在于掌握瞬態(tài)動力學的分析方法,特別是時間離散中時間步長的選取、初始條件的施加(本題為初位移)、時間歷程后處理等。進行模態(tài)分析計算固有頻率,其目的是為了時間步長的選取,時間步長一般取固有周期的1/40~1/20,采用此方法確定復(fù)雜模型的時間步長非常有用。

---

一、問題描述

一個單自由度系統(tǒng)的自由振動,質(zhì)量m= 0.5 kg,剛度k=5 kN/m,初始向下受拉力,使彈簧伸長0.2m,然后突然釋放,計算位移隨時間歷程響應(yīng)情況。

阻尼比取4種:

阻尼比1:ξ=0.0(無阻尼);

阻尼比2:ξ=0.02(欠阻尼);

阻尼比3:ξ=1.0(臨界阻尼);

阻尼比4:ξ=2.0(過阻尼)。

問題分析:

單自由度系統(tǒng)采用COMBIN40建模,實常數(shù)中輸入剛度k,阻尼系數(shù)c和質(zhì)量m。采用Newmark直接積分法,時間步長取固有周期1/30,約為0.002s。

阻尼比1:ξ=0(無阻尼)時,阻尼系數(shù)c=0;

阻尼比2:ξ=0.02(欠阻尼)時,阻尼系數(shù)c=2;

阻尼比3:ξ=1.0(臨界阻尼)時,阻尼系數(shù)c=100;

阻尼比4:ξ=2.0(過阻尼)時,阻尼系數(shù)c=200。

計算結(jié)果:ANSYS計算得出的固有頻率為15.915Hz,不同阻尼比條件下得出的位移隨時間變化曲線見下圖。

二、GUI步驟

---

(一)建立模型

1.進入ANSYS

程序→ ANSYS 15.0→ ANSYS Product Launcher→ 改變working directory到指定文件夾→ 在job name輸入:file。

2.設(shè)置計算類型

Main Menu> Preferences→選擇Structural→ OK。

3.定義單元屬性

(1)定義單元類型

Main Menu> Preprocessor>Element Type> Add/Edit/Delete→ Add→ 選擇COMBIN40單元,即在左列表框中選擇Combination,在右列表框中選擇Combination 40→ OK。

(2)單元選項設(shè)置:選中COMBIN40單元→ Option→ 在K3中選擇UY;在K4中選擇Closegap only;在K6中選擇Mass at node J→ OK→ Close。

4.定義實常數(shù)

Main Menu> Preprocessor> Real Constants> Add→Type 1→ OK→ Real Constant Set No.:1,彈簧剛度K1:5000,阻尼系數(shù)C:2,質(zhì)量M:0.5→ OK→ Close。

5.建立模型

(1)創(chuàng)建節(jié)點:Main Menu > Preprocessor> Create > Nodes> In Active CS→ NODE:1,X,Y,Z:0,0,0→ Apply→ NODE:2,XYZ輸入0,-0.5,0。

(2)創(chuàng)建單元:

①Main Menu> Preprocessor> Create> Elements>Auto Numbered> Thru Nodes→ 輸入節(jié)點1,回車,再輸入節(jié)點2→ OK。第1個節(jié)點是節(jié)點I,第二個節(jié)點是節(jié)點J。

②查看單元屬性:Utility Menu> List> Elements>Nodes + Attributes。

6.施加約束

Main Menu> Solution> Define loads> Apply> Structural>Displacement> On Nodes→ 拾取節(jié)點1→ OK→ Lab2:UY→ OK。

7.保存模型

Utility Menu > Files > Save as→ Save Database to輸入mesh.db→ OK。

(二)模態(tài)分析,計算固有頻率

1.模態(tài)分析求解

(1) 進入求解器:Main Menu> Solution。

(1)指定分析類型為模態(tài)分析:Main Menu > Solution >Analysis Type > New Analysis→ Modal→ OK。重要步驟!

(2)設(shè)置模態(tài)分析選項:Main Menu > Solution >Analysis Type > Analysis Options→ Mode extraction method:Block Lanczos(模態(tài)提取方法,默認);No. of modes to extract(模態(tài)提取階次):2;NMODE No. of modes to expand(模態(tài)擴展階次):2→ OK。

(3)求解:Main Menu > Solution >Solve > Current LS→ File > Close→ Solve Current Load Step→ OK→ Solutionis done→ Close。

2.后處理

列表查看各階固有頻率:Main Menu > GenenralPostproc > Results Summary→ 記錄各階固有頻率,對應(yīng)的是TIME/FREQ。單自由度只有第1階固有頻率,固有頻率為15.915Hz,固有周期1/15.915=0.0628s。

(三)瞬態(tài)動力學分析

1.調(diào)入模型

Utility Menu > Files > Resum from→ 選擇mesh.db→ OK。

2.瞬態(tài)動力學分析設(shè)置

(1) 進入求解器:Main Menu> Solution。

(2) 設(shè)置分析類型:Main Menu> Solution>Analysis Type> New Analysis→ 選擇Transient→ OK→ 求解方法選擇Full(完全積分法)→ OK。

3.第1載荷步求解設(shè)置

(1) 求解控制設(shè)置:Main Menu> Solution> AnalysisType> Sol’n Controls。

① Basic:時間控制Time Control中, Time at end of loadstep(載荷步結(jié)束時間):1e-6;Time increment(時間增量)中Time step size(時間步長):0.002;Frequency(輸出頻率)中選擇Write every Nth substep,where N=1。時間步長取固有周期0.0628s/30約為0.002s。

② Transient:Full Transient Option中選擇Transient effects(瞬態(tài)影響),選擇Stepped Loading(階躍加載)。

(2) 設(shè)置初值條件:Main Menu> Solution> DefineLoads> Apply> Initial Condit'n> Define→ 拾取節(jié)點2→ OK→ Lab:UY;VALUE:-0.2→ OK。

(3) 保存求解前文件:Utility Menu> Files>Save as→ 輸入Load.db→ OK。

(4) 求解:Main Menu> Solution>Solve> Current LS→ File> Close→ Solve Current Load Step→ OK→ Solution isdone→ Close。

4.第2載荷步求解設(shè)置

(1) 求解控制設(shè)置

Main Menu> Solution> Analysis Type> Sol’nControls→ Basic:時間控制Time Control中,Time at end of loadstep:4.5*0.0628。

(2) 求解:Main Menu> Solution>Solve> Current LS→ File> Close→ Solve Current Load Step→ OK→ Solution isdone→ Close。

5.時間歷程后處理

(1) 進入時間歷程后處理器:Main Menu> TimeHistPostpro。

(2) 定義位移隨時間變化變量:Main Menu> TimeHist Postpro> VariableViewer→ 選擇紅色方框中的按鈕→ DOF Solution→ Y-Component of displacement→ OK→ 輸入節(jié)點2→ OK。

(3)顯示位移隨時間變化曲線:Main Menu> TimeHistPostpro> Variable Viewer→ 選中UY_2→ 選擇紅色方框中的按鈕。

(4)列表顯示位移隨時間變化數(shù)據(jù):Main Menu> TimeHistPostpro> Variable Viewer→ 選中UY_2→ 選擇紅色方框中的按鈕。

---

三、ADPL步驟

k1=5e3 !剛度

m1=0.5 !質(zhì)量

pi=acos(-1)

fre1=(1/(2*pi))*(k1/m1)**0.5 !固有頻率,Hz

T1=1/fre1 !固有周期,s

Dt= T1/30

kesi=0.02 !阻尼比或相對阻尼系數(shù)0/0.02

c1=2* Kesi*( m1*k1)**0.5 !阻尼系數(shù)

/PREP7 !進入前處理器

ET,1,COMBIN40 !單元類型

KEYOPT,1,3,2 !單元的自由度UY

KEYOPT,1,4,1 !關(guān)閉間隙

KEYOPT,1,6,2 !質(zhì)量在J節(jié)點

R,1,k1,c1,m1 !彈簧剛度、阻尼系數(shù)和質(zhì)量

N,1,0,0,0 !創(chuàng)建節(jié)點1

N,2,0,-0.5,0 !創(chuàng)建節(jié)點2

E,1,2 !創(chuàng)建單元

D,1,UY !施加約束

FINISH !退出前處理器

SAVE,mesh,db, !保存模型

/SOLU !進入求解器

ANTYPE,2 !模態(tài)分析

MODOPT,LANB,2 !提取前2階

MXPAND,2, , ,0 !擴展前2階

SOLVE !求解

FINISH !退出求解器

/POST1 !通用后處理

SET,LIST !列固有頻率結(jié)果

RESUME, mesh,db, !讀入模型

/SOLU !進入求解器

ANTYPE,TRANS !瞬態(tài)動力學分析

TRNOPT,FULL !完全積分法

TIME,1E-6 !第1載荷步

DELTIM, T1/30 !積分時間步長

OUTRES,ALL,1 !輸出每步結(jié)果

KBC,1 !階躍加載

IC,2,UY,-0.2 !初始條件

SOLVE !求解

TIME,4.5*T1 !第2載荷步

SOLVE !求解

FINISH !退出求解器

/POST26 !進入時間歷程處理器

FILE,,RST !調(diào)入結(jié)果文件

NSOL,2,2,U,Y,1UY !節(jié)點2的Y向位移

NPRINT,1 !每1個點打印1次

PRVAR,2 !畫圖顯示曲線

/XRANGE,0,4.5* T1 !X軸范圍

PLVAR,2 !畫圖顯示曲線

相關(guān)標簽搜索：瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓