Bu çalışma, beynitik bir çeliğin ısıl işlem koşullarının dönüşüm kinetiğini, mikroyapı biçimini, dönüşüm kaynaklı iç gerilmeleri ve mekanik davranışı nasıl değiştirdiğini deneyler ile üç boyutlu faz-alanı simülasyonlarını birleştirerek incelemiştir. Fe–0,19C–1,48Si–2,38Mn ağırlık yüzdesi bileşimindeki laboratuvar çeliği, östenitlemenin ardından 673 K veya 723 K’de 45 dakika izotermal olarak bekletilmiştir. Deneylerde 673 K işlemi daha ince beynitik ferrit ve 765 ± 5 MPa akma dayanımı; 723 K işlemi ise daha kalın ferrit ve 667 ± 5 MPa akma dayanımı üretmiştir. Faz-alanı modeli bu sıcaklık eğilimini ve deneysel çekme eğrilerinin genel sıralamasını yeniden üretmiş; fakat beynitik ferrit kalınlığını her iki sıcaklıkta da fazla hesaplamıştır. Çok eksenli akma davranışı Barlat91 modeliyle kestirilmiş, ancak yalnızca x yönündeki tek eksenli çekme cevabı deneysel olarak sınanmıştır.

Model sonuçlarına göre düşük bekletme sıcaklığı daha büyük aşırı soğuma ve daha güçlü dönüşüm itici kuvveti sağlamaktadır. Isı çekme katsayısının 0,25 s⁻¹ değerinden 0,5 s⁻¹ değerine yükseltilmesi soğumayı ve dönüşümü hızlandırmakta, mekanik gevşeme için daha az zaman bırakmakta ve daha yüksek dönüşüm kaynaklı iç gerilmeler oluşturmaktadır. En büyük model akma yüzeyi 673 K ve 0,5 s⁻¹ koşulunda elde edilmiştir. Bununla birlikte 0,25 s⁻¹ kullanılan simülasyonlar deneysel çekme eğrileriyle daha yakın uyum göstermiştir. Bu nedenle en yüksek model dayanımı ile deneye en yakın ısıl yol aynı koşul değildir.

Türkiye açısından değerlendirme: Çalışmanın yaklaşımı; otomotiv sacları, demiryolu çelikleri, dişli ve rulman malzemeleri, ağır makine parçaları, aşınmaya dayanıklı çelikler ve ısıl işlem tasarımı açısından Türkiye’de uygulanabilir. Yerli çelik üreticileri ve araştırma merkezleri, farklı alaşım bileşimleri için faz-alanı modellerini dilatometre, SEM/EBSD, X-ışını kırınımı, sertlik ve çok yönlü mekanik deneylerle doğrulayarak daha az deneme üretimiyle ısıl işlem pencereleri geliştirebilir. Ancak mevcut sonuçlar tek bir laboratuvar bileşimi, iki izotermal sıcaklık ve tek eksenli çekme doğrulamasıyla sınırlıdır. Türkiye’de üretilen belirli bir çelik kalitesinin dayanımı, yorulma ömrü, kaynaklanabilirliği, üretim maliyeti veya endüstriyel fırın performansı bu çalışmadan doğrudan çıkarılamaz.

Araştırmanın temel problemi nedir?

Beynitik çelikler yüksek dayanım, kırılma tokluğu ve aşınma direncini aynı mikroyapıda birleştirebildikleri için otomotiv, demiryolu, havacılık ve ağır mühendislik uygulamalarında değerlidir. Ancak “beynit” tek biçimli bir faz değildir. Isıl işlem koşullarına bağlı olarak beynitik ferrit levhalarının:

• Kalınlığı,
• Uzunluğu,
• Yönelimi,
• Varyant dağılımı,
• Aralarındaki kalıntı östenit miktarı,
• Oluşturdukları yerel iç gerilmeler

değişebilir. Bu özellikler yalnızca tek eksenli akma dayanımını değil, farklı yönlerden uygulanan çekme, basma ve kayma yükleri altındaki plastik davranışı da etkiler.

Çalışmanın temel araştırma sorusu şöyle özetlenebilir: Bekletme sıcaklığı ve çevreye ısı aktarımının şiddeti beynitik dönüşümü nasıl değiştirir; oluşan üç boyutlu mikroyapı, iç gerilme ve kristal yönelimleri mekanik dayanım ile anizotropik akma yüzeyine nasıl yansır?

Çalışmanın önerdiği süreç–mikroyapı–özellik zinciri

Araştırmacılar aşağıdaki bağlantıyı kurmaya çalışmıştır:

