İçeriğe atla

Kayma Modülü

Teradigma sitesinden

Kayma Modülü: Katılarda Şekil Direncinin Atomik Kökenleri ve Nizam Boyutları

Katı cisimlerin hacimlerini değiştirmeksizin şekillerini direnç göstererek koruması, doğanın en dikkat çekici mekanik olgularından birini oluşturmaktadır. Kayma modülü olarak adlandırılan bu direnç parametresi — makaslama modülü ya da rijitlik modülü olarak da bilinir — bir malzemenin yüzeyine paralel uygulanan kuvvetler karşısındaki davranışını nicel olarak tanımlamakta; mühendislik hesaplamalarından yer kürenin iç yapısına, insan kıkırdağından ileri malzeme bilimine uzanan geniş bir yelpazeye yön vermektedir. Bir malzemenin Young modülü (E) çekme-basma karşısındaki davranışını, hacim modülü (K) ise hidrostatik basınç altındaki tepkisini ölçerken, kayma modülü (G) bu üç temel esneklik sabitinin tamamlayıcı öğesi olarak şekil değişimine özgü direnci ifade etmektedir.[1]


1. Temel Kavramlar ve Matematiksel Çerçeve

1.1 Kayma Gerilmesi ve Kayma Gerinimi

Bir katı cismin alt yüzeyi sabitlenirken üst yüzeyine paralel bir F kuvveti uygulandığında, kuvvetin etki ettiği yüzeyin alanı A ise, oluşan kayma gerilmesi (shear stress) şu şekilde tanımlanır:[2]

τ=FA

Bu gerilmenin etkisiyle cisim üst yüzeyinde Δx kadar yana doğru kayarken, alt ve üst yüzeyler arasındaki başlangıç yüksekliği L olmak üzere, boyutsuz bir büyüklük olan kayma gerinimi (shear strain) elde edilir:

γ=ΔxL=tanϕϕ

Burada ϕ, küçük deformasyonlar için kayma açısıdır. Hooke Yasası’nın kayma gerilmesi-gerinimi eşdeğerinde, iki büyüklük arasındaki orantılılık sabiti kayma modülü G’yi tanımlar:

τ=Gγ

Buradan kayma modülünün temel ifadesi elde edilir:

G=τγ=F/AΔx/L=FLAΔx

G’nin birimi Pascal (Pa) ya da N/m²’dir ve genellikle gigapaskal (GPa) mertebesinde ifade edilir.[3]

1.2 Üç Temel Esneklik Sabiti Arasındaki Bağlantı

İzotropik, lineer elastik ve homojen bir malzeme için dört temel esneklik parametresi — Young modülü E, kayma modülü G, hacim modülü K ve Poisson oranı ν — birbirinden bağımsız değildir. Bu parametrelerden yalnızca ikisi bağımsız seçilebilir; diğerleri analitik ilişkilerle hesaplanır.[4]

Young modülü ile kayma modülü arasındaki bağıntı, Poisson oranı üzerinden kurulur:

G=E2(1+ν)

Hacim modülü ile kayma modülü ise şu eşitlikle ilişkilendirilir:

K=EG3(3GE)

ya da eşdeğer biçimiyle:

K=2G(1+ν)3(12ν)

Bu bağıntılar yalnızca pratik bir kolaylık aracı değildir; aynı zamanda bir malzemenin tek bir türde yük altındaki davranışından hareketle diğer tüm yük türlerine karşı tepkisinin kestirilebileceğini göstermesi bakımından dikkat çekici bir nizam içermektedir. Metalik malzemelerde G genellikle E’nin 0,35–0,45 katı kadardır; bu oran, metalik bağların çekme-basma ile kesme yüklerine karşı gösterdikleri direncin tutarlı bir oransal ilişki sergilediğini ortaya koymaktadır.[5]

1.3 Torsiyonel Rijitlik ve Eksen Dönmesi

Kayma modülünün en doğrudan uygulaması, silindirik bir millerin burulma (torsiyon) analizidir. Uzunluğu L, yarıçapı r ve kayma modülü G olan bir mile T büyüklüğünde bir dönme momenti (torque) uygulandığında, milin iki ucu arasındaki burulma açısı θ şu şekilde hesaplanır:[6]

θ=TLGJ

