Kilitlenme ve Kum Saati: Ayrıcalığın Bittiği Yer

1. Neden ayrı bir motor?

FEM-01'in temel sonucu şuydu: kiriş elemanının kübik şekil fonksiyonları, eleman içinde yayılı yük yokken Euler–Bernoulli denkleminin tam çözümüdür. Bu yüzden ucundan yüklü bir konsol tek elemanla tam çözülür.

Sürekli ortamda — düzlem gerilme, düzlem şekil değiştirme, katı — bu ayrıcalık yoktur. Şekil fonksiyonları gerçek yerdeğiştirme alanını üretemez; yalnızca ona yaklaşır. Dolayısıyla ağ sıklaştırmak burada gerçekten gereklidir.

Ama mesele bundan ibaret değil. Yaklaşımın ne kadar kötü olabileceği, çoğu kişinin beklediğinden büyük.

2. Bir elemanın geçerliliği nasıl sınanır?

Yeni bir eleman formülasyonunun doğru olup olmadığı, bir ticari programın çıktısıyla karşılaştırılarak değil, elemanın kendi matematiksel özellikleriyle sınanır. İki temel ölçüt vardır:

• Rijit cisim hareketi. Eleman ötelenip döndürüldüğünde hiçbir şekil değiştirme üretmemelidir. Üretiyorsa formülasyon yanlıştır ve o eleman sahte iç kuvvet doğurur.
• Yama sınaması (patch test). Sabit şekil değiştirme alanı, keyfi bir ağla tam olarak üretilebilmelidir. Bu, elemanın yakınsayacağının güvencesidir.

Motorumuz her ikisini de makine hassasiyetinde geçiyor; sonuçlar test paketinde sabitlenmiştir.

3. Yama sınaması

Kurulum şöyledir: bir ağın bütün sınır düğümlerine, sabit tek eksenli gerilme durumuna karşılık gelen yerdeğiştirmeler atanır. İç düğümler serbest bırakılır. Eleman doğruysa çözülen iç gerilme tam olarak sabit çıkar.

4. Kilitlenme: CST neden bu kadar rijit?

Sabit şekil değiştirmeli üçgen (CST), adı üstünde, eleman içinde şekil değiştirmeyi sabit kabul eder. Yerdeğiştirme alanı doğrusaldır, türevi sabittir.

Şimdi eğilmeyi düşünün. Eğilen bir kirişte şekil değiştirme yükseklik boyunca doğrusal olarak değişir: üstte basınç, altta çekme, ortada sıfır. Sabit şekil değiştirmeli bir eleman bunu tanımı gereği temsil edemez.

5. Sayısal karşılaştırma

Konsol kiriş, düzlem gerilme. L = 1000, H = 100, t = 10 mm; uç yükü 1000 N. Referans, kesme katkısı dâhil Timoshenko çözümü: —.

FEM-01'deki kiriş elemanıyla karşılaştırın: orada tek eleman tam sonucu veriyordu. Burada dört CST elemanı %90 hata yapıyor. İki durum arasındaki fark, şekil fonksiyonlarının çözüm uzayını içerip içermemesinde.

6. Yakınsama daima aşağıdan

Tablodaki bütün değerler referansın altında. Bu bir rastlantı değil, varyasyonel bir sonuçtur.

Pratik değeri şudur: uyumlu elemanlarla hesaplanan yerdeğiştirme bir alt sınırdır. Ağınız kabaysa gerçek sehim hesapladığınızdan büyüktür — bu, sehim sınırı kontrolünde emniyetsiz tarafta olmak demektir. Gerilme için aynı basit ilişki geçerli değildir; orada yakınsama tek yönlü olmayabilir.

7. İlaç: indirgenmiş integrasyon

Kilitlenmenin kaynağı, sahte kayma enerjisinin sayısal integrasyonda fazla temsil edilmesidir. Q4 elemanında rijitlik matrisi normalde 2×2 Gauss noktasıyla integre edilir. Kayma terimini tek noktada (merkezde) integre etmek, o sahte enerjinin çoğunu hesaba katmaz.

Tablodaki üçüncü sütun bunun sonucudur: aynı ağda indirgenmiş integrasyon, tam integrasyondan belirgin biçimde doğru. 64×16 ağda hata —, tam integrasyonda ise —.

8. İlacın yan etkisi: kum saati modu

Şimdi tablonun ilk satırındaki üçüncü hücreye bakın.

Sebep şudur: rijitlik matrisi, integrasyon noktalarındaki şekil değiştirmelerden kurulur. Tek noktada integre edilen bir Q4 elemanında, merkez noktasında hiç şekil değiştirme üretmeyen ama sıfır olmayan yerdeğiştirme alanları vardır. Bu alanlar elemana hiç enerji harcatmaz — yani sahte sıfır enerjili modlardır. Her indirgenmiş Q4 elemanının böyle iki modu vardır.

Yeterince kısıtlanmış bir ağda komşu elemanlar birbirini engeller ve bu modlar uyanmaz. Kaba ve tek sıralı bir ağda ise engelleyecek kimse yoktur; sistem matrisi neredeyse tekil hâle gelir ve çözüm patlar.

9. Yakınsama mertebesini ölçmek

Bir modelin doğru yakınsayıp yakınsamadığını anlamak için tam çözümü bilmeniz gerekmez. Ardışık ağlardaki hatanın nasıl küçüldüğüne bakmak yeter — bu, FEM-01'de de kullandığımız araç.

Gözlenen mertebenin kararlı bir değere oturması, modelin düzgün yakınsadığının işaretidir. Oturmuyorsa — mertebe ağdan ağa zıplıyorsa — modelde bir sorun vardır: tekillik, bozuk eleman ya da yanlış tanımlanmış sınır koşulu.

10. Pratikte ne yapmalı

1. Eğilme baskınsa CST kullanmayın. Üçgen eleman ağ üretmesi kolay olduğu için caziptir; eğilmede kabul edilemez derecede rijittir. En azından ikinci dereceden (altı düğümlü) üçgen kullanın.
2. Kalınlık yönünde tek eleman kullanmayın. Tablodaki en kötü satır 4×1'dir. Eğilme gradyanını temsil edecek kadar eleman gerekir.
3. İndirgenmiş integrasyon kullanıyorsanız kum saatini kontrol edin. Deforme olmuş şekle bakın: komşu elemanlar zıt yönlerde "zikzak" yapıyorsa mod uyanmıştır.
4. Yerdeğiştirmenin alt sınır olduğunu unutmayın. Kaba ağda gerçek sehim hesapladığınızdan büyüktür.
5. Yakınsamayı ölçün, varsaymayın. En az üç ağ çözün ve hatanın nasıl küçüldüğüne bakın.
6. Yeni bir programa geçtiğinizde yama sınamasını koşturun. Sabit gerilme alanını tam üretemeyen bir eleman, hiçbir ağda güvenilir değildir.

11. Kendini sınama

1. CST elemanı eğilmede neden bu kadar rijittir? "Sabit şekil değiştirme" ifadesi cevabın neresinde?
2. Yerdeğiştirme temelli bir modelin gerçek yapıdan daima daha rijit olması hangi ilkeden çıkar?
3. İndirgenmiş integrasyon kilitlenmeyi neden azaltır? Hangi enerji terimi eksik hesaplanır?
4. Kum saati modu neden kaba ağda uyanır da sık ağda uyanmaz?
5. Bir ağı iki katına çıkardınız ve hata dörtte bire indi. Gözlenen mertebe kaçtır?
6. Yama sınamasını geçen bir eleman, her problemde doğru sonuç verir mi?