abaqus 常用 技巧 总结.doc

1、Abaqus 常用技巧总结本手册是由 simwe 等论坛上精华帖以及本人下载的其他资料整理，由于很多资料搜集已经很久，而且时间有限，都没有注明原作者，也没有最资料进行分类整理，见谅。建议阅读方式：目录ABAQUS 常用技巧总结 .1目录 .21.对 TIME INCREMENT 的根本理解 .42.ABAQUS 请问 MOMENT 的加载 .43.ABAQUS 计算时 C 盘的临时文件太大了，怎么改目录？ 54.CAE 中如何加预应力 .55.HYPERMESH 里面看到 ABAQUS 分析的结果 .56.X-Y PLOTS57.把上一次的分析结果作为下一次分析的初始条件该怎么做 .68.材料

2、方向与增量步 79.多个 INP 文件如何实现批处理 710.关于 ABAQUS 的任务管理 811.关于数据的输入输出 .912.后处理积分 .1013.接触分析激活杀死 1014.利用 QUEUE 的功能由本地机器向远程 UNIX 机器提交 ABAQUS 作业的方法 精华 .1015.利用命令进行计算时如何设置调用内存量 1316.清华大学 BBS 的 ABAQUS 精华 .1317.请问怎么实现双曲线 .3618.取消坐标系等的显示 .3719.如何在计算中修改材料特性 3720.输出计算过程中的总质量和总刚度矩阵 4021.先张预应力： .4022.用户子程序的使用 4123.怎样设定

3、用双 CPU 机器进行 ABAQUS 计算 .4124.中途停止正在运算的 JOB 4225.自适应网格技术 4226.ABAQUS 计算与内存 4227.质量缩放 .4328.ABAQUS 多处理器进行并行计算的效果研究 5429.YAHOO 讨论组摘录 -CONTACT+OVERCLOSURE.5530.原创：无限元建立方法，希望得到加分 7031.分享ABAQUS 使用问答 7632.转帖ABAQUS6.4 导入外来模型的几点小经验！ .8733.ABAQUS 的多图层绘图 .8834.子结构 8935.如何在不同的分析步改变材料的参数 .8936.模型的重启动分析RESTART 913

4、7. ABAQUS 的单位心得 .921.对 time increment 的根本理解abaqus 的 step 里有 maximum number of increment、initial increment、minimum increment 、maximum increment 四个量许多网友不知怎样设置合理，合理设置是建立在深刻理解基础上的。要理解这个问题，首先需要了解 abaqus 的计算过程和有限元计算收敛性问题， abaqus首先用 initial 值输入进行叠代计算，如果计算结果收敛，则继续以这个值代入计算下一步，如果不收敛，则自动减小时间步长（time increment）重

5、新计算直到收敛然后计算下一步。但是如果时间步长减小到最小值 minimum 时计算结果还是不收敛，则 abaqus 将停止计算，由此可知 maximum 值和 minimum 值分别是 abaqus 在收敛计算时时间步长的上下限，同时 total time=求和(time increment*number)，当时间步长很小时，需要计算的步数 number相应增大(电脑计算花的时间也随之增大 )，因此 number 一般要设置较大值。minimum 并不是越小越好，因为 1)number 即计算时间增大 2)abaqus 计算精度约在 10(-5),当时间步长小于这个值，计算结果已经没什么意义了

6、。有限元计算收敛性与（最小空间步长/时间步长）值有关，若 minimum 设为 10(-5),还是不收敛，可适当减小空间步长（即把网格画细点） ，当然还有一些其他办法，如果实在计算不了，也许是模型本身有点问题，或改为显示 explicit 计算总而言之，maximun number 要适当设置较大值，initial 可适当改小（如-2，-3 量级） ，minimum（-5 量级）不要修改， maximum 值影响不大，可不改2.moment 的加载请教各位大牛 我最近在算一个题目,一个大筒体上有三个接管端面固定,大筒体两端加载扭距,但是不知道如何加载.请各位高手能赐教 谢谢1. 将大筒体两端要