Isıl işlem koşulu → sıcaklık geçmişi → östenitten beynite dönüşüm kinetiği → beynitik ferrit morfolojisi ve kalıntı östenit → dönüşüm kaynaklı iç gerilmeler → çekme davranışı ve çok eksenli akma yüzeyi.

Bu zincirin ilk bölümü faz-alanı modeliyle, mekanik bölümü kristal plastisiteyle, makroskopik anizotropik akma davranışı ise Barlat91 fenomenolojik akma ölçütüyle temsil edilmiştir.

Deneysel çelik nasıl üretilmiştir?

Yaklaşık 80 kg ağırlığında bir laboratuvar külçesi vakum indüksiyon fırınında üretilmiştir. Külçenin başlangıç kesiti 140 × 140 mm²’dir. Malzeme 1200 °C’de homojenleştirilmiş, ardından 60 × 60 mm² kesitli kütüklere dövülmüştür. İkinci homojenleştirme beş saat sürmüş ve bunu fırında soğutma izlemiştir.

Bileşim optik emisyon spektroskopisiyle, karbon miktarı ise yanma analiziyle belirlenmiştir.

Isıl işlem nasıl gerçekleştirilmiştir?

Numuneler tuz banyosunda, A_c3 sıcaklığının 60 K üzerinde 300 saniye östenitlenmiştir. Kaynak kesin A_c3 değerini vermemektedir. Tam östenitik yapı elde edildikten sonra numuneler hızla ikinci bir tuz banyosuna aktarılmıştır.

İki izotermal işlem uygulanmıştır:

• 673 K’de 45 dakika bekletme,
• 723 K’de 45 dakika bekletme.

Bekletme sonrasında numuneler oda sıcaklığına su verilerek soğutulmuştur. Her iki işlemde de beynitik ferrit esaslı mikroyapı oluşmuştur. Deneysel görüntüler, 673 K’de oluşan ferrit levhalarının 723 K’de oluşanlardan daha ince olduğunu göstermektedir.

Deneysel mikroyapı nasıl incelenmiştir?

Numune yüzeyleri 1200 kum SiC kâğıda kadar zımparalanmış, 6 ve 1 µm elmas pasta ile parlatılmış ve %3 Nital ile dağlanmıştır. İkincil elektron görüntüleri alan emisyonlu Zeiss Sigma mikroskobunda:

• 30 µm açıklık,
• 15 kV hızlandırma gerilimi,
• 9 mm çalışma mesafesi

kullanılarak alınmıştır.

Çalışmanın Şekil 5’i 673 ve 723 K numuneleri için “ters kutup figürü haritaları” sunmaktadır. Bu tür yönelim haritaları normalde elektron geri saçılım kırınımı verisine dayanır; ancak bu sürümde EBSD cihazı, adım aralığı, indeksleme oranı ve veri temizleme işlemleri açıklanmamıştır.

Faz-alanı modeli neyi temsil etmektedir?

Faz-alanı yönteminde her faz veya kristal varyant, uzay ve zamanla değişen sürekli bir alan değişkeniyle gösterilir. Keskin bir sınırı açıkça takip etmek yerine, fazlar arasındaki sınır belirli bir kalınlığa sahip geçiş bölgesi olarak çözülür.

Faz alanlarının zamanla değişimi şu çok fazlı kinetik denklemle ifade edilmiştir:

\[ \dot{\phi}_{\alpha}(\mathbf{x},t) = -\frac{1}{N} \sum_{\beta\neq\alpha} M_{\alpha\beta} \left( \frac{\delta F}{\delta\phi_{\alpha}} - \frac{\delta F}{\delta\phi_{\beta}} \right). \]

Burada \(\phi_{\alpha}\) belirli bir faz veya varyantın alan değişkenini, \(M_{\alpha\beta}\) arayüz hareketliliğini, \(F\) ise toplam serbest enerjiyi göstermektedir.

Toplam serbest enerji:

\[ F=\int_{\Omega} \left( f_{\mathrm{chem}}+ f_{\mathrm{int}}+ f_{\mathrm{el}} \right)d\Omega \]

şeklinde tanımlanmıştır. Üç temel katkı şunlardır:

• Kimyasal enerji: Faz dönüşümünü sağlayan termodinamik itici kuvvet,
• Arayüz enerjisi: Faz ve varyant sınırlarının oluşma maliyeti,
• Elastik enerji: Dönüşüm şekil değiştirmesi ve kafes uyumsuzluğundan doğan mekanik enerji.

Karbon taşınımı nasıl ele alınmıştır?

Yerel karbon konsantrasyonu bir difüzyon denklemiyle çözülmüştür:

