SeyrüseferSeyir DefteriMühendislik → Takviyeli Konik Kabukların Doğrusal Burkulması

Takviyeli Konik Kabukların Doğrusal Burkulması

Seyir Defteri - Mühendislik
Cumartesi, 17 Kasım 2018

Takviyeli Koni BurkulmasıÖnceki bölümde Takviyesiz İnce Cidarlı Konik Kabuklar sonlu elemanlar yöntemi yaklaşımı ile kısaca ele alındıktan sonra denizaltı mukavim tekneleri üzerinde de yaygın olarak kullanılmasıyla önemli bir yapısal alt başlık olan kemerelerle desteklenmiş takviyeli koni konusuna sıra geldi.

Bu tür bir çalışmada denel doğrulama faaliyeti son derece önemli olduğu için konuya yine bu doğrultuda bir giriş yapmak uygun olacaktı ve kolayca ulaşılabileceği için [1] üzerinde yayınlanan deney ve hesaplama çalışmalarının temel alınması tercih edildi.

Üç farklı ölçüdeki kemere ile desteklenmiş eşdeğer simetrik metal konilerin harici hidrostatik yük altındaki plastik burkulmasının incelendiği söz konusu makalede yayınlanan denel ölçüm verileri ve hesaplama sonuçları ile bir karşılaştırma yapabilmek için öncelikle açık-kaynak Gmsh [2] yazılımı kullanılarak önişlem safhası gerçekleştirildi.

Bu tür bir yapıyı farklı şekillerde modellemek düşünülebilirdi, burada hem koni hem de kemereler için 2B cidar elemanları kullanılmak yönünde bir tercih yapıldı. Hesaplamaları için ise Calculix ccx kullanıldı. Aslına bakılırsa Calculix [3] gerek kiriş gerekse kabuk elemanları için de katı elemanlar kullanmakta, bir ve iki boyutlu elemanları da hesaplamadan önce üç boyutlu hâle dönüştürmektedir.


♦ gmsh:
  1. // Takviyeli İnce Cidarlı Koni Burkulma Doğrulama Çalışması İçin
  2. // Kaynak: Plastic buckling of ring-stiffened conical shells
  3. // under external hydrostatic pressure, 2004
  4. // Carl T.F. Ross - Andrew P.F. Little - Kehinde A. Adeniyi
  5. // Akma Mukavemeti= 250MPa / Esneklik Katsayısı= 193e09Pa / Poisson Oranı= 0,3
  6. // Üsküdar Mühendishanesi - http://uskudar.biz
  7.  
  8. // Konilerin Tanımlalanmaları
  9. // Koni10
  10. A = 0.0025;
  11. B = 0.0015;
  12. t = 0.0020;
  13.  
  14. // Koni11
  15. //A = 0.0025;
  16. //B = 0.0020;
  17. //t = 0.0020;
  18.  
  19. // Koni12
  20. //A = 0.0025;
  21. //B = 0.0025;
  22. //t = 0.0020;
  23.  
  24. // Temel Ölçüler
  25. r1 = 0.101675/2; // büyük yarıçap, m
  26. r2 = 0.060000/2; // küçük yarıçapı, m
  27. x0 = 0.116500; // uzunluk, m
  28. x1 = 0.013250; // kısa aralık
  29. x2 = 0.022500; // uzun aralık
  30.  
  31. Geometry.CopyMeshingMethod = 1;
  32.  
  33. Point(0) = {0, 0, 0}; // merkez
  34.  
  35. Point(10) = {0, r2+t, 0};
  36. Point(11) = {0.01325,0.032369930+t,0};
  37. Point(12) = {0.03575,0.036394340+t,0};
  38. Point(13) = {0.05825,0.040418750+t,0};
  39. Point(14) = {0.08075,0.044443160+t,0};
  40. Point(15) = {0.10325,0.048467570+t,0};
  41. Point(16) = {0.11650,r1+t,0};
  42.  
  43. Point(21) = {0.01325,0.032369930+B+t,0};
  44. Point(22) = {0.03575,0.036394340+B+t,0};
  45. Point(23) = {0.05825,0.040418750+B+t,0};
  46. Point(24) = {0.08075,0.044443160+B+t,0};
  47. Point(25) = {0.10325,0.048467570+B+t,0};
  48.  
  49. Line (30) = {10,11};
  50. Line (31) = {11,12};
  51. Line (32) = {12,13};
  52. Line (33) = {13,14};
  53. Line (34) = {14,15};
  54. Line (35) = {15,16};
  55.  
  56. Line (41) = {11,21};
  57. Line (42) = {12,22};
  58. Line (43) = {13,23};
  59. Line (44) = {14,24};
  60. Line (45) = {15,25};
  61.  
  62. Transfinite Line {30,35} = 5;
  63. Transfinite Line {31,32,33,34} = 8;
  64. Transfinite Line {41,42,43,44,45} = 3;
  65.  
  66. d1 = 24; // Çeyrek için düğüm sayısı
  67.  
  68. Extrude { {1,0,0} , {0,0,0} , -Pi/2} {
  69. Line{30,31,32,33,34,35,41,42,43,44,45};
  70. Layers{d1};
  71. Recombine;
  72. }
  73.  
  74. Extrude { {1,0,0} , {0,0,0} , -Pi/2} {
  75. Line{46,50,54,58,62,66,70,74,78,82,86};
  76. Layers{d1};
  77. Recombine;
  78. }
  79.  
  80. Extrude { {1,0,0} , {0,0,0} , -Pi/2} {
  81. Line{90,94,98,102,106,110,114,118,122,126,130};
  82. Layers{d1};
  83. Recombine;
  84. }
  85. Extrude { {1,0,0} , {0,0,0} , -Pi/2} {
  86. Line{134,138,142,146,150,154,158,162,166,170,174};
  87. Layers{d1};
  88. Recombine;
  89. }
  90.  
  91. // CalculiX için gerekli bâzı ayarlar
  92. Mesh.SecondOrderIncomplete=1; // 2.derece elemanları doğru ihrac etmek için gerekli
  93. Mesh.SaveGroupsOfNodes = 1; // düğüm gruplarının kaydedilmesi için
  94. Mesh.ElementOrder = 2; // ikinci derece elemanlar için 2
  95.  
  96. // CalculiX için adlandırmalar
  97. Physical Line("mesnetk") = {47,135,179,91};
  98. Physical Line("mesnetb") = {112,200,156,68};
  99. Physical Surface("cidar") = {93,97,101,105,109,113,181,185,189,193,197,201,137,141,145,149,153,157,49,53,57,61,65,69};
  100. Physical Surface("kemere") = {205,209,213,217,221,161,165,169,173,177,73,77,81,85,89,117,121,125,129,133};


