Yapısal mekanikte sonlu eleman yöntemi - Finite element method in structural mechanics

sonlu eleman yöntemi (FEM), başlangıçta karmaşık problemlerin sayısal çözümü için geliştirilmiş güçlü bir tekniktir. yapısal mekanik ve karmaşık sistemler için tercih edilen yöntem olmaya devam etmektedir. FEM'de, yapısal sistem bir dizi uygun sonlu elemanlar düğüm adı verilen ayrı noktalarda birbirine bağlı. Elemanların kalınlık gibi fiziksel özellikleri olabilir, termal Genleşme katsayısı, yoğunluk, Gencin modülü, kayma modülü ve Poisson oranı.

Tarih

Sonlu yöntemin kökeni, yapıların matris analizine kadar izlenebilir. ^[1]^[2] yer değiştirme veya sertlik matrisi yaklaşımı kavramının tanıtıldığı yer. 1950'lerde mühendislik yöntemlerine dayalı olarak sonlu eleman kavramları geliştirilmiştir. Sonlu elemanlar yöntemi gerçek ivmesini 1960'larda ve 1970'lerde John Argyris ve iş arkadaşları; -de Stuttgart Üniversitesi, tarafından Ray W. Clough; -de California Üniversitesi, Berkeley, tarafından Olgierd Zienkiewicz ve iş arkadaşları Ernest Hinton, Bruce Irons;^[3] -de Swansea Üniversitesi, tarafından Philippe G. Ciarlet; -de Paris Üniversitesi; -de Cornell Üniversitesi Richard Gallagher ve meslektaşları tarafından. Argyris gibi orijinal eserler ^[4] ve Clough ^[5] bugünün sonlu elemanlar yapısal analiz yöntemlerinin temeli oldu.

Eksenel, eğilme ve burulma sertlikleri gibi fiziksel özelliklere sahip düz veya eğimli tek boyutlu elemanlar. Bu tip eleman, kabloları, köşebentleri, kafes kirişleri, takviyeleri, ızgaraları ve çerçeveleri modellemek için uygundur. Düz elemanlar genellikle her iki uçta bir tane olmak üzere iki düğüme sahiptir, oysa eğimli elemanlar uç düğümler dahil en az üç düğüme ihtiyaç duyar. Elemanlar, merkez asıl üyelerin ekseni.

Membran etkisiyle yalnızca düzlem içi kuvvetlere direnç gösteren iki boyutlu elemanlar (düzlem stres, uçak Gerginlik ) ve enine kesme ve eğilme hareketi ile enine yüklere direnen plakalar (plakalar ve kabuklar ). Düz veya kavisli gibi çeşitli şekillere sahip olabilirler. üçgenler ve dörtgenler. Düğümler genellikle eleman köşelerine yerleştirilir ve daha yüksek doğruluk için gerekirse, eleman kenarları boyunca ve hatta elemanın içine ilave düğümler yerleştirilebilir. Elemanlar, gerçek katman kalınlığının orta yüzeyinde konumlandırılır.
Torus membranlar, kalın plakalar, kabuklar ve katılar gibi eksenel simetrik problemler için-şekilli elemanlar. Elemanların enine kesiti, daha önce açıklanan türlere benzer: ince plakalar ve kabuklar için tek boyutlu ve katılar, kalın plakalar ve kabuklar için iki boyutlu.
3 boyutlu katıları modellemek için üç boyutlu elemanlar makine bileşenler barajlar, setler veya toprak kütleleri. Ortak eleman şekilleri şunları içerir: tetrahedrallar ve altı yüzlüler. Düğümler, tepe noktalarına ve muhtemelen eleman yüzlerine veya elemanın içine yerleştirilir.

Eleman ara bağlantısı ve yer değiştirme

Öğeler yalnızca dış düğümlerde birbirine bağlıdır ve hepsi birlikte tüm alanı mümkün olduğunca doğru bir şekilde kapsamalıdır. Düğümler düğüme sahip olacak (vektör) yer değiştirmeler veya özgürlük derecesi çeviriler, rotasyonlar ve özel uygulamalar için üst düzey türevler yer değiştirmeler. Düğümler yer değiştirdiğinde, sürüklemek eleman formülasyonu tarafından belirlenen belirli bir şekilde birlikte elemanlar. Başka bir deyişle, elemandaki herhangi bir noktanın yer değiştirmesi olacaktır. enterpolasyonlu düğüm yer değiştirmelerinden ve bu, çözümün yaklaşık doğasının ana nedenidir.

