JasonMLC 发表于 2013-11-3 17:22:41

求解答--六面体八节点子程序UEL与ABAQUS提供的C3D8结果不同

Hello everyone，
我从九月份开始转向有限元以及ABAQUS。因为是初学者，所以在编写复杂的UEL程序之前，我想编写一个六面体八节点单元的UEL程序与ABAQUS提供的C3D8单元进行比较。我发现UEL算出的结果与利用C3D8的结果不一样，但又不知问题出在哪里。所以请求各位i大神给予指导。
所附文件为：
1.理论推导
2. SD11.inp 与Purely_Elastic_UEL.for
3. Job-1.inp
文件说明：
(1)理论推导是为了方便看UEL程序
(2)SD11.inp是Purely_Elastic_UEL.for的配套文件
(3) Job-1.inp是利用ABAQUS提供的C3D8单元而生成的inp文件。
我所计算的对象是一个六面体八节点单元受一个集中力作用。运行Job-1.inp可以输出刚度矩阵；而运行 SD11.inp 与Purely_Elastic_UEL.for可以在dat文件中看到AMATRX，RHS。两种方式得到的刚度矩阵不一样，不知是什么原因。
此外，在网上也有哪位仁兄提供的UEL程序，我的和他的结果相同。有一点需要注意的是，我在编写UEL时应力分量没有按照ABAQUS应力分量的储存约定(11,22,33,12,13,23)。我采用的是11,22,33,23,13,12。但这并不影响AMATRX与得到的节点位移，节点的力。
Waiting for your reply




amani 发表于 2013-11-4 16:06:52

本帖最后由 amani 于 2013-11-4 16:10 编辑

B-bar单元，斯坦福大学huges院士提出来的，采用应变增强的模式，非常经典的一种处理体积自锁现象的单元，针对不可压缩材料。

hillyuan 发表于 2013-11-3 17:59:24

C3D8 in Abaqus is in fact B-bar element, which decrease volume-locking of general hex element. You need modify your (0.1.8-1)

hillyuan 发表于 2013-11-3 18:00:18

C3D8 in Abaqus is in fact B-bar element, which decrease volume-locking of general hex element. You need modify your (0.1.8-1)

JasonMLC 发表于 2013-11-3 19:42:38

hillyuan 发表于 2013-11-3 17:59
C3D8 in Abaqus is in fact B-bar element, which decrease volume-locking of general hex element. You n ...

可以详细解释一下吗？我不太明白你讲的意思。

JasonMLC 发表于 2013-11-3 20:01:51

hillyuan 发表于 2013-11-3 17:59
C3D8 in Abaqus is in fact B-bar element, which decrease volume-locking of general hex element. You n ...

看了你的回复，我在网上也搜到一点内容，上面也讲了一些内容。
http://www.researchgate.net/post/Stiffness_matrix_of_abaqus_element_c3d8_and_by_my_8_noded_isoparametric_element_is_not_matching-can_anyone_help
http://web.mit.edu/calculix_v2.4/CalculiX/ccx_2.4/doc/ccx/node18.html
我想问一下，对于自己编的六面体八节点等参单元，它适用于哪些情形？

hillyuan 发表于 2013-11-4 11:38:58

JasonMLC 发表于 2013-11-3 20:01
看了你的回复，我在网上也搜到一点内容，上面也讲了一些内容。
http://www.researchgate.net/post/Stiffn ...

C3D8 in calculix is not B-bar element. You can refer to http://www.mycrisp.com/docs/BBar_explained.pdf, for example.

JasonMLC 发表于 2013-11-4 20:22:22

hillyuan 发表于 2013-11-4 11:38
C3D8 in calculix is not B-bar element. You can refer to http://www.mycrisp.com/docs/BBar_explained ...

I see. Thank you very much.

JasonMLC 发表于 2013-11-4 20:23:58

amani 发表于 2013-11-4 16:06
B-bar单元，斯坦福大学huges院士提出来的，采用应变增强的模式，非常经典的一种处理体积自锁现象的单元，针 ...

嗯嗯，非常感谢。

daofman 发表于 2015-2-7 01:37:45

学习下，发现做3维的比较少哦

Gatsbywong 发表于 2015-3-6 14:13:01

B-Bar单元确实很好用

旧日的轨迹 发表于 2016-8-23 09:44:03

你好，能麻烦你看一下帖子吗？
http://forum.simwe.com/forum.php?mod=viewthread&tid=1244996&extra=page%3D1%26filter%3Dtypeid%26typeid%3D127
我和你做了相同的工作，编了六面体8节点等参单元与C3D8单元进行了对比，计算结果变形量相差很大，是因为单元数目过多的原因吗？

kanglongwuhui 发表于 2018-12-19 20:36:37

amani 发表于 2013-11-4 16:06
B-bar单元，斯坦福大学huges院士提出来的，采用应变增强的模式，非常经典的一种处理体积自锁现象的单元，针 ...

老师，您好，在学习，很苦恼，能不能和您交流一下，qq2551945729，真心感谢，我没法给您留言，希望您能看得到

123456qq 发表于 2021-6-5 14:17:04

66666666666666

599705369 发表于 2021-6-5 21:30:26

怎么上传图片
页: [1]
查看完整版本: 求解答--六面体八节点子程序UEL与ABAQUS提供的C3D8结果不同