Matlab自编C3D8、C3D20、C3D10单元计算动荷载问题

概述：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的计算文件可以分享，如有需要，请私信我。