【ANSYS算例】7.3(1)帶有張拉的繩索的振動模態(tài)分析(GUI)一條每單位長度質(zhì)量為的懸索兩端固定，由于受拉具有初始應變，試確定此時懸索的應力以及所受的拉力，并計算帶有張拉的懸索的前3階自由振動頻率；對該系統(tǒng)進行模態(tài)分析。模型中的各項參數(shù)如表7-3所示，為與文獻結果進行比較，這里采用了英制單位。表7-3帶有張拉的繩索模型的參數(shù)材料參數(shù)幾何參數(shù)載荷彈性模量E=30×106psi繩索長度=100in初始應變=0.00543228密度ρ=0.00073lb·sc2/in4繩索截面積A=0.00306796in2解答：計算模型如圖7-3所示。圖7-3帶有張拉的繩索的計算模型這里選取足夠數(shù)量的單元以確保能夠充分體現(xiàn)動態(tài)分析的特征。首先進行靜力分析，采用3D的單向桿單元(tension-onlyorcompression-onlysparelements：\o"LINK10description"LINK10)。建模的要點：⑴定義分析類型，選取單元、實常數(shù)和材料參數(shù)；⑵依據(jù)算例定義節(jié)點和單元，設置輸出選項和位移約束；⑶采用初應變來施加張力，⑷為了在后續(xù)的模態(tài)分析中考慮預緊的初始張力，需在靜力分析與模態(tài)分析過程中，應同時將預應力效應打開<PSTRES,ON>；⑸求解預應力引起的應力狀態(tài)；⑹提取計算結果，進行模態(tài)分析；⑺用命令<*GET>提取結果，對于一些線型單元(桿、梁、板)，需要采用單元列表<ETABLE>來提取應力。最后將計算結果與參考文獻所給出的解析結果進行比較，見表7-4。表7-4ANSYS模型與文獻的解析結果的比較

