共享ABAQUS UMAT和VUMAT子程序gydF4y2Ba

阿米尔:< / p >

我同意你的观点，在研究界有许多UMAT、UEL和其他用户子程序。在过去的几年里，我一直在提倡我的子程序，并向公众免费提供它们，可以从gydF4y2Bahttp://web.utk.edu/~ygao7/publication.htmgydF4y2Ba

为更广泛的公众编制一份这些清单是值得的。也许我们可以让ABAQUS为我们背书。关于你的三个问题:< / p >

(1)内聚区元素:我的网页中有一个，可以很容易地推广到其他版本，如camachoo - ortiz。< / p >

(2)多晶塑性。这个我不知道。我只知道在Jeff Kysar的网页(哥伦比亚大学)上可以找到由Jeff Kysar修改的著名代码。< / p >

(3) gurson型模型:我不知道是否有人免费提供。无限大应变版本并不难实现。< / p >

我同意。大约15年前，在我编写了第一个UEL之后，我想知道ABAQUS是否有兴趣为各种用户定义的子例程建立一个信息交换中心或存储库。这将类似于数学/科学应用软件(如方程求解器、特征值提取例程等)的可访问库(如NIST)。随着时间的推移，我撰写了UMAT和UEL程序，用于渐进式失效分析，高阶板有限元，专用粘接单元，以及包含各种耦合黏聚区模型的退聚有限元，以连接1-D梁，2-D和3-D连续单元。< / p >

这些子程序只是简单地放在我个人电脑的文件夹中收集电子“灰尘”。最好将它们上传到带有简短的自述文件或文档文件的站点，以供其他人使用。< / p >

埃里克< / p >

似乎有两个不同的线程有相同的主题，都有有趣的评论。如果版主能把这两者合并就好了。gydF4y2Bahttp://www.万博manbetx平台m.limpotrade.com/node/2375gydF4y2Ba

你好;< / p >

在我的毕业设计中，老师希望我在ABAQUS中编写UMAT子程序。我在找一个能帮我解决这个问题的人。我真的很想完成这个项目，真的很需要别人的帮助。< / p >

从现在开始，我要感谢那些帮助过我的人。< / p >

我留下我的msn地址给那些想联系我的人…< / p >

alper_tastan@hotmail.comgydF4y2Ba

我完全同意!而且，这个问题在ABAQUS上并不存在。同样的想法也可以帮助ANSYS, Nastran等!< / p >

你好,gydF4y2Ba共享Abaqus UMAT和VUMAT子程序。gydF4y2BaSIMULIA有一个“Abaqus过程自动化门户”gydF4y2Bahttp://www.simulia.com/PAPortal/gydF4y2Ba该门户为SIMULIA (Abaqus)社区提供了一种机制，通过讨论论坛进行GUI定制和脚本技术的协作。用户可以上传、下载和共享解决方案和过程自动化组件。gydF4y2Ba

Abaqus过程自动化门户即将迎来5周年纪念，许多工程师发现它是一个有用的资源。目前已有近3000名注册会员，而且这个数字还在稳步增长。gydF4y2Ba

问候,RenegydF4y2Ba

工欲善其事，必先利其器！< / p >

大家好，我有一个问题，如果我有UEL，我可以用它来模拟动态问题吗?< / p >

利用混凝土损伤塑性模型，用abaqus对钢筋混凝土圆柱壳进行建模。我想在abaqus cae中模拟圆柱形壳中的子午和周向钢筋。我想知道我应该使用哪种钢筋几何形状(角或常数)来建模子午和周向钢筋。gydF4y2Ba在abaqus中可以绘制钢筋混凝土外壳的裂缝图吗?gydF4y2Ba请指引我gydF4y2Ba

你可以在下面的帖子中找到一个很好的UMAT教程< / p >

http://万博manbetx平台m.limpotrade.com/node/7576gydF4y2Ba

大家好！< / p >

我是一名博士生，正在学习ABAQUS很短的时间，现在我正在使用UMAT子程序来定义我自己的材料，但是当我用。f扩展名完成fortran文件的编码时，我不知道如何用ABAQUS连接我的子程序文件和输入文件(.inp)。< / p >

当我写C:\Temp> abaqus job=test user=…< / p >

输入文件:…< / p >

Abaqus错误:无法定位以下文件:< / p >

我在abaqus教程中找不到任何信息。< / p >

谢谢< / p >

嗨，伙计们，我正试图为塑料分级材料制作一个UMAT代码，有人有相同主题的经验吗?< / p >

你好< / p >

我是学数学的新手。之前我用的是Windows xp，但现在我用的是Windows 7和64位操作系统< / p >

我安装了abaqus 6.10, Intel visual FORTRAN 11.1.60和vs2010，但我无法配置它们，有人可以帮助并给出在FORTRAN中配置abaqus的过程吗< / p >

非常感谢< / p >

你好,!< / p >

在我用vumat进行的模拟中，温度会在几个积分点上高于它的熔化温度，在这种情况下我应该如何更新应力和应变?给它们零价值?我认为压力应该存在，但是如何计算呢?< / p >