Pratik hususlar

Uygulama açısından, sistemi şu şekilde modellemek önemlidir:

Modelin boyutunu küçültmek için simetri veya anti-simetri koşullarından yararlanılır.
Herhangi bir gerekli süreksizlik dahil olmak üzere yer değiştirme uyumluluğu, özellikle bitişik elemanlar farklı tipte, malzemede veya kalınlıkta olduğunda, düğümlerde ve tercihen eleman kenarları boyunca sağlanır. Birçok düğümün yer değiştirmelerinin uyumluluğu genellikle kısıtlama ilişkileri yoluyla empoze edilebilir.
Öğelerin davranışları, hem yerel hem de küresel olarak gerçek sistemin baskın eylemlerini yakalamalıdır.
Kabul edilebilir doğruluk sağlamak için eleman ağı yeterince ince olmalıdır. Doğruluğu değerlendirmek için mesh, önemli sonuçlar çok az değişiklik gösterene kadar rafine edilir. Daha yüksek doğruluk için en boy oranı elemanların mümkün olduğunca birliğe yakın olması gerekir ve daha yüksek gerilimli kısımlarda daha küçük elemanlar kullanılır gradyan.
Simetri eksenlerindeki düğümlere özel dikkat gösterilerek uygun destek kısıtlamaları uygulanır.

Büyük ölçekli ticari yazılım paketleri, genellikle ağı oluşturmak için kolaylıklar ve hem girdi verilerinin doğrulanmasını hem de sonuçların yorumlanmasını büyük ölçüde kolaylaştıran girdi ve çıktının grafiksel gösterimini sağlar.

FEM Yer Değiştirme Formülasyonuna teorik genel bakış: Elemanlardan sisteme ve çözüme

ZEE teorisi farklı perspektiflerde veya vurgularda sunulabilirken, gelişimi yapısal Analiz yoluyla daha geleneksel yaklaşımı izler sanal çalışma ilke veya minimum toplam potansiyel enerji prensibi. sanal çalışma ilke yaklaşımı, hem doğrusal hem de doğrusal olmayan malzeme davranışlarına uygulanabildiğinden daha geneldir. Sanal çalışma yöntemi, enerjinin korunumu: konservatif sistemler için, sisteme uygulanan bir dizi kuvvet tarafından eklenen iş, yapının bileşenlerinin gerinim enerjisi biçiminde sistemde depolanan enerjiye eşittir.

Prensibi sanal yer değiştirmeler yapısal sistem için dış ve iç sanal çalışmanın matematiksel kimliğini ifade eder:

{displaystyle {mbox {Harici sanal çalışma}} = int _ {V} delta {oldsymbol {epsilon}} ^ {T} {oldsymbol {sigma}}, dVqquad mathrm {(1)}}

Diğer bir deyişle, dış kuvvetler kümesi tarafından sistem üzerinde yapılan işin toplamı, sistemi oluşturan elemanlarda gerinim enerjisi olarak depolanan işe eşittir.

Yukarıdaki denklemin sağ tarafındaki sanal dahili çalışma, tek tek öğeler üzerinde yapılan sanal çalışmanın toplanmasıyla bulunabilir. İkincisi, her bir elemanın tepkisini tanımlayan kuvvet yer değiştirme fonksiyonlarının kullanılmasını gerektirir. Bu nedenle, yapının yer değiştirmesi, bireysel (ayrı) elemanların toplu olarak tepkisi ile tanımlanır. Denklemler, sistemin yanıtını bir bütün olarak (bir süreklilik) tanımlayan tek bir denklem yerine, yapının tek tek öğelerinin küçük alanı için yazılmıştır. İkincisi, çözülemez bir soruna, dolayısıyla sonlu elemanlar yönteminin faydasına neden olacaktır. Sonraki bölümlerde gösterildiği gibi, Denklem (1) sistem için aşağıdaki yönetim denge denklemine yol açar:

{displaystyle mathbf {R} = mathbf {Kr} + mathbf {R} ^ {o} qquad qquad qquad mathrm {(2)}}

nerede

{displaystyle mathbf {R}}

= sistemin düğümlerine uygulanan dış kuvvetleri temsil eden düğüm kuvvetlerinin vektörü.

{displaystyle mathbf {K}}

= bireyin kolektif etkisi olan sistem sertliği matrisi elemanların sertlik matrisleri :

{displaystyle mathbf {k} ^ {e}}

.

{displaystyle mathbf {r}}

= sistemin düğüm yer değiştirmelerinin vektörü.

{displaystyle mathbf {R} ^ {o}}

= Önceki düğüm kuvveti vektörüne zaten dahil edilmiş düğüm kuvvetleri dışındaki tüm dış etkileri temsil eden eşdeğer düğüm kuvvetlerinin vektörü R. Bu dış etkiler, dağıtılmış veya konsantre yüzey kuvvetlerini, vücut kuvvetlerini, termal etkileri, ilk gerilmeleri ve gerilmeleri içerebilir.

Desteklerin kısıtlamaları hesaba katıldığında, düğüm yer değiştirmeleri çözülerek bulunur. doğrusal denklem sistemi (2), sembolik olarak:

{displaystyle mathbf {r} = mathbf {K} ^ {- 1} (mathbf {R} -mathbf {R} ^ {o}) qquad qquad qquad mathrm {(3)}}

Daha sonra, tek tek elementlerdeki suşlar ve gerilmeler aşağıdaki gibi bulunabilir:

{displaystyle mathbf {epsilon} = mathbf {Bq} qquad qquad qquad qquad mathrm {(4)}}

{displaystyle mathbf {sigma} = mathbf {E} (mathbf {epsilon} -mathbf {epsilon} ^ {o}) + mathbf {sigma} ^ {o} = mathbf {E} (mathbf {Bq} -mathbf {epsilon} ^ {o}) + mathbf {sigma} ^ {o} qquad qquad qquad mathrm {(5)}}

nerede

{displaystyle mathbf {q}}

= düğüm yer değiştirmelerinin vektörü - sistem yer değiştirme vektörünün bir alt kümesi r söz konusu unsurlarla ilgilidir.

{displaystyle mathbf {B}}

= düğüm yer değiştirmelerini dönüştüren uzama yer değiştirme matrisi q elemanın herhangi bir noktasında suşlar.

{displaystyle mathbf {E}}

= Elemanın herhangi bir noktasında etkili gerilmeleri gerilmelere dönüştüren esneklik matrisi.

{displaystyle mathbf {epsilon} ^ {o}}

= elemanlardaki ilk suşların vektörü.

{displaystyle mathbf {sigma} ^ {o}}

= elemanlardaki başlangıç gerilmelerinin vektörü.

Uygulayarak sanal çalışma denklem (1) sisteme, eleman matrislerini kurabiliriz ${displaystyle mathbf {B}}$ , ${displaystyle mathbf {k} ^ {e}}$ ve sistem matrislerini birleştirme tekniği ${displaystyle mathbf {R} ^ {o}}$ ve ${displaystyle mathbf {K}}$ . Gibi diğer matrisler ${displaystyle mathbf {epsilon} ^ {o}}$ , ${displaystyle mathbf {sigma} ^ {o}}$ , ${displaystyle mathbf {R}}$ ve ${displaystyle mathbf {E}}$ bilinen değerlerdir ve doğrudan veri girişinden ayarlanabilir.

Enterpolasyon veya şekil fonksiyonları

İzin Vermek ${displaystyle mathbf {q}}$ tipik bir elemanın düğüm yer değiştirmelerinin vektörü. Öğenin diğer herhangi bir noktasındaki yer değiştirmeler, aşağıdakiler kullanılarak bulunabilir: interpolasyon sembolik olarak şu şekilde çalışır:

{displaystyle mathbf {u} = mathbf {N} mathbf {q} qquad qquad qquad mathrm {(6)}}

nerede

{displaystyle mathbf {u}}