工況結果的物理量及單位Reference7.3(1)的結果ANSYS結果兩種結果之比靜態(tài)Flb500.500.1.000Stress,psi162974.162974.1.000模態(tài)f1,Hz74.70875.0941.005f2,Hz149.42151.291.012f3,Hz224.12229.681.025Reference7.3(1)：.ThomsonWT.VibrationTheoryandApplications.2ndPrinting.NJ:Prentice-Hall,Inc.1965,264給出的基于圖形界面的交互式操作(stepbystep)過程如下。(1)進入ANSYS（設定工作目錄和工作文件）程序→ANSYS→ANSYSInteractive→Workingdirectory（設置工作目錄）→Initialjobname（設置工作文件名String）:String→Run→OK(2)設置計算類型ANSYSMainMenu：Preferences…→Structural→OK(3)定義單元類型ANSYSMainMenu：Preprocessor→ElementType→Add/Edit/Delete...→Add…→Link:3DBilinear10→OK（返回到ElementTypes窗口）→選擇Type1LINK10單擊Options…→K3：Tensiononly→OK→Close(4)定義實常數(shù)ANSYSMainMenu:Preprocessor→RealConstants…→Add/Edit/Delete...→Add…→選擇Type1LINK10→OK→RealConstantsSetNo.:1（第1號實常數(shù)）,AREA:306796E-8(繩索截面積),ISTRN:543248E-8→OK→Close(關閉RealConstants窗口)(5)定義材料參數(shù)ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic：EX:30e6→Note：確定→OK，Density(定義材料密度)→輸入DENS:73e-5,→OK→關閉材料定義窗口(6)構造懸索模型生成節(jié)點ANSYSMainMenu：Preprocessor→Modeling→Create→Nodes→InActiveCS→Nodenumber：1，X，Y，ZLocationinactiveCS：0，0，0Apply→Nodenumber：14，X，Y，ZLocationinactiveCS：100，0，0→OKANSYSMainMenu：Preprocessor→Modeling→Create→Nodes→FillbetweenNds點擊1，14號節(jié)點→OK→OK生成元素ANSYSMainMenu:Preprocessor→Modeling→Create→Elements→AutoNumbered→ThruNodes→點擊1，2號節(jié)點，生成第一個單元ANSYSMainMenu:Preprocessor→Modeling→Copy→Elements→AutoNumbered→選擇第1個元素→OK→Totalnumberofcopies:13→OK(7)模型加約束ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→DisplacementOnNodes→PickAll→OK→Lab2:AllDOF(全約束)→OK(8)求解預應力引起的應力狀態(tài)ANSYSMainMenu:Solution→AnalysisType→NewAnalysis→Static→OKANSYSMainMenu:Solution→AnalysisType→Sol’ncontrols在Basic標簽下選中Calculateprestresseffects→OKANSYSMainMenu：Solution→Solve→CurrentLS→OK→Yes(9)存儲應力計算數(shù)據(jù)ANSYSMainMenu：GeneralPostproc→ElementTable→DefineTable→Add…→Lab:STRS,Resultsdataitem:Bysequencenumber,右側框中選擇LS，下面框中輸入LS，1→OK→CloseANSYSUtilityMenu:Parameters→GetScalarData→分別選擇Resultdata,Elemtabledata→OK→Nameofparametertobedefined:STRSS;ElementnumberN:13→OKANSYSMainMenu:TimeHistPostpro→DefineVariables→Add→選擇Reactionforces→Listofitems:1→OK→Item,CompDataitem:StructforceFX→OK→closeANSYSMainMenu:TimeHistPostpro→StoreDataANSYSUtilityMenu:Parameters→GetScalarData→分別選擇Resultdata,Time-histvar’s→OK→Nameofparametertobedefined:FORCE;VariablenumberN:2,Datatoberetrieved:MaximumvalVmax→OK(10)模態(tài)分析計算ANSYSMainMenu:Solution→AnalysisType→NewAnalysis→Modal→OK→Close(warning)ANSYSMainMenu:Solution→AnalysisType→AnalysisOptions→[MODOPT]BlockLanczos,No.ofmodestoextract:3Expandmodeshapes:Yes,No.ofmodestoexpand:3[PSTRES]:Yes→OK→OKANSYSMainMenu:Solution→DefineLoads→delete→Structural→Displacement→OnNodes→選擇2~13號節(jié)點→OK→選擇UX→OK刪除2~13號節(jié)點上的UX約束ANSYSMainMenu:Solution→DefineLoads→delete→Structural→Displacement→OnNodes→選擇2~13號節(jié)點→OK→選擇UY→OK刪除2~13號節(jié)點上的UY約束ANSYSMainMenu：Solution→Solve→CurrentLS→OK→Yes(11)存儲模態(tài)計算數(shù)據(jù)并顯示全過程計算結果ANSYSUtilityMenu:Parameters→GetScalarData→分別選擇Resultdata,Modalresults→OK→Nameofparametertobedefined:FREQ1;VariablenumberN:1,Datatoberetrieved:FrequencyFREQ→OKANSYSUtilityMenu:Parameters→GetScalarData→分別選擇Resultdata,Modalresults→OK→Nameofparametertobedefined:FREQ2;VariablenumberN:2,Datatoberetrieved:FrequencyFREQ→OKANSYSUtilityMenu:Parameters→GetScalarData→分別選擇Resultdata,Modalresults→OK→Nameofparametertobedefined:FREQ3;VariablenumberN:3,Datatoberetrieved:FrequencyFREQ→OKANSYSUtilityMenu：List→Status→Parameters→AllParameters(12)退出系統(tǒng)ANSYSUtilityMenu：File→Exit→SaveEverything→OK【ANSYS算例】7.3(2)帶有張拉的繩索的振動模態(tài)分析(命令流)針對【ANSYS算例】7.3(1)的GUI操作，提供完整的命令流。解答給出的命令流如下。!%%%%[ANSYS算例]7.3(2)%%%%%begin%%%%%/PREP7!進入前處理ANTYPE,STATIC!設置分析類型為靜力結構分析PSTRES,ON!用于后面的模態(tài)分析中ET,1,LINK10!選取單元類型1(單向桿單元)KEYOPT,1,3,0!設置僅承受拉應力，KEYOPT(3)=0R,1,306796E-8,543248E-8!設置實常數(shù)，包括繩索截面積(306796E-8)，初始應變(543248E-8)MP,EX,1,30E6!定義材料的彈性模量(1號材料)MP,DENS,1,73E-5!定義材料的密度(1號材料)N,1!定義第1號節(jié)點N,14,100!定義第14號節(jié)點FILL!均分填滿第2~13號節(jié)點E,1,2!由節(jié)點1及節(jié)點2生成單元EGEN,13,1,1!依序復制生成13個單元D,ALL,ALL!對所有節(jié)點施加固定約束FINISH!前處理結束/SOLU!進入求解模塊，求解預應力引起的應力狀態(tài)SOLVE!求解FINISH!退出求解模塊/POST1!進入一般的后處理ETABLE,STRS,LS,1!針對LINK10單元,建立單元列表STRS,通過LS及特征號1來獲得單元的軸向應力*GET,STRSS,ELEM,13,ETAB,STRS!針對單元列表STRS,提取13號單元的應力FINISH!后處理結束/POST26!進入時間歷程后處理，處理支反力RFORCE,2,1,F,X!提取1號節(jié)點上的X方向支反力，并存儲到2號變量中STORE