有没有人编过壳单元的有限元程序呀？

小弟现在想对ansys的shell63单元进行编程验证，遇到了麻烦，请各位高手给我指点迷津：
1.shell63单元就是DKQ弯曲单元（由四个DKT弯曲单元组集而成）和平面应力单元的组合，那么dkq弯曲单元是什么样子的，如何由四个DKT弯曲单元组集而成？有没有相关的推导？
2.本人对于shell63单元的平面应力分析进行了推导，并运用matlab软件编写了有限元程序计算单元刚度矩阵和单元质量矩阵。在ansys中利用shell63单元的膜分析提取单元质量矩阵和刚度矩阵，结果对比较好！…………相关结果见附件。
3.由于没有相关dkq单元的资料，使用了四节点平板弯曲结果进行了推导，并运用matlab软件编写了有限元程序计算单元刚度矩阵和单元质量矩阵，其结果和ansys的结果相差较大。

我该怎么办呢？高手们多提宝贵建议！