下载文档原格式
用Abaqus 所遇到问题汇总（持续更新）2011-02-17 02:45:26 转载▼ 标签：分类：好好学习abaqus杂谈问题1：当Verification 所有都pass 的时候，仍然提示Problem during compilation - ifort.exe not found in PATH解决办法：找到ABAQUS 安装目录下的Commands 文件夹例如D:SIMULIAAbaqusCommands下的abq6101.bat，右键，编辑此文件，插入下面这行使之成为第一行：8 P O2 b W R4 y6 U5ucallquotX:yourdirIntelCompilerFortranversionIA32Binifortvar s.batquot，例如我的是：- k q V/ : EcallquotC:ProgramFilesIntelCompiler11.1070binia32ifortvars_ia32. batquot f9 G9 R C0 Lamp quot d/ Y d问题2：当使用UMAT 子程序是出现以下错误Error in job Job-line44: 630 elements have been defined with zero hourglass stiffness. You may use hourglass stiffness or change the elementtype. The elements have been identified in element setErrElemZeroHourGlassStiffness.解决办法：由于设置了减缩积分，所以出现沙漏现象，将其改成全积分或imcompatible 可解决，详细解析在《基于ABAQUS 的有限元分析和应用》的第510页。
问题3：提交作业后模型出现问题，standard.exe 停止工作，只生成dat 文件而没有找到msg 文件解决办法：黄色图标的文件即msg 文件，但文件类型显示为outlook，用记事本打开即可。
四大要务第一要务：学习如何发帖研读《ABAQUS版发帖必读》和《ABAQUS版版面小提示【新手必看】》，学会如何发帖。
“入乡随俗”是世间和谐的铁律。
所以你初来aba版第一件事应该是看这两贴。
有人用“紧急求助”“救命”为题，让人很是生气：别人光看题目也不知道你是什么问题，就算解决了也不能为后人搜索到。
所以，题目一定要显示问题之所在。
另外，论坛规则，不得讨论盗版问题，也请大家遵守。
第二要务：下载顶置的《ABAQUS常见问题汇总-2.0版》这是aba-aba斑竹根据常见问题而精心总结的，也是论坛里面各位高手集体智慧的结晶。
就我的经验，新手所提的大部分问题都能在其中找到解决方案，这样的话，节约了版面，也使问题解决的速度加快，利版利民！我以前也发过一些2.0里面有板书的帖子，后来仔细读了2.0，对它佩服的五体投地，真实我们新手的指明灯啊！有时，有些简单问题高手不高兴搭理，就是因为这类问题被问过太多次，比如如何给实体单元加弯矩啊（实体单元没有转动自由度得先建reference point，然后耦合）；如何施加随时间变化的载荷（amplitude）2.0中有板书。
第三要务：版内搜索和第二步一样，很多问题早已经得到解决，所以遇到问题（称之为你的“盲点”）可以先版内搜索，一样是加快问题解决速度。
而且，这样你还会发现一个“副产品”：由于同样是这个问题，别人提起的方式和角度不一样，各位高手解答的方式和角度也是各异（虽然异曲同工），你都浏览一次，能使你对这个盲点有全方位的认识！盲点快速变成两点。
比如你的帮助文件不能搜索了该怎么办？（我总结过一次）。
出现“应变速度大于波速”怎么办？----版上已经有Robert大侠的总结了，新手多半是第一条（单位问题）。
第四要务：学会看帮助文件Documentation其实帮助文件时最好的辅导教材，可惜是English版，很多人望而生畏（包括我），在论坛下载了好多中文的书籍和例子学习，后来经论坛里面的各位版主提醒就硬着头皮去看Documentation。
abaqus问答精华（ABAQUS question and answer essence）Q: pre tensioned prestressed reinforced how, please pointing.Q: I see in the document that you want to define a rebar in the InP file, but rebar can only be used for shell, membrane, and, solid, elements. I want to do now is a prestressed cable, not embedded in the shell, membrane, and solid the unit, a zipper but independent. The cable unit is going to use truss, but how do you use rebar on the truss? Please expert pointingAnother problem, I see someone else's InP file, is as follows: *rebar, element=continuum, material=rebar2, name=ubarTop1,1.005e-4,0.15,0.0,0.5,1The second line, the first is setname (Top1), and the second is the cross-sectional area (1.005e-4) of the rebar. What do the third, fourth and fifth mean? (0.15,0,0.5) the last one should be the direction, the 1 direction. Which expert pointed out the first three or four and five respectively?A: exerts prestressing force*initinial, conditions, type=stress, rebarElset, rebar name, the prestress value, the other prestress hold to keep the prestress value unchanged, my understanding is to prevent other members eat the prestress, caused by prestress loss. With this command, this loss is avoided, and the applied prestressing force is applied to the reinforcement.A: thank you, what you say should be put in prestressed rebar above, but I find that the truss unit cannot be defined as rebar, in fact, I think that truss would have when cable engineering prestressed steel strand tension is to play a role, tension, and in ABAQUS, truss itself is taut, without prestressA: I know simulation reinforced when the need to use rebar, but the bar you can do directly with truss to simulate, and LZ said actually just during the construction of prestressed tension is not the real meaning of prestress, such as prestressed bolt etc.. If it is a cable, it may be a mere personal view of the prestressing force.Q: ask: do a space steel frame structure, beam column with beam element, board with shell element, intended to use tie command (common node), but I don't know how to achieve?A: I think you can use *equation to implement the constraints of the shared node, which is defined by this command.A: because I order more, but I want to use CAE, enter the command editor in CAE, and then edit it, write the *equation command, specifying the constrained degrees of freedom (see the standard manual, written very clearly)Where is the difference between "Q:" Response spectrum analysis (response spectrum analysis) and "Modal dynamic analysis" (modal dynamic analysis)? Such as Response, spectrum, analysis, you can design the structure? But what are the Modal dynamic analysis for?A: as far as I know, modal, dynamic, analysis should be modal decomposition method for dynamic analysis. Decomposition of a single degree of freedom system and a finite combination reaction.How does Q:abaqus apply seismic loads?A: can refer to the example of ABAQUS 6.3, Seismic, Analysis, of, a, Concrete, Dam, Gravity, and so onCan use:1. *amplitude, name=amp, and input=seismicdata.dat input seismic waves2. *boundary, type=acceleration, amplitude=amp to apply the load.In <ABAQUS, Example, Problems, Manual>, 2.1.15, Seismic, analysis, of, a,, concrete, gravity, damA: that's a relative problem. You can deduce that kinetic equation,The result is that ANSYS is the relative time history that takes the foundation to fix the relative base of the structure, and the ABAQUS is applied to the boundary, and the absolute time history of the structure is solved. Thus, there is a difference between the applied acceleration and the time history.Applied on the boundary time history acceleration by explicit calculation speed (mainly depends on your CPU ABAQUS STANDERD) which is the explicit calculation, in addition to its EXPLICT module can also see specific examples, analysis of the of a Concrete Gravity Dam Seismic Analysis; ANSYS is used for solving implicit (mainly depends on you the hard disk speed and CPU), and the big storm (the default output should be a lot of data for it), ANSYS/LSDYNA can be used to explicitly calculate.I have worked together for time analysis, with a little experience.Q: my boss asked me to use the concrete damage plasticity model1. In nonlinear analysis, how do I know the damage to the structure and how to obtain the ultimate bearing capacity of the structure?2. How do you determine concrete, compression, damage, and concrete tension damage? None of the two seems to have been tested. If you can ignore compression damage, then the definition of tension damage is important because of the presence of the tension stiffening (the falling segment is relatively long). In ABAQUS's example of gravity dams, the tension, stiffening, and tension damage data seem to have some relationship.3. When defining concrete tension stiffening, I read some of the papers above about the Tesion stiffening specification, saying that if tensiong stiffening is defined too small, convergence will be difficult. However, the size of tensionstiffening is defined in *concrete tension stiffening, and the first two data are remaining, direct, stress, after, cracking, and direct cracking strain.4. Does the tension stiffenig have anything to do with the rate of reinforcement? If yes, what is the relationship like? Thank you all.I think A: damage plasticity model is mainly used in concrete under cyclic loading, the concrete compression damage and concrete tension damage only when tension stiffening seems to need to define, the definition of concrete in uniaxial tension constitutive relation, tension and damage should Never mind. It doesn't seem to have anything to do with the reinforcement ratio.塑性损伤模型既可以模拟循环拟静荷载的情况，也可计算结构在动力荷载（地震、风以及冲击荷载的作用），通过混凝土压缩损伤和混凝土受拉破坏考虑混凝土的拉压异性，张力加劲则可考虑拉伸屈服后的软还阶段本构行为，同样是混凝土损伤全过程的一部分，它与配筋率及网格划分密度的相互关系很大程度上影响到计算的收敛性。
ABAQUS常见问题汇总 - 1.0版作者：aba_aba2006.06.14在SIMWE论坛注册5个月以来，在ABAQUS版回帖860篇，知无不言，言无不尽。
在此整理一下，以方便大家查找。
下面绝大多数是我的回帖，有少部分是其他网友的回答，都注明了作者。
下面只摘录了帖子中的一些主要内容，有些地方可能上下文不太连贯，完整的讨论请大家根据相应链接去论坛上察看。
如果这些内容对你有帮助，希望你也能常上论坛来回答一下别人的问题。
众人拾柴，火焰才高。
目录点击小节标题，可以跳到相应的内容（有些WORD版本可能需要按住ctrl键）1. 论坛提问要诀 (4)2. 故障诊断基本方法 (4)3. 约束刚体位移 / NUMERICAL SINGULARITY (4)4．过约束（Overconstraint） (7)5．其他常见error和warning信息 (7)6．ABAQUS无法运行（安装、License、环境变量） (10)7．弹塑性分析 (14)7.1 塑性问题不收敛的常见现象 (14)7.2 接触问题和塑性材料不要用二阶单元 (14)7.3 不要在塑性材料上施加点载荷 (15)7.4 塑性材料参数 (15)7.5 屈服准则 (18)7.6 后处理 / 判断材料是否屈服 (20)8. 接触分析 (22)8.1 接触分析不收敛的常见现象和解决方法 (22)8.2 接触面上的网格密度 (22)8.3 接触面的法线方向 (23)8.4 过盈接触 (23)8.5 管土/桩土接触 (25)8.7 凹坑成型的接触问题 (34)8.8 刚体穿透 (37)8.9 接触力 (37)8.10 Explicit中的接触问题 (38)8.11 改变接触状态 (40)8.12 其他与接触有关的Error和Warning信息 (41)9. 各个量的单位 (42)9.1 各个量的单位要相互一致 (42)9.2 时间/增量步的含义 (44)10．导入 / 修改 / 求解inp文件 (45)10.1 把inp文件导入ABAQUS/CAE (45)10.2 运行inp文件 (46)10.3 修改inp文件 (47)11．用户子程序和FORTRAN (49)12．创建/修改part (50)13．在Assembly中定位 (51)14．划分网格 (52)15．选择单元类型 (54)16．截面属性和材料 (56)16.1 定义截面属性 (56)16.2 用OFFSET参数偏置shell的中性面 (57)16.3 材料属性 (57)17．载荷 / 速度场 (58)17.1 集中力 (58)17.2 线载荷 (59)17.3 面载荷 (60)17.4 Body Force / Gravity (60)17.5 随时间变化的载荷 / Amplitude (61)17.6 在实体单元上施加弯矩 / 旋转速度 (62)18．边界条件 (67)19．定义Set / 选取模型的局部 (67)20．弹簧 (68)21．约束（coupling, tie） (68)22．后处理 (69)22.1 显示局部坐标系上的结果 (69)22.2 绘制曲线（X–Y data） (69)22.3 field output / history output (70)22.4 显示结果 (71)22.5 输出shell积分点上的应力和应变 (73)23．多步骤分析 (74)23.1 重启动分析（Restart） (74)23.3 用*import 把计算结果传递到下一个过程中 (74)24．岩土分析 (75)24.1 初始地应力平衡 (75)24.2 基坑开挖 (77)24.3 钢筋混凝土（Rebar Layer和Embedded Element） (77)24.4 实体单元的截面力/弯矩/转角 (78)24.5 软土固结 (80)25．热分析 (81)26．动态分析 (83)26.1 Explicit分析 (83)26.2 分析结果与amplitude不一致 (84)26.3 固有频率 (84)26.4 *DYNAMIC分析 (85)27．屈曲分析 (86)28．裂纹 (87)29．连接单元（CONNECTOR） (88)30．惯性释放（INERTIA RELIEF） (88)31．非线性问题的算法 (89)32．ABAQUS和其它软件的比较 (89)33．使用ABAQUS帮助手册 (90)34．论坛的“全文搜索”功能 (90)35．关于找枪手 (91)1. 论坛提问要诀在论坛上的提问是否能够得到解答，在很大程度上取决于提问的方式。
abaqus问答精华Q:预拉钢筋怎样施加预应力，请各位指点~~~~Q:我在文档里看到要在inp文件定义一个rebar，但是rebar只能用于shell, membrane, and solid elements 。
我现在想做的是一个预应力拉索，不是镶嵌在shell, membrane, and solid 这些单元里的，而是独立的一根拉锁。
拉索单元打算用truss,但是怎样在truss上使用rebar啊？请高手指点还有个问题，我看到别人的inp文件，如下：*rebar,element=continuum,material=rebar2,name=ubartop1,1.005e-4,0.15,0.0,0.5,1第二行第一个是setname（top1),第二个是rebar的截面面积（1.005e-4)，那第三、第四、第五是指什么？（0.15,0, 0.5),最后一个应该是方向，是1方向。
哪位高人指点下第三、四、五项分别代表什么？A:施加预应力 *initinial conditions,type=stress,rebarelset,rebar name,所施加预应力的值 ,另prestress hold 为保持所施加的预应力的值不变，我的理解是防止别的构件吃掉所施加的预应力，造成所施加预应力的损失。
使用了这个命令之后就避免了这种损失，保证所施加的预应力都施加到了钢筋上。
A:谢谢指点，你所说的应该是把预应力加在rebar上面，但我发觉truss单元不能定义成rebar，其实是我多想了，tr uss本来就可以当拉索，实际工程中加预应力只是为了使钢绞线拉紧，起到张拉作用，而在abaqus里，truss本身就是拉紧的，不用施加预应力A:我知道模拟加强筋的时候需要用rebar，但钢筋确实可以直接用truss来模拟 ,而lz所说的预应力其实其实只是施工时的张力而并不是真正意义上的预应力，比如螺栓预应力之类的。
如果是索的话可能是要施加预应力的，仅个人看法。
【在 luxz (panda) 的大作中提到: 】: 哪位有什么好办法？: 感觉abaqus自己的前处理用得很不爽: patran可以么？我感觉abaqus/cae不错的呀,实体建模颇有pro/e的风格感觉patran也好不到哪去,太麻烦,也许patran功能更强一些(要不然也对不起它那硕大的身躯呀), 但对于不是很庞大的模型,根本体现不出来什么优势: air1大侠，本人给你作广告，为何不可？: 请问：弹塑性矩阵【D]与ddsdde有何联系，: 你用过板壳单元吗？stress=D*stran?d(stress)=ddsdde*d(stran)--那应该就是一样的，因为全量理论，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一段我觉得你不应该太急，第一步应该保证代码的正确性。
《ABAQUS有限元分析常见问答解答》常见问答汇总ABAQUS是一种常用的有限元分析软件，它可以用于模拟和分析实际结构的力学行为。
在使用ABAQUS进行有限元分析时，可能会遇到一些常见的问题。
以下是一些常见的问题和解答。
问题1：如何创建一个新模型？答：在ABAQUS中，可以通过两种方式创建一个新的模型。
一种是使用图形用户界面(GUI)创建模型，另一种是使用ABAQUS命令创建模型。
使用GUI创建模型时，可以使用预定义的几何图形和网格划分工具创建模型，然后定义各种物理特性和加载条件。
使用ABAQUS命令创建模型时，可以使用ABAQUS提供的命令语言输入各种命令以创建和定义模型。
问题2：如何定义加载条件？答：在ABAQUS中，可以通过两种方式定义加载条件。
一种是使用GUI的加载工具，在模型的几何图形上添加加载边界条件和约束。
通过加载工具可以选择各种加载类型，如力、压力、温度等，然后指定加载的大小和方向。
另一种方式是使用ABAQUS命令定义加载条件。
可以使用各种ABAQUS提供的命令语言来定义加载条件。
问题3：如何求解模型？答：在定义完模型和加载条件后，可以使用ABAQUS求解模型。
求解模型时，需要选择分析类型和求解选项。
ABAQUS提供了多种分析类型，如静力学、动力学、稳态热分析等。
在选择分析类型后，可以选择求解选项，如迭代求解方法、收敛准则等。
然后点击求解按钮，ABAQUS会自动计算模型的响应。
问题4：如何查看结果？答：在求解完成后，可以查看模型的结果。
ABAQUS提供了多种结果查看工具。
使用GUI时，可以使用图形查看器来查看结果。
图形查看器可以显示模型的几何图形、网格、应力应变分布等。
此外，还可以在图形查看器中绘制图表和动画来更直观地查看结果。
另外，也可以使用ABAQUS 提供的命令语言来查看结果。
可以使用命令语言输出结果数据以及绘制图表。
问题5：如何修改模型？问题6：如何提高模型的计算效率？答：ABAQUS计算模型的效率受到多种因素的影响。
ABAQUS常见问题汇总 -1.0版作者：aba_aba2006.06.14在SIMWE论坛注册5个月以来，在ABAQUS版回帖860篇，知无不言，言无不尽。
在此整理一下，以方便大家查找。
下面绝大多数是我的回帖，有少部分是其他网友的回答，都注明了作者。
下面只摘录了帖子中的一些主要内容，有些地方可能上下文不太连贯，完整的讨论请大家根据相应链接去论坛上察看。
如果这些内容对你有帮助，希望你也能常上论坛来回答一下别人的问题。
众人拾柴，火焰才高。
目录点击小节标题，可以跳到相应的内容（有些WORD版本可能需要按住ctrl键）1. 论坛提问要诀32. 故障诊断基本方法33. 约束刚体位移 / NUMERICAL SINGULARITY34．过约束（Overconstraint）55．其他常见error和warning信息56．ABAQUS无法运行（安装、License、环境变量）77．弹塑性分析107.1 塑性问题不收敛的常见现象107.2 接触问题和塑性材料不要用二阶单元107.3 不要在塑性材料上施加点载荷107.4 塑性材料参数107.5 屈服准则137.6 后处理 / 判断材料是否屈服148. 接触分析158.1 接触分析不收敛的常见现象和解决方法158.2 接触面上的网格密度168.3 接触面的法线方向168.4 过盈接触168.5 管土/桩土接触188.6 板料成形的接触问题218.7 凹坑成型的接触问题238.8 刚体穿透258.9 接触力258.10 Explicit中的接触问题268.11 改变接触状态278.12 其他与接触有关的Error和Warning信息289. 各个量的单位299.1 各个量的单位要相互一致299.2 时间/增量步的含义3110．导入 / 修改 / 求解inp文件3110.1 把inp文件导入ABAQUS/CAE3110.2 运行inp文件3110.3 修改inp文件3211．用户子程序和FORTRAN3412．创建/修改part3513．在Assembly中定位3514．划分网格3515．选择单元类型3716．截面属性和材料3816.1 定义截面属性3816.2 用OFFSET参数偏置shell的中性面3816.3 材料属性3917．载荷 / 速度场3917.1 集中力3917.2 线载荷4017.3 面载荷4117.4 Body Force / Gravity4117.5 随时间变化的载荷 / Amplitude4117.6 在实体单元上施加弯矩 / 旋转速度4218．边界条件4619．定义Set / 选取模型的局部4620．弹簧4621．约束（coupling, tie）4722．后处理4722.1 显示局部坐标系上的结果4722.2 绘制曲线（X–Y data）4722.3 field output / history output4822.4 显示结果4922.5 输出shell积分点上的应力和应变5023．多步骤分析5023.1 重启动分析（Restart）5023.2 子模型（Submodel）5123.3 用*import 把计算结果传递到下一个过程中51 24．岩土分析5124.1 初始地应力平衡5124.2 基坑开挖5324.3 钢筋混凝土（Rebar Layer和Embedded Element）53 24.4 实体单元的截面力/弯矩/转角5424.5 软土固结5525．热分析5626．动态分析5726.1 Explicit分析5726.2 分析结果与amplitude不一致5826.3 固有频率5826.4 *DYNAMIC分析5927．屈曲分析5928．裂纹6029．连接单元（CONNECTOR）6130．惯性释放（INERTIA RELIEF）6131．非线性问题的算法6132．ABAQUS和其它软件的比较6233．使用ABAQUS帮助手册6234．论坛的“全文搜索”功能6235．关于找枪手631. 论坛提问要诀在论坛上的提问是否能够得到解答，在很大程度上取决于提问的方式。
清华牛人的ABAQUS总结[5篇模版]第一篇：清华牛人的ABAQUS总结准确的说，应该是谈谈我对发在这里的帖子质量的看法。
因为专业的原因，我对有限元及其软件是有很深感情的。
又因为一直对清华深有好感，可惜因为一些原因最终没有来清华深造，所以比较关注清华的BBS。
写这篇文章要耗费我差不多一个完整的下午，但是我愿意。
我知道学有限元其实不是一件容易的事情，我把我的想法说出来，希望对初学者有所裨益。
坦率的说，我认为这里有限元板块的质量是不高的。
之所以如此，是因为在这里很多人问的问题是太简单而且对自己不负责任的。
这不是版主的错，是因为我们许多人还没有养成良好的专业素养和严谨的精神。
请不要轻易的否认我的这个评价，好吗？至少，现在请不要。
就是对我的话不屑一顾，也应该是在看完我的话之后吧。
我用我的思维方式来说话，并不是每一个人都会习惯，请见谅！我还要声明的是，我本人的水平一般般，自己也对自己有很多的不满，所以在这里说的很可能很幼稚或者有错误。
请大家指教！我们应该有一个良好的讨论气氛。
有限元对许多工科的人而言，其必要性和重要性不言而喻。
问题在于，应该怎样的学习它呢？学习它，至少不用它到处害人也害己的话，我觉得至少要在下面四个方面有些基本知识：1、有限元基本理论及其求解基本步骤（数学基础）；2、有限元专业英语（英语基础）；3、你自己所属专业的东东（专业基础）；4、几何造型及拓扑学知识（建模基础）。
这个排序是由重到轻的。
接下来，我首先说一说上面四个方面的意义和作用；之后谈一下为什么我认为在这里问的相当一部分问题是太简单而且对自己不负责任的。
1、做专业就要有做专业的样子。
咱们理工科的学生，没有辛苦的付出是不可能有真正收获的。
收获和付出在这里成正比。
常常有人觉得有限元的软件很难，不好学，不好用，很多东西搞不懂，一提就头痛。
其实这里面相当的一部分是有限元基本理论可以解决的问题，而不是软件的设计思想不好。
现在的商用有限元软件，比如我用过的abaqus，ansys，adina以及algor，应该说它们的界面已经很友好了，包括帮助文档等等都不错。
ABAQUS学习技巧总结（转帖）第一篇：ABAQUS学习技巧总结(转帖)ABAQUS学习总结1.ABAQUS中常用的单位制。
-（有用到密度的时候要特别注意）单位制错误会造成分析结果错误，甚至不收敛。
2.ABAQUS中的时间对于静力分析，时间没有实际意义（静力分析是长期累积的结果）。
对于动力分析，时间是有意义的，跟作用的时间相关。
3.更改工作路径4.对于ABAQUS/Standard分析，增大内存磁盘空间会大大缩短计算时间;对于ABAQUS/Explicit分析，生成的临时数据大部分是存储在内存中的关键数据，不写入磁盘，加快分析速度的主要方法是提高CPU的速度。
临时文件一般存储在磁盘比较大的盘符下提高虚拟内存5.壳单元被赋予厚度后，如何查看是否正确。
梁单元被赋予截面属性后，如休查看是否正确。
可以在VIEW的DISPLAY OPTION里面查看。
6.参考点对于离散刚体和解析刚体部件，参考点必须在PART模块里面定义。
而对于刚体约束，显示休约束，耦合约束可以在PART ,ASSEMBLY,INTERRACTION,LOAD等定义参考点.PART模块里面只能定义一个参考点，而其它的模块里面可以定义很多个参考点。
7.刚体部件（离散刚体和解析刚体），刚体约束，显示体约束离散刚体：可以是任意的形状，无需定义材料属性，要定义参考点，要划分网格。
解析刚体：只能是简单形状，无需定义材料属性，要定义参考点，不需要划分网格。
刚体约束的部件：要定义材料属性，要定义参考点，要划分网格。
显示体约束的部件:要定义材料属性，要定义参考点，不需要要划分网格(ABAQUS/CAE会自动为其要划分网格)。
刚体与变形体比较:刚体最大的优点是计算效率高，因为它在分析作业过程中不参与所在基于单元的计算，此外，在接触分析，如果主面是刚体的话，分析更容易收敛。
刚体约束和显示体约束与刚体部件的比较：刚体约束和显示体约束的优点是去除约束后，就可以立即变为变形体。
【在luxz (panda) 的大作中提到: 】: 哪位有什么好办法？: 感觉abaqus自己的前处理用得很不爽: patran可以么？我感觉abaqus/cae不错的呀,实体建模颇有pro/e的风格感觉patran也好不到哪去,太麻烦,也许patran功能更强一些(要不然也对不起它那硕大的身躯呀), 但对于不是很庞大的模型,根本体现不出来什么优势: air1大侠，本人给你作广告，为何不可？: 请问：弹塑性矩阵【D]与ddsdde有何联系，: 你用过板壳单元吗？stress=D*stran?d(stress)=ddsdde*d(stran)--那应该就是一样的，因为全量理论，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.exe error.(并且程序.dat.msg.中无任何错误信息，log说error see msg,but msg中无）。
现在我明白了，在用户程序中，WRITE(*,*)则写入log,write(6,*)则写入dat,等，在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，还是为上一步的值，因为我每次都为0，怀疑有误还望大侠提醒！: dat,等，: 在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，应该不是0，如果在每个增量步的最后你重新更新过他输入文件*.inp中，*depvar不要忘了写statev()的维数哦！: 还是为上一步的值，因为我每次都为0，怀疑有误: 还望大侠提醒！--你只有一个单元么？【在saturn (好好工作) 的大作中提到: 】: ~!~~~~~~~~~~~~~~~~~~~~~~~~~~~~这我以前还不知道: 应该不是0，如果在每个增量步的最后你重新更新过他: 输入文件*.inp中，*depvar不要忘了写statev()的维数哦！--【在air1 (无影如风※勤奋羊羊※课题在即) 的大作中提到: 】: 可是这个写的命令有用么: 对于输出而言: 你只有一个单元么？已经写了，定义了statev的维数，*depvar20,看来我的错误不在statev,因为我的输出显示在第三增量步中，statev不为0，但我还是有些不清楚uamt程序的执行过程，比如，一板一次壳元s4r，面内采用缩减积分，即1×1，厚度方向为5点simpson积分，按道理，没一增量步的迭代步，对每一积分点都会运行一次umat,即上问题，每一迭代该运行5次，那上面的statev,是每次迭代都会改变（不是每次运行）？请air1,beauty_life,saturn,......等各位相助看过书上例子，如果定义statev为（nblock，*）的数组，应该对每个积分点都进行迭代计算，从而更新变量值（写stateNew(nblock,*)请问Mikejwg，变量statev怎么输出？可以象应力应变那样进行绘图、动画等后处理吗？不好意思，我只在log文件中输出，还没有在odb中输出，想同element output等那样输出，只能在UV ARM中定义，但同statev是不同的，老兄你可以试试。
发信人: Mikejwg (mikejwg), 信区: FEA标题: 用abaqus中的UMAT的大侠看过来。
发信站: BBS 水木清华站(Tue Dec 25 15:36:27 2001)各位大侠，在加入自己弹塑性本构模型时，如何加入区分弹性部分与塑性部分（或非弹性部分）。
也就是说，在输出时，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 水木清华站(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中，statev是不能用在output中的，statev只是作为一个解的状态变量，说来惭愧，我是在umat中定义peeq,(peeq）的求解一般在弹塑性力学书上有（等效塑性应变），用write()写入一个临时文件，（（切记：这个文件unit号不要与abaqus中的重合，因为他有一些系统默认的文件号，））然后在uvarm中读取，以uvarm输出，因为uvarm可以以odb的形式输出，支持output,field,output,history不知谁还有更好的方法。