= elemanın herhangi bir {x, y, z} noktasındaki yer değiştirme vektörü.

{displaystyle mathbf {N}}

= matris şekil fonksiyonları olarak hizmet etmek interpolasyon fonksiyonlar.

Denklem (6), büyük ilgi gören diğer miktarlara yol açar:

Sanal düğüm yer değiştirmelerinin bir işlevi olan sanal yer değiştirmeler: ${displaystyle delta mathbf {u} = mathbf {N} delta mathbf {q} qquad qquad qquad mathrm {(6b)}}$
Eleman düğümlerinin yer değiştirmelerinden kaynaklanan elemanlardaki suşlar: ${displaystyle mathbf {epsilon} = mathbf {Du} = mathbf {DNq} qquad qquad qquad qquad mathrm {(7)}}$

nerede

{displaystyle mathbf {D}}

= matris diferansiyel operatörler kullanarak yer değiştirmeleri suşlara dönüştüren doğrusal esneklik teori. Denklem (7) matrisin B (4) içinde

{displaystyle mathbf {B} = mathbf {DN} qquad qquad qquad qquad mathrm {(8)}}

Öğenin sanal düğüm yer değiştirmeleriyle tutarlı sanal suşlar: ${displaystyle delta {oldsymbol {epsilon}} = mathbf {B} delta mathbf {q} qquad qquad qquad qquad mathrm {(9)}}$

Tipik bir unsurda dahili sanal çalışma

Tipik bir hacim öğesi için ${displaystyle V ^ {e}}$ , sanal yer değiştirmeler nedeniyle dahili sanal çalışma, (5) ve (9) 'un (1)' e değiştirilmesi ile elde edilir:

{displaystyle {mbox {Dahili sanal çalışma}} = int _ {V ^ {e}} delta {oldsymbol {epsilon}} ^ {T} {oldsymbol {sigma}}, dV ^ {e} = delta mathbf {q} ^ {T} int _ {V ^ {e}} mathbf {B} ^ {T} {ig {} mathbf {E} (mathbf {Bq} -mathbf {epsilon} ^ {o}) + mathbf {sigma} ^ { o} {ig}}, dV ^ {e} qquad mathrm {(10)}}

Eleman matrisleri

Öncelikle referans kolaylığı açısından, tipik bir elemanla ilgili aşağıdaki matrisler şimdi tanımlanabilir:

Eleman sertlik matrisi

{displaystyle mathbf {K} ^ {e} = int _ {V ^ {e}} mathbf {B} ^ {T} mathbf {E} mathbf {B}, dV ^ {e} qquad mathrm {(11)}}

Eşdeğer eleman yük vektörü

{displaystyle mathbf {Q} ^ {oe} = int _ {V ^ {e}} - mathbf {B} ^ {T} {ig (} mathbf {E} mathbf {epsilon} ^ {o} -mathbf {sigma} ^ {o} {ig)}, dV ^ {e} qquad mathrm {(12)}}

Bu matrisler genellikle sayısal olarak değerlendirilir. Gauss kuadratürü için Sayısal entegrasyon Onların kullanımı (10) 'u şu şekilde basitleştirir:

{displaystyle {mbox {Dahili sanal çalışma}} = delta mathbf {q} ^ {T} {ig (} mathbf {K} ^ {e} mathbf {q} + mathbf {Q} ^ {oe} {ig)} qquad mathrm {(13)}}

Sistem düğüm yer değiştirmeleri açısından eleman sanal çalışması

Düğüm yer değiştirme vektöründen beri q sistem düğüm yer değiştirmelerinin bir alt kümesidir r (bitişik elemanlarla uyumluluk için), değiştirebiliriz q ile r eleman matrislerinin boyutunu yeni sütunlarla ve sıfır satırlarıyla genişleterek:

{displaystyle {mbox {Dahili sanal çalışma}} = delta mathbf {r} ^ {T} {ig (} mathbf {K} ^ {e} mathbf {r} + mathbf {Q} ^ {oe} {ig)} qquad mathrm {(14)}}

basitlik için, artık genişletilmiş boyuta ve uygun şekilde yeniden düzenlenmiş satır ve sütunlara sahip olan eleman matrisleri için aynı sembolleri kullanıyoruz.

