ALGOR能否输出单元刚度阵

nemia

有一个非常简单的例子，在相同的网格下，自己写的程序和ALGOR的结果相对误差在3%左右，现在不知道哪里有问题。请问ALGOR能否输出单元刚度阵，我想比较一下。谢谢。

chenyq

目前没有发现这个功能! 可能是由于误差造成的!也可能是单元的基本参数有点差异造成的!呵呵,很不错,自己写程序!

nemia

其实就是一个很简单的平面应力程序。我再研究研究，看看问题出在哪里。

nemia

Algor可以输出刚度阵对角元素的最大最小及平均值。把这些值和我自己写的程度对比，我的程度得到的值都是偏大的。继续去研究问题所在，同时企盼有大侠支着。

nemia

找到输出整体刚度阵的地方了：Analysis--Parameters--Output，选中Stiffness matrix，分析之后，会生成一个*.mtx文件，刚度阵就存在这里。

[ 本帖最后由 nemia 于 2008-6-20 15:40 编辑 ]

haige

期待得到楼主的结果，我也很奇怪，为什么刚度矩阵会有差别。  楼主可以用第三方软件试试

nemia

发现了问题所在。在"Element Definition"中的"Compatibility"，程序默认是"Not Enforced"，这样Algor的结果和我的程序的结果就不同（Algor的刚度偏小，位移偏大）。如果将"Compatibility"设为"Enforced"，结果就相同了。手边没有其他商用软件了，所以未用其他软件验证。

tonnyw

Compatibility应该是用来定义单元是协调元还是非协调元。

prolee

都是高手啊 今天学到东西了

chuman

If the "Not Enforced" option is selected, gaps or overlaps will be allowed along inter-element boundaries. These elements are formulated using an assumed linear stress field. These elements are most effective as low aspect ratio rectangles. If the "Enforced" option is selected, overlaps or discontinuities will not be allowed along inter-element boundaries. These elements are formulated using an assumed linear displacement field. These elements can overestimate the stiffness of the structure. In general, a greater mesh density in the direction of the strain gradient is required to achieve the same level of accuracy as elements for which the "Not Enforced" option is selected.

nemia

赞楼上。这个问题当初学的时候根本没仔细学。看来要回炉另造了。

cavingao

原帖由 chuman 于 2008-7-3 23:27 发表
If the "Not Enforced" option is selected, gaps or overlaps will be allowed along inter-element boundaries. These elements are formulated using an assumed linear stress field. These elements are most e ...

谁能翻译一下这段解释文字？谢谢！

haige

参考王勖成的书176面 ，  一般自己写程序应该都是协调单元（由于差值函数的线性，单元边界不存在空隙）， 采用非协调单元后，将会增加内部节点的位移项，最后通过自由度的凝聚，将自由度还是缩减到节点上（过程实际上和子模型比较相近）。

采用非协调元 会增加计算量，但是精度会提高。尤其是对纯弯单元。