共享ABAQUS UMAT和VUMAT子例程gydF4y2Ba

阿米尔:< / p >

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

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

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

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

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

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

这些子程序只是简单地坐在文件夹在我的电脑收集电子“灰尘”。最好将它们上传到带有简短README或文档文件的站点，供其他人使用。< / 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子例程来定义我自己的材料，但当我完成编码fortran文件与.f扩展名，我不知道如何连接我的子例程文件和输入文件(.inp)与ABAQUS。< / p >

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

输入文件:…< / p >

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

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

谢谢< / p >

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

你好< / p >

我是abaqus umat的新手。以前我用的是Windows xp，现在我用的是Windows 7和64位操作系统< / p >

我安装了abaqus 6.10，英特尔visual FORTRAN 11.1.60和vs2010，但我不能配置它们，谁能帮助和给出[过程配置abaqus在FORTRAN< / 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, Intel Fortran编译器和Abaqus的困惑(安装，版本，路径，…)，我终于可以成功地运行UMAT子例程的分析。这次运行和其他一些运行都是从Abaqus文档中提取的。gydF4y2Ba
现在我要实现我自己的分析。所以我阅读了文档，并通过自己修改文档文件来编写UMAT文件。从算法的角度来看，该文件是正确的。但我想知道我如何可以链接的Fortran文件和Abaqus文件。我的意思是Abaqus如何实现Fortran文件是模型的兼容子例程?从逻辑上讲，子例程文件不应该与另一个abaqus文件兼容。(我不知道我是否能解释我的问题!)gydF4y2Ba

顺便说一下，我已经弄清楚其中一个链接是关于属性模块。我应该在“编辑材料”中定义一个USER MATERIAL，我也应该在那里定义与PROPS兼容的变量(我现在搞混了!!)< / p >

请引导我。< / p >

嗨朋友,< / p >

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

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

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

随着时间的推移，我已经编写了UMAT和UEL例程，用于渐进破坏分析，高阶板有限元，专用粘接单元，以及将各种耦合内聚区模型结合到连接1-D梁，以及2-D和3-D连续单元的decohesion有限元。gydF4y2Ba肾结石gydF4y2Ba

嗨< / p >

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

作业Base_subro中的错误:在连接Abaqus/标准用户子例程时出现问题。此错误可能是由于Abaqus用户子例程参数不匹配造成的。这些参数有时会在不同版本之间发生变化，因此可能需要调整与Abaqus以前版本一起使用的用户子例程。< / p >

我该怎么办?< / p >

谢谢< / p >

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

好运！< / p >

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

你好朋友< / p >

几天前，我写了一个具有不同杨氏模量的悬臂梁的子程序。对于2D光束是正确的，但对于3D光束写错误。< / p >

错误:为这个增量做了太多的尝试< / p >

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

我经常改变负载和步进增量，但这还不是真的。我不知道在哪里。< / p >

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

谢谢。< / p >

你好,< / p >

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

1.微软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

我需要umatgydF4y2Ba和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)和微分应力δ (σ)， ist我需要调用(KEY。Eq . 11)然后....< / p >

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

I = 1 ntens吗< / p >

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

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

Hy每个人都在工作关于一个项目的主动振动控制的FGM板和我使用ansys，但我有一个问题，如何做三维FGM板在ansys< / p >

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

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

任何帮助，谢谢< / p >

你好……< / p >

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

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

谢谢……gydF4y2Ba

为什么用核级石墨作为缓蚀剂比其他缓蚀剂好，为什么它比其他缓蚀剂经济。< / p >

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

有人能帮帮我吗?< / p >

谢谢< / p >

的问候。< / p >

你好,< / p >

我想学习关于Abacus主题的模拟程序的步骤(分裂霍普金森压杆压缩)，有和没有影响温度。该装置在两根长杆之间放置试样，然后用短杆撞击时，产生高应变率冲击。我们的数据收集系统包括一对应变计在每个棒中间。数据被记录并以初始脉冲w表示gydF4y2BaAve通过入射棒应变计，当波击中棒的末端时产生的反射脉冲反射回入射棒应变计，最后是一个透射脉冲，这是波通过样品后的剩余部分。在所有情况下，一旦波通过，应变计上的电压读数应归零。然而，有了新的条，信号在波通过一个测量器后仍然存在。gydF4y2Ba
我的版本是6.12gydF4y2Ba
任何人都可以帮助我gydF4y2Ba
问候gydF4y2Ba

< / p > < p > HI gydF4y2Ba

我是abaqus的新手，我正在为新胡克材料写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, François Hild。我研究了达索系统的“在ABAQUS中编写用户子程序”。我也读了例子。但是我的模型有一个新的参数，它与示例模型不相似。我想做Abaqus/Explicit分析。我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中有一个疑问，如果我在每一步都计算遵从性张量(如附件所示)，并将其逆得到本构张量，然后与strainInc相乘得到应力增量。我不会更新雅可比矩阵。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个道具，密度，strainInc, relSpinInc，gydF4y2Ba

4 tempOld, stretchOld, defgradOld, fieldOld，gydF4y2Ba

