1、ANSYS疲勞分析的例子ansys疲勞分析例子/units,si /title, Fatigue analysis of cylinder with flat head! *參數(shù)設定*Di=1000                                ! 筒體內(nèi)徑t=20          

2、60;                       ! 筒體厚度hc=nint(4*sqrt(Di/2*t)/10)*10          ! 模型中筒體長度tp=60                        &

3、#160;         ! 平板封頭厚度r1=10                                  ! 平板封頭外測過渡圓弧半徑r2=10                   

4、              ! 平板封頭內(nèi)側(cè)應力釋放槽圓弧半徑exx=2e5                                ! 材料彈性模量mu=0.3              &#

5、160;                  ! 材料泊松比p1=2                                   ! 最高工作壓力p3=2.88        

6、0;                       ! 水壓試驗壓力n1=2e4                                 ! 最高/最低壓力循環(huán)次數(shù)n2=5     &

7、#160;                             ! 水壓試驗次數(shù)! *前處理*/prep7et,1,82                               &#

8、160;! 設定單元類型keyopt,1,3,1                           ! 設定周對稱選項mp,ex,1,exx                            ! 定義材料彈性模量mp,nuxy,1,m

9、u                           ! 定義材料泊松比!  *  建立模型  *k,1,0,0                            

10、    ! 定義關鍵點k,2,Di/2+t,k,3,Di/2+t,-(tp+hc)k,4,Di/2,-(tp+hc)k,5,Di/2,-tpk,6,Di/2-r2,-tp                        ! 定義應力釋放槽圓弧中心關鍵點k,7,0,-tpl,1,2               &

11、#160;                  ! 生成線l,2,3l,3,4l,4,5l,5,7l,7,1LFILLT,1,2,r1                          ! 生成外測過渡圓弧   al,all        &

12、#160;                        ! 生成子午面 CYL4, kx(6),ky(6), r2,180              ! 生成應力釋放槽面域   ASBA,1,2             

13、60;                ! 面相減wprot,90                             ! 旋轉(zhuǎn)工作平面wpoff,kx(6)-3*r2               

14、      ! 移動工作平面 asbw,all                               ! 用工作平面切割子午面wprot,90                        &#

15、160;     ! 旋轉(zhuǎn)工作平面wpoff,tp+r2                          ! 移動工作平面 asbw,all                            

16、  ! 用工作平面切割子午面esize,5                                ! 設定單元尺寸MSHKEY,1                           

17、;    ! 設定映射剖分amesh,1                                ! 映射剖分面1amesh,3                          &#

18、160;     ! 映射剖分面3esize,2                                ! 設定單元尺寸MSHKEY,0                      

19、        ! 設定自由剖分amesh,4                                ! 自由剖分面4fini                      

20、            ! 退出前處理! *求解*/solu                                  ! 筒體端部施加軸向約束dl,3,uy             

21、60;                ! 筒體端部施加軸向約束dl,6,symm                             ! 平板封頭對稱面施加對稱約束time,1              &

22、#160;                  ! 載荷步1lsel,s,8                             ! 選擇內(nèi)表面各線段lsel,a,11,13lsel,a,15cm,lcom1,line      

23、0;                   ! 生成內(nèi)表面線組件SFL,all,PRES,p1,                       ! 內(nèi)表面施加內(nèi)壓alls                

24、                  ! 全選solve                                  ! 求解fini            

25、0;                      ! 退出求解器! *后處理*/post1                                 ! 進入后處理FTSIZE,1,2,2,    

26、;                      ! 設定ubcolor=#ff0000疲勞/color/b/u評定的位置數(shù)、事件數(shù)及載荷數(shù)FP,1,1e1,2e1,5e1,1e2,2e2,5e2           ! 根據(jù)ubcolor=#ff0000疲勞/color/b/u曲線輸入SN數(shù)據(jù)FP,7,1e3,2e3,5e3,1e4,2e4,5e4FP,13,1e5,2

27、e5,5e5,1e6, ,FP,19, ,FP,21,4000,2828,1897,1414,1069,724  FP,27,572,441,331,262,214,159   FP,33,138,114,93.1,86.2, ,  FP,39, ,! *  水壓試驗循環(huán) *fs,4760,1,1,1,0,0,0,0,0,0              ! 儲存節(jié)點4760對應其第一載荷的應力set,1,last

28、60;                            ! 讀入第一載荷步數(shù)據(jù)FSNODE,4760,1,2                        ! 儲存節(jié)點4760對應其第二載荷的應力fe,1,n2,p3/p1  

29、0;                       ! 設定事件循環(huán)次數(shù)及載荷比例系數(shù)! *  最高/最低壓力循環(huán) *fs,4760,2,1,1,0,0,0,0,0,0              ! 儲存節(jié)點4760對應其第一載荷的應力set,1,last       

30、                     ! 讀入第一載荷步數(shù)據(jù)FSNODE,4760,2,2                        ! 儲存節(jié)點4760對應其第二載荷的應力FE,2,n1,1,          

31、;                   ! 設定事件循環(huán)次數(shù)及載荷比例系數(shù) FTCALC,1                               ! 進行ubcolor=#ff0000疲勞/color/b/u計算（并記錄使用系數(shù)）fini2./prep7et,1

