清华大学bbs的abaqus精华.docx
- 文档编号:14820142
- 上传时间:2023-06-27
- 格式:DOCX
- 页数:22
- 大小:52.08KB
清华大学bbs的abaqus精华.docx
《清华大学bbs的abaqus精华.docx》由会员分享,可在线阅读,更多相关《清华大学bbs的abaqus精华.docx(22页珍藏版)》请在冰点文库上搜索。
清华大学bbs的abaqus精华
清华大学bbs的abaqus精华
大变形下此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说errorseemsg,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水木清华站(ThuDec2714:
19:
272001)
用uvar()勉强成功
标题:
在abaqus中,DSTRAN,STRAN为工程剪应变吗?
SatDec2909:
44:
032001)
同题,这涉及到ddsdde,peeq的计算
对,我也是最近才发现的
peeq怎么计算?
在umat中计算peeq,(因用umat后,系统缺省peeq=0)然后再uvarm输出
uvarm可以输出到哪里?
.odb可以么?
另外那个strav?
?
就是自己定义用于umat的那个数组里的数可以输出到.odb里么?
谢了,你这个成长中的大侠
在umat中,statev是不能用在output中的,statev只是作为一个解的
状态变量,说来惭愧,我是在umat中定义peeq,(peeq)的求解一般在
弹塑性力学书上有(等效塑性应变),用write()写入一个临时文件,
((切记:
这个文件unit号不要与abaqus中的重合,因为他有一些系统默认
的文件号,))
然后在uvarm中读取,以uvarm输出,因为uvarm可以以odb的形式输出,
支持output,field,output,history
不知谁还有更好的方法。
听说清华有人abaqus用的很熟了。
大家过奖了,我是没办法,要出论文,总得学习,上午版主还在qq中指点了我一番
YOUCANTRYASFOLLOWS:
*STEP,EXTRAPOLATION=NO,INC=2000000
*STATIC
0.001,500.0,0.00001,0.1
可以
并且在step里好像有选项选择
你是手写的程序么
新建一个step
然后定义边条
在上一step定义边条的节点上重新定义边条
而且,好像,是不是有个选项是new
这样就可以取消上一step的边条设置了
在以后的载荷步里重新定义边界条件,加参数new
比如在第二个载荷步里:
*STEP
......
......
*BOUNDARYCONDITIONS,OP=NEW
......
*ENDSTEP
Command方式:
abaqusjob=job_nameuser=source_file(或者obj_file);
CAE方式:
在最后submit的时候可以选择source_file(或者obj_file).
不太明白你的意思。
如果只是要得到所需要的如应力,应变等,那么在数据行请求
输出就行了啊,而且个人觉得abaqus的后处理功能基本上也够用了,直接输出到*.odb文件
中,用CAE或者Viewer后处理即可。
而*elprint输出的dat文件包含其他的信息,如版本号、处理内容等。
若想要直接输出到其他后处理软件的话,可以输出到resultfile文件*.fil,不过,
这可是二进制文件或者ASCII文件,还是需要处理的,或者直接用Hypermesh等后处理
程序打开做图。
好像是在tool->display->create->再选定一个单元集或者是节点集
然后操作是“取代”还是其他等等(就是里面的两个圆圈相交,
还是想减)。
你试试看看。
这个是要到后处理才行吧
非要计算以后形成*.odb吗
前处理里面怎么才可以啊
Python在abaqus中可以运行scrpit文件,你在CAEmanual中就
可以将其附录中的例子存为*.py文件,然后在CAE的菜单中
runscript就可以运行了。
BTW:
请问你的是abaqus的那个版本?
6.2or6.3?
6.3的有盗版了?
:
(1),请问如何得到M,C,K矩阵?
用什么命令?
试一试:
*elementmatrixoutput
:
(2),想要用Newmark方法求解.求解问题,在ABAQUS中如何实现?
或者
:
那里能找到类似的例子?
看看:
standardmanuals:
analysis的dynamicsanalysis和explicitmanuals
:
(3),ABAQUS中能实现这样的东西吗?
就是我需要平滑ABAQUS产生的位移场
:
还是这个平滑只能在其他环境中完成?
试试这个输出选项:
*Elprint,Position=averageatNodes
:
abaqus/CAE好像不能生成节点集,只能定义某个实体,后处理用这个实体上的节点。
:
看到CAE的例题都是这样做的。
NODESET是abaqus的基本功能,也是其方便使用之处,
在CAE中是支持的。
这个例子我做过,没有问题。
:
:
我用tool》set》create居然不能选择由CAE生成的节点(不是顶点)
这是因为你没有正确操作PartitionEdge这一步。
【在sunhaolan的大作中提到:
】
:
abaqus/CAE好像不能生成节点集,只能定义某个实体,后处理用这个实体上的节点。
ABAQUSdividesallofitsanalysisproceduresinto
twomaingroups:
linearperturbationandgeneral.
ABAQUSmakesaverycleardistinctionbetween
perturbationandgeneralanalysisproceduresbecause
loadingconditionsand"time"aredefineddifferently
forthetwocases.Furthermore,theresultsfromeach
typeofprocedureshouldbeinterpreteddifferently.
Theresponseofthemodelduringageneralanalysis
procedure,knownasageneralstep,maybeeither
nonlinearorlinear.Inastepthatusesaperturbation
procedure,whichiscalledaperturbationstep,
theresponsecanonlybelinear.ABAQUStreats
suchstepsasalinearperturbationaboutthepreloaded,
predeformedstate(knownasthebasestate)createdby
anypreviousgeneralsteps;therefore,itscapability
fordoinglinearsimulationsisrathermoregeneral
thanthatofapurelylinear
analysisprogram.
各位大侠,我在定义step的时候,前三个用的是standard中的相关命令.在
step4的时候,我想用explicit来求解,并且定义了amplitude,来模拟一个
位移随时间的简单变化,命令如下:
*stepname="~~~",nlgeom,inc=1000(好像explicit里没有inc这个参数)
*dynamic,explicit
2
*amplitudename=ramp
0,0,1,1,2,1
*boundaryamplitude=ramp
refpunch(刚体参照点)2,2,-0.03
*endstep
有问题吗?
另外standard的命令和explicit的命令能出现在同一个.inp文件中吗?
要转换的.
ABAQUS/StandardandABAQUS/Explicitareseparateprogrammoduleswith
differentdatastructures;hence,the
explicitdynamicsprocedurecannotbeusedinthesameanalysisasanyof
theproceduresinABAQUS/Standard.
However,ABAQUSprovidesacapabilitytoimportadeformedmeshandassociated
materialstatefrom
ABAQUS/ExplicitintoABAQUS/Standardandviceversa.Thisprocedureis
describedin``Transferringresultsbetween
ABAQUS/ExplicitandABAQUS/Standard,''Section7.6.2oftheABAQUS/Standard
User'sManual.
【在linyuwho(苍龙有悔)的大作中提到:
】
:
各位大侠,我在定义step的时候,前三个用的是standard中的相关命令.在
:
step4的时候,我想用explicit来求解,并且定义了amplitude,来模拟一个
:
位移随时间的简单变化,命令如下:
:
*stepname="~~~",nlgeom,inc=1000(好像explicit里没有inc这个参数)
:
*dynamic,explicit
:
2
:
*amplitudename=ramp
:
0,0,1,1,2,1
:
*boundaryamplitude=ramp
:
refpunch(刚体参照点)2,2,-0.03
:
*endstep
:
...................
feature菜单下有删除一项。
基准点是一种feature
:
工具栏不有个像铅笔一样的的东东么
我想把计算出来的应力,应变,位移等数据提取出来,在其他软件里处理,
请问该怎么办?
我在。
dat里没有找到这些数据阿?
。
res不是文本格式,
其他的好像也没有啊。
计算结果,没有文本格式,只能自己把数据倒出来。
在显示结果的模块中,名字忘了,
就是sketch前面那个,
在report菜单中,可以到处自己需要的文件。
文件格式是:
rpt
谢谢,但是这个为什么只有应力和应变呢?
计算出的位移的数据怎么会没有呢?
【在santong(supernova)的大作中提到:
】
:
计算结果,没有文本格式,只能自己把数据倒出来。
:
在显示结果的模块中,名字忘了,
:
就是sketch前面那个,
:
在report菜单中,可以到处自己需要的文件。
:
文件格式是:
rpt
发信人:
wjytj(赖子),信区:
FEA
你要在STEP模块fieldoutput和historyoutput改变默认
请求输出位移U
或者在*step中加上:
*NODEPRINT
U
计算完成后在*.DAT文件中就可以提取节点位移。
请问fieldoutput和historyoutput什么区别?
关键字*nodeoutput和*nodeprint有什么区别?
谢谢!
我的。
inp数据文件里有这么几行:
*Output,field,frequency=99999
*NodeOutput
U,RF
*ElementOutput
S,E,EE
*Output,history,variable=PRESELECT
*ElPrint,freq=999999
*NodePrint,freq=999999
其中在historyoutput中有*nodeprint关键字,但是没有指定输入U,
最后导出的数据没有U这一选项,请问是什么原因?
【在wjytj(赖子)的大作中提到:
】
:
你要在STEP模块fieldoutput和historyoutput改变默认
:
请求输出位移U
:
或者在*step中加上:
:
*NODEPRINT
:
U
:
计算完成后在*.DAT文件中就可以提取节点位移。
发信人:
Mikejwg(mikejwg),信区:
FEA
fieldoutput和historyoutput这在viewer模块中很明显。
fieldoutput是场量输出,historyoutput为历史数据输出,会记录场量随时间的曲线
至于*nodeoutput与*nodeprint的区别在于他们写入不同的后缀文件,
比如*.odb,*.dat等,在help中有详细介绍。
【在
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 清华大学 bbs abaqus 精华