有人能帮我解决这个问题吗?< / p >

谢谢< / p >

大家好,< / p >

我是德国波鸿鲁尔大学计算机工程专业的一名硕士学生。我正在实现一个新的压电元件作为我的硕士论文的一部分。我应该写一个UEL子程序并在abaqus中实现它。我would like to know whether anyone has written subroutines in Fortran 90 or further versions and implemented them in abaqus. All the Help documents in Abaqus have subroutines written in F77. If it possible in F90, then it will be very helpful if someone could point out the necessary modifications that one has to make to the F90 file to run it in Abaqus. I think many ppl have this question .I myself have been trying to figure out an answer from past 1 month but could not find a convincing solution anywhere.< / p >

我是一名硕士研究生，我的论文题目要求我建模一个分裂的霍普金森压力棒。我am modelling a SHPB assembly together with the specimen, and I am having trouble with placing the loads,creating the boundary constrain between the incident, transmitter bar and the specimen. I created the interaction between the assembly with a suppressed TIE, and after loading and submitting the job I kept on getting errors. The load did not transmit to up to a quarter of the incident bar.
我需要关于如何在组成组件的单独组件之间应用附件的建议，如何正确地定时步骤和加载顺序。gydF4y2Ba
我使用了ABAQUS Explicit< / p >

巴蒂尼Veluri< / p >

你好埃里克,< / p >

你能分享一下你的知识吗，因为我正在为疲劳载荷下的内聚区模型构建UEL。所以如果我参考你的知识库，它会帮助我得到一些好的观点。< / p >

电子邮件:gydF4y2Babadrinath.veluri@hotmail.comgydF4y2Ba

提前感谢< / p >

问候< / p >

巴蒂尼< / p >

嗨朋友,< / p >

我正在准备力学领域的博士论文，我想在Abaqus上模拟粉末的致密化，但我不能模拟热脱粘行为(聚合物重量随温度的损失)。有人知道吗?< / p >

我有comsol和Abauqs< / p >

谢谢你的回答!< / p >

belman@hotmail.frgydF4y2Ba

贝尔< / p >

嗨，亲爱的朋友们gydF4y2Ba
经过大量的努力和时间和混乱VS，英特尔Fortran编译器和Abaqus(安装，版本，路径，…)，我终于可以运行一个分析与UMAT子程序成功。这次运行和其他一些运行都是从Abaqus文档中提取的。gydF4y2Ba
现在我想自己实现一个分析。因此，我阅读了文档，并通过在文档文件中进行自己的更改来编写UMAT文件。从算法的角度来看，该文件是正确的。但是我想知道如何将Fortran文件和Abaqus文件链接起来。我的意思是Abaqus如何实现Fortran文件是模型的兼容子程序?从逻辑上讲，子例程文件不应该与另一个abaqus文件兼容。(我不知道我是否能解释我的问题!)gydF4y2Ba

顺便说一句，我发现其中一个链接是关于属性模块。我应该在“编辑材料”中定义一个用户材料，我也应该在那里定义与PROPS兼容的变量(我现在混淆了!!)。< / p >

请引导我。< / p >

嗨朋友,< / p >

我需要umat subrotine(Lemaitre损伤模型)< / p >

如果你有，请叫我来。谢谢你。< / p >

我的电子邮件:gydF4y2Baepeyman_hf@yahoo.comgydF4y2Ba

随着时间的推移，我撰写了UMAT和UEL程序，用于渐进式失效分析，高阶板有限元，专用粘接单元，以及包含各种耦合黏聚区模型的退聚有限元，以连接1-D梁，2-D和3-D连续单元。gydF4y2Ba肾结石gydF4y2Ba

嗨< / p >

我已经编译了一个fortran代码，并将其用作子程序，但我已经收到了这个错误:< / p >

工作Base_subro错误:链接期间的问题- Abaqus/标准用户子例程。此错误可能是由于Abaqus用户子例程参数不匹配造成的。这些参数有时会随着版本的不同而变化，因此可能需要调整与前一个Abaqus版本一起使用的用户子例程。< / p >

我该怎么办呢?< / p >

谢谢< / p >

我以前遇到过这个问题。也许在你的Fortran文件中有一些错误，你可以检查一下，特别是程序语法。< / p >

好运！< / p >

嗨，有人用Vumat subrotine做过各向异性硬化吗?如果你有任何各向异性硬化的示例代码请告诉我。< / p >

你好朋友< / p >

几天前，我写了一个各种杨氏模量的连续梁的子程序。对于二维波束是正确的，但是对于三维波束写入误差是正确的。< / p >

错误:对这个增量进行了太多尝试< / p >

Abaqus/Standard Analysis退出时出现错误-如果文件存在，请查看消息文件以获取可能的错误消息。< / p >

我经常改变负载和步长增量，但它还不是真的。我不知道哪里是不正确的。< / p >

请帮帮我.............< / p >

谢谢。< / p >

你好,< / p >

我最近购买了Intel Visual Fortran 11.1.072和Visual Studio 2010。我目前使用的是ABAQUS 6.10ef2。按照以下顺序安装后，< / p >