\[ \dot{c}= \nabla\cdot \left[ \sum_{\alpha}\phi_{\alpha}D_{\alpha}\nabla c_{\alpha} + \sum_{\alpha,\beta}J_{\alpha\beta} \right]. \]

\(D_{\alpha}\) her fazın difüzyon katsayısı matrisini, \(J_{\alpha\beta}\) ise yaygın arayüz yaklaşımında oluşabilecek yapay çözünen hapsini azaltmak için kullanılan karşı-hapsetme akısını temsil eder.

Faz-alanı ve difüzyon denklemleri düzenli bir ağda sonlu fark yöntemiyle ayrıklaştırılmış, zaman integrasyonu açık ileri Euler yöntemiyle yapılmıştır.

Dönüşüm kaynaklı gerilmeler nasıl hesaplanmıştır?

Model sonlu şekil değiştirme altında elastik, plastik ve dönüşüm bileşenlerini çarpımsal olarak ayırmıştır:

\[ \mathbf{F} = \mathbf{F}^{el} \mathbf{F}^{pl} \mathbf{F}^{tr}. \]

Burada:

• \(\mathbf{F}^{el}\), elastik deformasyon gradyanı,
• \(\mathbf{F}^{pl}\), plastik deformasyon gradyanı,
• \(\mathbf{F}^{tr}\), faz dönüşümünden doğan deformasyon gradyanıdır.

Elastik şekil değiştirme Green–Lagrange ölçüsü ve St. Venant–Kirchhoff hiperelastik yaklaşımıyla hesaplanmıştır. Mekanik denge problemi hızlı Fourier dönüşümü tabanlı spektral bir çözücüyle çözülmüştür.

Kristal plastisite modeli

Plastik deformasyon, kristaldeki etkin kayma sistemlerinin toplamıyla temsil edilmiştir:

\[ \mathbf{L}^{p} = \sum_{s=1}^{N} \dot{\gamma}^{s} \mathbf{m}^{s}\otimes\mathbf{n}^{s}. \]

\(\dot{\gamma}^{s}\) ilgili kayma sisteminin kayma hızını, \(\mathbf{m}^{s}\) kayma yönünü ve \(\mathbf{n}^{s}\) kayma düzlemi normalini göstermektedir.

Kayma hızı çözülmüş kayma gerilimine bağlı bir güç yasasıyla verilmiştir:

\[ \dot{\gamma}^{s} = \dot{\gamma}_{0} \left| \frac{\tau^{s}}{\tau_{c}^{s}} \right|^{n} \mathrm{sgn}(\tau^{s}). \]

Kayma direncinin plastik deformasyonla değişimi sertleşme yasasıyla hesaplanmıştır. Çalışma, kullanılan fenomenolojik kristal plastisite modelinin doğrudan sıcaklığa bağlı olmadığını açıkça belirtmektedir.

Isı çekme ve gizli ısı nasıl modellenmiştir?

Numune sıcaklığı Newton soğuma yasasının dönüşüm gizli ısısıyla genişletilmiş biçimiyle hesaplanmıştır:

\[ \dot{T}(t) = -r\left[T(t)-T_s\right] + \frac{Q}{\rho C_p}\dot{f}. \]

Burada:

• \(T\), numune sıcaklığı,
• \(T_s\), soğutma ortamının veya banyonun sıcaklığı,
• \(r\), ısı çekme katsayısı,
• \(Q\), beynitik dönüşümün gizli ısısı,
• \(\rho\), yoğunluk,
• \(C_p\), özgül ısı,
• \(\dot{f}\), dönüşen faz hacim oranının değişim hızıdır.

Dönüşüm sırasında açığa çıkan gizli ısı, soğumayı geçici olarak yavaşlatabilir veya sıcaklıkta küçük bir yeniden yükselme oluşturabilir. Şekil 3’te yaklaşık beşinci saniye çevresindeki plato veya yükseliş bu etkiyle açıklanmıştır.

Üç boyutlu simülasyonun kurulumu

Şekil 2’de her ısıl koşul için üç boyutlu mikroyapılar gösterilmektedir. Farklı renkler 24 Kurdjumov–Sachs beynitik ferrit varyantını, siyah bölgeler ise kalıntı östeniti temsil etmektedir.

Sıcaklık eğrileri ne göstermiştir?

Dört hesaplama da 900 K’den başlamıştır. Beklendiği gibi 673 K banyo koşulları daha düşük son sıcaklığa ulaşmıştır. Aynı banyo sıcaklığında:

