清华大学abaqus论坛精华帖

合集下载

用Abaqus所遇到问题汇总

用Abaqus所遇到问题汇总

用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,用记事本打开即可。

给simwe论坛新手的一些建议和经验

给simwe论坛新手的一些建议和经验

四大要务第一要务:学习如何发帖研读《ABAQUS版发帖必读》和《ABAQUS版版面小提示【新手必看】》,学会如何发帖。

“入乡随俗”是世间和谐的铁律。

所以你初来aba版第一件事应该是看这两贴。

有人用“紧急求助”“救命”为题,让人很是生气:别人光看题目也不知道你是什么问题,就算解决了也不能为后人搜索到。

所以,题目一定要显示问题之所在。

另外,论坛规则,不得讨论盗版问题,也请大家遵守。

第二要务:下载顶置的《ABAQUS常见问题汇总-2.0版》这是aba-aba斑竹根据常见问题而精心总结的,也是论坛里面各位高手集体智慧的结晶。

就我的经验,新手所提的大部分问题都能在其中找到解决方案,这样的话,节约了版面,也使问题解决的速度加快,利版利民!我以前也发过一些2.0里面有板书的帖子,后来仔细读了2.0,对它佩服的五体投地,真实我们新手的指明灯啊!有时,有些简单问题高手不高兴搭理,就是因为这类问题被问过太多次,比如如何给实体单元加弯矩啊(实体单元没有转动自由度得先建reference point,然后耦合);如何施加随时间变化的载荷(amplitude)2.0中有板书。

第三要务:版内搜索和第二步一样,很多问题早已经得到解决,所以遇到问题(称之为你的“盲点”)可以先版内搜索,一样是加快问题解决速度。

而且,这样你还会发现一个“副产品”:由于同样是这个问题,别人提起的方式和角度不一样,各位高手解答的方式和角度也是各异(虽然异曲同工),你都浏览一次,能使你对这个盲点有全方位的认识!盲点快速变成两点。

比如你的帮助文件不能搜索了该怎么办?(我总结过一次)。

出现“应变速度大于波速”怎么办?----版上已经有Robert大侠的总结了,新手多半是第一条(单位问题)。

第四要务:学会看帮助文件Documentation其实帮助文件时最好的辅导教材,可惜是English版,很多人望而生畏(包括我),在论坛下载了好多中文的书籍和例子学习,后来经论坛里面的各位版主提醒就硬着头皮去看Documentation。

abaqus问答精华(ABAQUSquestionandansweressence)

abaqus问答精华(ABAQUSquestionandansweressence)

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版

ABAQUS常见问题汇总 - 1.0版

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问答精华1

ABAQUS问答精华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所说的预应力其实其实只是施工时的张力而并不是真正意义上的预应力,比如螺栓预应力之类的。

如果是索的话可能是要施加预应力的,仅个人看法。

清华大学Abaqus论坛精华贴(simwe)

清华大学Abaqus论坛精华贴(simwe)

【在 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有限元分析常见问答解答》常见问答汇总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总结[5篇模版]

清华牛人的ABAQUS总结[5篇模版]第一篇:清华牛人的ABAQUS总结准确的说,应该是谈谈我对发在这里的帖子质量的看法。

因为专业的原因,我对有限元及其软件是有很深感情的。

又因为一直对清华深有好感,可惜因为一些原因最终没有来清华深造,所以比较关注清华的BBS。

写这篇文章要耗费我差不多一个完整的下午,但是我愿意。

我知道学有限元其实不是一件容易的事情,我把我的想法说出来,希望对初学者有所裨益。

坦率的说,我认为这里有限元板块的质量是不高的。

之所以如此,是因为在这里很多人问的问题是太简单而且对自己不负责任的。

这不是版主的错,是因为我们许多人还没有养成良好的专业素养和严谨的精神。

请不要轻易的否认我的这个评价,好吗?至少,现在请不要。

就是对我的话不屑一顾,也应该是在看完我的话之后吧。