Burada J=πr42 silindirik kesit için kutupsal atalet momentidir. GJ çarpımı “torsiyonel rijitlik” olarak adlandırılır ve bir millik uzunluk başına birim dönme açısına yol açmak için gerekli torku ifade eder.


Örnek 1: Çelik Mil Burulma Hesabı

Standart yapı çeliği için G80 GPa değeri alındığında, dış çapı 40 mm (r=0,02 m), uzunluğu L=1,5 m olan içi dolu bir mil T=500 N·m’lik bir dönme momentine maruz kalırsa:

J=π(0,02)42=π×1,6×10722,51×107 m4

Ayrıştırılamadı (sözdizim hatası): {\displaystyle \theta = \frac{500 \times 1{,}5}{80 \times 10^9 \times 2{,}51 \times 10^{-7}} = \frac{750}{20080} \approx 0{,}0373 \text{ rad} \approx 2{,}14° }

Bu değer, 1,5 metrelik bir çelik milin, standart endüstriyel yükler altında yalnızca yaklaşık 2° civarında bürüldüğünü göstermektedir. Aynı hesap alüminyum için (G26 GPa) yapıldığında açı yaklaşık 6,5°’ye çıkar — çelik milin yaklaşık üç katı. Bu mertebedeki fark, neden yüksek torsiyon altındaki mil ve aksların büyük ölçüde çelik ya da titanyum alaşımından tercih edildiğini doğrudan açıklamaktadır.[7]


Örnek 2: Kayma Modülü ile Poisson Oranı Hesabı

Bir malzemenin Young modülü E=200 GPa ve kayma modülü G=77 GPa olarak ölçüldüğünde, Poisson oranı şu şekilde belirlenir:

G=E2(1+ν)ν=E2G1=2002×771=20015411,2991=0,299

Bu sonuç tipik yapı çeliği değeri olan ν0,30 ile mükemmel örtüşmektedir. Dikkat çekici nokta şudur: G ve E değerleri ayrı ayrı ölçülmüş olmasına karşın, bu iki ölçümden bütünüyle farklı bir fiziksel büyüklük olan Poisson oranı — malzemenin yanal şişme eğilimi — tam olarak kestirilebilmektedir. Üç büyüklüğü aynı anda bağlayan bu bağıntı, izotropik esneklik teorisinin içsel tutarlılığını yansıtmaktadır.[8]


1.4 Atomik Kökenleri: Kayma Modülü Neden Maddenin İçinden Kaynaklanır?

Kayma modülü, makroskopik bir ölçüm olmasına karşın kökleri atomlar arası etkileşimlere dayanmaktadır. Kristal bir katıda atomlar denge konumlarında bulunurken, aralarındaki potansiyel enerji bir minimum sergiler. Yüzeye paralel bir kuvvet uygulandığında atomik düzlemler birbirlerine görece kayar; ancak bağlar kısmen bozulup kısmen gerildiğinden potansiyel enerji yükselir. İşte bu potansiyel enerji eğrisinin eğrilik katsayısı (ikinci türevi), kayma modülünü mikroskopik düzeyde belirler.[9]

Açısal bağımlı atomlar arası potansiyeller üzerine yapılan güncel çalışmalar (Zaccone, ETH Zurich), amorf camların kayma modülünü Cauchy-Born serbest enerji açılımı üzerinden türetmiş; merkezi çift-etkileşim potansiyellerinin yanı sıra bağ-bükme (bond-bending) etkileşimlerinin de kayma rijitliğine özgün katkılar yaptığını göstermiştir.[10] Bu bulgular, kayma modülünün atomik düzlem kayması ve uzun-menzilli korele hareketlerin birleşiminden doğduğunu ortaya koymaktadır — yani kayma modülü atomun “içinde” değil, atomların “arasındaki düzenlemenin” içinde saklıdır. Bu ayrım, ilerleyen kavramsal analizde merkezi bir önem taşıyacaktır.

Sıcaklık artışıyla beraber termal titreşimler (fononlar), atomar potansiyel kuyularının anharmonik katkılarından dolayı bağların efektif sertliğini zayıflatır. Bu nedenle G, mutlak sıfırdan erime noktasına kadar yaklaşık doğrusal bir düşüş sergiler; erime noktasında ise kritik kayma gerimine ulaşıldığında G hızla sıfıra yaklaşır. Lindemann erime kriteri bu çerçevede yeniden formüle edilmiş; Gc olarak gösterilen erime kristalinin ortalama kayma modülü, sıvı fazın farklı iç durum konfigürasyonlarına geçişi için kritik bir kayma gerinimi belirlediği gösterilmiştir.[11]