• \(r=0{,}5\ \mathrm{s}^{-1}\) daha hızlı başlangıç soğuması,
• \(r=0{,}25\ \mathrm{s}^{-1}\) daha uzun geçici soğuma süreci

oluşturmuştur.

Banyo sıcaklığı erişilecek son termal seviyeyi, ısı çekme katsayısı ise bu seviyeye ulaşma hızını belirleyen ana değişken olarak yorumlanmıştır.

Beynitik ferrit hacim oranı

Şekil 4’te beynitik ferrit hacim oranı sıcaklık azaldıkça yükselmiştir. 673 K koşulları, simülasyon süresi içinde 723 K koşullarından daha yüksek son beynitik ferrit oranına ulaşmıştır.

Isı çekme katsayısı son faz oranını sıcaklık kadar güçlü değiştirmemiş; fakat dönüşümün ne kadar hızlı gerçekleştiğini ve ortaya çıkan morfolojiyi belirgin biçimde etkilemiştir. 723 K ve 0,25 s⁻¹ koşulu en düşük son beynitik ferrit oranını vermiştir.

Deneysel ve hesaplanan ferrit kalınlığı

Deney ve simülasyon aynı eğilimi göstermiştir: bekletme sıcaklığı arttıkça beynitik ferrit kalınlaşmaktadır. Ancak hesaplanan değerler mutlak olarak büyüktür. Fark özellikle 673 K’de belirgindir.

Araştırmacılar bu sapmayı model çözünürlüğü ve görüntü tabanlı kalınlık belirleme yönteminin sınırlılığıyla ilişkilendirmiştir. Simülasyon istatistiklerinde ferrit kalınlığı yaklaşık 6–7 piksel ile temsil edildiğinden küçük farklılıklar ağ belirsizliği içinde kalabilmektedir.

Alan, uzunluk ve biçim istatistikleri

Şekil 6’da alan, en-boy oranı, uzunluk ve kalınlık dağılımları karşılaştırılmıştır. Kaynaktaki yorumlara göre:

• 673 K ve 0,5 s⁻¹, en büyük alan ve en uzun beynitik özellikleri üretmiştir.
• 673 K koşulları daha uzamış özelliklerle ilişkilendirilmiştir.
• 723 K koşulları daha kısa ve daha kompakt özellikler göstermiştir.
• 723 K ve 0,25 s⁻¹ koşulunda uzun termal maruziyet yanal büyüme ve kabalaşmaya daha fazla zaman vermiştir.
• Kalınlık dağılımlarındaki farklar ağ çözünürlüğü nedeniyle güvenilir biçimde ayrılamamıştır.

Burada “daha büyük alan veya uzunluk” ile “daha kalın ferrit” aynı ölçüt değildir. 673 K ve hızlı ısı çekme koşulunda özellikler daha uzun ve uzamış olabilirken levha kalınlığının fiziksel olarak daha ince olması beklenmektedir.

Alan ve uzunluk dağılımları yalnızca simülasyon görüntülerinden çıkarılmıştır. Bunların karşılaştırılabileceği deneysel alan ve uzunluk dağılımları çalışmada bulunmamaktadır.

Dönüşüm kaynaklı iç gerilmeler

Şekil 7, dört ısıl işlem koşulu sonrasındaki yerel von Mises gerilme dağılımlarını göstermektedir. Gerilmeler mikroyapı içinde homojen değildir; belirli ferrit levhaları, varyant sınırları ve mekanik uyumsuzluk bölgelerinde yoğunlaşmaktadır.

En yüksek yerel gerilmeler 673 K ve 0,5 s⁻¹ koşulunda hesaplanmıştır. Hızlı dönüşüm, dönüşüm şekil değiştirmesini kısa sürede meydana getirmekte ve plastik gevşemeye daha az zaman bırakmaktadır. Daha yavaş ısı çekme ise mekanik uyum için daha uzun süre sağlamaktadır.

Şekildeki ölçek yaklaşık 3,3 × 10⁷ ile 1,3 × 10⁹ Pa arasındadır. Bu değerler makroskopik çekme gerilimi değil, hesaplama hacmi içindeki yerel ve kısıtlanmış gerilmelerdir. Model çatlama, hasar ve bütün gerilme boşaltma mekanizmalarını içermediğinden mutlak tepe değerler doğrudan gerçek parçadaki ölçülebilir gerilme olarak yorumlanmamalıdır.

Çekme deneyleri

Tek eksenli çekme deneyleri B5×25 silindirik numunelerde gerçekleştirilmiştir:

• Ölçü çapı: 5 mm,
• Ölçü boyu: 25 mm,
• Deney sıcaklığı: oda sıcaklığı,
• Şekil değiştirme hızı: 0,001 s⁻¹,
• Cihaz: ZwickRoell Z100 elektromekanik üniversal test makinesi.

