请问叶老师: 模型参考《线性波浪加载下海底斜坡失稳机制的数值分析》如图所示。 海床:底边xy固定、两边x固定。底边长684m,左侧高30m,右侧高45m,边坡两边水平距离相等,坡长171m,坡高15m。 一阶stokes波:周期8.6s、波高25m、水深10m、静水位线采用边坡中点高度加水深(45+30)/2+10m 本构:摩尔库伦:粘聚力20kpa、摩擦角2度 分析步1:static 计算时间1s 不更新材料、刚度矩阵、坐标 分析步2:consolidation 计算时间600s 不更新材料、刚度矩阵、坐标 问题1:软件中力对节点压为负、拉为正。得出的结果所有节点最大主应力都为负如图。最小主应力如图,在坡肩右侧海床表面,对节点沿着海床水平方向的小主应力出现了正值即小主应力对节点为拉,我认为边坡节点大小主应力都应该为压应力。请问这个现象是正常的吗?边坡在波浪作用下会出现主应力为拉的情况吗? 问题2:我发现大小主应力一正一负的时候,min principal stress输出的是绝对值小的数值而不是代数值小的数值,结果如图。我这样理解对吗? 问题3:怎么找到海床在波浪作用下的临界状态的时间?取计算报错的前一秒或者看位移时程曲线突变的时刻吗?
针对问题一,我认为在坡肩部位出现小范围的水平方向拉应力是很正常的现象,之前的研究表明,坡肩部位出现水平拉应力是正确的,尤其是土比较软弱时。
针对问题二:小主应力应该是莫尔圆坐标系中横轴上最左边的那个值,和正负无关。关于这一点,软件研发团队需要进行确认并修复问题。
针对问题三:边坡的临界状态应该是坡体发生变形突变的状态,而不应该以计算不收敛的标志。