我用我的思维方式来说话,并不是每一个人都会习惯,请见谅!我还要声明的是,我本人的水平一般般,自己也对自己有很多的不满,所以在这里说的很可能很幼稚或者有错误。

请大家指教!我们应该有一个良好的讨论气氛。

有限元对许多工科的人而言,其必要性和重要性不言而喻。

问题在于,应该怎样的学习它呢?学习它,至少不用它到处害人也害己的话,我觉得至少要在下面四个方面有些基本知识:1、有限元基本理论及其求解基本步骤(数学基础);2、有限元专业英语(英语基础);3、你自己所属专业的东东(专业基础);4、几何造型及拓扑学知识(建模基础)。

这个排序是由重到轻的。

接下来,我首先说一说上面四个方面的意义和作用;之后谈一下为什么我认为在这里问的相当一部分问题是太简单而且对自己不负责任的。

1、做专业就要有做专业的样子。

咱们理工科的学生,没有辛苦的付出是不可能有真正收获的。

收获和付出在这里成正比。

常常有人觉得有限元的软件很难,不好学,不好用,很多东西搞不懂,一提就头痛。

其实这里面相当的一部分是有限元基本理论可以解决的问题,而不是软件的设计思想不好。

现在的商用有限元软件,比如我用过的abaqus,ansys,adina以及algor,应该说它们的界面已经很友好了,包括帮助文档等等都不错。

ABAQUS学习技巧总结(转帖)

ABAQUS学习技巧总结(转帖)

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会自动为其要划分网格)。

刚体与变形体比较:刚体最大的优点是计算效率高,因为它在分析作业过程中不参与所在基于单元的计算,此外,在接触分析,如果主面是刚体的话,分析更容易收敛。

刚体约束和显示体约束与刚体部件的比较:刚体约束和显示体约束的优点是去除约束后,就可以立即变为变形体。

hypermesh-abaqus论坛常见问题汇总

hypermesh-abaqus论坛常见问题汇总

HM——ABA接口问题简洁一些,引用小宝斑竹在接口问题中的体会:“关于hm-abaqus的接口补充说明经常用HM-abqus的人或者刚开始使用的人,对于软件的接口一般存在以下问题:1、INP文件导入abaqus出错。

2、在abaqus里选择加载面,设置材料属性不方便。

(也可以理解为没有几何元素的模型在网格上选取东东比较困难)除了以上的两种,暂时我还没发现其他的问题。

对于1,只要前处理没有除网格外的载荷信息,并且所有关键字名称都没有数字,那么恭喜你,它一定不会出错。

(当然,有些人说abaqus/cae有很多关键字不支持,但是hm支持。

我的建议是:有那功夫,或许INP文件都改完了)对于2,首先声明,所有在hm里建立的SET,都会出现在abaqusassembly 里的set里,所以在abaqus里加载的时候,都是可以调用的,你所要做的就是正确的建立nodesetor elementset。

很多人知道在abaqus part里也有个set,那个是干什么的呢?目前我就是用来操作材料施加的。

很多时候模型是各种复杂材料的混合,如果在abaqus里直接赋予,选取模型区域的手段有限(单个点,点到手抽筋;by angle,很多地方选不上,选到眼花;by set,在abaqus里建立partset的难度不下于手动点),我的建议:在hm里赋予一个空材料属性给相关的区域(hm里有几何元素,选起来简直就是小CASE),到了abaqus 里,软件会自动为你的每个材料区建立一个新的part set,这时候,空的材料属性要炸要炒随你便。

剩下的问题都不是问题。

”论坛问题汇总1、hypermesh导入abaqus有单元显示、无模型显示。

答:这个问题很常见,不仅在hypermesh_to_abaqus中有人问,在HM与其它软件接口也有人反复问。

首先要肯定的是导入过程没有任何问题。