Örnek 3: Seismolojide Kayma Hızı ve G İlişkisi

Yerin mantosu gibi büyük ölçekte katı davranan ortamlarda, cisimsel kayma dalgaları (S-dalgaları) vS hızıyla ilerler:

vS=Gρ

Üst manto için tipik değerler vS4,5 km/s ve ρ3300 kg/m3 alındığında:

G=ρvS2=3300×(4500)2=3300×2,025×1076,68×1010 Pa67 GPa

Bu değer, yüzey kayalarında ölçülen G20–40 GPa aralığının neredeyse iki katıdır. Sismolojik ölçümler sayesinde hiçbir zaman doğrudan erişilemeyen yerin iç katmanlarının mekanik özellikleri, yüzeyde ölçülen dalga hızlarından hesaplanmaktadır. Kilometre derinliklerdeki kayma modülü dağılımındaki ani değişimler, Mohorovičić süreksizliği gibi jeofizik sınırları belirlemekte; 2024 yılında Wyoming Kratonu kenarında saptanan kayma hızı laterel değişimleri, üst manto yapısındaki keskin sınırları litosfer-astenosfer etkileşimi bağlamında yorumlamak için kullanılmıştır.[12]


2. Güncel Araştırmalardan Bulgular

2.1 Kıkırdak Dokusunun Kayma Özellikleri ve Rijidite Perkolasyonu

2025 yılında yayımlanan kapsamlı bir derleme (Oliveira ve ark., Progress in Biomedical Engineering, IOP), eklem kıkırdağının dinamik kayma modülünün G* olarak sığır kıkırdağında 0,1–0,4 MPa aralığında ölçüldüğünü ve kayıp açısının (loss angle) ~10° mertebesinde kaldığını aktarmaktadır.[13] Bu büyüklükler, kıkırdağın kayma yükleri altındaki tepkisinin fibriler bir kollajen ağ ile agregkan proteoglikan matriksinin viskoelastik etkileşiminden doğduğunu ortaya koymaktadır.

Yapısal Kökenlerin Rijidite Perkolasyon Çerçevesinden Yorumlanması başlıklı bir çalışmada (PMC, Biophysical Journal), kıkırdak yüzeyine yakın katmanlarda kollajen ağının seyrekleşmesiyle kayma modülünün katının iç bölgelerine kıyasla 100 kata kadar düştüğü saptanmıştır.[14] Üstelik bu yüzey katmanı, kayma enerjisinin yaklaşık %90’ını dağıtarak adeta biyolojik bir şok emici olarak işlev görmektedir — bu özelliğin doku içinde derinlik boyunca hassas biçimde konumlandırılmış olması dikkat çekicidir.

2.2 Oksitik Yardımcı Metamalzemeler ve Sıfıra Yakın Kayma Modülü

2026 yılı başlarında yayımlanan EPJ Applied Metamaterials dergisindeki bir çalışma, negatif Poisson oranına sahip oksitik metamalzemelerin (auxetic metamaterials), geometrik birim hücre tasarımı sayesinde son derece yüksek kayma modülü, darbe direnci ve enerji absorpsiyon kapasitesi sergilediğini göstermiştir.[15] Klasik malzemelerin homojen ve izotropik kayma davranışının aksine, bu yapılarda kayma rijitliği şekil geometrisine bağlı olarak istenilen değerlere tasarlama yoluyla getirilebilmektedir. 2025 tarihli bir diğer derleme (Frontiers in Materials), bazı mekanik metamalzemelerin sıfıra yakın kayma modülü (vanishing shear modulus) sergilediğini bildirmekte; bu özelliğin şekil değişimi olmaksızın ses dalgalarının yönlendirilmesini mümkün kıldığını aktarmaktadır.[16]

2.3 Elastik Modüller ve Kimyasal Bağ Anharmonikliği