1.Microsoft Visual c++< / p >

2.英特尔Fortran< / p >

3.有限元分析< / p >

我的验证显示visual studio和FORTRAN的状态都是失败的。在这方面有人能帮我吗?< / p >

Mesbah< / p >

萨拉姆gydF4y2Ba
我需要gydF4y2Ba模型gydF4y2Ba的gydF4y2Ba股骨gydF4y2Ba有限元分析gydF4y2Ba以及如何gydF4y2Ba静载荷gydF4y2Ba在gydF4y2Ba皮质骨gydF4y2Ba
什么gydF4y2Ba是两者的关系gydF4y2Ba有限元分析gydF4y2Ba和gydF4y2BaumatgydF4y2Ba
谢谢你！gydF4y2Ba
我的gydF4y2Ba电子邮件:gydF4y2Basalah_cao@yahoo.frgydF4y2Ba

嗨朋友,gydF4y2Ba

我需要它。gydF4y2Ba和VmatgydF4y2Ba亚酪氨酸(GTN型)gydF4y2Ba

如果你有，请发给我。谢谢你。< / p >

我的电子邮件:gydF4y2BaggydF4y2Bauojunhang@163.comgydF4y2Ba

骨折;塑性< / p >

请gydF4y2Ba我需要gydF4y2Ba如何配置gydF4y2Ba有限元分析gydF4y2Ba与gydF4y2BafortrangydF4y2Ba和我gydF4y2Ba安装gydF4y2Ba有限元分析gydF4y2Ba6.11gydF4y2Ba和gydF4y2BaFORTRANgydF4y2Ba10gydF4y2Ba和gydF4y2Ba微软gydF4y2BavisuagydF4y2BalgydF4y2Ba10gydF4y2Ba程序gydF4y2Ba当gydF4y2Ba下运行gydF4y2Ba有限元分析gydF4y2Ba有一个gydF4y2Ba错误gydF4y2Ba：gydF4y2Ba问题gydF4y2Ba落在gydF4y2Ba编译gydF4y2Ba-gydF4y2Baifort.exegydF4y2Ba没有找到gydF4y2Ba在gydF4y2Ba路径gydF4y2Ba．gydF4y2Ba
将如何gydF4y2Ba是gydF4y2Ba请gydF4y2Ba
谢谢你！gydF4y2Ba

喂gydF4y2Ba
我需要使用umat Abaqus来模拟复合材料板的热力学效应gydF4y2Ba
请帮帮我< / p >

嗨……< / p >

如果方程是这样的< / p >

σ(n+1) = σ(n) + (σ)< / p >

如果我得到初始应力σ (n)和差应力δ (σ)，首先我需要调用(KEY)。Eq . 11)然后....< / p >

或者我需要做出完整的刚度矩阵，然后< / p >

I = 1吗< / p >

应力(i) = .......< / p >

非常感谢您的帮助。< / p >

为什么每个人都在做一个FGM板的主动振动控制项目，并使用ansys，但我有一个问题，如何在ansys中做一个三维FGM板< / p >

所以我需要一个apdl代码< / p >

我的电子邮件:gydF4y2Bakouider84@live.comgydF4y2Ba

有什么帮助吗，谢谢< / p >

你好……< / p >

我正在开始一个学术项目，我的目标是用abaqus模拟早期混凝土的开裂。< / p >

知道我刚刚开始，这是一个经典问题，我需要知道材料行为是否已经存在于Abaqus中?如果没有，有人知道如何创建相应的UMAT吗?< / p >

谢谢……gydF4y2Ba

为什么使用核级石墨作为慢化剂比其他慢化剂更好，为什么比其他慢化剂更经济。< / p >

大家好。我叫亚历克斯，是萨拉曼卡大学的博士生。我用的是轴对称模型，钢材料。我打算获得在很短的时间内(皮秒)击中主要负载模型(Tera帕斯卡)时产生的残余应力。本构模型是有效的:弹性各向同性，塑料约翰逊-库克，eos -up。首先使用Abaqus Explicit Abaqus，然后使用standard(我导入张力)，但Abaqus CAE / standard(版本6.13.2)不允许塑料或Johnson-cook或EOS。使弹簧回弹，从而产生残余应力。< / p >

请问有人能帮帮忙吗?< / p >

谢谢< / p >

的问候。< / p >

你好,< / p >

我想学习关于Abacus(分离式霍普金森压杆压缩)的模拟程序的步骤，该程序在温度和不影响温度的情况下进行，该装置在两根长杆之间放置试样，然后用短杆撞击，从而产生高应变率的影响。我们的数据收集系统由一对沿每个杆的中间应变计组成。记录数据并以初始脉冲w表示gydF4y2BaAve通过入射杆应变计，当波到达杆端并反射回入射杆应变计时产生的反射脉冲，最后是透射脉冲，即波通过样品后剩下的脉冲。在所有情况下，一旦波通过，应变片的电压读数应归零。然而，在新的条形图中，信号在波通过一个刻度后仍然存在。gydF4y2Ba
我的6.12版本gydF4y2Ba
任何人都可以帮助我gydF4y2Ba
问候gydF4y2Ba

