| 查看: 1559 | 回复: 7 | ||||
[求助]
ansys仿真过程中循环力加载的次数增加，同一点特定时间的应力值也发生改变已有1人参与
|
|
最近在做CT试样裂纹扩展仿真，我写 ansys 程序的时候遇到了两个问题。 我用一个形状为三角波的循环力去拉伸一个CT试样。 首先我写了一个 0 到 6 秒的循环，仿真成功结束，再次查看了第 6 秒最大应力值是66Mpa 我又写了一个 0 到 7秒的循环，程序在6.12秒就终止报错了，然后查看了第 6 秒的最大应力值，是295Mpa。 我又写了一个 0 到 11 秒的循环，仿真成功结束，再次查看了第 6 秒最大应力值是291Mpa。 我认为第 6 秒的应力值应该和循环次数没有关系，也就是说，循环次数只要大于 6 次，相同点第六秒的应力值都应该是保持不变的。 另外如果0-12秒能够循环，0-7秒也能够循环，不清楚为什么0-7秒的循环终止了。 三种情况的具体问题描述和相关ansys程序都在附件里了。 请大家帮忙解答一下我的疑惑 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ansys裂纹扩展仿真出现的问题描述.zip
2016-11-25 18:44:04, 290 K
» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
| 祝福 |
2楼2016-11-25 19:05:07
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
honorlyw: 金币+30, ★有帮助 2016-11-27 12:18:56
感谢参与，应助指数 +1
honorlyw: 金币+30, ★有帮助 2016-11-27 12:18:56
|
看了你的命令流，谈一点看法： 1 你的长度单位是什么？ 如果按照你的材料性能，长度单位似乎是 mm，而按照你定义关键点的数据，长度单位又似乎是 m，不知到底是什么？ 2 你的三段命令流，分别定义了不同的载荷工况，并使用 lswrite 写了相应的载荷工况文件，但是你的求解仍然使用 solve 命令，这样的话，将只是求解最后一个载荷步，而不是顺序求解各个载荷步。因此，6s 的情况，最后一步的压力是 1.76，结果很小；7s 和 11s 最后的载荷是 35.2，因此结果大得多。 正确的方法是使用 lssove, imin, imax, iinc (最小载荷步，最大载荷步及载荷步增量)。比如你的 6s 的情况，应该是： Lssove,1,6,1 3 你使用了理想塑性材料，屈服应力 345 MPa，但是 7s 和 11s 的最大应力都没有超过屈服应力，这种情况，最大应力与循环次数应该关系不大。 试算了一下，我的机器较差，费时太多，所以 7s 的情况没有算到底，是否不收敛或为什么不收敛，不太清楚。照理，改用 lssolve 后都应该收敛的。 另外，由于三种情况在 6s 时的压力不同，6s 时的结果还是不会一样的。 |
3楼2016-11-26 11:58:35
4楼2016-11-26 15:26:18
|
我按照你的建议改了一下程序，仿真在1s的时候报错了。 请你帮我看一下哪里出问题了。 !使用*do循环的0-11秒循环的程序段 finish /clear,start /COM, Structural /PREP7 ET,1,PLANE183 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,3.45e2,,,,, K,1,0,0.25,0, K,2,0,30,0, K,3,62.5,30,0, K,4,62.5,0,0, K,5,39.5,0,0, K,6,36.5,0,0, K,7,36.25,0.25,0, FLST,2,7,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 FITEM,2,5 FITEM,2,6 FITEM,2,7 A,P51X CYL4,12.5,16.25,6.25 ASBA,1,2 K,12,36.5,6,0, K,13,46.5,6,0, K,14,46.5,0,0, FLST,2,4,3 FITEM,2,6 FITEM,2,12 FITEM,2,13 FITEM,2,14 A,P51X aovlap,1,3 APLOT ESIZE,0.1,0 MSHKEY,0 AMESH,2 ESIZE,1,0 MSHKEY,0 AMESH,4 /SOL ANTYPE,4 TRNOPT,FULL LUMPM,0 NSEL,S,LOC,X,12.5 NSEL,R,LOC,Y,22.5 D,P51X, , , , , ,UX, , , , , allsel，all FLST,2,2,4,ORDE,2 FITEM,2,16 FITEM,2,-17 DL,P51X, ,UY, OUTRES,ALL,ALL, xc=5 *do,m,1,xc, TIME,2*m-1 AUTOTS,1 DELTIM,0.05,0.01,0.1,1 KBC,0 TSRES,ERASE FLST,2,2,4,ORDE,2 FITEM,2,8 FITEM,2,-9 SFL,P51X,PRES,35.2e0, LSWRITE,2*m-1, TIME,2*m AUTOTS,1 DELTIM,0.05,0.01,0.1,1 KBC,0 TSRES,ERASE FLST,2,2,4,ORDE,2 FITEM,2,8 FITEM,2,-9 SFL,P51X,PRES,1.76e0, LSWRITE,2*m, *enddo LSSOLVE,1,2*xc,1, FINISH |
5楼2016-11-26 16:42:26
6楼2016-11-26 17:02:54
7楼2016-11-26 17:46:10
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
【答案】应助回帖
|
试算了一下你的命令流，发现一些问题： 1 你的模型中没有对 Ux 的约束，你的命令流中有如下语句： NSEL,S,LOC,X,12.5 NSEL,R,LOC,Y,22.5 D,P51X, , , , , ,UX, , , , , 但是这里的 p51x 没有定义过，不能起作用，因而 Ux 的约束没有加上。若改为： D, 15955, , , , , ,UX, , , , , 则可以加上 Ux 的约束，但是计算时不到 Time=1，就不收敛了。 2 最关键的是：你的压力载荷 35.2 MPa (近 350 个大气压) 太大了，结构早早的进入了塑性。 删除你的命令流中材料的非线性部分，可以一直计算到底，得到的最大 Mises 应力达到 1741.0，远远超过了材料的屈服应力 345.0。 不知道你的实际载荷是多少？这方面需要仔细检查一下。 |
8楼2016-11-27 12:09:15