ANSYS作業(yè)四圓柱殼地非線性屈曲分析報(bào)告

1、實(shí)用標(biāo)準(zhǔn)文檔文案大全有限元原理與工程應(yīng)用題目圓柱殼的非線性屈曲分析年級(jí)與專業(yè)機(jī)械電子工程所在學(xué)院（系）機(jī)械工程學(xué)系梁結(jié)構(gòu)的瞬態(tài)完全法分析1要求用弧長(zhǎng)法進(jìn)行一個(gè)圓柱殼的非線性屈曲分析。一個(gè)對(duì)邊簡(jiǎn)支的圓柱殼，在其中心作用一個(gè)垂直的集中載荷。目的是分析當(dāng) 載荷大小為1000N時(shí)，A、B兩點(diǎn)的垂直位移（UY 。材料特性：EX=3.1x109Pa （楊氏模量）NUXY=0.3 （泊松比）幾何特性：R=2540 mm L=254 mmH=6.35 mm 9 =5.7 載荷：p=1000 N采用2D的shell單元定義了厚度，由于該圓柱殼結(jié)構(gòu)為對(duì)稱結(jié)構(gòu)，故實(shí)際 操作時(shí)選取1/4結(jié)構(gòu)對(duì)其進(jìn)行分析即可。2操作步

2、驟求解步驟（GUI方法）2.1定義工作文件名及工作標(biāo)題(1) 定義工作文件名：執(zhí)行 Utility Men u/File/Cha nge Job name命令，在彈 出的Cha nge Job name對(duì)話框中輸入文件名為“ buckle2. ”，同時(shí)勾上【Newlog and error files 】,單擊 川 按鈕。(2) 定義工作標(biāo)題：執(zhí)行Utility Menu/File/Cha nge Title命令，彈出ChangeTitle 對(duì)話框中輸入“ Analysis of the buckle ” ,然后單擊心I按鈕.。(3)重新顯示：執(zhí)行 Utility Menu/Plot/Repl

3、ot命令。2.2顯示工作平面(1)顯示工作平面，執(zhí)行 Utility Menu/WorkPlane/Display Working Plane命令。(2)關(guān)閉三角坐標(biāo)符號(hào)的顯示：執(zhí)行UtilityMenu/PlotCtrls/Wi ndowControls/Window Options 命令，彈出 Window Options 對(duì)話框。在 Location oftriad下拉列表框中選擇Not shown選項(xiàng)，單擊 按鈕。2.3建立圓柱殼結(jié)構(gòu)模型2.3.1定義單元類型(1)選擇單元類型：執(zhí)行 PreprocessorElement TypeAdd/Edit/Delete命令，彈出復(fù)選框【Ele

4、me nt Type，單擊“ Add”，彈出【Library of Eleme nt Types對(duì)話框中選擇“ Structural Shell ”和“Elastic 4node63，單擊 Apply 按鈕。(2) 定義單元實(shí)常數(shù)：執(zhí)行 PreprocessorReal C on sta ntsAdd /Edit/Delete 命令，彈出【Real Constants 復(fù)選框，單擊 “AdcT ，選擇 “Type1 Shell63 ”選項(xiàng)，單擊“ 0K，彈出如下圖所示對(duì)話框，將參數(shù)如圖3進(jìn)行設(shè)置，單擊 工按鈕。A 尺田 1 Cons-Unt Set Nunnber 1P for SHELL63

5、Element Type Reference No. 1Real Constant Set No.1Shell thickness at node 1 TKQ635at node J TKQ)at node K. TK(K)it node L TKCL)Elastic foundation stiffness EFSElement X-axi5 rotation THETABending mom of inertia ratio RM1Dist from mid surf to top CTOPDist from mid surf to bot CBOTAdded mass/unit area

6、A.DMSUAOKIApply 1匚an cel 1Help |VJ圖1 “ Real Constant”對(duì)話框(3) 定義材料參數(shù)：執(zhí)行 PreprocessorMaterial PropsMaterial Model 命令，彈出【Define Material Model Behavior】，在彈出的對(duì)話框中進(jìn)行設(shè)置， 依次點(diǎn)擊 StructuralLi nearElasticlsotropic，彈出【Lin ear IsotropicProperties for Material Number 1】，對(duì)話框中的參數(shù)設(shè)置如圖 9所示，單擊X 按鈕即可。A Lirear Isotropic

