二次开发的单元 后处理能显示图形码？

lishhit

询问的具体内容如题？

蓝牙

这个问题要分开回答：
1.任何自定义开发的单元都是可以采用tecplot软件，通过积分点处的应力应变信息绘制云图的
2.如果你开发的单元的自由度，积分点信息和单元库中某单元的一样，是有办法不用tecplot等第三方软件的。你可以详细看看单元手册，我目前因为时间的关系还没有具体琢磨这个事，我这里有篇文章，它就这样做的，你试试看

蓝牙

文章中有这样一句话，你搜索之
where the file “Out.elements.inp” is a copy of all the element definitions just with a high number(here 1000000) added to the element numbers. Thereby, it is possible to make visualization in astandard way in the ABAQUS/VIEWER

lishhit

谢谢 两位 大侠

kensington28

UEL定义的单元不能在后处理显示，但是可以用ABAQUS本身自带的单元覆盖你的自定义单元，ODB文件便可以在abaqus viewer显示。

蓝牙

5# kensington28
你能做个详细的例子吗？教程，我可以给技术积分

lr00011

大哥们，不用那么麻烦。ABAQUS也提供了后处理子程序，UVARM。这个程序还比较简单，如果你会用户单元子程序的话，编这个程序就是小意思。程序变量跟UMAT和uel差不多。
subroutine uvarm(uvar,direct,t,time,dtime,cmname,orname,
     1 nuvarm,noel,npt,layer,kspt,kstep,kinc,ndi,nshr,coord,
     2 jmac,jmatyp,matlayo,laccfla)
只要通过Abaqus utility routine GETVRM 获得积分点的数据就可以了。在ABAQUS培训资料上就有一个用这个子程序作用户单元自程序的后处理的。

kensington28

6# 蓝牙

我对帮助文档verification uel中的T3D2那个例子做了适当的修改。这里的用户单元的刚度矩阵是直接输入的，所以没有用到for文件。 本例子在用户单元的基础上，覆盖了T3D2单元，因此生成的odb文件可以在CAE里打开。

但是我发现在CAE里打开的odb计算结果属于ABAQUS本身自带的T3D2单元的，因此要真正对自定义单元的odb计算结果进行处理，估计还有困难。

我把inp文件，Python读取结果位移的文件传上来，希望大家能一起讨论，提出有效的方法。

蓝牙版主，我目前在做ABAQUS UEL单元的二次开发，想和你交个朋友，一起讨论，我的QQ：156516106， 方便的话，请加我，谢谢~~

蓝牙

7# lr00011

朋友，如果你能做一个4节点双线性的UEL，在abaqus中显示的例子
我给你加精，并给技术积分

4节点双线性的UEL在下面链接

http://forum.simwe.com/viewthread.php?tid=920967&extra=&page=1

lr00011

9# 蓝牙
大哥，我看了下你编程序，程序不知道对不对，细节的东西我没跟有限元公式比较。有些东西可能需要改进，比如*User element, nodes=4, type=U1001, properties=3, coordinates=2, variables=7改成variables=28就更好，因为你有4个积分点，4*7=28，对了你在程序里也说明了，但是却依然7，可能没改好。这个跟你最后输出应力应变有跟关系，你只输出总共7个变量，我没看明白，你到底输出那个点的应力应变，还是平均值，可能少一个do循环，我也不清楚。
DO I=1,3
DO J=1,8
SSTRAIN(I)=SSTRAIN(I)+B(I,J)*U(J)
ENDDO
ENDDO
DO I=1,3
DO J=1,3
SSTRESS(I)=SSTRESS(I)+D(I,J)*SSTRAIN(J)
ENDDO
ENDDO
用户单元子程序如果做变形方面后处理是很简单，直接skin一个标准单元就可以，然后材料属性设很小就可以。他的变形肯定跟你的开发单元子程序一样的。如果做场变量的话，你需要计算他的场值。比如应力应变等随便的场量，而这里显然你的应力应变没计算好，如果你把各个积分点的所关心的场值都算好了，我可以帮你写个显示场值的云图。当然基于现在程序也能作应力应变的云图，但是各个积分点的值都一样，也没什么意思吧，毕竟是错的。

蓝牙

是我疏忽了，上个衔接是我最早的线弹性的例子，当初主要是重点在于单元刚度矩阵和残余力，应变的计算没有问题，这个我验证过了，但是SVARS部分我做的不完整！


http://forum.simwe.com/thread-938271-1-1.html
这个是我前一段时间编的一个弹塑性的UEL，SVARS共有9*4个变量，采用的是平面应变单元，9个变量中，前4个为塑性应变分别是eps（11），eps（22），eps（33）eps（12）
4个流动应力sigma（11），sigm（22），sigm（33） sigm（12）和一个等效塑性应力
请用这个inp和for
还需要修改什么，你跟帖

蓝牙

欢迎其他会员也试着做下这个，这个问题可能是比较普遍的问题！

lr00011

改了一下，根据4节点双线性的UEL做了一个UVARM。好像uel应力计算错误，但是不影响显示，可能更能说明UVARM通用性。变量说明：
DO I=1,3
SSTRAIN(I)=ZERO
SSTRESS(I)=ZERO
ENDDO
DO I=1,3
DO J=1,8
SSTRAIN(I)=SSTRAIN(I)+B(I,J)*U(J)
ENDDO
ENDDO
CALL F_MATMUL(SSTRAIN,B,3,8,U,8,1)
CALL F_MATMUL(SSTRESS,D,3,3,SSTRAIN,3,1)
DO I=1,3
temporary_array(I)=SSTRAIN(I)
temporary_array(I+4)=SSTRESS(I)
ENDDO
temporary_array(4)=EMU*(SSTRESS(1)+SSTRESS(2))
DO I=1,nsvint
SVARS(nsvint*(Kintk-1)+I)= temporary_array(I)
user_variable(jelem,I,Kintk)=temporary_array(I)
END DO
有兴趣的朋友可以自己去运行玩玩。

蓝牙

13# lr00011

谢谢你的工作，先给你奖励，我还没有来得及尝试

你能不能重开一贴，把UVARM的基本原理，实现步骤，以及算例
做成完整的一个教程，如果你时间比较紧简单的介绍下也行


我相信你的努力是对论坛和部分使用UEL的会员非常是有价值的！

谢谢
工作完成后请跟帖给出链接！

lr00011

随便写了点。http://forum.simwe.com/viewthread.php?tid=941826&extra=