创建博客 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

神瑛侍者的博客

 
 
 

日志

 
 

abaqus FAQ  

2009-11-18 15:51:00|  分类: 默认分类 |  标签: |举报 |字号 订阅

Abaqus 使用FQA:

      Q: abaqus的图形如何copy?
      A: file>print>file格式为png,可以用Acdsee打开。

      Q: 用Abaqus能否计算[Dep]不对称的问题?
      A: 可以,并且在step里面的edit step对话框other里面的matrix solver有个选项。

      Q: 弹塑性矩阵【D】与ddsdde有何联系?
      A: stress=D*stran;d(stress)=ddsdde*d(stran)。

      Q: 在abaqus中,如果采用umat,利用自己的本构,如何让abaqus明白这种材料的弹塑性应
      变,也就是说,如何让程序返回弹性应变与塑性应变,好在output中输出,我曾想用最笨
      地方法,在uvarm中定义输出,利用getvrm获取材料点的值,但无法获取增量应力,材料常
      数等,研究了帮助中的例子,umatmst3.inp,umatmst3.for,他采用mises J2 流动理论,我
      在output history 显示他已进入塑性状态,但他的PE仍然为0!!?
      A: 用uvar( )勉强成功 。

      Q: 本人在用umat作本构模型时,
      *static,
      1,500,0.000001,0.1 此时要求的增量步很多,即每次增量要很小,
      *static
      1,500 时,在弹性向塑性过度时,出现错误,增量过大,出现尖点.?
      A: YOU CAN TRY AS FOLLOWS:
      *STEP,EXTRAPOLATION=NO,INC=2000000
      *STATIC
      0.001,500.0,0.00001,0.1。

      Q: 模型中存在两个物体的接触,计算过程中报错,怎么回事?
      A: 接触问题不收敛有两个方面不妨试试:
      一、在*CONTACT PAIR 里调试ADJUST参数;
      二、调一些模型参数,比如FRICTION等。 。

      Q: 在边界条件和加载时,总是有initial这个步,然后是我们自己定义的加载步,请问这
      个initial步,主要作用是什么?能不能去掉?
      A: 不能去掉,所有的分析都有,是默认的步。

      Q: A solid extrusion base feature 这句话是什么意思?
         extrusion、revolution等是什么意思?
      A: 这两的是三维建模时候,在画完二维图形,如何来生成三维图形,
      extrusion意思是你给定一个厚度,然后二维图形第三个方向上面伸展这么多形成三维图形
      revolution意思是你给定一个旋转轴,二维图形绕其旋转后形成三维轴对称图形 。

      Q: 偶在umat中调用求主应力函数
      CALL SPRINC(STRESS,PS,LSTR,NDI,NSHR)
      后,存储主应力得数组PS中
      各个主应力排列顺序是什么?
      PS1>PS2>PS3 ?
      PS1<PS2<PS3 ?
      PS1>PS3>PS2 ?
      A: 第二个 。

      Q: 在*USER MATERIAL的定义中,Δσ对Δε的偏倒数,即
      DDSDDE被称为:
      "material stiffness matrix";
      而在UMAT中,DDSDDE被称为:
      "material's Jacobian matrix"。
      请问DDSDDE和材料的切线刚度矩阵的关系是什么?二者是一个概念么?
      A: 一般说可以这样理解:
        σ=ddsdde*Δε
      有点像我们常说的弹塑性矩阵:Dep。

      Q: 请问field output和history output什么区别?
      关键字*node output和*node print有什么区别? ?
      A: field output和history output在viewer模块中很明显。
      field output是场量输出,history output为历史数据输出,会记录场量随时间的曲线
         
      至于*node output与*node print的区别在于他们写入不同的后缀文件,比如*.odb,*.dat
      等,在help中有详细介绍 *node output是给CAR或VIEWER做后处理用的,是二进制文件,
      而*node print则是写到DAT文件中,你可以自己看的文本文件。

      field output用来输出模型中较大部分(a large portion of the model)的那些输出频
      率较小的变量,如模型的等值线,变形图等;
      而history oupput用来输出模型中较小部分的那些输出频率较大的变量,如荷载作用点 处
      的荷载--位移图等。
      所以,应力,应变,位移,反力等既可以作为field output也可以作为history oupput 输
      出。 。

      Q: 我在学习ABAQUS时,遇到以下专有名词想请教各位。
      1、orphan mesh instance,
      2、self-contact,
      3、elastic foundation,
      4、convective interaction,
      5、amplitude,
      6、solution-dependent state variable,
      7、datum,
      8、thermal film condition ?
      A:
      1、独立网格实体;
      2、自接触;
      3、弹性基础,如弹性地基梁;
      4、对流相互作用;
      5、幅度,数值大小,如定义随时间变化的荷载;
      6、依赖于解的状态变量,在UMAT中经常用到,它的大小取决于某一增量步收敛解;
      7、辅助数据,在CAE中经常用于定义数据点、刚体参考点、辅助平面等;
      8、不知道:) 。

      Q: (1),请问如何得到M,C,K矩阵?用什么命令
      (2)ABAQUS中能实现这样的东西吗?就是我需要平滑ABAQUS产生的位移场 ,还是这个平滑
      只能在其他环境中完成 ?
      A: (1)试一试:*element matrix output
      (2)试试这个输出选项:
      *El print,Position=average at Nodes。

      Q: solution-dependent variable 和 time-dependent variable 这两种变量有什么差别
      ?
      A: solution-dependent variable
      从字面上理解,即为"与解答有关的变量",它的值与每一增量结束时的变量有关,常用于U
      MAT中,即SDV,可以由*DEPDV来定义其个数。
      time-dependent variable
      从字面上理解,即为"与时间有关的变量",常用于定义随时间变化的量如地震荷载等,可
      以由*APPLITUDE来定义。 。

      Q: S4R单元可以输出sth(就是section thickness单元厚度变化),但C3D8R却不能输出sth
      该怎么看单元层厚度方向的变化 ?
      A: 个人觉得
      定义了局部坐标方向,特别是对于各向异性问题,有限元在材料方向上计算应力,应
      变,输出也就在材料方向,假如第三方向为厚度方向,那第三方向的塑性应变就是他的厚
      度变化。 。

      Q: ABAQUS 在运行过程中,对于一些比较大的问题,经常会出现以下的信息:
      ABAQUS Info: License Timeout set to 70(可能是不同的值) minutes. ?
      A: 个人感觉是没有问题的,你放心吧 。

      Q: 在计算固结过程中,msg文件总是有一下信息 ***ERROR: TOO MANY ATTEMPTS MADE
      FOR TTHIS INCREMENT: ANALYSIS TERMINATED
      改变其最小时间变量也不行,太郁闷了,是不能收敛吗?各位大侠救命。请指教可能是什
      么原因?
      A: 很多原因
      比如边界条件不对 ;约束不够 ;接触定义不对;单元划分畸变等等 。

      Q: ?
      A: 个人觉得:
      umat实现自己的本构没有固定的方法,对于不同的本构有可能必须采用不同的方法。这要
      靠自己不断地摸索。有可能一种方法对于简单加载问题还行,但有可能对于复杂问题并不
      收敛。最重要一点,就是umat中采用的算法必须consistent.再就是ddsdde必须正确,(如
      果采用back_Euler 方法等一些算法,ddsdde错误有时不影响结果(对于简单加载问题没有
      影响,能收敛,),但对于复杂问题不收敛。
         
      uptonow,你这个算法对于Mises,hill,J2,J2d等一类的屈服函数是正确的,但具体的本构
      还要灵活运用,这我也正学习,正在摸索。

      有时,umat需要很强的有限元基础,并且对采用的本构要很熟悉,不要在一颗树上吊死才
      好。首先要确认自己的umat没有错误,如果没有,但就是不收敛(在不断减小加载步长的
      情况下,当然最好对步长不敏感,特别是对于粘弹性,粘塑性,内变量一类的材料,有的
      本构取决于背应力的计算)。 那就应该考虑换一种算法。

      一点体会,请大家探讨。

      Q: 第一次安装abaqus6.2.1,装完exceed6.2&3D后,开始装6.2.1的Product installation
      for network licensing,结果每到进度为75%时,弹出界面:
      An error ocurred during the move data process: -115
      Component:Complete Inslattation
      File Group: fg_common
      File:c:\abaqus\6.2-1\cae\External\ebt\adi3xcol.dll
      就停止安装了,重新安装依然如此,现在感到头都大了,哪位好心人帮帮我,不深感激 ?
      A: 没有关闭防火墙!

      Q: WARNING: THE SYSTEM MATRIX HAS 1 NEGATIVE EIGENVALUES.
      一般在什么情况下会发生系统矩阵出现负特征值?
      A: 很多情况
      比如边界条件定义不对;接触定义不对 。

      Q: 就是手册的umatmst3.for
      1)Mises各向同性的子程序中,调用硬化曲线函数中的Table(2,nvalue)这个二维数值
      保存是应力应变曲线的数据,可是ABAQUS并没有告诉UMAT中Table数值这些应力应变的数值
      啊,UMAT里面Table怎么知道的呢?
      2)RHS=SMISESS-EG3*DEQPL-SYIELD这个式子是什么意思?尤其是 3G*DEQPL是什么意
      思?
      A: 就是Constant=8中的从第三个开始的数值,数组名传递。 。

      Q: 我在计算时MSG文件出现如下错误,是不是由于节点数太多啦?该如何处理这个问题?

      ERROR: SPECIFIED STANDARD_MEMORY VALUE OF 8000000 IS TOO SMALL TO RUN THE
      NALYSIS. STANDARD_MEMORY MUST BE INCREASED. MINIMUM POSSIBLE VALUE IS
      23477555. LOOK AT MEMORY ESTIMATES SECTION OF .dat FILE FOR FURTHER
      INFORMATION ?
      A: 对abaqus_v6.env文件中的STANDARD_MEMORY的值进行修改。
      修改Site文件夹下的abaqus_v6.env中的配置,如下:
      #
      # System-Wide ABAQUS Environment File
      # -------------------------------------
      pre_memory = 33554432
      standard_memory = 33554432 ##*********修改这里******************
      #
      # NT specific settings 。

      Q: WARNING: THE SYSTEM MATRIX HAS 148 NEGATIVE EIGENVALUES
      1 ABAQUS VERSION 6.3-1 DATE 27-NOV-2002 TIME 22:08:00 PAGE..
      For use by None user license from HKS Inc. ..
      STEP 1 INCREMENT 1 STEP TIME 0.00
      STEP 1 S T E A D Y S T A T E S ..
      AUTOMATIC TIME CONTROL WITH -
      A SUGGESTED INITIAL TIME INCREMENT OF 0.300?
      A: 将你的INITIAL TIME INCREMENT 改小些试试看,不过,一般出现此类问题,多半是你
      的模型有问题,欠约束或者其它什么的。

      Q: ZERO PIVOT是什么意思 ?
      A: zero pivot可以理解为刚度矩阵出了问题,例如奇异。
          可能有不同的原因,如: 你所模拟的是软化性质材料,该点因破坏等原因而软化至不
      能吃劲; 模型有问题,如约束不够,或者是单元拓扑出错等;或者是你的UMAT中写的[dds
      dde]有错;。

      Q: abaqus-uamt的老问题,缺少 'ABA_PARAM.INC'文件 ?
      A: 在cvf6.5调试时,显示 缺少 'ABA_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.IN C.dp或
      ABA_PARAM.INC.sp改名为ABA_PARAM.INC。
      呵呵,他的意思是在Visual Fortran中调试其子程序,我觉得这是一个好办法,我当时也
      是这麽办的,毕竟在ABAQUS中调试是非常麻烦的,只有当你的UMAT没有语法或者明显的逻
      辑错误,你在ABAQUS中调试才能事半功倍。

      Q: abaqus中有没有布尔运算?能否在abaqus中将两个体粘成一个体?
      A: 可以,用constraint和interaction可以将两个part进行定位的。

      Q: abaqus里做静力分析时怎么加重力加速度?在inp文件里的card是什么?
      A: 如:
      *DLOAD
      ELSET, GRAV, 9.81, 0., 0., 1.
      有疑问查一些*DLOAD卡就清楚了。

      Q: 如何在ABAQUS CAE如何设定在某一节点上受集中力呢
      请问一下在ABAQUS CAE如何设定在某一节点上受集中力呢??
      我用CAE中的Set好象不能设定一个节点丫….
      请问在CAE要如何设定,或则祇能在inp下设定呢??
      A: 用set可以定义一个节点的。
      另外如果只是在一个节点上加集中力就没有必要用set了,在Load模块里creat a load,选
      择集中力,然后点中该节点(别忘了done)就可以输入力的数值和方向了
      设定一个节点当然可以了,在cae里一点的集中荷载还是挺容易加了。partition,
      datumn point, 都可以
      用CAE建的模型好象不能直接加在接点上,不过可以加在vertex上

      Q: 如何将两个part接在一起?
         几个人将一个工作分开来做,完成后要将造好型的part接在一起,再加边界条件和载荷
      。该如何把这几个part接起来?
      A: 先导入File->import->part,然后形成assembly,这样做试试看如何

      Q: abaqu cae中作梁的分析时可否定义不规则的profile形状?
      只提供了box,pipe,circle,角钢,方钢,T型钢,l型钢,工字钢,我可不可以定义像三角
      形空心钢管,或者梯形空心钢管这样的profile?
      A: 直接给出 I 和 G

      Q: 在CAE中建立PART时为什么只能先建立在X-Y平面上的模型
      在CAE中建立PART时为什么只能先建立在X-Y平面上的模型,可不可以直接建立空间的点,
      结点线,面,体呀
      A: 在图纸上画投影图当然是2维的了。
      这是它的建模思路,也是一种专利了。直接在空间里建立模型当然也是一种思路,哪个好
      呢?各有千秋。
      你要使用CAE,就要熟悉并掌握它的思路,把整体结果分解为各个小part,把每个part都简
      化到何以在平面上画出投影图并extrude或回转就得到part的空间实体。然后把这些part给
      instance,指定其section特性(就是其材料呀,厚度呀等在图形上显示不了的数据)装配
      起来构成复杂的模型。
      同一个part的草图可以有许多instance,你可以尽量利用装配功能来减少你建立part的工
      作量。
      对机械图比较熟悉的,可以想象为制作零件然后装配成机器的过程,对于土木类工程师可
      能开始会不太习惯,但是多做几个例子就会熟悉并喜欢上这种操作了。
      如果你实在是喜欢空间直接建模法,你可以用ANSYS建立模型划分单元后再用小程序变成AB
      AQUS的输入文件

      Q: abaqus的所有功能在cae中都能完成吗?还是需要在命令行中输入
      A: 大部分工作可以在CAE中实现,有些功能还需要在INPUT文件中修改,将CAE转化成*.INP
      ,就可。
      abaqus的前置确实不太好。 cae的功能一直没有长进。 如板料成形方面,回弹的分析就是
      不能在cae里完成,必须写inp文件来做。 写inp文件的效率可一直不高啊 因为一般成形我
      们最好用动态显示算法,而回弹最好静态隐式算法。
      不过在cae里就是不行。只要建立了动态算法的step,后面就不能建立静态算法的step。
      还有啊,abaqus的cae里面是先设置边界条件和输出,然后剖分网格,这个不适合大多数有
      限元前置处理的思路。很少前置软件是这样的。 比起其他的前置软件,如femap,patran
      ,还有工业用软件autoform,pam-stamp,差不少啊!

      Q: 如何给一个面的一部分加载
      A: Are you using /CAE? If so, please partition your surface and apply
      loading.
      对的,将矩形和圆区分开来,将圆定义为一个块就行了!
      将要加载的面上所有的节点设为一个nset就行了
      在矩形面上画个circle 然后就 在tool-》中partion

      Q: 怎么样使用Abaqus最好呢?
      1)用CAE
      2)直接填写Keyword文件
      3)用Python script
      4)用其它前处理及后处理(那么用什么前处理最好呢?Hypermesh还是?怎么样使用呢)
      A: 根据具体情况而定,不复杂几何模型模型最好在CAE中建模 这是最容易对模型控制的方
      式,无论是几何还是网格 写keyword会累死人的,除非很多相同的部分,简单的装配, 当
      年师兄用C++编了一个程序写斜拉桥input文件,input文件有N兆,佩服佩服 个人认为hype
      rmesh是最好的前处理软件,做复杂的大模型再好不过 python edit国外用的很多,当然用
      它来继续开发CAE的功能还是很方便!
      有限元要充分利用各种工具的特点和长处。
      比如造型,比如剖发网格,比如前置处理。
      我的习惯思路:
      ug 造型,导出模型 ;femap或者hypermesh 剖发网格 ;femap作前置,导出inp文件 ;手工
      检查inp文件,部分性的修改 ;最后用abaqus计算。 简单的模型可以在cae里作。 最近看
      着别人用patran,我也临时学点patran作前置,挺好

      Q: 想用patran调用abaqus,如何在ini文件中设置那
      A: 这个问题我也碰到过。 安装patran的时候,没有装abaqus的介质。 只要重新安装patr
      an选择用户自定义安装。注意abaqus的介质选项即可。

      Q: 在Abaqus/CAE中如何施加三角形线荷载?
      例如,一根悬臂梁上作用有三角形线荷载。
      只找到了均布荷载Uniform,还有就是User Subroutine了,
      这么简单的问题应该不用如此复杂吧?
      Inp文件中似乎也只有Dload。
      在Ansys中就很好处理此问题~~~~~~~~
      A: you can use either AMPLITUDE or TABULAR to define load values. If it is
      a
      triangle, just give a few data from the triangle load line. Then when you
      define load, be sure to specify the name of the AMPLITUDE or TABULAR.

      Q: Abaqus的背景只能黑色吗?怎么修改?
      A: 是个问题呀!我是把图片硬拷贝下来后用其他图形软件改变背景的。还有啊,图片怎么
      导出呢?导出是什么格式的呀?
      利用Abaqus的print功能可以将图导出到图形文件中,具体路径如下:
      File>Print,在弹出的菜单中,不选print viewport background,
      rendition可以是white&black,也可以是color.
      将Destination由printer改为file,自己定义file name,文件格式由PS改为
      Tiff图形格式。此文件可以用画笔打开,也可以插入到word中。可惜,图片文件质量太差
      ,还不如抓图得到的质量呢。
      抓图有个毛病,就是黑色北京需要反色,而反色后,"热"色变成"凉"色,看上去不大好。
      改变Contour的色谱,用相反排列色谱,反色后,应力最高的反而成了黄色,看着也不爽。
      大家有好办法的话,我就不用画笔程序再一点点的改了,虽然工作量不大,就是觉得不舒
      服!
      我试了一下,觉得还是直接保存为.png格式的图片比较好,而且文件也比较小,可以插入w
      ord。的确不错!不过需要注意,保存前最好把那些legend,title block和state block的
      字体颜色改为黑色,背景改为其他你喜欢的颜色,这样效果会更好,不至于看不清图上的
      文字。

      Q: ABAQUS建模问题:如何建立复合材料的模型呢?
      我看了CAE的文档,但是还是不知道如何建立下列模型。
      我要建立一个圆柱体,高度H,半径R=X1,其中x2>R>0为材料1,X1>R>x2为材料2,在CAE
      中如何实现?不考虑接触,因为两种材料性能相差不大。请高手指点
      A: 我想你应该首先建立两种材料的定义。
      然后,使用assign命令对于不同的模型部分,分配不同的材料。
      那要如何操作?首先建立一个圆柱PART,然后建立2种材料,那如何将2种材料分配给一个P
      ART INSTANCE的不同部分?或者如何将2个不同直径的PART组装为一个整体?请指教。
         would like to answer your question, although I have never played around

      Abaqus/cae. Since the geometry of the model is simple (a cylinder), my
      solution is you do not have to use /cae, and you manually write your input

      file or rewrite the input file generated by abaqus/cae. I guess I am kind
      of
      cheating in answering your question. However, I think some hints could be
      helpful.
      1) Again, take advantage of the symmetry if possible. If this is an
      axisymmetrical problem, you need to model only a plane. Maybe a quarter
      model?
      2) When you define the elements, use "elset" to group two different
      elements:
      such as:
      elset, elset=material1, generate
      1,100
      200,300
      elset, elset=material2, generate
      101,199
      which means the elements from #1 to #100, and from #200 to #300 are one
      material, while elements from #101 to #199 are the other.
      Then, you could define material properties by
      material, name=material1
      ......
      material, name=material2
      ......

给出了一个用户定义场USDFLD的例子,用来计算层状复合材料平板的损伤和失效.

用户子程序看

damagefailcomplate_cps4.f

输入文件:

damagefailcomplate_cps4.inp

damagefailcomplate_element.inp

damagefailcomplate_node.inp

运行

abaqus analysis job=damagefailcomplate_cps4 user=damagefailcomplate_cps4

  评论这张
 
阅读(589)| 评论(0)
推荐 转载

历史上的今天

最近读者

热度

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2014