清华大学bbs的abaqus精华Word文档格式.docx
- 文档编号:15246634
- 上传时间:2022-10-28
- 格式:DOCX
- 页数:16
- 大小:41.35KB
清华大学bbs的abaqus精华Word文档格式.docx
《清华大学bbs的abaqus精华Word文档格式.docx》由会员分享,可在线阅读,更多相关《清华大学bbs的abaqus精华Word文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
Sij=DijklEkl(满足张量求和约定)
即Stress=D*Strain;
而在增量理论中,
△S=D*△E(在有限变形中,△其实应该为应力的客观率)
--
似乎不对吧
大变形下此D非彼D ﻫ你看过黄克智的固体本构关系这本书么ﻫ如果你从全量理论和增量理论的角度上讲 ﻫ那似乎第一个Digkl就不对
你有第一个式子么
如果有,求导不久完了?
那应该就是一样的,因为全量理论,ﻫ:
Sij=DijklEkl(满足张量求和约定)ﻫ:
即Stress=D*Strain;
ﻫ:
而在增量理论中,
△S=D*△E(在有限变形中,△其实应该为应力的客观率)
是啊,大变形下的[D]与普通意义下的[D]在构型上是不一样的,ﻫ毕竟[D]大是变形历史的函数,而[D]小则不是,ﻫ我推导一种新的本构关系,
△Sij=Dijkl△Ekl(其中△为Jaumann率)ﻫ假设材料一开始就屈服(即屈服面为0)ﻫ想用壳单元,
:
似乎不对吧 ﻫ:
大变形下此D非彼D ﻫ:
你看过黄克智的固体本构关系这本书么
如果你从全量理论和增量理论的角度上讲 ﻫ:
那似乎第一个Digkl就不对
你有第一个式子么
如果有,求导不久完了?
是啊,大变形下的[D]与普通意义下的[D]在构型上是不一样的,
毕竟[D]大是变形历史的函数,而[D]小则不是,
我推导一种新的本构关系,ﻫ:
△Sij=Dijkl△Ekl(其中△为Jaumann率)ﻫ:
假设材料一开始就屈服(即屈服面为0)
想用壳单元,ﻫ唉,别提了,问题就出在,在abaqus中,明明写着可以考虑剪切效应,
可我打印出剪切力个数是,nshr=1,即只有S12,
那我的S13,S23就不知怎么计算,
(DDSDDE(5,5)无法计算,因为ntens=3,最多只能计算DDSDDE(3,3))ﻫ
你编umat编进去不久行了
他让用NDI,NSHR,NTENS表示变量,你就用这些表示变量
这样他就可以任意的计算了呀,而不在意实际计算的变量数
ﻫ:
唉,别提了,问题就出在,在abaqus中,明明写着可以考虑剪切效应,ﻫ:
可我打印出剪切力个数是,nshr=1,即只有S12,
那我的S13,S23就不知怎么计算,
(DDSDDE(5,5)无法计算,因为ntens=3,最多只能计算DDSDDE(3,3))ﻫ因为我的UMAT从abaqus传来的变量(当我选用壳单元时)ntens=3,ﻫndi=2,nshr=1,这样的话,STRESS为3个,STRESS(3),不可能计算
STRESS(4),STRESS(5),同理,STRAIN,DDSDDE也存在同样的情况
在中厚度板壳元中(MINDLIN)中,DDSDDE为5×
5 ﻫ
你的FOR文件收到,可惜研究的对象相差太大,看不懂!
不过你在文件中定义的变量有的根本就没有用到,
还有有的地方似乎是毫无用处的代码,比如:
SMISES一段
我觉得你不应该太急,第一步应该保证代码的正确性。
ﻫ另外看了你上面的讨论,我觉得你应该换一种单元试一试。
ﻫ比如用梁单元。
很有可能NDI和NSHR会不一样。
不知道你是不是一定要用到壳单元上?
如果必须用在壳单元上,可能S23,S31就必须自己在UAMT ﻫ中输出了。
具体的做法我也不清楚。
还有对DDSDDE,我的理解也是:
dS=DDSDDE*dE
ﻫ我的意思是不管他是几乘几的矩阵ﻫ一般的D矩阵就是6*6
你编的umat需要适应于任何单元
所以当然6*6的矩阵应该给完全
当算体的时候
如果是壳单元,他只有ntens=3
那就等于三呗,他用哪个就是哪个 ﻫ反正你的D6*6的矩阵放在哪了ﻫ如果是三维实体单元,ntens=6ﻫ所以你的D矩阵就全用上啦
记住你编的umat应该适应于所有的单元,所以不用考虑哪几个参数的值 ﻫ而且这个stress(3)也许是人家只用到了你D6*6的第一行,第二行,和的四行哪ﻫsigh ﻫ不知道我说明白了没有ﻫ:
因为我的UMAT从abaqus传来的变量(当我选用壳单元时)ntens=3, ﻫ:
ndi=2,nshr=1,这样的话,STRESS为3个,STRESS(3),不可能计算
STRESS(4),STRESS(5),同理,STRAIN,DDSDDE也存在同样的情况ﻫ:
在中厚度板壳元中(MINDLIN)中,DDSDDE为5×
5
ﻫ对于DDSDDE,其实就是给了一个应力应变的斜率,我感觉
真正计算的是利用更新的应力,DDSDDE我觉得就是为了确保收敛ﻫ真开心有这么多人讨论ABAQUS,特别是umatﻫ我国外的师兄说Umaterial在国外应用得也很广
对于DDSDDE,其实就是给了一个应力应变的斜率,我感觉 ﻫ:
真正计算的是利用更新的应力,DDSDDE我觉得就是为了确保收敛ﻫ:
真开心有这么多人讨论ABAQUS,特别是umat ﻫ:
我国外的师兄说Umaterial在国外应用得也很广ﻫﻫ各位大侠,看到我提出的uamt,及statev有这么多人感兴趣,实在很ﻫ高兴,希望大家能把摸索到的东西介绍给大家,现在我的程序能通了,
(因没有屈服面的概念),虽然结果是错误的,但总可以静下心来,不像以前
一运行自己的umat程序就standard.exeerror.(并且程序.dat.msg.中无
任何错误信息,log说error seemsg,butmsg中无)。
ﻫ现在我明白了,在用户程序中,WRITE(*,*)则写入log,write(6,*)则写入
dat,等,
在监控程序时,发现,statev(),是每次迭代和增量运行它,他就为0, ﻫ还是为上一步的值,因为我每次都为0,怀疑有误ﻫ还望大侠提醒!
dat,等,
在监控程序时,发现,statev(),是每次迭代和增量运行它,他就为0,
应该不是0,如果在每个增量步的最后你重新更新过他ﻫ输入文件*.inp中,*depvar不要忘了写statev()的维数哦!
还是为上一步的值,因为我每次都为0,怀疑有误
还望大侠提醒!
ﻫ--
你只有一个单元么?
ﻫ【在saturn(好好工作) 的大作中提到:
】
~!
~~~~~~~~~~~~~~~~~~~~~~~~~~~~这我以前还不知道ﻫ:
应该不是0,如果在每个增量步的最后你重新更新过他
输入文件*.inp中,*depvar不要忘了写statev()的维数哦!
--ﻫ【 在 air1(无影如风※勤奋羊羊※课题在即)的大作中提到:
】ﻫ:
可是这个写的命令有用么
对于输出而言ﻫ:
你只有一个单元么?
ﻫ已经写了,定义了statev的维数,ﻫ*depvarﻫ20,
ﻫ看来我的错误不在statev,
因为我的输出显示在第三增量步中,statev不为0,
但我还是有些不清楚uamt程序的执行过程,ﻫ比如,
一板一次壳元s4r,面内采用缩减积分,即1×
1,厚度方向为5点simpson
积分,按道理,
没一增量步的迭代步,对每一积分点都会运行一次umat,ﻫ即上问题,每一迭代该运行5次,ﻫ那上面的statev,是每次迭代都会改变(不是每次运行)?
请air1,beauty_life,saturn,......等各位相助
ﻫ看过书上例子,
如果定义statev为(nblock,*)的数组, ﻫ应该对每个积分点都进行迭代计算,从而更新变量值
(写stateNew(nblock,*)ﻫ请问Mikejwg,变量statev怎么输出?
可以象应力应变那样进行绘图、动画等后处理吗?
ﻫ不好意思,我只在log文件中输出,还没有在odb中输出,ﻫ想同elementoutput等那样输出,只能在UVARM中定义,但
同statev是不同的,老兄你可以试试。
ﻫ
发信人:
Mikejwg(mikejwg), 信区:
FEAﻫ标题:
用abaqus中的UMAT的大侠看过来。
。
ﻫ发信站:
BBS 水木清华站(TueDec2515:
36:
272001)
各位大侠,ﻫ在加入自己弹塑性本构模型时,如何加入区分弹性部分与塑性部分(或
非弹性部分)。
也就是说,在输出时,EE,LE,IE([PE]),因为在umat中,
我们只发现dstran,stress,除此并没有发现其它有关应力应变得变量。
ﻫ
你所说的区分弹性部分和塑性部分是指应变的弹塑性部分:
E=EE+PE吗?
还是指加载过程的弹塑性?
E=EE+PE
我觉得想从umat的输入参数里面区分是不可能的ﻫ可不可以你先算出弹性的,然后用stran减去弹性的,得到塑性的
我猜的哦ﻫ
我自己是明白的,但不知怎么让程序明白 ﻫ比如,在一定条件下可解出位移增量,abaqus得到应变增量dstran,stran, ﻫ而弹性部分,dEEij=dSij/(2G),dEPij=dstranij-dEEij
但怎么控制程序了?
ﻫ默认得情况又是怎样计算的
谢谢ﻫﻫ我不清楚你的本构,
我大概觉得你如果是用屈服面的话ﻫ你是不是可以用statev()记录每一步的EE,EP
然后通过应力偏量来计算DEE,DEP,然后自己更新EE,EP,就不与stran扯上关系ﻫ我以前见过用Umat算J2流动理论的时候,是这么做的
首先我的本构不能采用屈服面的概念(一种新型的弹塑性本构模型), ﻫ我也看了那个采用j2流动法则的例子,我可以在umat中的statev重定义
ee,pe,但statev不能像EE,PE,IE那样在odb中输出。
再说在umat中的ee,pe与inp中的ee,pe是不同的。
谢谢 ﻫ
看来你只能试试uvarm了ﻫ呵呵, ﻫ
Mikejwg(mikejwg),信区:
FEA
标题:
Re:
ABAQUS中的弹塑性应变的计算问题?
?
发信站:
BBS水木清华站(ThuDec 2714:
19:
27 2001) ﻫ用uvar()勉强成功ﻫ
标 题:
在abaqus中,DSTRAN,STRAN为工程剪应变吗?
ﻫSatDec29 09:
44:
032001)
同题,这涉及到ddsdde,peeq的计算ﻫﻫ对,我也是最近才发现的 ﻫ
peeq怎么计算?
ﻫﻫ在umat中计算peeq,(因用umat后,系统缺省peeq=0)然后再uvarm输出ﻫ
uvarm可以输出到哪里?
.odb可以么?
ﻫ另外那个strav?
?
就是自己定义用于umat的那个数组里的数可以输出到.odb里么?
ﻫ谢了,你这个成长中的大侠
在umat中,statev是不能用在output中的,statev只是作为一个解的
状态变量,说来惭愧,我是在umat中定义peeq,(peeq)的求解一般在 ﻫ弹塑性力学书上有(等效塑性应变),用write()写入一
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 清华大学 bbs abaqus 精华
![提示](https://static.bdocx.com/images/bang_tan.gif)