清华大学bbs的abaqus精华

清华大学清华大学 bbs 的的 abaqus 精华精华 大家看看吧 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 的维数 depvar 20 看来我的错误不在 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 等那样输出 只能在 UVARM 中定义 但 同 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 不知谁还有更好的方法 听说清华有人 abaqus 用的很熟了 大家过奖了 我是没办法 要出论文 总得学习 上午版主还在 qq 中指点了我一番 YOU CAN TRY AS FOLLOWS STEP EXTRAPOLATION NO INC 2000000 STATIC 0 001 500 0 0 00001 0 1 可以 并且在 step 里好像有选项选择 你是手写的程序么 新建一个 step 然后定义边条 在上一 step 定义边条的节点上重新定义边条 而且 好像 是不是有个选项是 new 这样就可以取消上一 step 的边条设置了 在以后的载荷步里重新定义边界条件 加参数 new 比如在第二个载荷步里 STEP BOUNDARY CONDITIONS OP NEW END STEP Command 方式 abaqus job job name user source file 或者 obj file CAE 方式 在最后 submit 的时候可以选择 source file 或者 obj file 不太明白你的意思 如果只是要得到所需要的如应力 应变等 那么在数据行请求 输出就行了啊 而且个人觉得 abaqus 的后处理功能基本上也够用了 直接输出到 odb 文 件 中 用 CAE 或者 Viewer 后处理即可 而 el print 输出的 dat 文件包含其他的信息 如版本号 处理内容等 若想要直接输出到其他后处理软件的话 可以输出到 result file 文件 fil 不过 这可是二进制文件或者 ASCII 文件 还是需要处理的 或者直接用 Hypermesh 等后处理 程序打开做图 好像是在 tool display create 再选定一个单元集或者是节点集 然后操作是 取代 还是其他等等 就是里面的两个圆圈相交 还是想减 你试试看看 这个是要到后处理才行吧 非要计算以后形成 odb 吗 前处理里面怎么才可以啊 Python 在 abaqus 中可以运行 scrpit 文件 你在 CAE manual 中就 可以将其附录中的例子存为 py 文件 然后在 CAE 的菜单中 run script 就可以运行了 BTW 请问你的是 abaqus 的那个版本 6 2 or 6 3 6 3的有盗版了 1 请问如何得到 M C K 矩阵 用什么命令 试一试 element matrix output 2 想要用 Newmark 方法求解 求解问题 在 ABAQUS 中如何实现 或者 那里能找到类似的例子 看看 standard manuals analysis 的 dynamics analysis 和 explicit manuals 3 ABAQUS 中能实现这样的东西吗 就是我需要平滑 ABAQUS 产生的位移场 还是这个平滑只能在其他环境中完成 试试这个输出选项 El print Position average at Nodes abaqus CAE 好像不能生成节点集 只能定义某个实体 后处理用这个实体上的节点 看到 CAE 的例题都是这样做的 NODE SET 是 abaqus 的基本功能 也是其方便使用之处 在 CAE 中是支持的 这个例子我做过 没有问题 我用 tool set create 居然不能选择由 CAE 生成的节点 不是顶点 这是因为你没有正确操作 Partition Edge 这一步 在 sunhaolan 的大作中提到 abaqus CAE 好像不能生成节点集 只能定义某个实体 后处理用这个实体上的节点 ABAQUS divides all of its analysis procedures into two main groups linear perturbation and general ABAQUS makes a very clear distinction between perturbation and general analysis procedures because loading conditions and time are defined differently for the two cases Furthermore the results from each type of procedure should be interpreted differently The response of the model during a general analysis procedure known as a general step may be either nonlinear or linear In a step that uses a perturbation procedure which is called a perturbation step the response can only be linear ABAQUS treats such steps as a linear perturbation about the preloaded predeformed state known as the base state created by any previous general steps therefore its capability for doing linear simulations is rather more general than that of a purely linear analysis program 各位大侠 我在定义 step 的时候 前三个用的是 standard 中的相关命令 在 step 4的时候 我想用 explicit 来求解 并且定义了 amplitude 来模拟一个 位移随时间的简单变化 命令如下 step name nlgeom inc 1000 好像 explicit 里没有 inc 这个参数 dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 刚体参照点 2 2 0 03 end step 有问题吗 另外 standard 的命令和 explicit 的命令能出现在同一个 inp 文件中吗 要转换的 ABAQUS Standard and ABAQUS Explicit are separate program modules with different data structures hence the explicit dynamics procedure cannot be used in the same analysis as any of the procedures in ABAQUS Standard However ABAQUS provides a capability to import a deformed mesh and associated material state from ABAQUS Explicit into ABAQUS Standard and vice versa This procedure is described in Transferring results between ABAQUS Explicit and ABAQUS Standard Section 7 6 2 of the ABAQUS Standard User s Manual 在 linyuwho 苍龙有悔 的大作中提到 各位大侠 我在定义 step 的时候 前三个用的是 standard 中的相关命令 在 step 4的时候 我想用 explicit 来求解 并且定义了 amplitude 来模拟一个 位移随时间的简单变化 命令如下 step name nlgeom inc 1000 好像 explicit 里没有 inc 这个参数 dynamic explicit 2 amplitude name ramp 0 0 1 1 2 1 boundary amplitude ramp refpunch 刚体参照点 2 2 0 03 end step feature 菜单下有删除一项 基准点是一种 feature 工具栏不有个像铅笔一样的的东东么 我想把计算出来的应力 应变 位移等数据提取出来 在其他软件里处理 请问该怎么办 我在 dat 里没有找到这些数据阿 res 不是文本格式 其他的好像也没有啊 计算结果 没有文本格式 只能自己把数据倒出来 在显示结果的模块中 名字忘了 就是 sketch 前面那个 在 report 菜单中 可以到处自己需要的文件 文件格式是 rpt 谢谢 但是这个为什么只有应力和应变呢 计算出的位移的数据怎么会没有呢 在 santong supernova 的大作中提到 计算结果 没有文本格式 只能自己把数据倒出来 在显示结果的模块中 名字忘了 就是 sketch 前面那个 在 report 菜单中 可以到处自己需要的文件 文件格式是 rpt 发信人 wjytj 赖子 信区 FEA 你要在 STEP 模块 field output 和 history output 改变默认 请求输出位移 U 或者在 step 中加上 NODE PRINT U 计算完成后在 DAT 文件中就可以提取节点位移 请问 field output 和 history output 什么区别 关键字 node output 和 node print 有什么区别 谢谢 我的 inp 数据文件里有这么几行 Output field frequency 99999 Node Output U RF Element Output S E EE Output history variable PRESELECT El Print freq 999999 Node Print freq 999999 其中在 history output 中有 node print 关键字 但是没有指定输入 U 最后导出的数据没有 U 这一选项 请问是什么原因 在 wjytj 赖子 的大作中提到 你要在 STEP 模块 field output 和 history output 改变默认 请求输出位移 U 或者在 step 中加上 NODE PRINT U 计算完成后在 DAT 文件中就可以提取节点位移 发信人 Mikejwg mikejwg 信区 FEA field output 和 history output 这在 viewer 模块中很明显 field output 是场量输出 history output 为历史数据输出 会记录场量随时间的曲线 至于 node output 与 node print 的区别在于他们写入不同的后缀文件 比如 odb dat 等 在 help 中有详细介绍 在 wsun 的大作中提到 请问 field output 和 history output 什么区别 关键字 node output 和 node print 有什么区别 谢谢 我的 inp 数据文件里有这么几行 Output field frequency 99999 Node Output U RF Element Output S E EE Output history variable PRESELECT El Print freq 999999 发信人 wjytj 赖子 信区 FEA 然也 补充一句 node output 是给 CAR 或 VIEWER 做后处理用的 是二进制文件 而 node print 则是写到 DAT 文件中 你可以自己看的文本文件 在 Mikejwg 的大作中提到 field output 和 history output 这在 viewer 模块中很明显 field output 是场量输出 history output 为历史数据输出 会记录场量随时间的曲 至于 node output 与 node print 的区别在于他们写入不同的后缀文件 比如 odb dat 等 在 help 中有详细介绍 在 wsun 的大作中提到 请问 field output 和 history output 什么区别 关键字 node output 和 node print 有什么区别 谢谢 我的 inp 数据文件里有这么几行 Output field frequency 99999 Node Output 发信人 octopuszy 猴哥 信区 FEA node output 是在 odb 文件中输出 br node print 是在 dat 文件中输出 在 Mikejwg mikejwg 的大作中提到 field output 和 history output 这在 viewer 模块中很明显 field output 是场量输出 history output 为历史数据输出 会记录场量随时间的曲 至于 node output 与 node print 的区别在于他们写入不同的后缀文件 比如 odb dat 等 在 help 中有详细介绍 file print file 文件类型选择 png 或 tiff 在 Dshmy 的大作中提到 谢谢 一般说可以这样理解 ddsdde 有点像我们常说的弹塑性矩阵 Dep 在 wjytj 的大作中提到 在 USER MATERIAL 的定义中 对 的偏倒数 即 DDSDDE 被称为 material stiffness matrix 而在 UMAT 中 DDSDDE 被称为 material s Jacobian matrix 请问 DDSDDE 和材料的切线刚度矩阵的关系是什么 二者是一个概念么 第二个 在 octopuszy 猴哥 的大作中提到 偶在 umat 中调用求主应力函数 CALL SPRINC STRESS PS LSTR NDI NSHR 后 存储主应力得数组 PS 中 各个主应力排列顺序是什么 PS1 S2 S3 PS1 S2S3 PS2 谢谢 abaqus job jobname parametercheck abaqus job jobname datacheck abaqus job jobname continue 在 gentlewtj gentlewtj 的大作中提到 那有没有检查错误的方法 在 uptonow 到目前为止 的大作中提到 1 在边界条件和加载时 总是有 initial 这个步 然后是我们自己定义的 加载步 请问这个 initial 步 主要作用是什么 能不能去掉 不能去掉 所有的分析都有 是默认的步 2 A solid extrusion base feature 这句话是什么意思 extrusion revolution 等是什么意思 这两的是三维建模时候 在画完二维图形 如何来生成三维图形 extrusion 意思是你给定一个厚度 然后二维图形第三个方向上面伸展这么多形成三维图形 revolution 意思是你给定一个旋转轴 二维图形绕其旋转后形成三维轴对称图形 谢谢 发信人 fangxj fxj 信区 FEA 标 题 Re abaqus 中接触问题请教 发信站 BBS 水木清华站 Fri Jun 14 15 51 20 2002 接触问题不收敛有两个方面不妨试试 一 在 CONTACT PAIR 里调试 ADJUST 参数 二 调一些模型参数 比如 FRICTION 等 在 octopuszy 猴哥 的大作中提到 偶的模型中存在两个物体的接触 计算过程中报错 WARNING OVERCLOSURE OF CONTACT SURFACES FOUNDATIONOUTSIDE and BUCKETOUTSIDE IS TOO SEVERE CUTBACK WILL RESULT YOU MAY WANT TO CHANGE THE VALUE OF HCRIT ON THE CONTACT PAIR OPTION WARNING CONVERGENCE JUDGED UNLIKELY INCREMENT WILL BE ATTEMPTED AGAIN WITH A TIME INCREMENT OF 9 76563E 05 ERROR TOO MANY ATTEMPTS MADE FOR THIS INCREMENT ANALYSIS TERMINATED 不知该如何解决 standard 使用隐式作动力分析 explicit 用显式作动力分析 你该用 explicit 在 uptonow 到目前为止 的大作中提到 请问这两本书有什么不同 我做的是圆形杯子的深冲的模拟过程 我这里只有 standard 一书 那位有 explicit 书的大虾请留个言 我想复印 谢谢 一般说来 你那样是正确的 如果不行 你可以在 dat log msg 中找到中断的增量步 假如为60 不妨可以 在 res pulse 中 restart read step 1 inc 59 在 ypm 红塔山 的大作中提到 我的一个 abaqus 程序 结构动力学 用的 standard 模块 假设名 为 pulse inp 一个 step inc 1000 restart 命令的参数为 restart write frequency 10 在算了一半之后人为停掉 然 后我想从中止的地方继续运算 我定义一个 res pulse inp 其中就一个命令 restart read 然后在 command 窗口 输入 abaqus job res pulse oldjob pulse 结果程序从起始点开始运算 而不是从中断的地方 请问是什么 原因 刚学不久 希望大侠们指点指点 thanks 发信人 Mikejwg mikejwg 信区 FEA uptonow 一直叫我在水木清华上写一点东西 我一直忙于毕业论文 再说也是困难一个一个 出现 永远有解不完的问题 就本人体会 写一点东西 在 umat 中 我总结三点 1 求试算应力 stress stress De dstran 其中 De 为弹性矩阵 2 据应力势函数 屈服函数 求解满足条件的应力 应变 在此我们要注意试算应力与真实应力的关系 sigma true Sigma tri De dpstran dpstr an 为塑性应变增量 据此满足屈服条件 F Sigmaij 0 求得真实的应力 更新应力 stress stress De dpstran 3 与2联合 弹塑性区域的转换 判断是否屈服 4 组集 ddsdde 与个人的本构有关 这4步都要在 UMAT 实现 采用上一迭代步构造 ddsdde 然后 dsigma ddsdde dstran 不可取 以我个人经验 收敛很 慢 并且敏感与步长 希望大家交流 发信人 uptonow 为什么受伤的总是我 信区 FEA 衷心感谢 Mikejwg 大侠的该出手时 就出手 希望 Mikejwg 大侠以后经常出手 造福于 UMAT 的万民 在 Mikejwg mikejwg 的大作中提到 uptonow 一直叫我在水木清华上写一点东西 我一直忙于毕业论文 再说也是困难一 出现 永远有解不完的问题 就本人体会 写一点东西 在 umat 中 我总结三点 1 求试算应力 stress stress De dstran 其中 De 为弹性矩阵 2 据应力势函数 屈服函数 求解满足条件的应力 应变 在此我们要注意试算应力与真实应力的关系 sigma true Sigma tri De dpstran d an 为塑性应变增量 据此满足屈服条件 F Sigmaij 0 求得真实的应力 更新应力 stress stress De dpstran 3 与2联合 弹塑性区域的转换 判断是否屈服 4 组集 ddsdde 与个人的本构有关 发信人 wjytj 赖子 信区 FEA 精华区有一些帖子 很有帮助 大家感兴趣可以去看看 在 Mikejwg 的大作中提到 uptonow 一直叫我在水木清华上写一点东西 我一直忙于毕业论文 再说也是困难一 出现 永远有解不完的问题 就本人体会 写一点东西 在 umat 中 我总结三点 1 求试算应力 stress stress De dstran 其中 De 为弹性矩阵 2 据应力势函数 屈服函数 求解满足条件的应力 应变 在此我们要注意试算应力与真实应力的关系 sigma true Sigma tri De dpstran d an 为塑性应变增量 据此满足屈服条件 F Sigmaij 0 求得真实的应力 更新应力 stress stress De dpstran 3 与2联合 弹塑性区域的转换 判断是否屈服 发信人 uptonow 为什么受伤的总是我 信区 FEA 精华区的帖子太少了 还需要你补充呢 不要推辞了 推辞就是逃避责任呵呵 开玩笑的 在 wjytj 赖子 的大作中提到 精华区有一些帖子 很有帮助 大家感兴趣可以去看看 发信人 edodo 大漠孤沙 信区 FEA 老哥可以挺身而出呀 哈哈 期待中 在 uptonow 为什么受伤的总是我 的大作中提到 精华区的帖子太少了 还需要你补充呢 不要推辞了 推辞就是逃避责任呵呵 开玩笑的 发信人 air1 无影如风 勤奋36小时 永不言败 信区 FEA 我曾经犯的错误 1 发现456方向的应变是工程应变 2 全量法可以直接根据更新后的应变求应力 单独求 ddsdde 矩阵 3 ddsdde 何时更新对于收敛性很有影响 在 Mikejwg mikejwg 的大作中提到 uptonow 一直叫我在水木清华上写一点东西 我一直忙于毕业论文 再说也是困难一 出现 永远有解不完的问题 就本人体会 写一点东西 在 umat 中 我总结三点 1 求试算应力 stress stress De dstran 其中 De 为弹性矩阵 2 据应力势函数 屈服函数 求解满足条件的应力 应变 在此我们要注意试算应力与真实应力的关系 sigma true Sigma tri De dpstran d an 为塑性应变增量 据此满足屈服条件 F Sigmaij 0 求得真实的应力 更新应力 stress stress De dpstran 3 与2联合 弹塑性区域的转换 判断是否屈服 4 组集 ddsdde 与个人的本构有关 发信人 wjytj 赖子 信区 FEA 1 发现456方向的应变是工程应变 说白了 因为在 UMAT 中应变是以列向量的形式存储的 而初始刚度 DDSDDE 的形式为 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 所以 应变是工程应变 而如果 DDSDDE 的形式为 直接用张量表示 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 2 0 0 0 0 0 0 2 那么 应变就应该是应变张量中的应变了 也就是说 采用工程应变只是为了方便应用而 已 比如 求弹性应变能 可以直接将应力列向量的转置乘以应变列向量 显然 这在编 程中很容易实现 2 全量法可以直接根据更新后的应变求应力 单独求 ddsdde 矩阵 同感 应该可以减少不少工作量 3 ddsdde 何时更新对于收敛性很有影响 希望 air1能够说得更详细些 比如在开始是可以采用初始刚度 在后来可以采用切线刚度 可是如何控制 编程如何实现呢 在 air1 的大作中提到 我曾经犯的错误 1 发现456方向的应变是工程应变 2 全量法可以直接根据更新后的应变求应力 单独求 ddsdde 矩阵 3 ddsdde 何时更新对于收敛性很有影响 在 Mikejwg mikejwg 的大作中提到 uptonow 一直叫我在水木清华上写一点东西 我一直忙于毕业论文 再说也是困难 出现 永远有解不完的问题 就本人体会 写一点东西 在 umat 中 我总结三点 1 求试算应力 stress stress De dstran 其中 De 为弹性矩阵 发信人 wjytj 赖子 信区 FEA 2 全量法可以直接根据更新后的应变求应力 单独求 ddsdde 矩阵 在我的本构中 DDSDDE 的确定需要知道 STRESS 中的某些量 那么 请问 按照你的意思和经验 在 UMAT 中是否可以先得到 STRESS 比如 利用全量法 然后再得到 DDSDDE 如果可以的话 就应该可以减少迭代求解的次数了 谢谢 假定采用 mises 屈服准则 1 在调用 UMAT 之前 ABAQUS 传递给 UMAT 本次增量开始时的应力 sigma 0 总应变 E 应变增 量 delta E 状态变量保存 弹性应变 塑性应变 等效塑性应变 请问 状态变量保存的弹性应变 塑性应变是否等于 ABAQUS 传递给 UMAT 的总应变 2 然后在 UMAT 中利用上述的 ABAQUS 传递的量和状态变量得到 DDSDDE 矩阵 然后返 回给 ABAQUS ABAQUS 根据 delta sigma ddsdde delta E 并且得到本次增量结束时的应力 sigma 1 sigma 0 delta sigma sigma 0 ddsdde delta E 3 然后更新本次增量结束时的状态变量 弹性应变 塑性应变 等效塑性应变以供下次调 用 UMAT 请问手册上 UMAT 必须更新应力 可是根据上述我的理解好像是 ABAQUS 根据 UMAT 提 供的本 次增量的 DDSDDE 在 ABAQUS 中更新 请问到底是怎么回事 谢谢 个人觉得 在 Abaqus 中更新应力是必须的 而且必须准确 这是你的本构的核心 DDSDDE 则可以是不准确的 比如 你可以保持 DDSDDE 一直为初始刚度矩阵 当然 这样做的问 题 会是收敛速度较慢或者根本不收敛 在 Manual 中也说过 DDSDDE 不影响结果 前提是 只要能得到结果 与大家共同讨论 在 uptonow 的大作中提到 假定采用 mises 屈服准则 1 在调用 UMAT 之前 ABAQUS 传递给 UMAT 本次增量开始时的应力 sigma 0 总应变 E 应 量 delta E 状态变量保存 弹性应变 塑性应变 等效塑性应变 请问 状态变量保存的弹性应变 塑性应变是否等于 ABAQUS 传递给 UMAT 的总应变 2 然后在 UMAT 中利用上述的 ABAQUS 传递的量和状态变量得到 DDSDDE 矩阵 然后返 回 ABAQUS ABAQUS 根据 delta sigma ddsdde delta E 并且得到本次增量结束时的应 sigma 1 sigma 0 delta sigma sigma 0 ddsdde delta E 3 然后更新本次增量结束时的状态变量 弹性应变 塑性应变 等效塑性应变以供下 期盼了好久的 abaqus6 3终于安装成功 寄希望于此来成功 import 我的 igs 模型 unigraphics v18生成的 但是 结果却不尽如人意 现将我的感受说一下 1 6 3在 import 进模型后 它的形状基本上与真实形状差不多 比6 2更像原始模型 但是 依然是 invalid geometry 2 经过 delete invalid geometry 后 它确实成了 imprecise geometry 意味着你 以划分网格了 但是此时的形状已经面目全非了 删除了太多的东西了 得到的 imprecise geometry 就毫无意义了 3