版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、Good is good, but better carries it.精益求精,善益求善。abaqus问答浙江大学Cambridge-bbs浙江大学bbs采用abaqus的cae进行力学问题的分析,其对模型的处理存在很多的技巧,对abaqus的一些分析技巧进行一些概述,希望对大家有所帮助。1.abaqus的多图层绘图abaqus的cae默认一个视区仅仅绘出一个图形,譬如contor图,变形图,xy曲线图等,其实在abaqus里面存在一个类似于origin里面的图层的概念,对于每个当前视区里面的图形都可以建立一个图层,并且可以将多个图层合并在一个图形里面,称之为OverlayPlot,譬如你可
2、以在同一副图中,左边绘出contor图,右边绘出xy图等等,并且在abaqus里面的操作也是很简单的。1.首先进入可视化模块,当然要先打开你的模型数据文件(。odb)2.第一步要先创建好你的图形,譬如变形图等等3.进入view里面的overlayplot,点击creat,创建一个图层,现在在viewportlayer里出现了你创建的图层了4.注意你创建的图层,可以看到在visible下面有个选择的标记,表示在视区里面你的图层是否可见,和autocad里面是一样,取消则不可见current表示是否是当前图层,有些操作只能对当前图层操作有效,同cadname是你建立图层的名称,其他的属性值和你的模
3、型数据库及图形的类型有关,一般不能改动的。5.重复2-4步就可以创建多个图层了6.创建好之后就可以选择plot/apply,则在视区显示出所有的可见的图层1.什么是子结构子结构也叫超单元的(两者还是有点区别的,文后会谈到),子结构并不是abaqus里面的新东东,而是有限元里面的一个概念,所谓子结构就是将一组单元组合为一个单元(称为超单元),注意是一个单元,这个单元和你用的其他任何一种类型的单元一样使用。2.为什么要用子结构使用子结构并不是为了好玩,凡是建过大型有限元模型的兄弟们都可能碰到过计算一个问题要花几个小时,一两天甚至由于单元太多无法求解的情况,子结构正是针对这类问题的一种解决方法,所以
4、子结构肯定是对一个大型的有限元模型的,譬如在求解非线性问题的时候,因为对于一个非线性问题,系统往往经过多次迭代,每次这个系统的刚度矩阵都会被重新计算,而一般来说一个大型问题往往有很大一部分的变形是很小的,把这部分作为一个子结构,其刚度矩阵仅要计算一次,大大节约了计算时间。3.哪些情况可以使用子结构前面提到的非线性问题,包括了很小变形的或者线弹性部分可以使用子结构,特别是当模型中有很多相同的部分时,提到的最多的一个例子就是桌子的四条腿,四条腿作为子结构(因为基本时弹性变形)可以包括了很多的实体单元,可以大大提高效率再一个就是问题确实太大,只有采用子结构将问题分成很多块,计算出结果后再次采用子结构
5、分块计算,一直到能对每块单独计算为止。4.abaqus中子结构的特点及要注意的问题子结构是一组单元的集合,但是在子结构中仅仅只有你指定的那些节点的自由度会保留下来而其他节点的自由度都被消除了,其他的节点均是通过线性插值的方式获得求解;子结构是通过你指定的节点与其他的单元建立联系的;在abaqus的6.4版本中只有子结构这个概念没有超单元了,其区别就是子结构可以求得单元内部准确的解而超单元不行;当你定义子结构的时候不要包含太多的单元,因为单元的刚度矩阵集成的时候会花掉太多的时间,可以用更多的含有较少单元的子结构代替。5.在abaqus中子结构的用法一般包括如下部分,首先你要利用substruct
6、uregenerate和retainednodaldofs定义一个子结构,然后你可以定义子结构的内部荷载,边界条件一旦你已经定义了子结构以后你就可以象利用一般的单元一样使用子结构了,譬如输出请求等等。因为内容太多只能另外的文章再述了,大家也可以直接看看手册,要注意的是,abaquscae是不支持子结构的。如何在不同的分析步改变材料的参数我所了解的大概有三种不同的方法:1.最强大的当然是采用umat的方式,不过需要有深厚的有限元基础,一般人不推荐使用2.采用场变量,不过功能相对简单3.采用abaqus的import命令将前面分析的结果传递到新的分析之中这里介绍下第二种方法*什么是场变量所谓场变量
7、,我的理解就是一个环境变量,它建立了一个与材料参数之间的中介,虽然不能直接指定材料参数在不同的分析步具有不同的值,但是通过场变量,间接的达到了目的。*怎样使用场变量其实场变量用的较多的实在热力学和流体力学的分析种,这里介绍的仅仅是在固体力学中的用法1.定义场变量*你可以在initial中指定场变量的值,格式如下initialconditions,typefield,variablen(场变量的编号)Set1(你定义的结点集),1.0(场变量的值)场变量是通过编号来识别的,一次只能定义一个场变量*你也可以直接在分析步中指定场变量的值,格式如下*field,variable=1Set-1,1当然也
8、可以同时使用initial和field,当你指定的场变量改变时,默认材料的参数是在增量步间线性变化的。2.建立材料参数和场变量之间的联系如果你用cae,在prop模块里面的材料参数一般都有Numberoffieldvariables,场变量都是从1开始的,你也可以选择多个场变量。填入场变量的值和材料参数间的关系,譬如杨是模量泊松比field1200.E9,0.3,1.180.E9,0.3,2.如果你用的是命令格式,则在inp文件里面键入:*ELASTIC,DEPENDENCIES=1200.E9,0.3,1.180.E9,0.3,2.第四个参数表示场变量的值3.注意场变量在不同分析步中的值有你
9、在不同的分析步中指定,如果没有指定,材料参数默认微是场变量1的值,例如*STEP,name=step1*STATIC*FIELD,VARIABLE=1NALL,1.*ENDSTEP*STEP,name=step2*STATIC*FIELD,VARIABLE=1NALL,2.*ENDSTEPcae步支持场变量,所以你必须自己更改inp文件,4.技巧如果材料的参数变化比较复杂,一般是利用幅值曲线来定义场变量值的变化,*FIELD,VARIABLE=1,amplitude?对幅值曲线步清楚的自己可以看手册*模型的重启动分析restart按理说restart不应该算是一个分析的技巧,而是一个常识,不过
10、呢可能有很多朋友没有建过大型模型导致restart也用的较少,所以也介绍下1.什么是restart你的job可能包含多个step,可是如果你的模型很大,可能会有这样一种情况,当你花了几天几夜,终于分析好的时候,你发现thefirststep的边界条件设置的有问题,这对于你真是晴天霹雳,于是你只好重新来过,可是低二天你发现你的电脑restart,这时的你可能只能问上帝了,howcanido?*restart,就是将一个复杂的模型分析过程分成很多的阶段,甚至是一个increatmentstep一个阶段,你可以对每个阶段的结果进行检验,然后进入下一个阶段进行分析。2.重启动需要那些文件对于stand
11、ard来说,.res,.mdl,.stt,.prt,.odb,这些文件是用于重启动的,explict是.abq,.stt,.prt,.odb.3.如何在一个分析中设置重启动来生成以上文件。这里只介绍下在standard的用法,其实很简单?inp文件里面加入*RESTART,WRITE,FREQUENCY=N就可以了cae默认加入了重启选项,不过可以在stepoutput-restartrequest里面设置输出的频率,也就是frequency。*技巧:因为res文件包含了模型的几乎全部信息,所以非常大,你可以设置overlay参数使后面的数据覆盖吊前面的数据,不过restart的话你也只能从最
12、后一个增量步开始4.如何重启你要指定一个重启点,inp文件里面加上*RESTART,READ,STEP=step,INC=increment就可以了cae中更简单,首先在modeleditattribute里面选择restart,指定前面分析的job名和你想重启动的开始分析步和增量步就可以了,然后在job里面指定重新创建的工作类型,restart,thatsall.5.注意重启动不能改变你的原始分析中的任何参数,也就是说,你的启动点的模型必须和原始分析中的模型完全一致的,所以不要企图采用restart的方法来改变边界条件,材料参数或者网格的密度等等。这些需要另外的技巧来实现。Q:第一次安装ab
13、aqus6.2.1,装完exceed6.2&3D后,开始装6.2.1的Productinstallationfornetworklicensing,结果每到进度为75%时,弹出界面:Anerrorocurredduringthemovedataprocess:-115Component:CompleteInslattationFileGroup:fg_commonFile:c:abaqus6.2-1caeExternalebtadi3xcol.dll就停止安装了,重新安装依然如此,现在感到头都大了,哪位好心人帮帮我,不深感激?A:没有关闭防火墙!Q:WARNING:THESYSTEMMATRI
14、XHAS1NEGATIVEEIGENVALUES.一般在什么情况下会发生系统矩阵出现负特征值?A:很多情况。比如边界条件定义不对;接触定义不对。Q:就是手册的umatmst3.for1)Mises各向同性的子程序中,调用硬化曲线函数中的Table(2,nvalue)这个二维数值保存是应力应变曲线的数据,可是ABAQUS并没有告诉UMAT中Table数值这些应力应变的数值啊,UMAT里面Table怎么知道的呢?2)RHSSMISESS-EG3*DEQPL-SYIELD这个式子是什么意思?尤其是3G*DEQPL是什么意思?A:就是Constant8中的从第三个开始的数值,数组名传递。Q:我在计算时
15、MSG文件出现如下错误,是不是由于节点数太多啦?该如何处理这个问题?ERROR:SPECIFIEDSTANDARD_MEMORYVALUEOF8000000ISTOOSMALLTORUNTHENALYSIS.STANDARD_MEMORYMUSTBEINCREASED.MINIMUMPOSSIBLEVALUEIS23477555.LOOKATMEMORYESTIMATESSECTIONOF.datFILEFORFURTHERINFORMATION?A:对abaqus_v6.env文件中的STANDARD_MEMORY的值进行修改。修改Site文件夹下的abaqus_v6.env中的配置,如下:
16、#System-WideABAQUSEnvironmentFile#pre_memory=33554432standard_memory=33554432#*修改这里*#NTspecificsettings。Q:WARNING:THESYSTEMMATRIXHAS148NEGATIVEEIGENVALUES1ABAQUSVERSION6.3-1DATE27-NOV-2002TIME22:08:00PAGE.ForusebyNoneuserlicensefromHKSInc.STEP1INCREMENT1STEPTIME0.00STEP1STEADYSTATES.AUTOMATICTIMECON
17、TROLWITH-ASUGGESTEDINITIALTIMEINCREMENTOF0.300?A:将你的INITIALTIMEINCREMENT改小些试试看,不过,一般出现此类问题,多半是你的模型有问题,欠约束或者其它什么的。Q:ZEROPIVOT是什么意思?A:zeropivot可以理解为刚度矩阵出了问题,例如奇异。可能有不同的原因,如:你所模拟的是软化性质材料,该点因破坏等原因而软化至不能吃劲;模型有问题,如约束不够,或者是单元拓扑出错等;或者是你的UMAT中写的ddsdde有错;。Q:abaqus-uamt的老问题,缺少ABA_PARAM.INC文件?A:在cvf6.5调试时,显示缺少A
18、BA_PARAM.INC文件!这个没有任何关系的,这个错误将在ABAQUS调用UMAT的时候自动会找到,仅仅有这个错误将没有任何影响的。也就是说,ABAQUS中调用的时候,实际并不存在这个错误。FT,忘了说一句了,你把ABA_PARAM.INC.dp或ABA_PARAM.INC.sp拷到你的程序工作空间后,应该将把ABA_PARAM.INC.dp或ABA_PARAM.INC.sp的后缀.sp或.dp去掉,即将ABA_PARAM.INC.dp或ABA_PARAM.INC.sp改名为ABA_PARAM.INC。呵呵,他的意思是在VisualFortran中调试其子程序,我觉得这是一个好办法,我当时
19、也是这麽办的,毕竟在ABAQUS中调试是非常麻烦的,只有当你的UMAT没有语法或者明显的逻辑错误,你在ABAQUS中调试才能事半功倍。Q:abaqus中有没有布尔运算?能否在abaqus中将两个体粘成一个体?A:可以,用constraint和interaction可以将两个part进行定位的。Q:abaqus里做静力分析时怎么加重力加速度?在inp文件里的card是什么?A:如:*DLOADELSET,GRAV,9.81,0.,0.,1.有疑问查一些*DLOAD卡就清楚了。Q:如何在ABAQUSCAE如何设定在某一节点上受集中力呢请问一下在ABAQUSCAE如何设定在某一节点上受集中力呢?我用
20、CAE中的Set好象不能设定一个节点丫.请问在CAE要如何设定,或则祇能在inp下设定呢?A:用set可以定义一个节点的。另外如果只是在一个节点上加集中力就没有必要用set了,在Load模块里creataload,选择集中力,然后点中该节点(别忘了done)就可以输入力的数值和方向了设定一个节点当然可以了,在cae里一点的集中荷载还是挺容易加了。partition,datumnpoint,都可以。用CAE建的模型好象不能直接加在接点上,不过可以加在vertex上。Q:如何将两个part接在一起?几个人将一个工作分开来做,完成后要将造好型的part接在一起,再加边界条件和载荷。该如何把这几个pa
21、rt接起来?A:先导入File-import-part,然后形成assembly,这样做试试看如何Q:abaqucae中作梁的分析时可否定义不规则的profile形状?只提供了box,pipe,circle,角钢,方钢,T型钢,l型钢,工字钢,我可不可以定义像三角形空心钢管,或者梯形空心钢管这样的profile?A:直接给出I和GQ:在CAE中建立PART时为什么只能先建立在X-Y平面上的模型在CAE中建立PART时为什么只能先建立在X-Y平面上的模型,可不可以直接建立空间的点,结点线,面,体呀A:在图纸上画投影图当然是2维的了。这是它的建模思路,也是一种专利了。直接在空间里建立模型当然也是一
22、种思路,哪个好呢?各有千秋。你要使用CAE,就要熟悉并掌握它的思路,把整体结果分解为各个小part,把每个part都简化到何以在平面上画出投影图并extrude或回转就得到part的空间实体。然后把这些part给instance,指定其section特性(就是其材料呀,厚度呀等在图形上显示不了的数据)装配起来构成复杂的模型。同一个part的草图可以有许多instance,你可以尽量利用装配功能来减少你建立part的工作量。对机械图比较熟悉的,可以想象为制作零件然后装配成机器的过程,对于土木类工程师可能开始会不太习惯,但是多做几个例子就会熟悉并喜欢上这种操作了。如果你实在是喜欢空间直接建模法,你
23、可以用ANSYS建立模型划分单元后再用小程序变成ABAQUS的输入文件。Q:abaqus的所有功能在cae中都能完成吗?还是需要在命令行中输入A:大部分工作可以在CAE中实现,有些功能还需要在INPUT文件中修改,将CAE转化成*.INP,就可。abaqus的前置确实不太好。cae的功能一直没有长进。如板料成形方面,回弹的分析就是不能在cae里完成,必须写inp文件来做。写inp文件的效率可一直不高啊因为一般成形我们最好用动态显示算法,而回弹最好静态隐式算法。不过在cae里就是不行。只要建立了动态算法的step,后面就不能建立静态算法的step。还有啊,abaqus的cae里面是先设置边界条件
24、和输出,然后剖分网格,这个不适合大多数有限元前置处理的思路。很少前置软件是这样的。比起其他的前置软件,如femap,patran,还有工业用软件autoform,pamstamp,差不少啊!Q:如何给一个面的一部分加载A:Areyouusing/CAE?Ifso,pleasepartitionyoursurfaceandapplyloading.对的,将矩形和圆区分开来,将圆定义为一个块就行了!将要加载的面上所有的节点设为一个nset就行了在矩形面上画个circle然后就在tool中partionQ:怎么样使用Abaqus最好呢?1)用CAE2)直接填写Keyword文件3)用Pythonsc
25、ript4)用其它前处理及后处理(那么用什么前处理最好呢?Hypermesh还是?怎么样使用呢)A:根据具体情况而定,不复杂几何模型模型最好在CAE中建模这是最容易对模型控制的方式,无论是几何还是网格写keyword会累死人的,除非很多相同的部分,简单的装配,当年师兄用C+编了一个程序写斜拉桥input文件,input文件有N兆,佩服佩服个人认为hypermesh是最好的前处理软件,做复杂的大模型再好不过pythonedit国外用的很多,当然用它来继续开发CAE的功能还是很方便!有限元要充分利用各种工具的特点和长处。比如造型,比如剖发网格,比如前置处理。我的习惯思路:ug造型,导出模型;fem
26、ap或者hypermesh剖发网格;femap作前置,导出inp文件;手工检查inp文件,部分性的修改;最后用abaqus计算。简单的模型可以在cae里作。最近看着别人用patran,我也临时学点patran作前置,挺好。Q:想用patran调用abaqus,如何在ini文件中设置那?A:这个问题我也碰到过。安装patran的时候,没有装abaqus的介质。只要重新安装patran选择用户自定义安装。注意abaqus的介质选项即可。Q:在Abaqus/CAE中如何施加三角形线荷载?例如,一根悬臂梁上作用有三角形线荷载。只找到了均布荷载Uniform,还有就是UserSubroutine了,这么
27、简单的问题应该不用如此复杂吧?Inp文件中似乎也只有Dload。在Ansys中就很好处理此问题A:youcanuseeitherAMPLITUDEorTABULARtodefineloadvalues.Ifitisatriangle,justgiveafewdatafromthetriangleloadline.Thenwhenyoudefineload,besuretospecifythenameoftheAMPLITUDEorTABULAR.Q:Abaqus的背景只能黑色吗?怎么修改?A:是个问题呀!我是把图片硬拷贝下来后用其他图形软件改变背景的。还有啊,图片怎么导出呢?导出是什么格式的呀
28、?利用Abaqus的print功能可以将图导出到图形文件中,具体路径如下:FilePrint,在弹出的菜单中,不选printviewportbackground,rendition可以是white&black,也可以是color.将Destination由printer改为file,自己定义filename,文件格式由PS改为Tiff图形格式。此文件可以用画笔打开,也可以插入到word中。可惜,图片文件质量太差,还不如抓图得到的质量呢。抓图有个毛病,就是黑色北京需要反色,而反色后,热色变成凉色,看上去不大好。改变Contour的色谱,用相反排列色谱,反色后,应力最高的反而成了黄色,看着也不爽。
29、大家有好办法的话,我就不用画笔程序再一点点的改了,虽然工作量不大,就是觉得不舒服!我试了一下,觉得还是直接保存为.png格式的图片比较好,而且文件也比较小,可以插入word。的确不错!不过需要注意,保存前最好把那些legend,titleblock和stateblock的字体颜色改为黑色,背景改为其他你喜欢的颜色,这样效果会更好,不至于看不清图上的文字。Q:ABAQUS建模问题:如何建立复合材料的模型呢?我看了CAE的文档,但是还是不知道如何建立下列模型。我要建立一个圆柱体,高度H,半径RX1,其中x2R0为材料1,X1Rx2为材料2,在CAE中如何实现?不考虑接触,因为两种材料性能相差不大。
30、请高手指点。A:我想你应该首先建立两种材料的定义。然后,使用assign命令对于不同的模型部分,分配不同的材料。那要如何操作?首先建立一个圆柱PART,然后建立2种材料,那如何将2种材料分配给一个PARTINSTANCE的不同部分?或者如何将2个不同直径的PART组装为一个整体?请指教。wouldliketoansweryourquestion,althoughIhaveneverplayedaroundAbaqus/cae.Sincethegeometryofthemodelissimple(acylinder),mysolutionisyoudonothavetouse/cae,andyo
31、umanuallywriteyourinputfileorrewritetheinputfilegeneratedbyabaqus/cae.IguessIamkindofcheatinginansweringyourquestion.However,Ithinksomehintscouldbehelpful.1)Again,takeadvantageofthesymmetryifpossible.Ifthisisanaxisymmetricalproblem,youneedtomodelonlyaplane.Maybeaquartermodel?2)Whenyoudefinetheelemen
32、ts,useelsettogrouptwodifferentelements:suchas:elset,elset=material1,generate1,100200,300elset,elset=material2,generate101,199whichmeanstheelementsfrom#1to#100,andfrom#200to#300areonematerial,whileelementsfrom#101to#199aretheother.Then,youcoulddefinematerialpropertiesbymaterial,name=material1material
33、,name=material2Q:abaqus的图形如何copy?A:fileprintfile格式为png,可以用Acdsee打开。Q:用Abaqus能否计算Dep不对称的问题?A:可以,并且在step里面的editstep对话框other里面的matrixsolver有个选项。Q:弹塑性矩阵【D】与ddsdde有何联系?A:stress=D*stran;d(stress)=ddsdde*d(stran)。Q:在abaqus中,如果采用umat,利用自己的本构,如何让abaqus明白这种材料的弹塑性应变,也就是说,如何让程序返回弹性应变与塑性应变,好在output中输出,我曾想用最笨地方法,
34、在uvarm中定义输出,利用getvrm获取材料点的值,但无法获取增量应力,材料常数等,研究了帮助中的例子,umatmst3.inp,umatmst3.for,他采用misesJ2流动理论,我在outputhistory显示他已进入塑性状态,但他的PE仍然为0!?A:用uvar()勉强成功。Q:本人在用umat作本构模型时,*static,1,500,0.000001,0.1此时要求的增量步很多,即每次增量要很小,*static1,500时,在弹性向塑性过度时,出现错误,增量过大,出现尖点?A:这样尝试:*STEP,EXTRAPOLATION=NO,INC=2000000*STATIC0.00
35、1,500.0,0.00001,0.1。Q:模型中存在两个物体的接触,计算过程中报错,怎么回事?A:接触问题不收敛有两个方面不妨试试:一、在*CONTACTPAIR里调试ADJUST参数;二、调一些模型参数,比如FRICTION等。Q:在边界条件和加载时,总是有initial这个步,然后是我们自己定义的加载步,请问这个initial步,主要作用是什么?能不能去掉?A:不能去掉,所有的分析都有,是默认的步。Q:Asolidextrusionbasefeature这句话是什么意思?extrusion、revolution等是什么意思?A:这两的是三维建模时候,在画完二维图形,如何来生成三维图形,e
36、xtrusion意思是你给定一个厚度,然后二维图形第三个方向上面伸展这么多形成三维图形revolution意思是你给定一个旋转轴,二维图形绕其旋转后形成三维轴对称图形。Q:偶在umat中调用求主应力函数CALLSPRINC(STRESS,PS,LSTR,NDI,NSHR)后,存储主应力得数组PS中各个主应力排列顺序是什么?PS1PS2PS3?PS1PS1PS3PS2?A:第二个。Q:在*USERMATERIAL的定义中,对的偏倒数,即DDSDDE被称为:materialstiffnessmatrix;而在UMAT中,DDSDDE被称为:materialsJacobianmatrix。请问DDS
37、DDE和材料的切线刚度矩阵的关系是什么?二者是一个概念么?A:一般说可以这样理解:ddsdde*有点像我们常说的弹塑性矩阵:Dep。Q:请问fieldoutput和historyoutput什么区别?关键字*nodeoutput和*nodeprint有什么区别??A:fieldoutput和historyoutput在viewer模块中很明显。fieldoutput是场量输出,historyoutput为历史数据输出,会记录场量随时间的曲线。至于*nodeoutput与*nodeprint的区别在于他们写入不同的后缀文件,比如*.odb,*.dat等,在help中有详细介绍*nodeoutpu
38、t是给CAR或VIEWER做后处理用的,是二进制文件,而*nodeprint则是写到DAT文件中,你可以自己看的文本文件。fieldoutput用来输出模型中较大部分(alargeportionofthemodel)的那些输出频率较小的变量,如模型的等值线,变形图等;而historyoupput用来输出模型中较小部分的那些输出频率较大的变量,如荷载作用点处的荷载-位移图等。所以,应力,应变,位移,反力等既可以作为fieldoutput也可以作为historyoupput输出。Q:我在学习ABAQUS时,遇到以下专有名词想请教各位。1、orphanmeshinstance,2、self-cont
39、act,3、elasticfoundation,4、convectiveinteraction,5、amplitude,6、solution-dependentstatevariable,7、datum,8、thermalfilmcondition?A:1、独立网格实体;2、自接触;3、弹性基础,如弹性地基梁;4、对流相互作用;5、幅度,数值大小,如定义随时间变化的荷载;6、依赖于解的状态变量,在UMAT中经常用到,它的大小取决于某一增量步收敛解;7、辅助数据,在CAE中经常用于定义数据点、刚体参考点、辅助平面等;8、Q:(1),请问如何得到M,C,K矩阵?用什么命令?(2)ABAQUS中能实
40、现这样的东西吗?就是我需要平滑ABAQUS产生的位移场,还是这个平滑只能在其他环境中完成?A:(1)试一试:*elementmatrixoutput(2)试试这个输出选项:*Elprint,Position=averageatNodes。Q:solution-dependentvariable和time-dependentvariable这两种变量有什么差别?A:solution-dependentvariable从字面上理解,即为与解答有关的变量,它的值与每一增量结束时的变量有关,常用于UMAT中,即SDV,可以由*DEPDV来定义其个数。time-dependentvariable从字面上
41、理解,即为与时间有关的变量,常用于定义随时间变化的量如地震荷载等,可以由*APPLITUDE来定义。Q:S4R单元可以输出sth(就是sectionthickness单元厚度变化),但C3D8R却不能输出sth。该怎么看单元层厚度方向的变化?A:个人觉得定义了局部坐标方向,特别是对于各向异性问题,有限元在材料方向上计算应力,应变,输出也就在材料方向,假如第三方向为厚度方向,那第三方向的塑性应变就是他的厚度变化。Q:ABAQUS在运行过程中,对于一些比较大的问题,经常会出现以下的信息:ABAQUSInfo:LicenseTimeoutsetto70(可能是不同的值)minutes.?A:个人感觉
42、是没有问题的,你放心吧。Q:在计算固结过程中,msg文件总是有一下信息*ERROR:TOOMANYATTEMPTSMADEFORTTHISINCREMENT:ANALYSISTERMINATED改变其最小时间变量也不行,太郁闷了,是不能收敛吗?请指教可能是什么原因?A:很多原因,比如边界条件不对;约束不够;接触定义不对;单元划分畸变等等。Q:A:个人觉得:umat实现自己的本构没有固定的方法,对于不同的本构有可能必须采用不同的方法。这要靠自己不断地摸索。有可能一种方法对于简单加载问题还行,但有可能对于复杂问题并不收敛。最重要一点,就是umat中采用的算法必须consistent.再就是ddsd
43、de必须正确,(如果采用back_Euler方法等一些算法,ddsdde错误有时不影响结果(对于简单加载问题没有影响,能收敛,),但对于复杂问题不收敛。uptonow,你这个算法对于Mises,hill,J2,J2d等一类的屈服函数是正确的,但具体的本构还要灵活运用,这我也正学习,正在摸索。有时,umat需要很强的有限元基础,并且对采用的本构要很熟悉,不要在一颗树上吊死才好。首先要确认自己的umat没有错误,如果没有,但就是不收敛(在不断减小加载步长的情况下,当然最好对步长不敏感,特别是对于粘弹性,粘塑性,内变量一类的材料,有的本构取决于背应力的计算)。那就应该考虑换一种算法。下面续贴本人摘自
44、剑桥大学BBS上的FAQ:1.GeneralQuestionsQ1.1:HowmuchdiskspacedoIneed?Minimumof10MBytestorunsimpleexamplesifusingPATRANforPreorPostprocessing.Probablyabout20MBytesformediumsizedproblems.Theamountofdiskspacerequireddependsonthenumbersofnodes/elementspresentinthemeshandtheoutputfrequencyrequestedforoutputtingda
45、tatothe*.fil,*.datand*.resfiles.ThefollowingABAQUSfilesarecreatedintheusershomedirectory.*.dat,*.fil,*.res,*.msg,*.sta,*.log,*.job,*.023(deletedatthecompletionoftheABAQUSjob)._Q1.2:WhattypeofanalysiscanIdo?TheproceduresavailableinABAQUSarelistedbelow(inalphabeticalorder):BUCKLECOUPLEDTEMPERATURE-DIS
46、PLACEMENT(steadystateandtransient)COUPLEDTHERMAL-ELECTRICAL(steadystateandtransient)DYNAMICFREQUENCYGEOSTATICHEATTRANSFER(steadystateandtransient)MASSDIFFUSION(steadystateandtransient)MODALDYNAMICRANDOMRESPONSERESPONSESPECTRUMSOILS,CONSOLIDATIONSOILS,STEADYSTATESTATICSTEADYSTATEDYNAMICSVISCO_Q1.3:Wh
47、atoutputoptions(hardcopy)arethere?WithABAQUS5.7anysizeplots(includingA3andA4)canbeobtained.IfusingPATRANofFEMGVhardcopiesofanysize(includngA4andA3)canbeobtained._Q1.4:Whatpreandpost-processingprogramsareavailableotherthanABAQUS/Post?PATRANandFEMGVareavailableintheteachingsystem.Boththeseprogramsread
48、theresultsfile*.filforpostprocessing.IfyouarepostprocessingusingFEMGVandhavewrittenresultsatelementGausspointsthenyouneedtowritethegausspointco-ordinatesaswell.IncludethefollowinglinesintheABAQUSinputfile:*ELFILE,POSITION=INTEGRATIONPOINTSSECOORDOtherwiseitwillnotbepossibletousetheFEMGVprogramtopost
49、processtheelementgausspointresultsieFEMGVwouldnotknowthelocationofthegausspoints._Q1.5:IsABAQUS/Preavailable?No.ABAQUS/Preisnotavailableintheteachingsystem.HoweverPATRANwhichprovidesthesamefunctionalityandwhichhasthesamelookandfeelofABAQUS/Preisavailableintheteachingsystem.TypepatrantorunthePATRANpr
50、ogram.Q1.6:IsABAQUS/Explicitavailable?ABAQUS/Explicitisonlylicencedontw900.Youneedtorlogintotw900torunthisversionofABAQUS.ABAQUS/Standardisalsoavailableontw900._Q1.7:Whaton-linedocumentationareavailableandhowdoIaccessit?Typeabaqus57docintheCUEDteachingsystemtoaccessthefullsetoftheABAQUSUsersmanual(V
51、olumesI,IIandIII).Atpresentthisistheonlyon-linedocumentationavailableexceptforABAQUSReleaseNotes(100pages)andABAQUSSiteGuide(80pages)._Q1.8:WhatunitsareusedinABAQUS?ThereisnoinherentsetofunitsusedinABAQUS.Itisuptotheusertodecideonaconsistentsetofunitsandusethatunits.Typicalsetsofunits:12Length-metre
52、smmForce-NewtonsNewtonsTime-secondsecondMass-Kgtonne(*)Density-Kg/m3tonne/mm3Stress-N/m2N/mm2(=MPa)YoungsModulus-N/m2N/mm2(=MPa)*1tonne=1000kilogramsDecideontheunitsbeforeyoustartpreparingyourdata.Thisiscritical.Ifyoustarttypinginthenodalco-ordinatesthatmeansyouhavealreadydecidedonwhatunitstouserfor
53、theLengthparameter.ThisonlyleavesthechoicefortheunitsofForce.Itisnotagoodideatochoosemmforlength,NewtonforforceandthenspecifytheYoungsmodulusinKN/m2.2.JobsQ2.1:HowdoIrunsmalljobs?Usethefollowingcommandabaqusjob=job-idinteractiveExample:abaqusjob=plate_Q2.2:HowdoIgoaboutrunningmanysmalljobs?Createaba
54、tchfile(sayaba.run)withonelineperanalysisasshownbelow:abaqusjob=analysis-ainteractiveabaqusjob=analysis-binteractiveabaqusjob=analysis-cinteractiveabaqusjob=analysis-dinteractiveThenmakethefileanexecutableusingthefollowingunixcommand:chmodu+xaba.runThentypeaba.runtoexecutetheABAQUSjobsoneatatimewhil
55、eyouareloggedON.Thisisonlysuitableforsmalljobswhichonlytakeafewminutestorun.Thesejobswillrunoneatatimeandinsequence.Thisispreferabletosubmittingallthejobsatthesametime(forexampletypingtheabovecommandsdirectlyattheterminalwithouttheinteractiveparameter).Thiswillputastrainontheserveranditsresourcesand
56、inconveniencetheotherusersaswell.FormediumtolargejobsusethebatchcommandavailableintheCUEDteachingsystem._Q2.3:HowdoIrunlargejobsusingbatch?IntheCUEDteachingsystemusethebatchcommand.Torunthejobontw500ortw900serversuse:batch-QX-mbaoabaqusjob=job-idinteractiveTorunajobspecificallyonthetw900serveruse:ba
57、tch-QN-mbaoabaqusjob=job-idinteractiveExample:batch-QN-mbaoabaqusjob=plateinteractiveTheprogressofthesubmittedjobcanbemonitoredusingthebatchqcommand.Usethebatchrmcommandtodeleteanybatchjobsyouhadsubmittedbeforethesearerun,ifyouchangeyourmind.Seethemanpagesonbatchq,batchrmformoredetails.Example:Typem
58、anbatchq._Q2.4:HowdoIrunABAQUS/Post?Usingthefollowingcommand:abaquspostjob=job-idExample:abaquspostjob=cantilever_Q2.5:HowdoIrunABAQUS/Plot?Usingthefollowingcommand:abaqusplotjob=job-iddevice=cpsorhglorx11Usedevice=x11toviewtheplotsonthescreen.Useoneoftheotheroptions(cpsforcolourpostscriptorhglforhp
59、gl)tocreateahardcopyoftheplot.Example:abaquspostdevice=cpsjob=cantileverQ2.6:HowdoIgetcopiesoftheABAQUSexamplesinputdatafiles(*.inp)?TheABAQUSdatafilesusedintheexamplesmanualcanbefoundinthe/export/abaqus/samples/exastddirectory.Similarlythedatafilesintheverificationmanualcanbefoundin/export/abaqus/s
60、amples/verstddirectory.Thesefileswillhavetheextensionnameinp.Allfilenamesareinlowercase.usethecpcommandtocopytherelevantfile.Example:cp/export/abaqus/samples/exastd/1010101.inp.Thiswillcopythe1010101.inpfiletothecurrentdirectory.Ifthesedirectoriesdontexistthencanusetheabaqusfetchcommand.Example:abaq
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论