673 K işlemi, 723 K’ye göre yaklaşık 98 MPa daha yüksek akma dayanımı sağlamıştır. Araştırmacılar bu farkı daha yüksek beynitik ferrit oranı, daha ince ferrit yapısı ve düşük sıcaklığın oluşturduğu daha güçlü dönüşüm itici kuvvetiyle ilişkilendirmiştir.

Deney tekrar sayısı, ±5 MPa değerinin standart sapma mı başka bir belirsizlik ölçüsü mü olduğu ve akma dayanımının hangi ofset yöntemiyle belirlendiği kaynakta açıklanmamıştır.

Deneysel ve hesaplanan çekme eğrileri

Şekil 8, deneysel eğrilerle dört simülasyon eğrisini karşılaştırmaktadır. Model her iki ısı çekme katsayısında da 673 K numunesinin 723 K numunesinden daha yüksek gerilme taşıdığını yeniden üretmiştir.

\(r=0{,}25\ \mathrm{s}^{-1}\) kullanılan eğriler deneylere en yakın sonuçları vermiştir. Buna karşılık \(r=0{,}5\ \mathrm{s}^{-1}\) kullanılan simülasyonlar daha yüksek dayanım öngörmüştür. Kaynak bunu:

• Daha hızlı dönüşüm,
• Daha ince ferrit oluşumu,
• Plastik akışa karşı daha fazla mikroyapısal engel,
• Daha yüksek kalıntı iç gerilme

ile açıklamaktadır.

Barlat91 akma yüzeyi neyi göstermektedir?

Tek eksenli çekmede akma tek bir gerilme değeriyle ifade edilebilir. Çok eksenli gerilme durumlarında ise elastik ve plastik bölgeleri ayıran bir akma yüzeyi gerekir.

Çalışmada altı gerilme bileşenini ve anizotropik davranışı temsil edebilen Barlat91 ölçütü kullanılmıştır. Modelin genel biçimi:

\[ \Phi=2\sigma_y^n \]

eşitliğiyle tanımlanmış; ikinci ve üçüncü gerilme değişmezleri ile altı anizotropi parametresi kullanılmıştır. Kaynak BCC malzemeler için \(n=6\), FCC malzemeler için \(n=8\) değerini belirtmektedir. Çok fazlı beynitik ferrit–kalıntı östenit yapısı için nihai Barlat91 uyumunda hangi \(n\) değerinin kullanıldığı açıkça yazılmamıştır.

Önce x yönündeki çekme davranışı deneysel eğriye göre değerlendirilmiştir. Daha sonra aynı kristal plastisite parametreleriyle çekme–basma ve kayma ağırlıklı ek yük yolları sayısal olarak uygulanmıştır. Bu sayısal akma noktaları Barlat91 yüzeyinin belirlenmesinde kullanılmıştır.

Akma yüzeylerinin temel sonucu

Şekil 11’deki yüzeyler dairesel değildir. Bu durum, mikroyapının izotropik von Mises ölçütüyle tam olarak temsil edilemeyecek yön bağımlı bir plastik tepki verdiğini göstermektedir.

Akma yüzeyi büyüklüğünün sıralaması genel olarak şöyledir:

1. 673 K, 0,5 s⁻¹: en büyük yüzey ve en yüksek model akma direnci,
2. 723 K, 0,5 s⁻¹,
3. 673 K, 0,25 s⁻¹,
4. 723 K, 0,25 s⁻¹: en küçük yüzey.

Düşük sıcaklık ve hızlı ısı çekme, daha yüksek dönüşüm itici kuvveti, daha hızlı mikroyapı oluşumu ve daha yüksek iç gerilmeler yoluyla akma yüzeyini dışarı doğru genişletmiştir.

Neden akma davranışı anizotropiktir?

Çalışma anizotropiyi dört mikroyapısal etkene bağlamaktadır:

• Beynitik ferrit levhalarının yönlü ve uzamış biçimi,
• 24 Kurdjumov–Sachs varyantının eşit olmayan hacim oranları,
• Varyantlara bağlı dönüşüm şekil değiştirmeleri,
• Mikroyapı içindeki heterojen kalıntı gerilmeler.

Yük ferrit levhalarına paralel, dik veya kayma yönünde uygulandığında farklı kayma sistemleri ve faz sınırı kısıtları etkinleşebilir. Bu nedenle akma direnci yalnızca toplam ferrit oranına değil, yükleme yönüne de bağlıdır.