2023 yılında Angewandte Chemie International Edition’da yayımlanan bir çalışma (Isotta ve ark.), termoelektrik malzemelerde kayma modülünün atomlar arası potansiyel enerjisinin anharmonik katkılarına son derece duyarlı olduğunu ortaya koymuştur.[17] Farklı kristalografik yönlerde bağların gerilmesi ya da kayması, Cij esneklik tensörünün her bir elemanına sıcaklık bağımlılığı kazandırmaktadır. Saf harmonik potansiyelde sabit kalması beklenen bu elemanlar, gerçek katılardaki anharmonik terimler nedeniyle sıcaklıkla sistematik biçimde azalmaktadır.

Demir ve ağır geçiş metallerinden oluşan alaşımlar üzerinde yürütülen moleküler dinamik simülasyonları (2023–2024, Journal of Materials Research and Technology), kayma modülü yumuşamasının erime süreci için Lindemann kriteri ile doğrudan örtüştüğünü teyit etmiştir.[18] Bu bağlantı, dislokasyon dinamiğini ve kırılma mekaniklerini anlama açısından kritik veriler sunmaktadır.


3. Kavramsal Analiz

3.1 Nizam, Gaye ve Sanat Analizi

Kayma modülünün sayısal değeri maddeden maddeye muazzam biçimde değişmektedir: Yumuşak biyolojik dokularda G104105 Pa düzeyindeyken, çelik için G80×109 Pa’ya, elmas için ise G500×109 Pa’ya ulaşmaktadır. Bu, beş ila altı büyüklük mertebesini kapsayan bir aralıktır. Her bir malzemenin kendine özgü G değeri o malzemeyi oluşturan atomların kimyasal kimliği, kristal örgü geometrisi ve bağ türünün son derece hassas bir bileşimiyle belirlenmektedir; bu mertebedeki bir çeşitlilik, tek bir parametrenin ne denli çok boyutlu bilgiyi kodladığını ortaya koymaktadır.

Özellikle dikkat çekici olan şudur: Eklem kıkırdağı için saptanan G*0,1–0,4 MPa değeri, yüzeyden dokuya doğru kademeli olarak artmakta ve üstelik enerji dağıtımının %90’ı yüzeydeki en düşük G değerine sahip bölgeye yüklenmektedir.[19] Bu, israfı en aza indirip yükleri en uygun biçimde dağıtan bir işlevsel hiyerarşi örneğidir. Eğer kıkırdak boyunca G tekdüze dağılsaydı, darbe enerjisi tek bir noktada yoğunlaşacak ve doku kalıcı hasara uğrayacaktı; ne var ki derinliğe bağlı bu kademeli G dağılımı, onlarca yıl süren günlük mekanik yüklere karşın dokunun bütünlüğünü korumaktadır. Bir ihtiyaca göre işlev dağıtma ilkesinin somut bir tezahürü olarak dikkat çekmektedir.

Aynı nizam, çok farklı bir ölçekte, yerkürenin katmanlı yapısında da karşımıza çıkmaktadır. Sismik dalgaların iç kattaki kayma modülü profilini gözler önüne sermesi, manto dışında —hiçbir zaman erişilemeyen derinliklerde— hangi malzemelerin hangi mekanik özelliklerle var olduğunu açıkça ortaya koymaktadır. 67 GPa’lık üst manto kayma modülü değeri yüzeyden derinleşerek artar; bu artış, üst manto kayaçlarının tektoniğin yavaş hareketlerine izin verirken ani deprem sarsıntılarına karşı katı gibi davranmasını sağlayan kritik bir özellik olarak göze çarpmaktadır. Yerkürenin yalnızca dıştan gözlemlenebilir sismik dalga hızlarıyla iç yapısı bu denli hassas biçimde tespit edilebiliyorsa, bu sonuç üç büyüklüğün — G, ρ ve vS — tutarlı bir nizam içinde ilişkilendirildiğini göstermektedir.

Meseleye “bu bize ne anlatıyor?” diye sorulduğunda, alışkanlık perdesinin kaldırıldığı her anda şaşırtıcı bir gerçekle karşılaşılmaktadır: G, E ve ν gibi üç parametre aynı anda bilinirse, geri kalan tüm esneklik özelliklerinin eksiksiz türetilebilmesi — bir malzemenin bütün boyutlardaki mekanik davranışının yalnızca iki bağımsız sabite indirgenmesi — fiziksel dünyanın katı cisimler için gereksiz özgürlük dereceleri içermediğini gösterir. Bu yoğun-bütüncül kodlama, yapının bilgi taşıyan bir sistem olarak örgütlendiğine işaret eder.