32、,plane42MPTEMP,1,0  mpdata,ex,1,2e8mpdata,prxy,1,0.3rectng,0,200,0,100,cyl4,100,50,25asba,1,2smrtsize,3amesh,allfinish/solvnsel,s,loc,x,0d,all,uxd,1,uysfl,2,pres,0,31allsel,allsolvefinish/post1plnsol,s,eqv,0,1FP,1,100,200,500,1e3,1500,2e3   FP,7,1e4,1e5,1e6,2e6,3e6,5e6FP,13,6e6,7

33、e6,8e6,9e6,10e6,11e6 FP,19,12e6,15e6 FP,21,150,120,110,100,95,90 FP,27,85,80,75,70,65,60 FP,33,55,50,45,40,35,30 FP,39,29,25 *set,n_num,node(100,75,0)fl,1,n_numfsnode,n_num,1,1,fs,n_num,1,2,1,fe,1,10000,2,evelftcalc,1這有個疲勞的實例，不過沒見他加周期載荷，結(jié)果也看不懂商業(yè)軟件msc.fatigue我以前做過疲勞分析方面的，那時候是用ANSYS做的，覺得效果很不好，感覺ANSYS對

34、材料SN曲線主要進行的是線性插值法處理的，誤差蠻大的msc.fatigue這個軟件做疲勞壽命分析還是滿不錯的，ABAQUS做壽命分析還沒接觸過/Ha-Q1d 不知道目前有沒有人研究剩余壽命，感覺目前這塊兒還是個空白要做好疲勞分析，首先是要做好靜強度分析，)W I8DQ.YJW7e y!j$"qP1s V/_a 一般疲勞分析都是以靜強度計算結(jié)果得到應力應變，然后以此為基礎，進行計算疲勞的！%E$X j3?-!r y !P5cfJz 發(fā)動機方面目前用的比較多的都是FEMFATMSC_Fatigue疲勞分析應用與實例結(jié)構(gòu)疲勞壽命分析姚衛(wèi)星：我可以在MSC.FATIGUE中定義自己的材料常數(shù)

35、嗎？  & N% d/ , h! 6 ! Y2 n/ y( Q* g7 N：可以。利用所附的材料管理模快PFMAT就可輸入您個人使用的材料常數(shù)，例如楊氏系數(shù)(Young's modulus, E)、最大拉力(UTS)以及相關的疲勞常數(shù)等等。細節(jié)可參考User's Manual Vol.1的第三章。  " ' W  + k4 c6 N  + Q8 p#   U' 6 N+ ：MSC.FATIGUE可以與哪些FEA的軟件配合？  

36、( F9 G+ z. ?4 . M( H, a1 ?：因MSC.FATIGUE早已與MSC/PATRAN完全整合，所以只要是PATRAN可以接受的有限元分析結(jié)果，F(xiàn)ATIGUE都可以通過讀取PATRAN的*.db文件讀到應力或應變的結(jié)果進行疲勞計算；例如，MSC.NASATRAN、ABAQUS、MARC、ANSYS、MSC/PATRAN's FEA、MSC/AFEA等等。  / a* S, u. 0 X! Z. L$ t0 j/ V  9 R3 ?/ g/ Z1 O* ：MSC.FATIGUE適用于2D或3D模型？  ) X&

37、#39; m: I% ?8 i5 U8 y：MSC.FATIGUE本身并沒有限制模型的種類，因為它只讀取經(jīng)過有限元分析計算后的應力或應變張量(tensor)，若從后處理的角度來看則只支持2D與3D(shell & solid)的模型。  9 l" 3 d* P- % F. v  8 y/ L" M& k& * w0 L: ：如何建立正確的有限元分析模型以進行疲勞壽命評估？  $ M' D* i* X7 T, . |* P+ f：  " P &#

38、160;1 c( T& m# b  n3   : t% E7 x2 r$ p* x) - C" B$ N幾何與約束條件必須準確描述。  1 i& _! V  s. M: m& V應用板單元時得確定在厚度方向的變形可以忽略。  - v/ 1 _0 G) c. c, L7 y/ 使用的單元種類必須考慮到計算后的節(jié)點應力與應變能充分提供疲勞分析所要的參數(shù)，必要時得采用較高階的單元。  % W9 M4 s& N. & # G在應力集中

39、、幾何不連續(xù)或疲勞可能發(fā)生位置，單元的密度較高，主要是希望能正確算出這一區(qū)的應力。  + i2 C& X, A* C. B4 O( j盡量少用三角形(triangular)與楔形(wedge)單元，若已使用則須注意其單元品質(zhì)。在不同單元相接位置或不同厚度的板元相接時，須注意到其產(chǎn)生的應力不連續(xù)(fictitious stress raisers)要盡量降低。  ( A" 2 d  Q8 X' ; d' V6 c) R) i! E若可能的話，可將有限單元分析的應變與應變儀測得的結(jié)果作比較。  : b% N: _. V1 K" b/ x( w! ; b- M  9 4 n) c! h$ T：應取哪一個應力來做疲勞分析較好，例如：Von Mises, Max. Abs. Principal等？  . W/ v( K$ s5 p9 ：以MSC.NAST