在此,引用老向版主的一段话来回答“版上不停的有人问为什么HM不能导出几何.看的人都烦了.为什么要导出几何呢?* H/ _/ m' j; C? @不同的软件,对几何的理解是不一样的,所以才有这么多的不同的几何格式.; E0 H- x8 ?0 m5 D k如果要导出几何,HM还得去研究你abaqus/ansys/patran内部是如何理解几何的,这是个浩大的工程.- M) S0 M! \( \你应该知道,对于求解器来说,它只需要知道节点,单元,材料,载荷等信息就能够求解了.要几何干什么呢?' X- q3 w G) A6 H8 A5 j" i, d: \* i$ c3 E& ~( C6 x4 n" V# R2 I几何模型的作用仅仅是为了得到节点,网格..一旦有节点,网格有了,几何模型就可以扔掉了.后处理程序本身也是基于有限元模型的,而不是基于几何模型的.! D6 K6 C' ?7 r9 j8 g你既然打算用HM做前处理,就干脆一点,把所有的东西都在HM里面做好,然后提交给abaqus/nastran计算就行了.p4 l9 W, t! u9 X( }HM对abaqus求解器的支持一点不比abaqus/cad本身差,只有及少数特殊的场合,需要手工添加一些语句.”2、从HM中把一个装配件划分好的网格生成inp文件后,导入abaqus中去,为什么是一个零件,怎样才能是一个装配件?答:part, instance是abaqus独有的概念。

对于ABAQUS初学者的一些经验分享

对于ABAQUS初学者的一些经验分享

首先声明:个人原创,初学参考,高手拍砖!学习ABAQUS有一段时间了,这个软件真是让人又爱又恨(爱它的功能强大,恨自己不懂运用),想必初学者也会有和我一样的想法吧。

很多人苦恼自己建的模型不收敛,作业提交不上,尝试过几次就崩溃了。

这就是计算机软件,它只懂它的言语,稍有不慎,哪怕是个很低级的bug,也会让人纠结很久,但一旦解决,心中成就感油然而生,这也就是ABAQUS的魅力所在。

废话不多说,下面根据自己的经验,为初学者提出下面三点建议:1、多尝试、不怕失败ABA令人苦恼的莫过于他的错误信息(警告信息有时可忽略,不太重要),因为一旦出现,作业提交失败,就可能意味着前面的工作竹篮打水。

很多时候你自己也不知道到底是那个地方设置有问题,ABA 里面有那么多模块不说,每个模块里的菜单、子菜单不说,有时仅仅是一个单薄的对话框就能让新手纠结个半死,为什么呢?因为选项太多,不知道要改变哪个,又不知道要改成什么样的才合适。

自己拿着鼠标点了几百次甚至几千次(CAE建模)好不容易建成的模型就这样走不通了,没有任何结果。

走不下去了,怎么办?多尝试、再试试。

对于最开始接触ABA的新手,一般来说会照着自己已有资料中参考的步骤去完成一个模型,这是一个快速熟悉ABA这种操作性强的软件的好方法。

可是很多时候发现照着做也还是行不通,回想自己看着书盯着电脑一个一个选项选、一个一个参数输,还是没有结果,心中的挫败感会油然而生。

这个时候会有两种选择:第一、放弃这个例子,去完成另外一个,不必在一个歪脖子树上吊死;第二、抱着那颗歪脖子树不放。

我一开始学习ABA是按照庄拙先生的那本书上的例子操作的,做到第二个例子的时候就走不通了,后来我选择放弃并重新开始另一个例子,结果还是失败。

“欲速则不达”,一句古话一针见血。

当第三个例子继续“倒下”时,我开始反思自己是不是过于急躁,因为我觉得自己能把它做出来,不会在“傻瓜操作”上爬不起来。

这样下去是没有任何意义的,不仅浪费时间而且没有什么收获。

abaqus常用技巧总结

abaqus常用技巧总结