Kayma ağırlıklı yüklemede sıcaklık etkisi neden azalmıştır?

Akma yüzeyinin bazı kayma ağırlıklı bölgelerinde 673 ve 723 K eğrileri birbirine yaklaşmıştır. Araştırmacılar bunu, kayma yüklemesinin daha geniş bir kristal yönelimi grubunu etkinleştirmesi ve tek bir baskın morfolojik yönün etkisini azaltmasıyla açıklamaktadır.

Özellikle 0,5 s⁻¹ koşullarında 2, 4, 9, 10, 13, 15, 20 ve 23 numaralı bazı Kurdjumov–Sachs varyantlarının benzer hacim oranlarına sahip olması, kayma yönündeki toplu cevabın birbirine yaklaşmasına katkı sağlayabilir. Bu açıklama model tabanlı bir yorumdur ve çok eksenli deneyle sınanmamıştır.

Çalışmanın güçlü yönleri

• Laboratuvar çeliği üretilmiş ve kimyasal bileşimi ölçülmüştür.
• İki kontrollü izotermal beynit işlemi uygulanmıştır.
• Mikroyapı ve tek eksenli çekme davranışı deneysel olarak incelenmiştir.
• Faz dönüşümü, karbon difüzyonu, gizli ısı ve mekanik etkileşim aynı çerçevede birleştirilmiştir.
• Üç boyutlu mikroyapıda 24 kristal varyantı ve kalıntı östenit temsil edilmiştir.
• Simülasyonun ferrit kalınlığını fazla hesapladığı açıkça raporlanmıştır.
• Yerel iç gerilmelerin makroskopik çekme gerilimi olmadığı belirtilmiştir.
• Tek eksenli doğrulamadan çok eksenli tahmine geçişin model tabanlı olduğu açıklanmıştır.
• Isıl işlem, mikroyapı ve akma yüzeyi arasında bütünleşik bir bağlantı kurulmuştur.

Çalışmanın sınırlılıkları

• Çalışma hakem değerlendirmesinden geçmemiştir.
• Yalnızca tek bir çelik bileşimi incelenmiştir.
• Yalnızca 673 ve 723 K izotermal bekletme sıcaklıkları deneysel olarak karşılaştırılmıştır.
• Simülasyon beynitik ferrit kalınlığını nicel olarak fazla hesaplamıştır.
• Simülasyon alan ve uzunluk dağılımları deneysel olarak doğrulanmamıştır.
• İç gerilmeler X-ışını, nötron kırınımı veya mikrogerilme yöntemiyle ölçülmemiştir.
• Hasar, çatlama ve kırılma mekanizmaları modele dâhil edilmemiştir.
• Yalnızca x yönündeki tek eksenli çekme cevabı deneysel olarak değerlendirilmiştir.
• Çok eksenli akma yüzeyleri ve kayma cevabı deneysel olarak doğrulanmamıştır.
• Çekme deneyi tekrar sayısı ve belirsizlik hesabı açıklanmamıştır.
• EBSD veya ters kutup figürü haritası üretim ayrıntıları verilmemiştir.
• Barlat91 parametrelerinin tam sayısal değerleri ve ters tanımlama ayrıntıları sunulmamıştır.
• Tablo 5’te gerilme ve sertleşme parametrelerinin birimleri tutarsız görünmektedir.
• Faz-alanı hesaplarının zaman adımı, toplam hesaplama süresi, donanımı ve yakınsama incelemesi bildirilmemiştir.
• Veri ve kod için açık depo bağlantısı sunulmamıştır.

Çalışma neyi desteklemektedir?

• Düşük izotermal sıcaklık bu çelikte daha ince beynitik ferrit ve daha yüksek akma dayanımıyla ilişkilidir.
• Faz-alanı modeli sıcaklığa bağlı ferrit kalınlığı eğilimini nitel olarak yakalayabilir.
• Isı çekme katsayısı dönüşüm hızını ve iç gerilme gelişimini etkileyebilir.
• Hızlı dönüşüm daha yüksek yerel dönüşüm gerilmeleri oluşturabilir.
• 673 K işlemi 723 K işleminden daha yüksek deneysel akma dayanımı sağlamıştır.
• Faz-alanı/kristal plastisite modeli deneysel tek eksenli dayanım sıralamasını yeniden üretmiştir.
• Simüle edilen beynitik mikroyapı yön bağımlı akma davranışı göstermektedir.
• Barlat91, modelden elde edilen anizotropik akma noktalarını makroskopik bir yüzeyle temsil edebilir.

Çalışma neyi kanıtlamamaktadır?

