瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析
2017-04-26 by:CAE仿真在線 來源:互聯網
瞬態(tài)動力學分析(亦稱時間歷程分析)是用于確定承受任意的隨時間變化載荷結構的動力學響應的一種方法。可以用瞬態(tài)動力學分析確定結構在穩(wěn)態(tài)載荷、瞬態(tài)載荷和簡諧載荷的隨意組合作用下的隨時間變化的位移、應變、應力及力。載荷和時間的相關性使得慣性力和阻尼作用比較重要。如果慣性力和阻尼作用不重要,就可以用靜力學分析代替瞬態(tài)分析。
瞬態(tài)動力分析用在以下的設計中:
-
承受各種沖擊載荷的結構,如:汽車中的門和緩沖器、建筑框架以及懸掛系統(tǒng)等。
-
承受各種隨時間變化載荷的結構,如:橋梁、地面移動裝置以及其它機器部件。
-
承受撞擊和顛簸的家庭和辦公設備,如:便攜式電話、筆記本電腦和真空吸塵器等。
下面將展現最簡單的單自由度系統(tǒng)的自由振動的瞬態(tài)動力學ANSYS分析,重點在于掌握瞬態(tài)動力學的分析方法,特別是時間離散中時間步長的選取、初始條件的施加(本題為初位移)、時間歷程后處理等。進行模態(tài)分析計算固有頻率,其目的是為了時間步長的選取,時間步長一般取固有周期的1/40~1/20,采用此方法確定復雜模型的時間步長非常有用。
一、問題描述
一個單自由度系統(tǒng)的自由振動,質量m= 0.5 kg,剛度k=5 kN/m,初始向下受拉力,使彈簧伸長0.2m,然后突然釋放,計算位移隨時間歷程響應情況。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys仿真分析圖片1](http://m.09bigdata.com/i/g/a6/a6034d4dab3fc49f122b279e6247a3dfr.png)
阻尼比取4種:
阻尼比1:ξ=0.0(無阻尼);
阻尼比2:ξ=0.02(欠阻尼);
阻尼比3:ξ=1.0(臨界阻尼);
阻尼比4:ξ=2.0(過阻尼)。
問題分析:
單自由度系統(tǒng)采用COMBIN40建模,實常數中輸入剛度k,阻尼系數c和質量m。采用Newmark直接積分法,時間步長取固有周期1/30,約為0.002s。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys仿真分析圖片2](http://m.09bigdata.com/i/g/01/01a503a4612da9ce3bb38de95bcbc7e5r.png)
阻尼比1:ξ=0(無阻尼)時,阻尼系數c=0;
阻尼比2:ξ=0.02(欠阻尼)時,阻尼系數c=2;
阻尼比3:ξ=1.0(臨界阻尼)時,阻尼系數c=100;
阻尼比4:ξ=2.0(過阻尼)時,阻尼系數c=200。
計算結果:ANSYS計算得出的固有頻率為15.915Hz,不同阻尼比條件下得出的位移隨時間變化曲線見下圖。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys結構分析圖片3](http://m.09bigdata.com/i/g/b8/b88613b66a666a0613cc67bad58fbd0dr.png)
二、GUI步驟
(一)建立模型
1.進入ANSYS
程序→ ANSYS 15.0→ ANSYS Product Launcher→ 改變working directory到指定文件夾→ 在job name輸入:file。
2.設置計算類型
Main Menu> Preferences→選擇Structural→ OK。
3.定義單元屬性
(1)定義單元類型
Main Menu> Preprocessor>Element Type> Add/Edit/Delete→ Add→ 選擇COMBIN40單元,即在左列表框中選擇Combination,在右列表框中選擇Combination 40→ OK。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys結構分析圖片4](http://m.09bigdata.com/i/g/3b/3b64d480e75d8950649361fc7ccd089br.jpg)
(2)單元選項設置:選中COMBIN40單元→ Option→ 在K3中選擇UY;在K4中選擇Closegap only;在K6中選擇Mass at node J→ OK→ Close。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析圖片6](http://m.09bigdata.com/i/g/35/35d4538302dc9a0cd3ec98aa634fc292r.jpg)
4.定義實常數
Main Menu> Preprocessor> Real Constants> Add→Type 1→ OK→ Real Constant Set No.:1,彈簧剛度K1:5000,阻尼系數C:2,質量M:0.5→ OK→ Close。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析圖片7](http://m.09bigdata.com/i/g/a5/a5a3d814db562c4e4e492a02dc777bf5r.jpg)
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。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析圖片8](http://m.09bigdata.com/i/g/48/48ca22300de604618520f3e319a9e7e7r.jpg)
(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。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys培訓的效果圖片9](http://m.09bigdata.com/i/g/c6/c6e637ea756b4eab231fc15b786e065br.jpg)
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。重要步驟!
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys培訓的效果圖片10](http://m.09bigdata.com/i/g/48/48075e17a188141c957dc88afd3789der.jpg)
(2)設置模態(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。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys培訓的效果圖片11](http://m.09bigdata.com/i/g/8a/8aaecdfa48eb9a33476e3f44730dc71br.jpg)
(3)求解:Main Menu > Solution >Solve > Current LS→ File > Close→ Solve Current Load Step→ OK→ Solutionis done→ Close。
2.后處理
列表查看各階固有頻率:Main Menu > GenenralPostproc > Results Summary→ 記錄各階固有頻率,對應的是TIME/FREQ。單自由度只有第1階固有頻率,固有頻率為15.915Hz,固有周期1/15.915=0.0628s。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys仿真分析圖片12](http://m.09bigdata.com/i/g/c9/c955e526bca529b8f9dfc7a25bfd957cr.jpg)
(三)瞬態(tài)動力學分析
1.調入模型
Utility Menu > Files > Resum from→ 選擇mesh.db→ OK。
2.瞬態(tài)動力學分析設置
(1) 進入求解器:Main Menu> Solution。
(2) 設置分析類型:Main Menu> Solution>Analysis Type> New Analysis→ 選擇Transient→ OK→ 求解方法選擇Full(完全積分法)→ OK。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys仿真分析圖片13](http://m.09bigdata.com/i/g/07/077ea0fd7e8bbeb25ca4ae3f421941e8r.jpg)
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析圖片14](http://m.09bigdata.com/i/g/8b/8b6504847779d46d474f0698c251dadbr.jpg)
3.第1載荷步求解設置
(1) 求解控制設置:Main Menu> Solution> AnalysisType> Sol’n Controls。
① Basic:時間控制Time Control中, Time at end of loadstep(載荷步結束時間):1e-6;Time increment(時間增量)中Time step size(時間步長):0.002;Frequency(輸出頻率)中選擇Write every Nth substep,where N=1。時間步長取固有周期0.0628s/30約為0.002s。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析圖片15](http://m.09bigdata.com/i/g/22/228c658c3154e0616c314161182b55bdr.jpg)
② Transient:Full Transient Option中選擇Transient effects(瞬態(tài)影響),選擇Stepped Loading(階躍加載)。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys圖片16](http://m.09bigdata.com/i/g/7d/7def512fcfe83e78b81c8fb477082339r.jpg)
(2) 設置初值條件:Main Menu> Solution> DefineLoads> Apply> Initial Condit'n> Define→ 拾取節(jié)點2→ OK→ Lab:UY;VALUE:-0.2→ OK。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys圖片17](http://m.09bigdata.com/i/g/f4/f441033a334d21c1403cfc5be027b05fr.jpg)
(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載荷步求解設置
(1) 求解控制設置
Main Menu> Solution> Analysis Type> Sol’nControls→ Basic:時間控制Time Control中,Time at end of loadstep:4.5*0.0628。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys培訓的效果圖片18](http://m.09bigdata.com/i/g/3c/3c50d5a962e0cfa61bc881dbd642ded3r.jpg)
(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。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys培訓的效果圖片19](http://m.09bigdata.com/i/g/04/042e69cfd58695bbf467df53eff6eae1r.jpg)
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析案例圖片20](http://m.09bigdata.com/i/g/8d/8debd8c8151b3f1a98bc4c78bbf85437r.jpg)
(3)顯示位移隨時間變化曲線:Main Menu> TimeHistPostpro> Variable Viewer→ 選中UY_2→ 選擇紅色方框中的按鈕。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析案例圖片21](http://m.09bigdata.com/i/g/c2/c2c131727acc68124c5c13cc82fa2884r.jpg)
(4)列表顯示位移隨時間變化數據:Main Menu> TimeHistPostpro> Variable Viewer→ 選中UY_2→ 選擇紅色方框中的按鈕。
![瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析ansys分析案例圖片22](http://m.09bigdata.com/i/g/13/134c0d8b7d333d274968f98b28078cb5r.jpg)
三、ADPL步驟
k1=5e3 !剛度
m1=0.5 !質量
pi=acos(-1)
fre1=(1/(2*pi))*(k1/m1)**0.5 !固有頻率,Hz
T1=1/fre1 !固有周期,s
Dt= T1/30
kesi=0.02 !阻尼比或相對阻尼系數0/0.02
c1=2* Kesi*( m1*k1)**0.5 !阻尼系數
/PREP7 !進入前處理器
ET,1,COMBIN40 !單元類型
KEYOPT,1,3,2 !單元的自由度UY
KEYOPT,1,4,1 !關閉間隙
KEYOPT,1,6,2 !質量在J節(jié)點
R,1,k1,c1,m1 !彈簧剛度、阻尼系數和質量
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 !列固有頻率結果
RESUME, mesh,db, !讀入模型
/SOLU !進入求解器
ANTYPE,TRANS !瞬態(tài)動力學分析
TRNOPT,FULL !完全積分法
TIME,1E-6 !第1載荷步
DELTIM, T1/30 !積分時間步長
OUTRES,ALL,1 !輸出每步結果
KBC,1 !階躍加載
IC,2,UY,-0.2 !初始條件
SOLVE !求解
TIME,4.5*T1 !第2載荷步
SOLVE !求解
FINISH !退出求解器
/POST26 !進入時間歷程處理器
FILE,,RST !調入結果文件
NSOL,2,2,U,Y,1UY !節(jié)點2的Y向位移
NPRINT,1 !每1個點打印1次
PRVAR,2 !畫圖顯示曲線
/XRANGE,0,4.5* T1 !X軸范圍
PLVAR,2 !畫圖顯示曲線
相關標簽搜索:瞬態(tài)動力學專題-單自由度系統(tǒng)自由振動ANSYS分析 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