7、Properties for Material Number 1Linear kotropic Material Properties for Material Number 1T1Temper atmre5EXAdd Temperature Delete TemperatureGraphKCancel |H elpJ圖2“ Density”對(duì)話框2.3.2生成幾何模型生成點(diǎn)執(zhí)行 PreprocessorModelingCreateKeypoints In Active CS 命令，彈 出【Create Keypoints In Active Coordinate System對(duì)話框，依次在

8、“ 1,0,0 ”，“2,0,2000 ”，“3,2540,0 ”，“4,2540,254 ”，最后生成的結(jié)果如圖 11 所示。POINTETYPE NUM生成線L-K4ANJAJJ 17 2 Ji 314:37:114WYk_3Analyaxs pthe buckle生成曲面廠ALines about AxhAROTATLints 誦bat 心話ARCArc length in degrees5.7NSEjNq. of rea segmerrtsOKApplyCa n celHelpL/TYPE NUMANJAN 17 201319：3D：4L圖3生成矩形劃分網(wǎng)格ELEi-ftrrsANJA

9、M 17 20131*33;056,8,1ELZtHArialysis the buckleANJAH 17 20130VJlLOTTMIFAnalysis of the buckleANJAM 17 201319:53s322.4步驟二：恢復(fù)數(shù)據(jù)庫(kù)文件“ buckle2.db ”Utility Menu FileResume from2.5 步驟三：進(jìn)入求解器。Main Menu soluti on2.6 步驟四：定義分析類型和選項(xiàng)1、 選擇菜單路徑 Main MenuSolution-Analysis Type-New Analysis.單 擊“ Static ”來選中它然后單擊OK2、擇

10、菜單路徑 Main MenuSolutionUnabridged MenuAnalysis Options。 出現(xiàn)對(duì)話框。3、 單擊Large deform effects option（大變型效應(yīng)選項(xiàng)）使之為ON，然后 單擊OK2.7 步驟五：在結(jié)點(diǎn)1的Y方向施加一個(gè)大小為-250的力Mai n menu Solutio n -Load -Apply Force/Mome nt On Nodes2.8步驟六：設(shè)置載荷步選項(xiàng)1、擇菜單路徑 Main Menu SolutionUnabridged MenuLoad stepopts-Time/Freque nc Time and substps

11、，時(shí)間和時(shí)間步選項(xiàng)對(duì)話框出現(xiàn)。2、Number of substeps （ 子步數(shù)）鍵入 30。3、擇菜單路徑 Main Menu SolutionUnabridged Menu Load step opts-Output ctrls DB/Resuls File ，對(duì)數(shù)據(jù)庫(kù)和結(jié)果文件寫入的控制對(duì)話框 出現(xiàn)。4、在“ Item”中，選擇“ all ”5、對(duì)“ FREQ ,選擇“ Every Substep ”2.9 步驟七：選擇弧長(zhǎng)法1、選擇菜單路徑 Main Menu SolutionUnabridged MenuLoad stepopts-NonlinearArc-length opts。

12、Arc-length opts T （時(shí)間和時(shí)間步選項(xiàng)）對(duì)話框出現(xiàn)。2、將“ KEY （Arc-length method on/off） 設(shè)置為“ ON3、對(duì)“ MAXARC（Maximum multiplier）, 輸入 44、單擊OK2.10 步驟八：求解問題1、 選擇菜單路徑 Main MenuSolution-Solve-Current LS。2、檢閱狀態(tài)窗口中的信息然后單擊 close o3、單擊Solve Current Load Step（ 求解當(dāng)前載荷步）對(duì)話框中的OK開始求 解。2.11 步驟九：進(jìn)行所需要的后處理。3求解步驟（命令流方法）FINI/CLE/PREP7smr

13、t,offET,1,SHELL63,1R,1,6.350 ! SHELL THICKNESSMP,EX,1,3102.75MP,NUXY,1,0.3!CREATE FINITE ELEMENT MODELR1 = 2540 ! SHELL MID-SURFACE RADIUSL = 254 ! HALF THE LENGTHPI = 4*ATAN(1) ! VALUE OF PI COMPUTEDTHETA = 0.1*180/PI ! 0.1 RADIANS CONVERTED TO DEGREESCSYS,1 ! CYLINDRICAL CO-ORDINATE SYSTEMN,1,R1,9