Çelik için G/E0,39, alüminyum için G/E0,38, bakır için G/E0,38, titanyum için G/E0,38 değerleri elde edilmektedir. Farklı atom kütlesi, kristal yapısı ve elektronik konfigürasyona sahip bu metallerin G/E oranlarının bu denli yakın çıkması, ilk bakışta tesadüfmüş gibi görünse de metalik bağın doğasından — elektron gazının oluşturduğu perdeleme düzeninin bağ türü ne olursa olsun benzer mekanik oranlar doğurmasından — kaynaklanmaktadır. Farklı hammaddelerin aynı oransal nizama tabi kılınması, düzenin hammaddeye değil, daha üst bir ilkeye ait olduğuna işaret etmektedir.

3.2 İndirgemeci Yaklaşımların ve Fail-Meful Hatasının Eleştirisi

Popüler anlatılarda kayma modülü çoğu kez şu türden ifadelerle aktarılmaktadır: “Çelik kayma kuvvetlerine karşı koyar”, “Malzeme şeklini korumaya çalışır”, “Atomlar bağları sağlamlaştırır”. Bu dilbilimsel kısayollar, bir aktarım kolaylığı sağlasa da nedenselliği araçtan faile yanlış biçimde devreder.

Gerçekte çelik ne “karşı koyar” ne de “çalışır”. Demir atomlarının metalik bağlardan oluşan kristal örgüsünde, atomlar arası potansiyel enerjinin yüzeye-paralel deformasyona karşı gösterdiği eğrilik katsayısı belirli bir değer almaktadır; bu değerin ölçüm sonucu G80 GPa olarak saptanmakta, deformasyona karşı oluşan gerilme bu sabit üzerinden hesaplanmaktadır. Yani “direnç”, atomun kastına değil, atomların belirli bir kristal örgü içinde düzenlenmesinin ürünüdür. Bir piyano akordatörünün teli gererek sesi düzeltmesi gibi: telin kendisi “müziği kastetmez”; ses, telin titreşim frekansından, o frekansı belirleyen gerilim ve kütlesel dağılımından ortaya çıkar.

Benzer bir hata, “kayma modülü yüksek olduğu için çelik tercih edilir” ifadesinde de saklıdır. Sanki malzeme, kendi yüksek G değerine sahip olduğu için bir amaca hizmet edecek şekilde var olmuş izlenimi uyanmaktadır. Oysa fiziksel açıklamanın yönü tersine işlemektedir: Mühendis, belirli bir işlev için gereken kayma rijitliğini hesaplar, bu rijitliği sağlayacak atomik yapıyı ve kristal örgüyü araştırır, ardından bu yapıyı bünyesinde barındıran malzemeyi seçer. İşlev bir amacı, amaç bir seçimi, seçim ise bir iradeyi gerektirmektedir.

Ayrıca indirgemeci bir hataya dikkat çekmek gerekmektedir: Kayma modülünü yalnızca “atomlar arası potansiyelin ikinci türevinin kristal oralaması” biçiminde tanımlamak, büyüklüğü doğru verir ancak mekanizmanın tamamını yakalamaz. Çünkü G, dislokasyon yoğunluğuna, kristal tane sınırlarına, aşılama elementlerine ve hatta işleme tarihine bağlıdır; yani aynı demir atomlarından oluşan ama farklı biçimde işlenmiş iki çelik, farklı G değerleri sergileyebilir. Hammadde (atomlar) sabitken “sanat” (düzenleme biçimi) değişmekte ve fiziksel özellik buna bağlı olarak değişmektedir. Bu “kayma modülü atomlarda mı, yoksa atomların tertibinde mi?” sorusu, gerçek açıklamanın nerede aranması gerektiğini belirlemektedir.[20]

3.3 Hammadde ve Sanat Ayrımı Analizi

Kayma modülü için hammadde ayrımını en çarpıcı biçimde gösteren örnek oksitik (auxetic) metamalzemelerdir. Bu yapılarda negatif Poisson oranı ve son derece yüksek kayma modülü, atomların kimyasından değil, birimlerin geometrik düzenlemesinden kaynaklanmaktadır.[21] Reentrant (içe çekik) hücre geometrisine sahip bir bal peteği yapısı, çekildikçe genişleyen bir malzeme sergiler; oysa aynı madde farklı bir geometrik düzenlemede çekilince daralır. Hammadde — polimer, metal ya da seramik — değişmemiştir; yalnızca düzenleme değişmiştir ve bu düzenleme tamamen farklı kayma özellikleri doğurmuştur.