< / p > < p > HI gydF4y2Ba

我是abaqus的新手，我正在为新hookean材料编写VUMAT。早些时候我做了UMAT，它运行成功了 

gydF4y2Ba

但是当运行VUMAT时，它给出错误“无法打开包含文件”vaba_param。公司”。我正在正确地发出命令 

gydF4y2Ba

还有其他可能吗?输入文件错误会导致这种错误吗?< / p >< / p >

在这里加上我的两点意见:< / p >

http://万博manbetx平台m.limpotrade.com/node/17545gydF4y2Ba

谁能上传非线性粘弹性材料的Schapery材料模型(1-D或3-D)的UMAT ?< / p >

喂,< / p >

谁能解释一下ABAQUS是如何计算弹性-塑性问题的PEEQ、PENER和SENER的?其次，我想知道ABAQUS是如何计算这种问题的等效应变的。这将是非常有帮助的，如果有人帮助我或给我一些链接。< / p >

感谢你< / p >

Braj< / p >

亲爱的所有,< / p >

点头,< / p >

我将为我的模型使用两个不同的umat。第一个定义了超弹性，第二个定义了我的模型中的骨重塑过程。我的问题是，如何将它们结合起来同时使用?< / p >

在输入文件中，您应该有两个具有不同名称的不同用户材料。变量cmname保存材料名称。< / p >

在fortran子程序文件中，您的代码应该大致如下:< / p >

子程序umat ()< / p >

变量声明< / p >

如果cmname .eq.“matl1。然后< / p >

——本构方程< / p >

如果cmname .eq.matl2然后< / p >

——本构方程< / p >

结束< / p >

我不记得具体的语法了。请阅读abaqus手册。< / p >

祝你好运< / p >

Brunda< / p >

大家好！< / p >

我想实现陶瓷在弹道冲击下的多尺度模型。我的参考论文:脆性固体的动态破碎:一个多尺度模型，Christophe Denoual, franois Hild。我学习了达索系统的“在ABAQUS中编写用户子程序”。我也读了一些例子。但是我的模型有一个新的参数，和例子模型不一样。我想用Abaqus/显式分析。我f any one can help me regarding this it would be great. Please read the paper and give me some idea of how to implement in Abaqus using Subroutine VUMAT.If you have the subroutine it would be great help if you please send the subroutine to my email Id:kosuribs99@gmail.comgydF4y2Ba

亲爱的所有人gydF4y2Ba
我想改变一个现有的材料模型，使它的一些参数可以作为应变率的函数而改变。例如，我想在abaqus的简单塑性材料模型中替换流动应力。应该注意的是，我有一个应变率函数，它与应变率有不正常的关系。< / p >

使用选项PLASTIC, DEPENDENCIES应该是可行的< / p >

否则:< / p >

ABAQUS邮件列表gydF4y2Ba
https://groups.yahoo.com/group/ABAQUSgydF4y2Ba

你好,< / p >

我为“脆性材料动态破碎的概率损伤模型”编写了一个子程序，《应用力学进展》，2010年。我在UMAT子程序中所遵循的程序如下:< / p >

1.我calculate the principal stresses.
2.根据2010年论文式(84)随机生成一个破坏应力。gydF4y2Ba
3.对于每个主应力，我计算修正缺陷生长密度，如2010年论文式(83)所示。gydF4y2Ba
4.每个主应力对应的损伤值根据2010年论文式(81)给出的条件求出(注:D1、D2、D3三个正交方向均采用各向异性损伤模型)gydF4y2Ba
5.我update the damage tensor/compliance tensor using (33)(JPEG file is attached)as provided in multi-scale model paper(2002).
6.我update the stress and repeat the steps:1-6.< / p >

注意:1。材料是各向同性弹性只需要E、泊松比、材料参数。给出了其他材料参数，如威布尔参数。参考所提供的文件进行澄清。< / p >

2.我am dealing with Impact problems so I need VUMAT< / p >

如果我在每一步计算柔度张量(如附件所示)并将其逆以获得本构张量，然后与应变相乘以获得应力增量，我对VUMAT有疑问。我不会更新雅可比矩阵。Can I do this in VUMAT to calculate stressNew?.Is my subroutine right?.Please check the VUMAT.txt file.Am I updating the stress correctly?< / p >

