概述：matlab自編程序求解C3D8（一階六面體）、C3D20（二階六面體）、C3D10（二階四面體）單元的剛度矩陣、質量矩陣和阻尼矩陣，并采用newmark求解動荷載問題，計算結果與abaqus保持一致。

視頻首先采用abaqus計算懸臂梁受動荷載算例，與此同時修改inp文件，并導出abaqus計算的總體剛剛度矩陣、質量矩陣和阻尼矩陣。然后運行matlab讀取inp文件的模型信息，自己求解單元的剛度矩陣、質量矩陣和阻尼矩陣與abaqus計算的結果進行對比。

結論：

1、C3D8（一階六面體）采用B-Bar修正避免剪切自鎖，將協調質量矩陣處理為集中質量矩陣。

2、C3D20（二階六面體）單元的質量為協調質量矩陣。

3、C3D10（二階四面體）單元質量為協調質量矩陣，特別注意，質量矩陣采用15積分點的錘子積分計算。

4、自編newmark積分程序求解動荷載問題。

5，最終計算的剛度矩陣、質量矩陣和阻尼矩陣均與abaqus保持一致。求解動荷載問題的位移、速度和加速度與abaqus保持一致。

注：

1、最近感冒，視頻中有兩處（或者三處）吸鼻涕聲音，請諒解。

2、視頻里面abaqus的計算文件可以分享，如有需要，請私信我。