Benzer bir örnek kıkırdak dokusundaki hiyerarşik ölçekte karşımıza çıkmaktadır: Kollajen molekülleri (hammadde) tek başına G1–10 MPa mertebesinde sertliğe sahiptir; ancak eklem yüzeyine yakın bölgede bu liflerin gevşek bir örgü oluşturduğu katmanda, ağırlıklı olarak agregkan matriksiyle birleşerek oluşturulan düzenleme, doku düzeyinde G0,1 MPa seviyesine düşmekte, bu bölge ise enerji dağıtıcı bir şok emici olarak işlev görmektedir.[22] Kollajen lifin kendi G değeri değişmemiştir; yoğunluk ve yönlenme düzeni değişmiş, bu düzenleme değişikliği ise tamamen farklı bir işlev doğurmuştur.

Yere bir kutu karbon lifi (hammadde) dökülse, bu lifler kendi kendine örgülenip derinliğe bağlı farklılaşmış kayma modülü profiline sahip, yaşam boyu aşınmaya dirençli bir kıkırdak oluşturabilir mi? Kollajen ve agregkan bu yapının lifleri ve matrikidir; ama dokunun işlevsel planı — hangi katmanda ne kadar G olacağı, nasıl bir hiyerarşik geçiş kurulacağı — hangi kaynak tarafından belirlenmektedir?

Bu soruyu, matematiğin bir yanıtı olarak şu gözlem pekiştirmektedir: İzotropik bir malzeme için iki bağımsız esneklik sabiti bilinince, malzemenin tüm olası yük yönlerindeki mekanik tepkisi tam anlamıyla belirlenmektedir. Hammaddenin kendisi belirsizken, bu ikili kodun tamamını taşıyan düzenleme son derece belirlidir. Hammadde sanatın malzemesidir; ancak iki sabit içindeki bu yoğun mekanik bilgi, malzemenin kendi başına sahip olabileceği bir şey değildir.


4. Sonuç

Kayma modülü, ilk bakışta mühendislik hesaplamalarının teknik bir parametresi olarak görünmektedir. Ancak derinlemesine incelendiğinde, bu büyüklük birçok boyutu bir arada barındıran bir kavramsal kesişim noktasına dönüşmektedir: Atomik bağ potansiyelinin ikinci türevi olan G, aynı zamanda izotropik malzemelerin bütün mekanik davranışını iki bağımsız sabitlik içinde sıkıştıran bir kodlama şemasının parçasıdır. Kıkırdağın milyon yıllık mekanik uyuma izin veren derinliğe bağlı G dağılımı, oksitik metamalzemelerde geometrinin kimyadan bağımsız olarak G değerini belirlemesi ve yerin iç yapısının kayma dalgaları üzerinden okunabilmesi, bu parametre etrafında çarpıcı bir nizam tablosunu ortaya koymaktadır.

Kayma modülünün atomik kökeni, kristalografik bağımlılıkları ve sıcaklık duyarlılığı, malzemedeki bilginin yalnızca hammaddede değil, hammaddenin tertip biçiminde saklı olduğunu açıkça ortaya koymaktadır. Üç temel esneklik sabitinin yalnızca iki bağımsız parametre üzerine kurulu olması, fiziksel dünyada bir ekonomi ve nizam anlayışına işaret etmektedir. Oksitik metamalzemelerde irade-yüklü geometrik seçimlerle G’nin istenilen değerlere getirilebilmesi ise tersine mühendislikte dahi irade ve bilgi olmadan istenen sonucun elde edilemeyeceğini göstermektedir. Deliller ışığında, kayma modülünün atomik, yapısal ve sistematik boyutlarında gözlemlenen bu dizaynsal hassasiyet ve nizam; nihai karar okuyucunun kendi aklına ve vicdanına bırakılmaktadır.


