結論： 1、C3D8（一階六面體）采用B-Bar修正避免剪切自鎖，將協調質量矩陣處理為集中質量矩陣。 2、C3D20（二階六面體）單元的質量為協調質量矩陣。 3、C3D10（二階四面體）單元質量為協調質量矩陣，特別注意，質量矩陣采用15積分點的錘子積分計算。 4、自編newmark積分程序求解動荷載問題。 5，最終計算的剛度矩陣、質量矩陣和阻尼矩陣均與abaqus保持一致。

θ′≠β′的物理來源 第八自由度：截面畸變（箱形截面輪廓變形）與橫向雙力矩 橫隔板/墜板加強體系：離散力法求解與連續化等效 數學同構：隔板連續化方程與彈性地基梁方程的對比及邊界層衰減效應 第三章 空間穩定、幾何剛度與有限元落地邏輯 偏心荷載下的壓彎扭耦合與幾何剛度矩陣構建 變系數特征值問題與 Galerkin 近似解法 有限元軟件底層架構：7自由度單元的形函數基底與剪切自鎖假設