7、施加扭矩的节点分别定义为两个 Nset:left, right. 2. 分别在大筒体两端的圆心处定义两个 reference node: rp-left,rp-right. 3. 用如下命令将两个节点集绕 3 轴旋转的自由度与参考点耦合起来,其他自由度度是否耦合根据具体问题而定: *KINEMATIC COUPLING, REF NODE=rp-left left, 6,6 *KINEMATIC COUPLING, REF NODE=rp-right right, 6,6 4. 在两个参考点上施加绕 3 轴旋转的弯矩. 友情提醒:reference node 也有自由度,注意相应的边界条件 .

8、3.abaqus 计算时 c 盘的临时文件太大了，怎么改目录？临时目录是 Windows 自己定义的 ,可以在系统环境变量中修改.4.CAE 中如何加预应力具体没作过，看看*PRESTRESS HOLD 和*INITIAL CONDITIONS, TYPE=SOLUTION, REBAR这两个命令以及 ABAQUS Analysis Users Manual “Defining reinforcement,” Section 2.2.3 “Defining rebar as an element property,” Section 2.2.4 这两节吧。5.hypermesh 里面看 aba

9、qus 分析的结果、你在 abaqus 中计算完成后，将结果文件输出到*.fil.2、利用 hyperworks 提供的 hmabaqus.exe(在安装目录下的 Altairhw7.0translators 中)3、在控制台下运行 hmabaqus *.fil *res，执行完成后就生成了相应的 res 文件4、在 hyperview 中打开你的模型文件*.inp 和结果文件*.res，就可以查看你的结果了gook luck!6.X-Y Plots* STEP: pre-load*Step, name=pre-load, nlgeompre-loading*Static0.01, 1., 1

10、e-05, 0.1* LOADS* Name: pt-load Type: Concentrated force*Cload_G5, 2, -200.E6*Output, history, frequency=1*node output, nset=_G5CF2,U2*element output, elset=_G5E22, S22*monitor, node=_G5, dof=27.如何把上一次分析结果作为下一次分析的初始条件Q 如果想把上一次的分析结果作为下一次分析的初始条件该怎么做 A 使用 LDREAD 命令 首先需要注意下面两个问题 1 每一个 ANSYS 的实体模型的面或体都要定

11、义对应的单元类型编号材料属性编号实常数编号单元坐标系编号这些参数在整个分析过程中保持不变而这些编号对应的属性在各个步骤中是不同的 2 网格划分要满足所有步骤的要求单元类型必须兼容 步骤 1 创建实体模型2 创建多个物理环境 设定一个物理环境中的单元类型材料属性实常数坐标系等 将这些参数的编号赋给实体模型的面或体 施加基本物理载荷和边界条件 设定求解选项 选择一个标题使用 PHYSICS, WRITE 命令将物理环境存入文件中 3 清楚当前的物理环境命令是 PHYSICS, CLEAR 4 重复第二步准备下一个物理环境 8.材料方向与增量步材料方向：针对各向异性材料（如板金材料、复合材料等）变形

12、体，材料方向定义材料的某一特定方向如纤维方向。该方向随着变形体的移动而移动，旋转而旋转。输出的场变量值都以材料方向为参考，有时便于数据处理，分析计算结果等。abaqus 中增量步abaqus 中把所有载荷按一定的要求分成若干载荷步 step，每一步 step 根据 abaqus 自动载荷增量，分成若干增量 increments，每一增量施加一定的载荷，然后每一增量通过若干迭代步 iteration 进行迭代，当系统达到平衡时，迭代结束，完成一个增量。当所有的增量都完成后，计算结束；反之，计算可能出现发散。这时，可以通过采用多钟方法（如调整放大质量系数，单元网格优化等）调整增量大小，使计算继续进

13、行9.多个 inp 文件如何实现批处理我喜欢用 fortran写成这样：! nodedeform.f90! FUNCTIONS:! nodedeform - Entry point of console application.! Example of displaying Hello World at execution time.!*! PROGRAM: nodedeform! PURPOSE: Entry point for Hello World sample console application.!*program nodedeformUSE DFLIB! implicit non