Kaynakça

  1. Britannica. (2024). Shear modulus. Encyclopedia Britannica. https://britannica.com/science/shear-modulus
  2. Vedantu Physics. (2024). Shear modulus (modulus of rigidity): Definition, formula & relation to elasticity. https://www.vedantu.com/physics/shear-modulus-elastic-moduli
  3. Vedantu Physics. (2024). Shear modulus (modulus of rigidity): Definition, formula & relation to elasticity. https://www.vedantu.com/physics/shear-modulus-elastic-moduli
  4. Vedantu Physics. (2024). Relation between elastic constants: Formula, derivation & Poisson’s ratio. https://www.vedantu.com/physics/relation-between-elastic-constants
  5. SAM Materials. (2025). Shear modulus (modulus of rigidity). Stanford Advanced Materials. https://www.samaterials.com/content/shear-modulus-(modulus-of-rigidity).html
  6. FIRGELLI Auto. (2026). Shear modulus interactive calculator. https://www.firgelliauto.com/blogs/engineering-calculators/shear-modulus-calculator
  7. FIRGELLI Auto. (2026). Shear modulus interactive calculator. https://www.firgelliauto.com/blogs/engineering-calculators/shear-modulus-calculator
  8. Vedantu Physics. (2024). Relation between elastic constants: Formula, derivation & Poisson’s ratio. https://www.vedantu.com/physics/relation-between-elastic-constants
  9. Isotta, E., et al. (2023). Elastic moduli: A tool for understanding chemical bonding and thermal transport in thermoelectric materials. Angewandte Chemie International Edition, 62(9), e202213649. https://doi.org/10.1002/anie.202213649
  10. Zaccone, A. (2008/2023). The shear modulus of metastable amorphous solids with strong central and bond-bending interactions. arXiv:0807.3656. https://arxiv.org/pdf/0807.3656
  11. Buchenau, U., & Zorn, R. (2013). A new view of the Lindemann criterion. arXiv:1311.4739. https://arxiv.org/abs/1311.4739
  12. Noise.earth/Utah. (2024). Upper mantle earthquakes along the edge of the Wyoming craton. https://www.noise.earth.utah.edu/grl70276.pdf
  13. Oliveira, S., et al. (2025). A guide to articular cartilage functioning: A comprehensive review, current challenges and mechanobiological solutions. Progress in Biomedical Engineering, 7(3), 032009. https://doi.org/10.1088/2516-1091/ade83a
  14. Ghadiali, S., et al. (2022). Structure-function relations and rigidity percolation in the shear properties of articular cartilage. PMC – Biophysical Journal. https://pmc.ncbi.nlm.nih.gov/articles/PMC8836800/
  15. EPJ Applied Metamaterials. (2026). The structure design and application of metamaterials with negative Poisson’s ratio. EPJ Applied Metamaterials, 13, 5. https://doi.org/10.1051/epjam/epjam250035
  16. Frontiers in Materials. (2024). A review on the auxetic mechanical metamaterials and their applications in the field of applied engineering. Frontiers in Materials. https://www.frontiersin.org/journals/materials/articles/10.3389/fmats.2024.1453905/full
  17. Isotta, E., et al. (2023). Elastic moduli: A tool for understanding chemical bonding and thermal transport in thermoelectric materials. Angewandte Chemie International Edition, 62(9), e202213649. https://doi.org/10.1002/anie.202213649
  18. Daw, M., & Foiles, S. (2023). Angular-dependent interatomic potential for large-scale atomistic simulation of W-Mo-Nb ternary alloys. Journal of Materials Research and Technology. https://www.sciencedirect.com/science/article/abs/pii/S0927025623007280
  19. Ghadiali, S., et al. (2022). Structure-function relations and rigidity percolation in the shear properties of articular cartilage. PMC – Biophysical Journal. https://pmc.ncbi.nlm.nih.gov/articles/PMC8836800/
  20. Zaccone, A. (2008/2023). The shear modulus of metastable amorphous solids with strong central and bond-bending interactions. arXiv:0807.3656. https://arxiv.org/pdf/0807.3656
  21. EPJ Applied Metamaterials. (2026). The structure design and application of metamaterials with negative Poisson’s ratio. EPJ Applied Metamaterials, 13, 5. https://doi.org/10.1051/epjam/epjam250035
  22. Ghadiali, S., et al. (2022). Structure-function relations and rigidity percolation in the shear properties of articular cartilage. PMC – Biophysical Journal. https://pmc.ncbi.nlm.nih.gov/articles/PMC8836800/