18.取消坐标系等的显示.................................................................................. 37 19.如何在计算中修改材料特性 ....................................................................... 37 20.输出计算过程中的总质量和总刚度矩阵 ..................................................... 40 21.先张预应力: ............................................................................................ 40 22.用户子程序的使用 ..................................................................................... 41 23.怎样设定用双CPU机器进行ABAQUS计算 ................................................ 41 24.中途停止正在运算的JOB .......................................................................... 42 25.自适应网格技术 ......................................................................................... 42 26.ABAQUS计算与内存 ................................................................................ 42 27.质量缩放 ................................................................................................... 43 28.ABAQUS多处理器进行并行计算的效果研究............................................. 54 29.YAHOO讨论组摘录--CONTACT+OVERCLOSURE .................................. 55 30.原创:无限元建立方法,希望得到加分 ..................................................... 70 31.[分享]ABAQUS 使用问答 ......................................................................... 76 32.[转帖]ABAQUS6.4 导入外来模型的几点小经验! .................................... 87 33.ABAQUS的多图层绘图 ............................................................................. 88 34.子结构 ....................................................................................................... 89 35.如何在不同的分析步改变材料的参数 ........................................................ 89 36.模型的重启动分析-RESTART ................................................................. 91 37. ABAQUS的单位心得................................................................................ 92

XFEMtutor问题汇总

XFEMtutor问题汇总

《ABAQUS6.9版本XFEM(扩展有限元)例子的详细图解step by step》帖子的问题汇总已做出解答部分1、问一下为什么要添加赋予材料取向,谢谢!答:①我也试了,材料方向可以不设置。

同时Damage Stabilization也必须设置,否则不收敛。

Damage Stabilization是为了增加收敛性的,如果一个问题的非线性程度比较低,比如损伤应力比较低,相对应需要的能量也比较少的情况下,很容易收敛,所以就不需要Damage Stabilization,如果问题属于高度非线性,不设置这些参数,问题就可能发散。

兄台感觉如何?②发现这里有个多余的步骤,设置方向,对各向同性材料截面设置中选择homogeous 后,不用再设方向!2、Material模块中的操作的“3.赋予材料取向”时看不到“在part Plate中创建的4个集合:all,bottom,top和fixZ”,为何?答:能看到all这个集合,看不到bottom, top和fixZ,个人感觉后三个集合只是面或者点,而包含材料的集合只能是一个体集合,所以在赋予材料属性的时候将后面三个集合过滤掉了。

不知道是否正确?回复:果真如此!3、集合bdisp是只包含db这个点,还是包含db这个点和装配体?答:应该是bdisp这个集合只包含db这个点,因为在接下来的定义接触方程时,用到bdsip集合作为第二行,必须是只包含点。

如下图的提示中所述。

另外:对于定义这个方程的作用以及对于计算结果的后处理,希望楼主能再给出些讲解。

(P2)4、关于参考点的问题①定义这个方程的作用是将一个面的运动与一个点的运动相联系。

具体到此例,即为:底面bottom在x(自由度1)方向上与点bdisp的运动一致。

因为在x方向上的载荷是施加在点bdisp点上。

②回复楼主:谢谢楼主的讲解!另外,不用参考点加载,去掉参考点与底部的方程约束直接将底部的x方向载荷加载在底面上,计算出来的结果是一样的,不知道设置这个参考点的意义何在?或许是考虑到收敛性的问题吗?③个人认为加这个参考点的作用是为了以后输出加载点的位移和反力用的。

abaqus自己遇到问题汇总

abaqus自己遇到问题汇总

Too many attempts made for this incrementThrere may be many factors that you should check. Some of them you might want to check predominantly are,1_Material properties and equivalency of units 2_ Mesh size and type 3_Boundary conditions 4_constraints such as rigid body motion 5_Step size and no of increments Also look in Abaqus documentation for Common problems in convergence of solution.自己的问题是由于模型没有出现几何非线性,不需要将Nlgeom设置为ON,我自己的做法大大增加了分析收敛难度。