14、eresult = SYSTEMQQ (abaqus job=3d05ctJ interactive)result = SYSTEMQQ (abaqus job=3d05ctJmodel-2 interactive)! result = SYSTEMQQ (abaqus job=test-1 interactive)! result = SYSTEMQQ (abaqus job=test-2 interactive)! result = SYSTEMQQ (abaqus job=test-3 interactive)! result = SYSTEMQQ (abaqus job=test-4

15、interactive)! result = SYSTEMQQ (abaqus job=test-5 interactive)! result = SYSTEMQQ (abaqus job=test-6 interactive)! result = SYSTEMQQ (abaqus job=test-7 interactive)! result = SYSTEMQQ (abaqus job=test-8 interactive)end program nodedeform10.关于 ABAQUS 的任务管理关于 ABAQUS 的任务管理 精华 chibiko 发帖: 26 积分: 1 于 20

16、04-07-03 03:56 - 前些日子的计算量很大，自己的机器又不是很快，一个任务要算上几天。一开始，在算题的时候由于 CPU 占用很大，基本上做不了别的。捉摸了好几天，现在终于有点掌握ABAQUS 的任务管理方法了，主要有两点，这里和大家共享。我的系统是 WindowsXP，使用的是 ABAQUS6.4。 1。ABAQUS6.4 新增了任务管理的命令，可以暂停、恢复、和终止一个正在背景运行的任务，方法如下（在命令行输入并运行）： 任务暂停： abaqus suspend job=job-name 任务恢复： abaqus resume job=job-name 任务终止： abaqus

17、terminate job=job-name 其中任务暂停（suspend）的时候， windows 任务管理栏中仍会保留 standard/explicit 的计算线程，只是不再使用 CPU 资源，当任务恢复（resume）的时候继续工作。 任务终止则就像 CAE 中提交的任务的 KILL 功能类似，直接 cut 掉正在运行的任务，不可恢复。其实也就和在 windows 任务管理栏中强行终止差不多，但属于合法操作。 2。就是利用 windos 的任务管理栏了，使用它的 processes 管理功能时，不仅仅可以强行关掉一个正在运行的任务，还可以通过改变它的优先级别来改变其对 CPU 的占用程

18、度。共有Realtime，High，AboveNormal ，Normal，BelowNormal ， Low 五个可选等级，windows 默认的等级是 Normal，此时所有的任务都在随机抢占 CPU 资源，一般 ABAQUS 在运行时想要运行别的程序就比较困难了，特别是一些大程序。在觉得机器使用时有明显的延迟时，就可以把 ABAQUS 任务的优先级别设的低一些，就可以腾出 CPU 资源给别的级别高的任务了，不用机器的时候再把 ABAQUS 任务的级别调上去，这样就可以娱乐、工作两不误了。 You may try this method which was posted at simwe

19、forum, but not my original idea. Anyway, I will feel very happy even if it probably just helps you just a little bit.11.关于数据的输入输出1.输出数据到 dat 文件：*NODE PRINT, NSET=nset_name，FREQ=1COORD得到的是变形前的坐标还是变形变形后的坐标？偶在前面已发文问过，偶有时得到的是变形前的，有时得到是变形后的，一直没有弄明白这个问题。偶想得到的是变形后的坐标2. 其实 abaqus 自己就带有相关的功能:abaqus job=job-1

20、 suspend.可以将计算挂起.如果你需要重新进行运算输入 abaqus job=job-1 resume.3.如何输出大量节点的时间历程曲线？在环境文件 abaqus_v6.env 中添加一句max_history_requests=0 即可。 12.后处理积分CAE 自带此功能，比如对速度时程积分得到位移时程，可以这样操作：利用 XY Data-ODB History Output 将速度时程存为 V，然后 XY Data-Operate on XY Data-Operations看到那个 integrate (X)这个可能还是只是实现对时间的积分，如果对 dx，dy 的积分如何实现呢？