子例程vumat (gydF4y2Ba

1 . nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal，gydF4y2Ba

2 stepTime, totalTime, dt, cmname, coordMp, charLength，gydF4y2Ba

3 . props, density, strain, relspin，gydF4y2Ba

旧的，伸展的，退化的，旧的，field旧的，gydF4y2Ba

3 . statold, enerinterold, enerInelasOld，gydF4y2Ba

tempNew, stretchNew, defgradNew, fieldNew，gydF4y2Ba

5 . statennew, enerinternew, enerInelasNewgydF4y2Ba

CgydF4y2Ba

包括“VABA_PARAM。公司的gydF4y2Ba

CgydF4y2Ba

CgydF4y2Ba

此算法中不使用所有以(*)为维度的数组gydF4y2Ba

尺寸道具(nprops)，密度(nblock)，gydF4y2Ba

1 coordMp (nblock *),gydF4y2Ba

2 charLength(*)， strain (nblock,ndir+nshr)，gydF4y2Ba

3 relSpinInc(*)， tempOld(*)，gydF4y2Ba

4 stretchOld(*)， defgradOld(*)，gydF4y2Ba

5 . fieldOld(*)， stressOld(nblock,ndir+nshr)，gydF4y2Ba

6 statold (nblock,nstatev)， enerinterold (nblock)，gydF4y2Ba

7 . enerInelasOld(nblock)， tempNew(*)，gydF4y2Ba

8 stretchNew(*)， defgradNew(*)， fieldNew(*)，gydF4y2Ba

9 . statennew (nblock,ndir+nshr)， statennew (nblock,nstatev)，gydF4y2Ba

1 enerinternew (nblock)， enerInelasNew(nblock)gydF4y2Ba

CgydF4y2Ba

字符* 80 cmnamegydF4y2Ba

CgydF4y2Ba

CgydF4y2Ba

***********************************************************************************gydF4y2Ba

** vumat, abaqus / standard中包含**的概率损伤模型gydF4y2Ba

**弹性性能和威布尔参数**gydF4y2Ba

* * * *gydF4y2Ba

***********************************************************************************gydF4y2Ba

***********************************************************************************gydF4y2Ba

**gydF4y2Ba

**gydF4y2Ba

**gydF4y2Ba

*用户子程序gydF4y2Ba

参数(M = 3, N = 3, 0 = 0. d0,一个= 1. d0, 2 = 2. d0, 3 = 3. d0,gydF4y2Ba

+ 6 = 6。D0, 9 = 9. D0, XALPHA = 0.31, = 4.18)gydF4y2Ba

CgydF4y2Ba

维XIDEN (M, N),我(6,6),x (6,6),gydF4y2Ba

+ STR (M, N), XDAM (nblock, M), AJAC (6,6), STRN (M, N),gydF4y2Ba

+ XPRIN (nblock, M), XDEF (nblock, M), XPRINN (nblock, M),gydF4y2Ba

+ XFAIL (nblock), XPRINR (nblock, M)gydF4y2Ba

维度eigVal(nblock,M)， eigVec(nblock,M)gydF4y2Ba

真正的V, B,, RKO RKT, RKTH, RKF RKgydF4y2Ba

整数的种子gydF4y2Ba

种子= 7654321gydF4y2Ba

CgydF4y2Ba

C xdam =各主方向对应的损伤值，xprint =主应力gydF4y2Ba

C xdef =各主应力对应的修正生长缺陷密度gydF4y2Ba

xprincipal - principal stress new, xprincipal - principal stress oldgydF4y2Ba

C -主应力速率gydF4y2Ba

CgydF4y2Ba

C --------------------------------------------------------------------gydF4y2Ba

CgydF4y2Ba

C指定材料属性gydF4y2Ba

所有尺寸均以米为单位gydF4y2Ba

CgydF4y2Ba

E =道具(1)gydF4y2Ba

XNUE =道具(2)gydF4y2Ba

XDENSITY =道具(3)gydF4y2Ba

XPOROSITY =道具(4)gydF4y2Ba

XM =道具(5)gydF4y2Ba

XMEANFS =道具(6)gydF4y2Ba

XEFFVOL =道具(7)gydF4y2Ba

ALAMBDA = E * XNUE / (1 + XNUE) /(两人* XNUE)gydF4y2Ba

阿姆河= E / (1 + XNUE) / 2gydF4y2Ba

V =√E / XDENSITYgydF4y2Ba

我= 1,MgydF4y2Ba

J = 1, NgydF4y2Ba

STR (I, J) = 0gydF4y2Ba

最后做gydF4y2Ba

最后做gydF4y2Ba

nten = ndir + nshrgydF4y2Ba

nten, I = 1吗gydF4y2Ba

nten, J = 1吗gydF4y2Ba

AJAC (I, J) = 0gydF4y2Ba

ENDDOgydF4y2Ba

ENDDOgydF4y2Ba

C状态变量:1 ~ 3表示损伤值在0 ~ 1之间变化gydF4y2Ba

C状态变量:4-失效应力gydF4y2Ba

C状态变量:5、6、7主应力gydF4y2Ba

C状态变量:8、9、10缺陷密度对应每个主应力gydF4y2Ba

if (stepTime.eq.zero) thengydF4y2Ba

101 I = 1,nblock吗gydF4y2Ba

B = ((XEFFVOL) * ((XMEANFS) * * XM))gydF4y2Ba

PK =跑(种子)gydF4y2Ba

Q = 1 / (ONE-PK)gydF4y2Ba

XFAIL (i) = (XMEANFS) *((日志(Q)) * * (1 / XM))gydF4y2Ba

AJAC (1, - 1) = (ALAMBDA + 2 * 12)gydF4y2Ba

AJAC (2, 2) = (ALAMBDA + 2 * 12)gydF4y2Ba

AJAC (3) = (ALAMBDA + 2 * 12)gydF4y2Ba

AJAC (4, 4) = 12gydF4y2Ba

AJAC (5,5) = 12gydF4y2Ba

AJAC (6,6) = 12gydF4y2Ba

AJAC(1、2)= ALAMBDAgydF4y2Ba

AJAC(1、3)= ALAMBDAgydF4y2Ba

AJAC (2, 1) = ALAMBDAgydF4y2Ba

AJAC (2,3) = ALAMBDAgydF4y2Ba

AJAC (3,1) = ALAMBDAgydF4y2Ba

AJAC (3,2) = ALAMBDAgydF4y2Ba

print*，“失败压力”gydF4y2Ba

打印*,XFAIL(我)gydF4y2Ba

nten, j = 1吗gydF4y2Ba

stressNew (i, j) = stressOld (i, j) +gydF4y2Ba

2 (DOT_PRODUCT(AJAC(j，:)， strain (i，:))))gydF4y2Ba