下面这个问题我在SIMWE论坛里,看到,然后在百思论坛也看到, 发帖人是07年12月求助,到08年1月中旬都没人帮他解决啊.本人看到,于是粘贴下来.今天谈谈我的看法.问题:关于HETVAL子程序flux定义,通过自己编写的程序(fortran)盘算了每个单元的生热率,盼望将该盘算成果调入abaqus子程序中,并将该值赋给HETVAL子程序中的flux1变量,作为传热分析中每个单元的内热源,请教研发网的朋友怎么样给每个单元分辨赋给内热源呢?帮助文档的例子中好像是对某个材料赋给内热源的(*heat generation, HETVAL),要给每个单元的生热分辨赋值该怎么做呢?查下载说通过UEXTERNALDB子程序可以将外部成果调进来,查赞助文档没有找到有关可以把外部成果调进来的例子,不知哪位有这方面的经验或者例子, 马上要交论文了,可是温度场还没算出来,请大家帮帮忙!!!!!多谢了!论坛回答:(1)要给每个单元的生热分辨赋值就对每个单元定义material,在abq中,*heat generation被视为材料的特点,故*heat generation不同,就得定义多种材料.我的回答: simwe仿真论坛的回答固然没有错,但是谁会去为为一个单元定义一个集合,再去给只有一个单元的集合赋予材料属性呢??论坛上高手众多,我不知为什么就没人去回答这个问题,是他太简单了吗?未必吧.就算简单,但是在别人看来,却是困难,我们应该尊重别人的困难,虽然可能在你看来不是什么了不起的事.我在这里有个笨方法,一定可以成功的.第一,生热率的云图可以用UVARM输出.怎么将每个单元的生热率赋给指定的单元呢?请看:SUBROUTINE UVARM(UVAR,DIRECT,T,TIME,DTIME,CMNAME,ORNAME,1 NUVARM,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,NDI,NSHR,COORD,2 JMAC,JMATYP,MATLAYO,LACCFLA)NOEL为单元编号,即模型中单元的标签号,你可以这样if(NOEL==269) thenuvar(1)=160else if......elseuvar(1)=80end if这样就可以指定生热率给指定的单元.我这只是个简单的示范,其实可以while循环最佳啊生热率具体数值其实也可在这里生成,具体的程序还没有.熟悉FORTRAN的人都可以做出来.要用到输入输出代码.第二,还有一个是温度分析.hetval.其FLUX定义,也可借鉴上面的方法.SUBROUTINE HETVAL(CMNAME,TEMP,TIME,DTIME,STATEV,FLUX,1 PREDEF,DPRED)HETVAL是没有NOEL这个参数的.但是要注意查看他的文档.can be used in conjunction with user subroutine USDFLD if it is desired to redefine any field variables before they are passed in; and它是可以和USDFLD连用的,如果变量在传入HETVAL之前,需要用USDFLD进行重定义.而USDFLD函数定义如下:SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)这里有NOEL.而且它还可以调用辅助函数GETVRM.功能强大.第三,生热率的应力应变从文件中直接输出.在CAE中用report也可,在inp中用命令EL PRINT输出到DAT文件中也可以.经过一些处理(线性插值,傅里叶变换等),然后用FORTRAN程序去读取.在这里有个技巧就是单元数据有多个数据项,一个单元号,六个应力值,六个应变值,最好用TYPE和MODULE将其封装起来(结构体中最好把其应力应变的关系函数的值也定义一个量).调用起来才不乱.这是我最近查阅ABAQUS文档,学习FORTRAN,再加上论坛上晃荡,和思考所得的一些东西.敬请受用.希望对大家能有所助.关于子程序UVARM,HETVAL,USDFLD以及辅助函数GETVRM具体细节和应用,等过段时间, 有时间我用具体例子,再细细说来.敬请期待.Hetval 子程序丢失的问题:由于自己材料设置有误多设置了一个HEAT GENERATIONTime increment required is less than the minimum specified??????????出现这个错误可能是多方面原因造成的,你可以采用排除法逐个因素去找1.可以把初始增量步最小增量步调小,最大增量步的数目调大。

清华BBSabaqus经典问题解答