21、比方说，Q=v*A已经知道流速的矢量分布，怎么得到流量的分布了！这也是可以的，关键在于你的 x-y data 数据的两轴数据是什么？其积分就是 y 数据值对 x值的积分！看看：CAE user manual ：33.4.4Overview of X Y data operations13.接触分析激活杀死Usage: Use both of the following options:*MODEL CHANGE, TYPE=CONTACT PAIR, ADDsurface_1, surface_2*CONTACT INTERFERENCEsurface_1, surface_2, Usage:

22、 *MODEL CHANGE, TYPE=CONTACT PAIR, REMOVE14.利用 queue 的功能由本地机器向远程 unix 机器提交abaqus 作业的方法 精华利用 queue 的功能由本地机器向远程 unix 机器提交 abaqus 作业的方法假设：1. 远程 unix 机器的安装目录为 /abaqus，2. 远程计算机 IP 地址为 192.168.0.166，主机名为 ibmlinux3. 本地计算机 IP 地址为 192.168.0.18，主机名为 training下面的讨论涉及的内容相应改变Step 1 分别设置本地计算机和远程计算机之间的主机名和 IP 地址对应。

23、1. 修改本地 hosts 文件，windows 上一般位于 c:windowssystem32driversetc 下，unix 上位于/etc 下，加入远程计算机主机名和 IP 地址对应行，如：192.168.0.166 ibmlinux2. 修改远程 hosts 文件，位于/etc 下，加入本地计算机主机名和 IP 地址对应行，如：192.168.0.18 trainingStep 2 设置本地计算机对远程计算机运行 rsh 和 rcp 命令正确。1. 两台计算机登陆名和密码一致。举例而言，若远程计算机登陆名和密码为 root/root，则设置本地计算机也用同样的用户名和密码登陆。2.

24、修改远程计算机.rhosts 文件，该文件位于对应用户名登陆后的主目录下，如/home/root，加入本地计算机的 IP 地址使得本地计算机可以访问到远程计算机。3. 在本地计算机的提示符下键入rsh ibmlinux ls测试 rsh 是否能够正常使用，如果可以列出远程计算机主目录下文件，代表 rsh 和 rcp 工作正常。Step 3 设置远程计算机该用户的默认登陆程序为 csh，修改 passwd 文件，位于/etc 下，如：root:Ty91eFGzybEE2:0:3:/:/usr/bin/cshStep 4 修改远程计算机的 abaqus 环境文件 abaqus_v6.env在环境文