ENDDOgydF4y2Ba

stateNew (4) = XFAIL(我)gydF4y2Ba

种子=种子+ 1gydF4y2Ba

101年继续gydF4y2Ba

endifgydF4y2Ba

100 I = 1,nblock吗gydF4y2Ba

CgydF4y2Ba

C将应力写成矩阵形式gydF4y2Ba

CgydF4y2Ba

j = 1,3gydF4y2Ba

STR(j,j) = stressOld(i,j)gydF4y2Ba

最后做gydF4y2Ba

STR(1,2) = stressOld(1, 4)gydF4y2Ba

STR(2,1) = stressOld(1,4)gydF4y2Ba

STR(1,3) = stressOld(1,5)gydF4y2Ba

STR(3,1) = stressOld(1,5)gydF4y2Ba

STR(2,3) = stressOld(1,6)gydF4y2Ba

STR(3,2) = stressOld(1,6)gydF4y2Ba

求应力矩阵的特征值gydF4y2Ba

调用vspring (NBLOCK, STR, eigVal, ndir, nshr)gydF4y2Ba

C是主重音gydF4y2Ba

3, j = 1吗gydF4y2Ba

XPRIN (i, j) = eigVal (i, j)gydF4y2Ba

ENDDOgydF4y2Ba

Print *， 'old principal stress'gydF4y2Ba

打印*,XPRIN(我,:)gydF4y2Ba

CgydF4y2Ba

C将应力写成矩阵形式gydF4y2Ba

CgydF4y2Ba

j = 1,3gydF4y2Ba

STRN(j,j) = stressNew(i,j)gydF4y2Ba

最后做gydF4y2Ba

STRN(1,2) = stressNew(1, 4)gydF4y2Ba

STRN(2,1) = stressNew(1,4)gydF4y2Ba

STRN(1,3) = stressNew(1, 5)gydF4y2Ba

STRN(3,1) = stressNew(1,5)gydF4y2Ba

STRN(2,3) = stressNew(1,6)gydF4y2Ba

STRN(3,2) = stressNew(1,6)gydF4y2Ba

求应力新矩阵的特征值gydF4y2Ba

调用vspring (NBLOCK, STRN, eigVal, ndir, nshr)gydF4y2Ba

C是主重音gydF4y2Ba

3, j = 1吗gydF4y2Ba

XPRINN (i, j) = eigVal (i, j)gydF4y2Ba

ENDDOgydF4y2Ba

Print *， 'new principal stress'gydF4y2Ba

打印*,XPRINN(我,:)gydF4y2Ba

C求主应力速率gydF4y2Ba

3, j = 1吗gydF4y2Ba

XPRINR (i, j) = (XPRINN (i, j) -XPRIN (i, j)) / dtgydF4y2Ba

ENDDOgydF4y2Ba

打印*，“主应力率”gydF4y2Ba

打印*,XPRINR(我,:)gydF4y2Ba

C求各主应力对应的缺陷密度gydF4y2Ba

C求各主应力对应的损伤值gydF4y2Ba

C参考方程83gydF4y2Ba

3, j = 1吗gydF4y2Ba

IF (XPRINR(i,j).LE.ZERO) THENgydF4y2Ba

stateNew (i, j) = stateOld (i, j)gydF4y2Ba

其他的gydF4y2Ba

如果((XPRIN (i, j) .GT.ZERO)。。(XPRINR (i, j) .GT.ZERO))gydF4y2Ba

IF (xprint (i,j).LE.XFAIL(i)) THENgydF4y2Ba

stateNew (i, j + 7) = 0gydF4y2Ba

stateNew (i, j) = 0gydF4y2Ba

其他的gydF4y2Ba

stateNew (i, j + 7) = K * ((XPRIN (i, j)) * * XM)gydF4y2Ba

cx是方程81右边的值gydF4y2Ba

XBETA = 6 * * ((0.33 * V) * * 3) * stateNew (i, j + 7)gydF4y2Ba