清华BBSabaqus经典问题解答
谢谢
我不清楚你的本构,
我大概觉得你如果是用屈服面的话
你是不是可以用statev()记录每一步的EE,EP
然后通过应力偏量来计算DEE,DEP,然后自己更新EE,EP,就不与stran扯上关系
我以前见过用Umat算J2流动理论的时候,是这么做的
首先我的本构不能采用屈服面的概念(一种新型的弹塑性本构模型),
dat,等,
在监控程序时,发现,statev(),是每次迭代和增量运行它,他就为0,
还是为上一步的值,因为我每次都为0,怀疑有误
还望大侠提醒!
: dat,等,
: 在监控程序时,发现,statev(),是每次迭代和增量运行它,他就为0,
应该不是0,如果在每个增量步的最后你重新更新过他
输入文件*.inp中,*depvar不要忘了写statev()的维数哦!
另外看了你上面的讨论,我觉得你应该换一种单元试一试。
比如用梁单元。很有可能NDI和NSHR会不一样。
不知道你是不是一定要用到壳单元上?
如果必须用在壳单元上,可能S23,S31就必须自己在UAMT
中输出了。具体的做法我也不清楚。
还有对DDSDDE,我的理解也是: dS=DDSDDE*dE
我的意思是不管他是几乘几的矩阵
我也看了那个采用j2流动法则的例子,我可以在umat中的statev重定义
ee,pe,但statev不能像EE,PE,IE那样在odb中输出。
再说在umat中的ee,pe与inp中的ee,pe是不同的。
谢谢
看来你只能试试uvarm了
呵呵,
发信人: Mikejwg (mikejwg), 信区: FEA
同statev是不同的,老兄你可以试试。
发信人: Mikejwg (mikejwg), 信区: FEA

清华大学Abaqus论坛精华贴(simwe)

清华大学Abaqus论坛精华贴(simwe)

【在 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用户子程序

Home浅谈ABAQUS用户子程序李青清华大学工程力学系摘要本文首先概要介绍了ABAQUS的用户子程序和应用程序,然后从参数,功能两方面详细论述了DLOAD, UEXTERNALDB, URDFIL三个用户子程序和GETENVVAR,POSFIL,DBFILE三个应用程序,并详细介绍了ABAQUS的结果文件(.FIL)存储格式。

关键字ABAQUS,用户子程序,应用程序,结果文件一、前言:ABAQUS为用户提供了强大而又灵活的用户子程序接口(USER SUBROUTINE)和应用程序接口(UTILITY ROUTINE)。

ABAQUS 6.2.5一共有42个用户子程序接口,13个应用程序接口,用户可以定义包括边界条件、荷载条件、接触条件、材料特性以及利用用户子程序和其它应用软件进行数据交换等等。

这些用户子程序接口使用户解决一些问题时有很大的灵活性,同时大大的扩充了ABAQUS的功能。

例如:如果荷载条件是时间的函数,这在ABAQUS/CAE 和INPUT 文件中是难以实现的,但在用户子程序DLOAD中就很容易实现。

二.在ABAQUS中使用用户子程序ABAQUS的用户子程序是根据ABAQUS提供的相应接口,按照FORTRAN语法用户自己编写的代码。

在一个算例中,用户可以用到多个用户子程序,但必须把它们放在一个以.FOR为扩展名的文件中。

运行带有用户子程序的算例时有两种方法,一是在CAE中运行,在EDIT JOB菜单的GENERAL子菜单的USER SUBROUTINE FILE对话框中选择用户子程序所在的文件即可;另外是在ABABQUS COMMAND用运行,语法如下:ABAQUS JOB=[JOB] USER¡[.FOR]¡C用户在编写用户子程序时,要注意以下几点:1.用户子程序不能嵌套。

即任何用户子程序都不能调用任何其他用户子程Home序,但可以调用用户自己编写的FORTRAN子程序和ABAQUS应用程序。

相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

【在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不知谁还有更好的方法。

相关文档
最新文档