25、件最后加入queue_name =“nice_queue“after_prefix = “-b“queue_prefix = “-q“nice_queue = “/bin/sh -c nice ./%S 1%L 2 而在增量理论中， S=D*E(在有限变形中，其实应该为应力的客观率） - 似乎不对吧 大变形下此 D 非彼 D 你看过黄克智的固体本构关系这本书么 如果你从全量理论和增量理论的角度上讲 那似乎第一个 Digkl 就不对 你有第一个式子么 如果有，求导不久完了？ : 那应该就是一样的，因为全量理论， : Sij=DijklEkl（满足张量求和约定） : 即 Stress=D*Stra

26、in; : 而在增量理论中， : S=D*E(在有限变形中，其实应该为应力的客观率） 是啊，大变形下的D 与普通意义下的 D在构型上是不一样的， 毕竟D大是变形历史的函数，而D 小则不是， 我推导一种新的本构关系， Sij=DijklEkl (其中为 Jaumann 率） 假设材料一开始就屈服（即屈服面为 0） 想用壳单元， : 似乎不对吧 : 大变形下此 D 非彼 D : 你看过黄克智的固体本构关系这本书么 : 如果你从全量理论和增量理论的角度上讲 : 那似乎第一个 Digkl 就不对 : 你有第一个式子么 : 如果有，求导不久完了？ : 是啊，大变形下的D与普通意义下的D在构型上是不一样的

27、， : 毕竟D大是变形历史的函数，而D小则不是， : 我推导一种新的本构关系， : Sij=DijklEkl (其中为 Jaumann 率） : 假设材料一开始就屈服（即屈服面为 0） : 想用壳单元， 唉，别提了，问题就出在，在 abaqus 中，明明写着可以考虑剪切效应， 可我打印出剪切力个数是，nshr=1,即只有 S12, 那我的 S13,S23 就不知怎么计算， （DDSDDE(5,5)无法计算，因为 ntens=3,最多只能计算 DDSDDE(3,3) 你编 umat 编进去不久行了 他让用 NDI,NSHR,NTENS 表示变量，你就用这些表示变量 这样他就可以任意的计算了呀，而

28、不在意实际计算的变量数 : 唉，别提了，问题就出在，在 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 也存在同样的情况 在中厚度板壳元中(

29、MINDLIN)中，DDSDDE 为 55 你的 FOR 文件收到，可惜研究的对象相差太大，看不懂！ 不过你在文件中定义的变量有的根本就没有用到， 还有有的地方似乎是毫无用处的代码，比如：SMISES 一段 我觉得你不应该太急，第一步应该保证代码的正确性。 另外看了你上面的讨论，我觉得你应该换一种单元试一试。 比如用梁单元。很有可能 NDI 和 NSHR 会不一样。 不知道你是不是一定要用到壳单元上？ 如果必须用在壳单元上，可能 S23，S31 就必须自己在 UAMT 中输出了。具体的做法我也不清楚。 还有对 DDSDDE，我的理解也是: dS=DDSDDE*dE 我的意思是不管他是几乘几的矩

30、阵 一般的 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,nsh

31、r=1,这样的话，STRESS 为 3 个，STRESS(3),不可能计算 : STRESS(4),STRESS(5),同理，STRAIN,DDSDDE 也存在同样的情况 : 在中厚度板壳元中(MINDLIN) 中，DDSDDE 为 55 对于 DDSDDE，其实就是给了一个应力应变的斜率，我感觉 真正计算的是利用更新的应力，DDSDDE 我觉得就是为了确保收敛 真开心有这么多人讨论 ABAQUS，特别是 umat 我国外的师兄说 Umaterial 在国外应用得也很广 : 对于 DDSDDE，其实就是给了一个应力应变的斜率，我感觉 : 真正计算的是利用更新的应力， DDSDDE 我觉得就是为

32、了确保收敛 : 真开心有这么多人讨论 ABAQUS，特别是 umat : 我国外的师兄说 Umaterial 在国外应用得也很广 各位大侠，看到我提出的 uamt,及 statev 有这么多人感兴趣，实在很 高兴，希望大家能把摸索到的东西介绍给大家，现在我的程序能通了， （因没有屈服面的概念),虽然结果是错误的，但总可以静下心来，不像以前 一运行自己的 umat 程序就 standard.exe error.(并且程序.dat.msg.中无 任何错误信息，log 说 error see msg,but msg 中无） 。 现在我明白了，在用户程序中，WRITE(*,*)则写入 log,writ

33、e(6,*)则写入 dat,等， 在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为 0， 还是为上一步的值，因为我每次都为 0，怀疑有误 还望大侠提醒！ : dat,等， : 在监控程序时，发现，statev() ，是每次迭代和增量运行它，他就为 0， 应该不是 0，如果在每个增量步的最后你重新更新过他 输入文件*.inp 中，*depvar 不要忘了写 statev()的维数哦！ : 还是为上一步的值，因为我每次都为 0，怀疑有误 : 还望大侠提醒！ - 你只有一个单元么？ 【 在 saturn (好好工作) 的大作中提到 : 】 : !这我以前还不知道 : 应该不是 0

34、，如果在每个增量步的最后你重新更新过他 : 输入文件*.inp 中，*depvar 不要忘了写 statev()的维数哦！ - 【 在 air1 (无影如风勤奋羊羊课题在即) 的大作中提到: 】 : 可是这个写的命令有用么 : 对于输出而言 : 你只有一个单元么？ 已经写了，定义了 statev 的维数， *depvar 20, 看来我的错误不在 statev, 因为我的输出显示在第三增量步中，statev 不为 0， 但我还是有些不清楚 uamt 程序的执行过程， 比如， 一板一次壳元 s4r，面内采用缩减积分，即 11，厚度方向为 5 点 simpson 积分，按道理， 没一增量步的迭代步

35、，对每一积分点都会运行一次 umat, 即上问题，每一迭代该运行 5 次， 那上面的 statev,是每次迭代都会改变（不是每次运行）？ 请 air1,beauty_life,saturn,等各位相助 看过书上例子， 如果定义 statev 为（nblock，*）的数组， 应该对每个积分点都进行迭代计算，从而更新变量值 （写 stateNew(nblock,*) 请问 Mikejwg，变量 statev 怎么输出？ 可以象应力应变那样进行绘图、动画等后处理吗？ 不好意思，我只在 log 文件中输出，还没有在 odb 中输出， 想同 element output 等那样输出，只能在 UVARM

36、中定义，但 同 statev 是不同的，老兄你可以试试。 发信人: Mikejwg (mikejwg), 信区: FEA 标 题: 用 abaqus 中的 UMAT 的大侠看过来。 。 。 发信站: BBS 水木清华站 (Tue Dec 25 15:36:27 2001) 各位大侠， 在加入自己弹塑性本构模型时，如何加入区分弹性部分与塑性部分（或 非弹性部分） 。也就是说，在输出时，EE,LE,IE(PE)，因为在 umat 中， 我们只发现 dstran,stress,除此并没有发现其它有关应力应变得变量。 你所说的区分弹性部分和塑性部分是指应变的弹塑性部分：E=EE+PE 吗？ 还是指加载

37、过程的弹塑性？ E=EE+PE 我觉得想从 umat 的输入参数里面区分是不可能的 可不可以你先算出弹性的，然后用 stran 减去弹性的，得到塑性的 我猜的哦 我自己是明白的，但不知怎么让程序明白 比如，在一定条件下可解出位移增量，abaqus 得到应变增量 dstran,stran, 而弹性部分，dEEij=dSij/(2G),dEPij=dstranij-dEEij 但怎么控制程序了？ 默认得情况又是怎样计算的 谢谢 我不清楚你的本构， 我大概觉得你如果是用屈服面的话 你是不是可以用 statev（）记录每一步的 EE，EP 然后通过应力偏量来计算 DEE,DEP，然后自己更新 EE，E

38、P，就不与 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 水

39、木清华站 (Thu Dec 27 14:19:27 2001) 用 uvar()勉强成功 标 题: 在 abaqus 中，DSTRAN,STRAN 为工程剪应变吗？Sat Dec 29 09:44:03 2001) 同题，这涉及到 ddsdde,peeq 的计算 对，我也是最近才发现的 peeq 怎么计算？ 在 umat 中计算 peeq,(因用 umat 后，系统缺省 peeq=0)然后再 uvarm 输出 uvarm 可以输出到哪里？.odb 可以么？ 另外那个 strav？就是自己定义用于 umat 的那个数组里的数可以输出到.odb 里么？ 谢了，你这个成长中的大侠 在 umat 中，

40、statev 是不能用在 output 中的，statev 只是作为一个解的 状态变量，说来惭愧，我是在 umat 中定义 peeq,(peeq）的求解一般在 弹塑性力学书上有（等效塑性应变） ，用 write()写入一个临时文件， （切记：这个文件 unit 号不要与 abaqus 中的重合，因为他有一些系统默认 的文件号， ） ） 然后在 uvarm 中读取，以 uvarm 输出，因为 uvarm 可以以 odb 的形式输出， 支持 output,field,output,history 不知谁还有更好的方法。 听说清华有人 abaqus 用的很熟了。 大家过奖了，我是没办法，要出论文，总

41、得学习，上午版主还在 qq 中指点了我一番 YOU CAN TRY AS FOLLOWS: *STEP,EXTRAPOLATION=NO,INC=2000000 *STATIC 0.001,500.0,0.00001,0.1 可以 并且在 step 里好像有选项选择 你是手写的程序么 新建一个 step 然后定义边条 在上一 step 定义边条的节点上重新定义边条 而且，好像，是不是有个选项是 new 这样就可以取消上一 step 的边条设置了 在以后的载荷步里重新定义边界条件，加参数 new 比如在第二个载荷步里： *STEP *BOUNDARY CONDITIONS,OP=NEW *END

42、 STEP Command 方式：abaqus job=job_name user=source_file（或者 obj_file); CAE 方式：在最后 submit 的时候可以选择 source_file（或者 obj_file). 不太明白你的意思。 如果只是要得到所需要的如应力，应变等，那么在数据行请求 输出就行了啊，而且个人觉得 abaqus 的后处理功能基本上也够用了，直接输出到 *.odb 文件 中，用 CAE 或者 Viewer 后处理即可。 而*el print 输出的 dat 文件包含其他的信息，如版本号、处理内容等。 若想要直接输出到其他后处理软件的话，可以输出到 re

43、sult file 文件*.fil，不过， 这可是二进制文件或者 ASCII 文件，还是需要处理的，或者直接用 Hypermesh 等后处理 程序打开做图。 好像是在 tool-display-create再选定一个单元集或者是节点集 然后操作是“取代” 还是其他等等（就是里面的两个圆圈相交， 还是想减） 。你试试看看。 这个是要到后处理才行吧 非要计算以后形成*.odb 吗 前处理里面怎么才可以啊 Python 在 abaqus 中可以运行 scrpit 文件，你在 CAE manual 中就 可以将其附录中的例子存为*.py 文件，然后在 CAE 的菜单中 run script 就可以运行

44、了。 BTW:请问你的是 abaqus 的那个版本？ 6.2 or 6.3？6.3 的有盗版了？ : (1),请问如何得到 M,C,K 矩阵?用什么命令? 试一试：*element matrix output : (2),想要用 Newmark 方法求解.求解问题,在 ABAQUS 中如何实现?或者 : 那里能找到类似的例子? 看看：standard manuals:analysis 的 dynamics analysis 和 explicit manuals : (3),ABAQUS 中能实现这样的东西吗?就是我需要平滑 ABAQUS 产生的位移场 : 还是这个平滑只能在其他环境中完成 ?

45、试试这个输出选项： *El print,Position=average at Nodes : abaqus/CAE 好像不能生成节点集，只能定义某个实体，后处理用这个实体上的节点。 : 看到 CAE 的例题都是这样做的。 NODE SET 是 abaqus 的基本功能，也是其方便使用之处， 在 CAE 中是支持的。 这个例子我做过，没有问题。 : : 我用 toolsetcreate 居然不能选择由 CAE 生成的节点（不是顶点） 这是因为你没有正确操作 Partition Edge 这一步。 【 在 sunhaolan 的大作中提到: 】 : abaqus/CAE 好像不能生成节点集，只能

46、定义某个实体，后处理用这个实体上的节点。 ABAQUS divides all of its analysis procedures into two main groups: linear perturbation and general. ABAQUS makes a very clear distinction between perturbation and general analysis procedures because loading conditions and “time“ are defined differently for the two cases. Furthe

47、rmore, the results from each type of procedure should be interpreted differently. The response of the model during a general analysis procedure, known as a general step, may be either nonlinear or linear. In a step that uses a perturbation procedure, which is called a perturbation step, the response

48、 can only be linear. ABAQUS treats such steps as a linear perturbation about the preloaded, predeformed state (known as the base state) created by any previous general steps; therefore, its capability for doing linear simulations is rather more general than that of a purely linear analysis program. 各位大侠，我在定义 step 的时候，前三个用的是 standard 中的相关命令.在 step 4 的时候，我想用 explicit 来求解，并且定义了 amplitude,来模拟一个 位移随时间的简单变化，命令如下： *step name=“，nlgeom,inc=1000(好像 explicit 里没有 inc 这个参数） *dynamic,explicit ,2 *am