結構勁度矩陣縮減問題

[kassimali]

加入土木人討論區粉絲團

以下題目的如果矩陣沒經過縮減,算出來的位移應該是正確的
但(c)小題指定考慮對稱性縮減勁度矩陣
而學生與同學討論,參考您力學電子報95年雲科大碩班考題解法
重新推導縮減勁度矩陣的型式
結果發現最後以這縮減後的勁度矩陣下去算位移,竟然跟沒縮減有相當的落差...
請問老師,我在(c)小題的推導過程是否有哪裡出了問題?
或是有簡便的作法?
矩陣演算部分我已用MATLAB算過,代大矩陣元素的步驟也有檢查過...
不好意思麻煩老師解惑了,謝謝

P.S.以下照片基於相機陽春,可能有些微糢糊........

[紫煌]

路過給個建議....
反矩陣那邊好像少除了一個det(k)
2階反矩陣會等於原矩陣主對角對調，非對角項差負號，再除上det(K)

[kassimali]

inv([0.2605,-0.0202;-0.0202,0.0991])
這是我用MATLAB進行的驗算的指令......
驗算結果確實有除行列式值

但神奇的是, 我試著照紫煌大您講的除了之後, 答案真的對了...
問題在於我的反矩陣已經除過了= =怎麼想都覺得很怪啊

不好意思麻煩您了,謝謝

[kassimali]

那個......怎麼目前為止只有紫煌大哥一個人回應......?= =

我在猜想, 是不是這題打從一開始r編號順序使我無法用勁度濃縮下去算?
因為r1對應外力0, r2為-p, r3為0, r4為-p......
還是說根本沒影響?

懇請幫忙了,謝謝