- 积分
- 0
- 注册时间
- 2013-10-29
- 仿真币
-
- 最后登录
- 1970-1-1
|
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
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|