matlab有限元求解问题（调用ansys的单刚）

kingpaul0411

小弟最近一直在做matlab求解有限元问题，发现求出的结果和ansys对应不上，向各各位达人求救。
说说我大约的问题：
问题是3*3m的一个悬臂版，一侧固支，另一侧中点加Z方向1000N的力，分成10*10的网格，弹性模量是靠近固定端的一半为6.9e10*0.0217，靠近载荷的一半为6.9e10，泊松比0.3，我用的是壳单元，做法是把在ansys建立模型，导出单元刚度阵和对号数组，然后利用他的单刚和对号数组来组装我的总刚，在这一过程中，我对比过单刚一致，总刚由于太大我对比了一部分的行列差距都很小。
但是结果差距非常大（最大位移我0.06，ansys0.15），但是当两侧杨氏模量相同时（都为6.9e10），差距并不大（最大位移我0.040，ansys0.041），
反复比对，仍没能找到错误希望高手指点，m文件，两个data文件，以及可能用到的ansys的apdl代码都在附件里了，请不吝赐教感激不尽！！！！