14、0 ! NODES 1 AND 2 ARE CREATED AT POINTSN,2,R1,90,L ! A AND B RESPECTIVELY.K,1,R1,90K,2,R1,(90-THETA)K,3,R1,90,LK,4,R1,(90-THETA ),LESIZE,2 ! TWO DIVISION ALONG THE REGION BOUNDARYA,1,3,4,2AMESH,1NUMMRG,NODE!APPLY BOUNDARY CONDITIONSNSEL,S, LOC,Z,0DSYM,SYMM,ZNSEL,S, LOC,Y,90DSYM,SYMM,XNSEL,S, LOC,Y,(

15、90-THETA)D,ALL,UX,UY,UZNSEL,ALLFINISHSAVE,BUCKLE2,DBRESUME,BUCKLE2,DB/SOLUTIONANTYPE,STATICNLGEOM,ON ! LARGE DEFLECTION TURNED ONOUTRES,1 ! WRITE SOLUTION ON RESULTS FILE FOR EVERY SUBSTEP F,1,FY,-250 ! 1/4 TH OF THE TOTAL LOAD APPLIED DUE TO SYMMETRY NSUBST,30 ! BEGIN WITH 30 SUBSTEPSARCLEN,ON,4SOL

16、VEFINISH/POST26NSOL,2,1,U,Y ! STORE UY DISPLACEMENT OF NODE 1NSOL,3,2,U,Y ! STORE UY DISPLACEMENT OF NODE 2 PROD,4,1,LOAD,4*250 ! TOTAL LOAD IS 4*250 DUE TO SYMMETRYPROD,5,2,-1 ! CHANGE SIGNS OF THE DISPLACEMENT VALUES PROD,6,3,-1*GET,UY1,VARI,2,EXTREM,VMIN*GET,UY2,VARI,3,EXTREM,VMINPRVAR,2,3,4 ! PR

17、INT STORED INFORMATION/AXLAB,X, DEFLECTION (MM)/AXLAB,Y, TOTAL LOAD (N)/GRID,1/XRANGE,0,35/YRANGE,-500,1050XVAR,5PLVAR,4 ! PLOT LOAD WITH RESPECT TO -UY OFXVAL,6PLVAR,4FINISH變形情況（取的是中間一半，所以最邊緣的下去是最多的）S7EF= 3TJR =4：TniE=2IMC =31.67嚴(yán)-fNCOAL OLUSICM3TTF=r3TJR =41TIME：=1SX恆VRSYSOHQ=31.674SWtJAN17 2 JI 3

18、20:00=51-15,554-8H0L2,59213.66424.737-H .01-2 r 953, l?a13,20L30,273ANJAN 17 20131959:20DEFLECTION 柯）* ANSYS POST26 VARIABLE LISTING *TIMEUY1 UY2 UY4 PRODUYLOAD0.32910E-01-0.3691880.253726E-01 32.91010.98776E-01-1.166760.794300E-01 98.77610.21372-2.783370.179024213.7180.31393-4.479230.252211313.9260

19、.40097-6.204270.266634400.9670.47620-7.915110.186102476.1950.53934-9.57903-0.277909E-01 539.3420.58825-11.1678-0.416425588.2530.61938-12.6474-1.02361619.3790.62881-13.9696-1.88783628.8090.62455-14.5513-2.42280624.5540.60237-15.4301-3.52204602.3690.55743-16.1745-4.96401557.4260.49699-16.6688-6.537674

20、96.9890.42664-16.9572-8.17302426.6430.35056-17.0791-9.82553350.5600.27176-17.0549-11.4714271.7570.19281-16.8847-13.0985192.8070.11640-16.5532-14.6988116.4020.45547E-01-16.0402-16.262145.5466* ANSYS POST26 VARIABLE LISTING *TIME 1 UY 2 UY 4 PRODUYUYLOAD-0.17078E-01-15.3379-17.7744-17.0779-0.71163E-01-14.4770-19.2219-71.1634-0.12159-13.5555-20.6035-121.589-0.17848-12.7847-21.9226-178.484-0.24620-12.51