Sistem sanal çalışması

Tüm öğeler için dahili sanal çalışmayı (14) toplamak, (1) 'in sağ tarafını verir:

{displaystyle {mbox {Sistem dahili sanal çalışması}} = toplam _ {e} delta mathbf {r} ^ {T} {ig (} mathbf {k} ^ {e} mathbf {r} + mathbf {Q} ^ {oe } {ig)} = delta mathbf {r} ^ {T} {ig (} sum _ {e} mathbf {k} ^ {e} {ig)} mathbf {r} + delta mathbf {r} ^ {T} toplam _ {e} mathbf {Q} ^ {oe} qquad mathrm {(15)}}

Şimdi (1) 'in sol tarafı düşünüldüğünde, sistem harici sanal çalışması şunlardan oluşur:

Düğüm kuvvetleri tarafından yapılan iş R: ${displaystyle delta mathbf {r} ^ {T} mathbf {R} qquad mathrm {(16)}}$
Dış güçler tarafından yapılan iş ${displaystyle mathbf {T} ^ {e}}$ bölümünde ${displaystyle mathbf {S} ^ {e}}$ elementlerin kenarları veya yüzeyleri ve vücut kuvvetleri tarafından ${displaystyle mathbf {f} ^ {e}}$

{displaystyle sum _ {e} int _ {S ^ {e}} delta mathbf {u} ^ {T} mathbf {T} ^ {e}, dS ^ {e} + sum _ {e} int _ {V ^ {e}} delta mathbf {u} ^ {T} mathbf {f} ^ {e}, dV ^ {e}}

(6b) 'nin değiştirilmesi şunu verir:

{displaystyle delta mathbf {q} ^ {T} sum _ {e} int _ {S ^ {e}} mathbf {N} ^ {T} mathbf {T} ^ {e}, dS ^ {e} + delta mathbf {q} ^ {T} toplamı _ {e} int _ {V ^ {e}} mathbf {N} ^ {T} mathbf {f} ^ {e}, dV ^ {e}}

veya

{displaystyle -delta mathbf {q} ^ {T} toplamı _ {e} (mathbf {Q} ^ {te} + mathbf {Q} ^ {fe}) qquad mathrm {(17a)}}

aşağıda tanımlanan ek eleman matrislerini sunduk:

{displaystyle mathbf {Q} ^ {te} = - int _ {S ^ {e}} mathbf {N} ^ {T} mathbf {T} ^ {e}, dS ^ {e} qquad mathrm {(18a)} }

{displaystyle mathbf {Q} ^ {fe} = - int _ {V ^ {e}} mathbf {N} ^ {T} mathbf {f} ^ {e}, dV ^ {e} qquad mathrm {(18b)} }

Tekrar, Sayısal entegrasyon değerlendirmeleri için uygundur. Benzer bir yedek q (17a) içinde r vektörleri yeniden düzenleyip genişlettikten sonra verir

{displaystyle mathbf {Q} ^ {te}, mathbf {Q} ^ {fe}}

:

{displaystyle -delta mathbf {r} ^ {T} toplamı _ {e} (mathbf {Q} ^ {te} + mathbf {Q} ^ {fe}) qquad mathrm {(17b)}}

Sistem matrislerinin montajı

(16), (17b) 'yi toplamak ve toplamı (15)' e eşitlemek şunu verir: ${displaystyle delta mathbf {r} ^ {T} mathbf {R} -delta mathbf {r} ^ {T} toplamı _ {e} (mathbf {Q} ^ {te} + mathbf {Q} ^ {fe}) = delta mathbf {r} ^ {T} {ig (} toplam _ {e} mathbf {k} ^ {e} {ig)} mathbf {r} + delta mathbf {r} ^ {T} toplamı _ {e} mathbf {Q} ^ {oe}}$

Sanal yer değiştirmelerden beri ${displaystyle delta mathbf {r}}$ keyfi ise, önceki eşitlik şu şekilde azalır:

${displaystyle mathbf {R} = {ig (} toplam _ {e} mathbf {k} ^ {e} {ig)} mathbf {r} + toplam _ {e} {ig (} mathbf {Q} ^ {oe} + mathbf {Q} ^ {te} + mathbf {Q} ^ {fe} {ig)}}$