• 673 K ve 0,5 s⁻¹ koşulunun gerçek üretimde her açıdan en iyi ısıl işlem olduğunu kanıtlamamaktadır.
• Hesaplanan yerel iç gerilmelerin gerçek çelikte aynı büyüklükte olduğunu göstermemektedir.
• Çok eksenli akma yüzeylerini deneysel olarak doğrulamamaktadır.
• Modelin başka çelik bileşimlerine doğrudan aktarılabileceğini kanıtlamamaktadır.
• Daha yüksek çekme dayanımının daha yüksek tokluk, yorulma ömrü veya aşınma direnci sağlayacağını göstermemektedir.
• Kalıntı östenitin deformasyon sırasında martenzite dönüşümünü deneysel olarak izlememektedir.
• Endüstriyel ölçekteki sıcaklık gradyanlarını, parça geometrisini veya fırın değişkenliğini çözmemektedir.
• Barlat91 modelinin von Mises, Hill48 veya başka ölçütlerden deneysel olarak üstün olduğunu kanıtlamamaktadır.

Mühendislik açısından anlamı

Çalışmanın en önemli katkısı, ısıl işlem sıcaklığını yalnızca nihai sertlik veya çekme dayanımıyla ilişkilendirmek yerine, aradaki mikroyapısal mekanizmaları çözmeye çalışmasıdır. Böyle bir çerçeve geliştirilebildiğinde üretici, her olası ısıl işlemi fiziksel olarak denemeden önce:

• Beynit dönüşüm hızını,
• Ferrit kalınlığını ve yönelimini,
• Kalıntı östenit dağılımını,
• Yerel iç gerilme bölgelerini,
• Farklı yük yönlerindeki akma direncini

sayısal olarak tarayabilir.

Ancak bu çalışma henüz tam anlamıyla endüstriyel bir “sanal malzeme laboratuvarı” değildir. Mutlak mikroyapı boyutlarının, çok eksenli mekanik cevabın ve yerel gerilmelerin daha kapsamlı deneylerle doğrulanması gerekmektedir.

Çalışmanın Yöntemi ve Bulguları

Araştırma tasarımının teknik özeti

Ana nicel bulgular

• 673 K’de deneysel beynitik ferrit kalınlığı 0,15 ± 0,05 µm’dir.
• 723 K’de deneysel beynitik ferrit kalınlığı 0,26 ± 0,08 µm’dir.
• 673 K faz-alanı sonucu 0,341 µm’dir.
• 723 K faz-alanı sonucu 0,382 µm’dir.
• 673 K’de deneysel akma dayanımı 765 ± 5 MPa’dır.
• 723 K’de deneysel akma dayanımı 667 ± 5 MPa’dır.
• Düşük sıcaklık işlemi yaklaşık 98 MPa daha yüksek deneysel akma dayanımı sağlamıştır.
• 0,25 s⁻¹ kullanılan simülasyonlar deneysel çekme eğrilerine en yakın sonuçları vermiştir.
• 0,5 s⁻¹ kullanılan simülasyonlar daha yüksek model dayanımı ve daha yüksek iç gerilme üretmiştir.
• En büyük hesaplanan akma yüzeyi 673 K ve 0,5 s⁻¹ koşuluna aittir.
• İç gerilme haritalarındaki üst ölçek yaklaşık 1,3 GPa’dır; bu değer yerel model gerilmesidir.

Görsellerin bilimsel işlevi

Kaynak içindeki teknik tutarsızlıklar

• Tablo 5’te gerilme parametrelerinin birimleri, verilen sayısal değerlerle birlikte fiziksel olarak olağan dışı büyüklükler üretmektedir.
• Sertleşme denklemlerinde \(a\) üssü kullanılırken Tablo 5 “Hardening Index m” vermektedir; iki sembol arasındaki ilişki açıklanmamıştır.
• Barlat91 için BCC ve FCC’ye ayrı \(n\) değerleri verilmiş, ancak iki fazlı mikroyapının makroskopik uyumunda seçilen değer bildirilmemiştir.
• Ters kutup figürü haritaları için EBSD deney koşulları sunulmamıştır.
• Başlıkta deneysel doğrulama vurgusu bulunmasına rağmen çok eksenli akma davranışı deneysel değil, model tabanlıdır.

Bulguların dikkatli yorumu

Çalışma, ısıl işlem sıcaklığının deneysel mikroyapı kalınlığı ve tek eksenli akma dayanımı üzerindeki etkisini güvenilir bir eğilim düzeyinde göstermektedir. Faz-alanı modeli bu eğilimi yeniden üretmektedir; fakat mutlak ferrit kalınlığı açısından henüz nicel doğruluğa ulaşmamıştır.