(RKO = ((1-XDAM (i, j)) * (XBETA * (totalTime * * 2))) / 2gydF4y2Ba

RKT = ((1 - (XDAM (i, j) + ((dt * (RKO) / 2)))gydF4y2Ba

1 * (XBETA * ((totalTime + dt / 2) * * 2)) / 2)gydF4y2Ba

RKTH = (1 - (XDAM (i, j) + (dt * RKT) / 2))gydF4y2Ba

1 * (XBETA * ((totalTime + dt / 2) * * 2)) / 2gydF4y2Ba

RKF = (1 - (XDAM (i, j) + (dt * RKTH)))gydF4y2Ba

1 * (XBETA * ((totalTime + dt) * * 2)) / 2gydF4y2Ba

RK = + 2 * RKT RKO + 2 * RKTH + RKFgydF4y2Ba

stateNew (i, j) = stateOld (i, j) + ((dt / 6) * RK)gydF4y2Ba

ENDIFgydF4y2Ba

其他的gydF4y2Ba

stateNew (i, j + 7) = 0gydF4y2Ba

stateNew (i, j) = 0gydF4y2Ba

ENDIFgydF4y2Ba

ENDIFgydF4y2Ba

ENDDOgydF4y2Ba

打印*，“损坏值”gydF4y2Ba

打印*,stateNew(我,1)gydF4y2Ba

打印*,stateNew(我,2)gydF4y2Ba

打印*,stateNew(我,3)gydF4y2Ba

C状态变量:1 ~ 3表示损伤值在0 ~ 1之间变化gydF4y2Ba

C状态变量:4-失效应力gydF4y2Ba

stateNew (4) = XFAIL(我)gydF4y2Ba

stateNew(我,5)= XPRIN(我,1)gydF4y2Ba

stateNew (6) = XPRIN(我,2)gydF4y2Ba

stateNew (7) = XPRIN(我,3)gydF4y2Ba

C更新应力矩阵，更新柔度张量gydF4y2Ba

C矩阵和逆顺应张量得到本构张量gydF4y2Ba

cxc =柔度张量，xs =柔度张量的逆gydF4y2Ba

XC (1, - 1) = 1 / E / (1-stateNew(我,1))gydF4y2Ba

XC(1、2)= (-XNUE) / EgydF4y2Ba

XC(1、3)= (-XNUE) / EgydF4y2Ba

XC(1、4)= 0.0gydF4y2Ba

XC(1、5)= 0.0gydF4y2Ba

XC (1,6) = 0.0gydF4y2Ba

XC (2, 1) = (-XNUE) / EgydF4y2Ba

XC (2, 2) = 1 / E / (1-stateNew(我,2))gydF4y2Ba

XC(2、3)= (-XNUE) / EgydF4y2Ba

XC(2、4)= 0.0gydF4y2Ba

XC(2、5)= 0.0gydF4y2Ba

XC (2,6) = 0.0gydF4y2Ba

XC (3,1) = (-XNUE) / EgydF4y2Ba

XC (3 2) = (-XNUE) / EgydF4y2Ba

XC (3) = 1 / E / (1-stateNew(我,3))gydF4y2Ba

XC(3、4)= 0.0gydF4y2Ba

XC(3、5)= 0.0gydF4y2Ba

XC(3、6)= 0.0gydF4y2Ba

我(4,1)= 0.0gydF4y2Ba

我(4,2)= 0.0gydF4y2Ba

我(4,3)= 0.0gydF4y2Ba

我(4,4)= (1 + XNUE) / ((ONE-stateNew(我,2))* * (XALPHA))gydF4y2Ba

3 / ((ONE-stateNew(我,3))* * (XALPHA)) / EgydF4y2Ba

我(4、5)= 0.0gydF4y2Ba

我(4、6)= 0.0gydF4y2Ba

我(5、1)= 0.0gydF4y2Ba

XC (5, 2) = 0.0gydF4y2Ba

我(5,3)= 0.0gydF4y2Ba

XC (5,4) = 0.0gydF4y2Ba

XC (5,5) = (1 + XNUE) / ((ONE-stateNew(我,3))* * (XALPHA))gydF4y2Ba

4 / ((ONE-stateNew(我,1))* * (XALPHA)) / EgydF4y2Ba

我(5、6)= 0.0gydF4y2Ba

我(6,1)= 0.0gydF4y2Ba

XC (2) = 0.0gydF4y2Ba

我(6,3)= 0.0gydF4y2Ba

XC (4) = 0.0gydF4y2Ba

我(6,5)= 0.0gydF4y2Ba

XC (6,6) = (1 + XNUE) / ((ONE-stateNew(我,1))* * XALPHA)gydF4y2Ba

5 / ((ONE-stateNew(我,2))* * XALPHA) / EgydF4y2Ba

用高斯乔丹法求x =柔度张量的逆gydF4y2Ba

调用INVERSEA (XC XS)gydF4y2Ba

打印*，'顺应张量的逆'gydF4y2Ba

打印*,XS (:,)gydF4y2Ba

nten, j = 1吗gydF4y2Ba

stressNew (i, j) = stressOld (i, j) + DOT_PRODUCT (XS (j,:), strainInc(我:))gydF4y2Ba

ENDDOgydF4y2Ba

100年继续gydF4y2Ba

返回gydF4y2Ba

结束gydF4y2Ba

******************** 年底VUMAT子例程 ******************************gydF4y2Ba

*************************************************************************gydF4y2Ba

用高斯约当法求6x6矩阵的C **逆gydF4y2Ba

*************************************************************************gydF4y2Ba

用高斯约当法求6x6矩阵的C **逆gydF4y2Ba

子例程INVERSEA (A, C)gydF4y2Ba

CgydF4y2Ba

包括“VABA_PARAM。公司的gydF4y2Ba

CgydF4y2Ba

参数(M = 6, N = 6)gydF4y2Ba

维度(M, N)、B(6、12),C (M, N)gydF4y2Ba

整我,J, NgydF4y2Ba

C **构建增广矩阵bgydF4y2Ba

40j =1 n吗gydF4y2Ba

30 i =1 n吗gydF4y2Ba

B (I, J) = (I, J)gydF4y2Ba

30日继续gydF4y2Ba

40继续gydF4y2Ba

我= 1,NgydF4y2Ba

J = N + 1, 2 N *gydF4y2Ba

如果((j-i) .eq.n)则gydF4y2Ba

B (I, J) = 1.0gydF4y2Ba

其他的gydF4y2Ba

B (I, J) = 0.0gydF4y2Ba

ENDIFgydF4y2Ba

ENDDOgydF4y2Ba

ENDDOgydF4y2Ba

160 i =1 n吗gydF4y2Ba

(我:)= B(我,:)/ B(我)gydF4y2Ba

150j =1 n吗gydF4y2Ba

如果(i . j .)那么gydF4y2Ba

B: B (J) = (J:) B - B (J,我)*(我,:)gydF4y2Ba

ENDIFgydF4y2Ba

150年继续gydF4y2Ba

160年继续gydF4y2Ba

我= 1,NgydF4y2Ba

J = N + 1, 2 N *gydF4y2Ba

C(我J-6) = B (I, J)gydF4y2Ba

ENDDOgydF4y2Ba

ENDDOgydF4y2Ba

打印*，' c矩阵第一行'gydF4y2Ba

:打印*,C (1)gydF4y2Ba

返回gydF4y2Ba

结束gydF4y2Ba

有谁能告诉我如何在Ababqus CAE上模拟试样在受拉伸载荷和试样中间45度处断裂或剪切的情况吗?如果你有模拟的程序，请通过我的邮箱发给我< / p >

bogale.sintayehu@gmail.comgydF4y2Ba

非常感谢!!< / p >

ABAQUS邮件列表gydF4y2Ba
https://groups.yahoo.com/group/ABAQUSgydF4y2Ba

大家好! !< / p >

我是abaqus和fortran的新手。谁能告诉我，如果我可以使用USDFLD子程序来改变摩擦角相对于塑性应变根据Vermeer 1984< / p >

我是博士，也是UMAT的新用户。有没有人可以帮助我调试我的UMAT建模< / p >

聚合物的记忆效果如何?有人有很好的经验吗?请在这方面给我很大的帮助。提前感谢。我的邮箱是gydF4y2Baliqinan1986@gmail.comgydF4y2Ba

Hy,< / p >

我是博士，也是子程序和abaqus的新用户。我希望有人能帮助我，我有一个流体动力学问题，我模拟密封(弹性体-超弹性材料)与流体接触，由雷诺兹方程控制(我在子程序中求解)。密封变形大，停止模拟不收敛。< / p >

你知道如何解决这个问题，如何优化计算时间吗?< / p >

提前感谢< / p >

我的邮箱:gydF4y2Baabdelhak.azzi@univ-poitiers.frgydF4y2Ba

你好< / p >

我想在VUMAT中建模FG材质。为此，我需要每个质点的x和y分量。但是Abaqus只给我一个x和y。而每个nblock由多个材质点组成。< / p >

请引导我。< / p >

谢谢,< / p >

Reza Pourshab< / p >

订阅并寻求帮助gydF4y2Ba
ABAQUS邮件列表gydF4y2Ba
https://groups.yahoo.com/group/ABAQUSgydF4y2Ba

在发布前搜索存档。< / p >

有关子例程的介绍，请获取该文件gydF4y2Ba
http://万博manbetx平台m.limpotrade.com/files/WritinggydF4y2Ba用户子程序与ABAQUS.pdf< / p >

祝你好运< / p >

弗兰克< / p >

大家好!< / p >

我想用UMAT模拟非饱和土的裂缝形成，有人能帮我写代码吗?< / p >

谢谢< / p >

即使忽略初始价格，也需要一个Intel XE编译器。编译器需要Visual studio。还有英特尔GPU。< / p >

大家好，我已经实现了一个UMAT子程序，并为一个塑料增量调试了它。有一个问题是，当一个增量结束时，应力和所有状态变量都保存为一个增量，在下一个增量开始时，我希望我有相同的应力和状态变量，而事实并非如此(最奇怪的是，在一个塑性增量之后，我的状态变量被更新(不为零)，对于下一个增量，它将所有这些设置为等于零)。谁能告诉我是什么问题，我应该如何解决这个问题?问候,Faezeh< / p >