从nastran中输出单元刚度矩阵到PCH文件中，利用此代码可以从文件中提取出单元刚度矩阵，方便后面的分析与计算

各向异性材料的刚度矩阵计算M文件，适合初学者学习和运用

六、单元刚度矩阵 五、应力 四、应变 七、等效结点力

三单元刚度矩阵，用计算三单元刚度矩阵计算，输入坐标即得到矩阵

函数 K=Matrix_Assembly（元素） 输入： --------- Elements：一个结构包含 Elements{i}.K 和 Elements{i}.DOFs（第 i 个元素的刚度和自由度）。 输出： --------- K：稀疏全局刚度矩阵

基于C语言的钻柱有限元单元刚度矩阵计算.pdf

函数 K=Fast_Matrix_Assembly（元素） 输入： --------- Elements：一个结构包含 Elements{i}.K 和 Elements{i}.DOFs（第 i 个元素的刚度和自由度）。 输出： --------- K：稀疏全局刚度矩阵 参考----------------- Cuvelier，François，Caroline Japhet和Gilles Scarella。 “在 Matlab 和 Octave 中执行有限元矩阵组装的有效方法。” arXiv 预印本 arXiv:1305.3122 (2013)。

在CAE领域，从学校、实验室的自研算法到实现真正的商业化软件是一条无比漫长的道路。本文不研究有限元的新方法、新理论，本文只是研究商用有限元软件的实现方式。有限元的理论发展了几十年已经相当成熟，商用有限元软件同样也是采用这些成熟的有限元理论，只是在实际应用过程中，商用软件在这些传统的理论基础上会做相应的修正以解决工程中遇到的不同问题，且各家软件的修正方法都不一样，每个主流商用软件手册中都会注明各个单元的理论采用了哪种理论公式，但都只是提一下用什么方法修正，很多没有具体的实现公式。 我们关注CAE中的结构有限元，所以选择了商用结构有限元软件中文档相对较完备的Abaqus来研究内部实现方式。一方面我们查阅Abaqus软件手册得到修正方法的说明，另一方面我们自己编程实现简单的结构有限元求解器，通过自研求解器和Abaqus的结果比较结合理论手册如同管中窥豹一般来研究Abaqus的修正方法，从而猜测商用有限元软件的内部计算方法。在研究的同时，准备将自己的研究成果记录下来写成一个系列文章，希望对那些不仅仅满足使用软件，而想了解软件内部实现方法甚至是做自己的软件的朋友有些帮助。由于水平有限，里面可能有许多错误，欢迎交流讨论。

更多详细内容可以下载《有限元分析》曾攀的代码附件

matlab编写的桁架有限元输入，通过输入单元的长度，弹性模量等，可以生成刚度矩阵