3 stress sold, stateOld, enerInternOld, enerInelasOld，gydF4y2Ba

6 tempNew, stretchNew, defgradNew, fieldNew，gydF4y2Ba

5 . stressNew, statennew, enerInternNew, enerInelasNew)gydF4y2Ba

CgydF4y2Ba

包括“VABA_PARAM。公司的gydF4y2Ba

CgydF4y2Ba

CgydF4y2Ba

C在此算法中不使用所有以(*)表示维度的数组gydF4y2Ba

维度道具(nprops)，密度(nblock)，gydF4y2Ba

1 coordMp (nblock *),gydF4y2Ba

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

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

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

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

6 stateOld(nblock,nstatev)， enerInternOld(nblock)，gydF4y2Ba

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

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

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

1 enerInternNew(nblock)， enerInelasNew(nblock)gydF4y2Ba

CgydF4y2Ba

字符* 80 cmnamegydF4y2Ba

CgydF4y2Ba

CgydF4y2Ba

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

** vumat，概率损伤模型的abaqus /标准纳入**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 =各主方向对应的损伤值，xprin =主应力gydF4y2Ba

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

C xprinn-principalstress是新的，xprinn-principalstress是旧的gydF4y2Ba

C xprinr -主应力速率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,10gydF4y2Ba

if (steptime .eq. 0) thengydF4y2Ba

101 I = 1,nblockgydF4y2Ba

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

打印*，“失败压力”gydF4y2Ba

打印*,XFAIL(我)gydF4y2Ba

nten, j = 1吗gydF4y2Ba

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

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

ENDDOgydF4y2Ba

stateNew (4) = XFAIL(我)gydF4y2Ba

种子=种子+ 1gydF4y2Ba

101年继续gydF4y2Ba

endifgydF4y2Ba

做100i = 1,nblockgydF4y2Ba

CgydF4y2Ba

C把应力写成矩阵形式gydF4y2Ba

CgydF4y2Ba

j = 1,3吗gydF4y2Ba

STR(j,j) = stress sold (i,j)gydF4y2Ba

最后做gydF4y2Ba

STR(1,2) = stress sold (i,4)gydF4y2Ba

STR(2,1) = stress sold (i,4)gydF4y2Ba

STR(1,3) = stress sold (i,5)gydF4y2Ba

STR(3,1) = stress sold (i,5)gydF4y2Ba

STR(2,3) = stress sold (i,6)gydF4y2Ba

STR(3,2) = stress sold (i,6)gydF4y2Ba

求应力矩阵的特征值gydF4y2Ba

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

cevr是主应力gydF4y2Ba

3, j = 1吗gydF4y2Ba

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

ENDDOgydF4y2Ba

打印*，'旧主重音'gydF4y2Ba

打印*,XPRIN(我,:)gydF4y2Ba

CgydF4y2Ba

C把应力写成矩阵形式gydF4y2Ba

CgydF4y2Ba

j = 1,3吗gydF4y2Ba

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

最后做gydF4y2Ba

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

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

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

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

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

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

求应力新矩阵的特征值gydF4y2Ba

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

cevr是主应力gydF4y2Ba

3, j = 1吗gydF4y2Ba

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

ENDDOgydF4y2Ba

打印*，'新主应力'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 . 0) THENgydF4y2Ba

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

其他的gydF4y2Ba

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

IF (XPRIN(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

C x是等式右边的值81gydF4y2Ba

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

C xc =顺应张量，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

用guass Jordan方法求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 ngydF4y2Ba

30i =1, ngydF4y2Ba

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

30日继续gydF4y2Ba

40继续gydF4y2Ba

我= 1,NgydF4y2Ba

J = N + 1, 2 N *gydF4y2Ba

If ((j-i) .eq.n)则gydF4y2Ba

B (I, J) = 1.0gydF4y2Ba

其他的gydF4y2Ba

B (I, J) = 0.0gydF4y2Ba

ENDIFgydF4y2Ba

ENDDOgydF4y2Ba

ENDDOgydF4y2Ba

160 i =1, ngydF4y2Ba

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

150 j = 1,ngydF4y2Ba

If (i.ne.j) thengydF4y2Ba

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 >

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

聚合物的Shope记忆效应?谁有很好的经验?请在这方面给我一个很大的帮助。提前谢谢你。我的邮箱是gydF4y2Baliqinan1986@gmail.comgydF4y2Ba

Hy,< / p >

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

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

提前谢谢你< / p >

我的电子邮件:gydF4y2Baabdelhak.azzi@univ-poitiers.frgydF4y2Ba

你好< / p >

我想在VUMAT中建模FG材料。为此，我需要每个材料点的x和y分量。但是Abaqus在每个nblock中只给我一个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/WritinggydF4y2Bapdf的用户子例程< / p >

祝你好运< / p >

弗兰克< / p >

大家好!< / p >

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

谢谢< / p >

即使你忽略了最初的价格，你也需要一个Intel XE编译器。编译器需要Visual studio。还有英特尔的图形处理器。< / p >

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