İç gerilme ve çok eksenli akma yüzeyi sonuçları, fiziksel mekanizmaları araştırmak ve yeni deneyler tasarlamak açısından değerlidir. Ancak doğrudan ölçülmedikleri için doğrulanmış malzeme verisi değil, deneyle sınanması gereken model kestirimleri olarak ele alınmalıdır.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: Experimentally validated process–microstructure–property relations of bainitic steels derived from phase-field simulations

Yazarlar: Dhanunjaya Kumar Nerella, Muhammad Adil Ali, Oguz Gulbay, Oleg Shchyglo ve Ingo Steinbach.

Yazar sıralaması: Kaynakta verildiği biçimde korunmuştur.

Eş birinci yazar: Eşit katkı veya eş birinci yazarlık bilgisi yer almamaktadır.

Sorumlu yazar: Dhanunjaya Kumar Nerella.

İletişim adresi: Yüklenen metinde e-posta adresi yer almamaktadır. Ruhr Üniversitesi Bochum’un resmî araştırmacı kaydı `dhanunjaya.nerella@rub.de` adresini vermektedir.

Kurum 1: Interdisciplinary Centre for Advanced Materials Simulation, Ruhr University Bochum, Universitätsstraße 150, 44801 Bochum, Kuzey Ren-Vestfalya, Almanya.

Kurum 2: Steel Institute, RWTH Aachen University, Intzestraße 1, 52072 Aachen, Kuzey Ren-Vestfalya, Almanya.

DOI:10.2139/ssrn.7176521

Resmî kaynak bağlantısı:SSRN resmî kayıt sayfası

Yayın platformu: SSRN.

Gönderim tarihi: 3 Temmuz 2026.

Yayın yılı: 2026.

Dergi: Belirli bir dergi adı veya kabul bilgisi bu sürümde yer almamaktadır.

Yayınevi: Metin “Preprint submitted to Elsevier” ifadesini taşımaktadır; ancak belirli bir Elsevier dergisi veya kabul edilmiş yayın kaydı belirtilmemiştir.

Kaynak türü: Deneysel ısıl işlem, mikroyapı karakterizasyonu ve çekme testlerini üç boyutlu faz-alanı, kristal plastisite ve fenomenolojik akma yüzeyi modellemesiyle birleştiren preprint araştırma makalesi.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir.

Yazar katkıları: Bu sürümde ayrı bir CRediT veya görev bazlı yazar katkı beyanı bulunmamaktadır.

Finansman: Araştırma Almanya Federal Ministry of Research, Technology and Space tarafından DiStEL projesi, proje numarası 13XP5226E kapsamında desteklenmiştir.

Çıkar çatışması: Ayrı bir çıkar çatışması veya rekabet eden menfaat beyanı bu sürümde yer almamaktadır.

Veri erişimi: Faz-alanı simülasyon verilerinin makul talep üzerine yazarlardan sağlanabileceği belirtilmiştir. Açık veri deposu bağlantısı verilmemiştir.

Kod erişimi: Kaynak kodu, kullanılan yazılım sürümü, çalışma betikleri veya yeniden üretim paketi paylaşılmamıştır.

Deneysel doğrulama sınırı: Deneysel karşılaştırma mikroyapı kalınlığı eğilimi ve x yönündeki tek eksenli çekme cevabıyla sınırlıdır. İç gerilmeler, diğer yük yönleri ve Barlat91 çok eksenli yüzeyleri model tabanlı kestirimlerdir.

Morfoloji doğrulama sınırı: Faz-alanı modeli sıcaklığa bağlı kalınlaşma eğilimini yakalamış; ancak ferrit kalınlığını özellikle 673 K’de fazla hesaplamıştır. Simüle edilen alan ve uzunluk dağılımları deneysel verilerle karşılaştırılmamıştır.

Kaynak tutarlılığı uyarısı: Kristal plastisite parametre tablosundaki gerilme ve sertleşme birimleri açıklama gerektiren büyüklükler vermektedir. Barlat91 üs seçimi ve EBSD ölçüm koşulları da açık değildir. Bu noktalar kaynak sessizce düzeltilmeden belirtilmiştir.

Bu Türkçe içerik, yüklenen çalışmanın deneysel yöntemleri, denklemleri, tabloları, mikroyapı görüntüleri, gerilme haritaları ve mekanik sonuçları temel alınarak hazırlanmıştır. Çalışmada bulunmayan bir endüstriyel üretim başarısı, yorulma performansı, tokluk artışı, parça güvenliği veya Türkiye’ye özgü mekanik özellik iddiası eklenmemiştir.