Velhâsıl yukarıdaki Gmsh betiği ile [1]'de tanımlanan üç adet koni geometrisi (10, 11 ve 12) için gereken hesaplama örgüsü kolayca oluşturulabilir ve istenilen çözücüye uygun olarak ihraç edilebilir ki burada çözücü CCX olacağı için hazırlık da buna uygun olarak yapılmış ve S8 türü elamanlara elde edebilmek amacıyla düzenli örgü meydana getirilmiştir.

Bu safhada yukarıdaki önişlem verileri kullanılarak, sadece doğrusal burkulma yaklaşımı ile elde edilen sonuçların [1] üzerinde verilen deney ve hesaplama sonuçları ile karşılaştırılması aşağıdaki çizelgede görülebilir.


MPa Yaklaşım Koni 10 Koni 11 Koni 12
CCX Doğrusal 63,30 73,55 84,48
Ansys Doğrusal 62,35 74,22 83,15
Pc Deney 14,30 15,80 18,80

Çizelge.1) Calculix (2.13) ile hesaplanan, hidrostatik basınç altındaki burkulma yüklerinin (MPa) [1]'de beyan edilen ticarî Ansys yazılımı ve Deney sonuçlarıyla karşılaştırılması.


Takviyeli Koni için CCX ile hesaplanan ilk doğrusal burkulma modu

Resim.1) Takviyeli Koni 12 için CCX ile hesaplanan, hidrostatik yük altındaki ilk doğrusal burkulma modu.


Takip edecek bölümde hem takviyeli hem de takviyesiz konik kabukların, eldeki mevcut deney sonuçları ile birlikte doğrusal-olmayan burkulma davranışları ele alınmaya çalışılacak...

♦ Kaynaklar

1. Plastic Buckling of Ring-Stiffened Conical Shells Under External Hydrostatic Pressure, 2004, Carl T.F. Ross - Andrew P.F. Little - Kehinde A. Adeniyi
2. Gmsh genelağ sitesi- http://gmsh.info
3. Calculix genelağ sitesi - http://calculix.de
 







Telif Hakkı © 1997-2020 [uskudar.biz] - sürüm 5.5.1 - Bütün Hakları Saklıdır. Kullanım şartları için tıklayın!
Joomla! GNU/GPL lisansı altında özgür bir yazılımdır.