孔隙弹性材料固结过程中临界初始时间步长的问题。

亲爱的所有,

我有一个问题，关于临界初始时间步在孔弹性材料在固结过程中。

我们知道，对于隐式积分法，它是无条件稳定的。但对于孔隙弹性材料的固结，

有一个由Vermeer和Verruijt(1981)导出的临界时间步长t_critical，它也被包含在abaqus软件(http://129.25.22.58:2180/v6.7/books/exa/default.htm)。

据我所知，使用初始时间步长δ ta_t< t_critical，在压力结果中会有固结。

但我读了一些Massimiliano Ferronato等人(2001)的论文“有限元孔隙弹性方程的病态调节”，

他们提到，“t_critical的概念基本上与找到解决方案的难度有关，而不是与解决方案本身的准确性有关，如果愿意支付足够高的计算成本，那么对于小于t_critical的时间步长，解决方案也可能证明相当准确”。

从上面的陈述来看，初始时间步长delta

我的问题比这更复杂，我的目标是模拟脑脊液循环系统，使用脑脊液填充区域(如脑室和蛛网膜下腔)和脑组织的不同渗透率。模型脑脊液为纯达西氏层，脑组织为多孔弹性材料。

从结果来看，我在压力溶液中没有看到任何溶解。但我应该相信压力的结果吗?我感兴趣的是200秒内的压强，但是根据t_critical方程，估计的t_critical可能已经大于200秒了，这是相当矛盾的。这使得初始时间步长不可能小于t_critical。

根据Massimiliano Ferronato等人(2001)的说法，它似乎是可信的，但根据Abaqus手册，一个人应该总是有一个大于t_critical?

如果能得到一些建议，我将不胜感激。

非常感谢。

BR,

X李