<p>子程序不會有問題，對照的各種線彈性的例子都是一樣的</p>

各位老師，我用UMAT寫的線彈性的子程序計算結果和ABAQUS不用子程序的結果不一樣是怎么回事？

瀏覽：1579 回答：7

子程序不會有問題，對照的各種線彈性的例子都是一樣的

邀請回答我來回答

全部回答

(4)

默認最新

AnnaP1893767

你好，我抽中的話費是20，不是10元哦

2017年7月31日

評論 2 點贊

藍牙

這要看哪些結果不一樣了，一般載荷位移曲線如果對的上，應力分布和應變分布稍微有區(qū)別不能說你是錯的，我給你舉個例子

我在做UEL中就發(fā)現(xiàn)自己做的單元剛度矩陣和abaqus自帶的不同。后來我找到了原因

For fully integrated first-order isoparametric elements (4-node elements in two dimensions and 8-node
elements in three dimensions) the actual volume changes at the Gauss points are replaced by the average volume change of the element. This is also known as the selectively reduced-integration technique, because the order of integration is reduced in selected terms, or as the technique, because the order of integration is reduced in selected terms, or as the technique, since the strain-displacement relation (β-matrix) is modified. This technique helps to prevent mesh locking and, thus, provides accurate solutions in incompressible or nearly incompressible cases: see Nagtegaal et al. (1974). In addition, Abaqus uses the average strain in the third (out-of-plane) direction for axisymmetric and generalized plane strain problems. Hence, in the two-dimensional elements only the in-plane terms need to be modified. In the three-dimensional elements the complete volumetric terms are modified. This may cause slightly different behavior between plane strain elements and three-dimensional elements for which a plane strain condition is enforced by boundary conditions.

也就是說在ABAQUS中為了避免單元的自鎖現(xiàn)象，它采用了用平均體積變形代替真實體積變形，對標準的平面應變進行了修正！這就是為什么推導出來的單元剛度矩陣和UEL中的剛度矩陣不同的原因了！
也可以得到這樣的結論，在涉及到單元受到彎矩作用的情況,abaqus 的一階4節(jié)點的等參單元要比ansys的plan42 option(2)=1的結果更接近真實情況！

在umat中也可能出現(xiàn)這樣的情況，比如塑性更新，abaqus為了提供收斂性，也許不是直接采用徑向返回方法，而是采用了一些處理方法，這可能對后期應力和應變的外推磨平等有影響！

不管怎么樣，宏觀的載荷位移曲線應該要對的上

2017年1月7日

已采納評論 1 點贊