错误警告：inistate data cannot be edited after the first load step

xzyzqs

执行下列命令流，出现了rt的错误信息：
（前处理省略）
/solu
time,1
autots,on
nlgeom,on
pred,on
NROPT,UNSYM     !打开非对称求解器
outres,all,all
allsel
nsel,s,loc,x,0  
dsym,symm,y    !左侧边界对称边界
nsel,s,loc,x,4.5
d,all,ux
nsel,s,loc,y,0
d,all,all
nsel,s,loc,y,9
d,all,all
allsel
acel,0,9.8
inistate,write,1,,,,,s   !写入初始应力文件，其中1表示生成初始文件，s表示Output stresses
save
solve

time,2
inistate,read,epm,,'F:\numerical simulation\Energy piles\normal pile\MT'   !读入初始应力文件
solve

time,3
sfl,9,pres,6e3       !600kPa
allsel
solve
save
finish

czthisme

论坛真的好。。

htbbzzg

命令流的以下部分：
time,1
.
.
.
inistate,write,1,,,,,s   !写入初始应力文件，其中1表示生成初始文件，s表示Output stresses
save
solve

time,2
inistate,read,epm,,'F:\numerical simulation\Energy piles\normal pile\MT'   !读入初始应力文件
solve

第一步定义初应力是可以的，但是求解后，应力状态已经发生了变化，第二步又定义初应力就不对了，错误信息的意思就是指出这一点。不知道你这样做是要模拟什么情况？

xzyzqs

我模拟桩基中荷载分布：先地应力平衡，然后在桩顶施加一个力荷载。模拟的方法是：第一步：将施加重力荷载计算初始应力场，然后采用inistate,write,1,,,,,s 写入初始应力文件；第二步（地应力平衡），我采用inistate,readd读取初始应力场计算，获取地应力平衡后的位移场（效果）；第三步:施加桩顶荷载，得到桩内荷载分布。之前参考了师访博士的书《ANSYS二次开发及应用实例详解》第四章中采矿过程应力分析命令流（其中也存在地应力平衡的过程），将前处理、求解、后处理分别写成宏文件，然后采用重启动的方式实现。但是我直接写成一个inp文件，之前是将各步分别/solu，solve,finish，出现的问题是我只能看到最后一步的结果（time3）;借鉴一些专家的意见，将第二步和第三步的/solu删除了，然后再求解，出现了最上面错误信息，按楼上的理解应该是不能在第二步定义初始应力（但是两次定义初始应力，其实一次是写入，一次是读取，会矛盾？）。
那么我现在再试试重启动分析，如果成功了，我再来说下！

htbbzzg

    这种情况，第二步不需要读取第一步计算的结果，软件会自动用第一步的结果作为初始条件，再考虑第二步的载荷，然后继续计算。

xzyzqs

htbbzzg 发表于 2015-3-10 16:14
这种情况，第二步不需要读取第一步计算的结果，软件会自动用第一步的结果作为初始条件，再考虑第二步的 ...

很高兴你能参与我的讨论，我一二部的目的是为了地应力平衡，即：得到重力作用下的应力场，又消除了因为重力而产生的位移场，然后在这样的基础上再去求解桩顶力荷载荷载诱发的应力。在所有涉及到地应力平衡似乎都用到这样一对命令：inistate,wirte(提取重力场文件)  和 inistate,read（抵消位移场）。

不知道htbbzzg是否做过地应力场平衡（只需要重力引发的应力场（作为初始条件），而不需要位移。）

至于上面的问题，我重新考虑了重启动分析，已经可以得到结果了！

xzyzqs

很高兴你能参与我的讨论，我一二部的目的是为了地应力平衡，即：得到重力作用下的应力场，又消除了因为重力而产生的位移场，然后在这样的基础上再去求解桩顶力荷载荷载诱发的应力。在所有涉及到地应力平衡似乎都用到这样一对命令：inistate,wirte(提取重力场文件)  和 inistate,read（抵消位移场）。

不知道htbbzzg是否做过地应力场平衡（只需要重力引发的应力场（作为初始条件），而不需要位移。）

至于上面的问题，我重新考虑了重启动分析，已经可以得到结果了！

htbbzzg

  如果是这种情况，试试以下做法：
  1  第一步：做重力分析，将应力结果写出到文件，作为下一次分析的初应力。然后退出 /solu
  2  再次进入 /solu，开始另一次分析，使用第一次计算存储的 (重力分析的) 应力结果文件作为初应力，然后施加其它载荷继续计算。
  这里的关键是：两次分析都是独立的，即第二步计算并不是从第一步的结果继续进行的，原因是你只使用了第一步 (单纯重力分析) 结果中的应力结果，而没有使用位移结果，因此并不是第一步的继续。
  这种做法有一个问题：由于两次分析是互相独立的，因此最终的结果文件中只有第二次分析的结果，第一次的分析结果不再存在。解决办法是：在完成第一步分析后，修改工作文件名，这样，第二步的结果将会保存在另一个结果文件中，第一步的结果还是在最初的结果文件中。  如果采用重启动，那么第二步的结果中仍会包含第一步的位移结果。

xzyzqs

我在分析完第一步后添加/filname,epm1来修改工作名，部分命令流如下：
。。。
solve
finish

/filname,epm1      !前面的文件名是epm
/solu
antype,,restart
time,2
....
solve
save
finish


但是出现了下面的错误信息：
File epm1.esav not found. Previous solution(required for restarting) cannot be recovered.

htbbzzg

这种情况，你需要将原来以 epm 命名的文件都改名为以 epm1 命名的的文件，或者拷贝为以 epm1 命名的文件。