(2) ile yapılan karşılaştırma şunu göstermektedir:

Sistem sertlik matrisi, elemanların sertlik matrislerinin toplanmasıyla elde edilir:

{displaystyle mathbf {K} = toplam _ {e} mathbf {k} ^ {e}}

Eşdeğer düğüm kuvvetlerinin vektörü, elemanların yük vektörlerinin toplanmasıyla elde edilir:

{displaystyle mathbf {R} ^ {o} = toplam _ {e} {ig (} mathbf {Q} ^ {oe} + mathbf {Q} ^ {te} + mathbf {Q} ^ {fe} {ig)} }

Pratikte, eleman matrisleri ne genişletilir ne de yeniden düzenlenir. Bunun yerine, sistem sertlik matrisi ${displaystyle mathbf {K}}$ bireysel katsayılar eklenerek birleştirilir ${displaystyle {k} _ {ij} ^ {e}}$ -e ${displaystyle {K} _ {kl}}$ alt simgeler ij, kl, elemanın düğüm yer değiştirmelerinin ${displaystyle {q} _ {i} ^ {e}, {q} _ {j} ^ {e}}$ sırasıyla sistemin düğüm yer değiştirmeleri ile eşleşir ${displaystyle {r} _ {k}, {r} _ {l}}$ . Benzer şekilde, ${displaystyle mathbf {R} ^ {o}}$ bireysel katsayılar eklenerek birleştirilir ${displaystyle {Q} _ {i} ^ {e}}$ -e ${displaystyle {R} _ {k} ^ {o}}$ nerede ${displaystyle {q} _ {i} ^ {e}}$ maçlar ${displaystyle {r} _ {k}}$ . Bu doğrudan ilavesi ${displaystyle {k} _ {ij} ^ {e}}$ içine ${displaystyle {K} _ {kl}}$ prosedüre adını verir Doğrudan Sertlik Yöntemi.

Ayrıca bakınız

Referanslar

^ Çerçeveli Yapıların Matris Analizi, Jr.William Weaver, 3. Baskı, James M. Gere, Springer-Verlag New York, LLC, ISBN 978-0-412-07861-3, 1966
^ Matris Yapısal Analiz Teorisi, J. S. Przemieniecki, McGraw-Hill Kitap Şirketi, New York, 1968
^ Hinton, Ernest; Irons, Bruce (Temmuz 1968). "Deneysel verilerin sonlu elemanlar kullanılarak düzleştirilmesi en küçük kareler". Gerginlik. 4 (3): 24–27. doi:10.1111 / j.1475-1305.1968.tb01368.x.
^ Argyris, J.H ve Kelsey, S. Enerji teoremleri ve Yapısal Analiz Butterworth Scientific yayınları, Londra, 1954
^ Clough, R.W, "Düzlem Gerilme Analizinde Sonlu Eleman." Bildiriler, 2. ASCE Elektronik Hesaplamalar Konferansı, Pittsburgh, Eylül 1960

[1] Çerçeveli Yapıların Matris Analizi, Jr.William Weaver, 3. Baskı, James M. Gere, Springer-Verlag New York, LLC, ISBN 978-0-412-07861-3, 1966

[2] Matris Yapısal Analiz Teorisi, J. S. Przemieniecki, McGraw-Hill Kitap Şirketi, New York, 1968

[3] Hinton, Ernest; Irons, Bruce (Temmuz 1968). "Deneysel verilerin sonlu elemanlar kullanılarak düzleştirilmesi en küçük kareler". Gerginlik. 4 (3): 24–27. doi:10.1111 / j.1475-1305.1968.tb01368.x.

[4] Argyris, J.H ve Kelsey, S. Enerji teoremleri ve Yapısal Analiz Butterworth Scientific yayınları, Londra, 1954

[5] Clough, R.W, "Düzlem Gerilme Analizinde Sonlu Eleman." Bildiriler, 2. ASCE Elektronik Hesaplamalar Konferansı, Pittsburgh, Eylül 1960

[1]

[2]

[3]

[4]

[5]