巷道围岩流变模拟讨论

cumtmbqq1

模型简介：埋深800，测压系数1.5，地应力平衡后，水平应力30MPA，垂直应力20MPA,采用cvisc(基于摩尔库伦准则的博格斯粘塑性流变模型)模型，计算过程中出现以下两种错误：
   （1）设置set large 时，算一会就会出现  非法几何 ，即网格严重扭曲，但是我查看后发现，只有
        个别单元变形超级变态，不合乎常理，而且是个别单元突然变形特别离谱。
   （2）不用set large 时，计算时突然所有模型全部进入塑性状态，接着出现错误，数值溢出。

采用如下方法进行了调试，但问题均没有解决：
    （1）怀疑是几何模型的问题，于是生成了很简单的几何模型，用同样的参数计算，出现同样错误。
    （2）怀疑是地应力过大导致单元变形离谱，尝试了，把地应力改为之前的1/3，计算成功，改为1/2，又计算失败，尝试增加粘聚力，结果的确收敛，但是塑性区范围很小，达不到我的目的。
    （3）最后还尝试用摩尔库伦准则算，还是出现同样的问题。
命令文件：

各位同仁是否遇到这样的问题，有什么比较好的解决方法分享一下吧，小弟感激不尽。

cheryboy

学习一下

cid2012

自重计算开大变形干什么，建议你看看彭文斌那本书第一个例子。

cumtmbqq1

cid2012 发表于 2011-12-31 20:37
自重计算开大变形干什么，建议你看看彭文斌那本书第一个例子。

自重计算为什么不能考虑大变形，小变形也是一种特殊的大变形，道理就跟：一次函数是特殊的二次函数一样。。。。之所以自重考虑大变形，主要是因为大变形考虑了网格的变形，自重作用下，本身是会产生网格变形的。
个人觉得，现在很多数值计算相关书籍，都是抄袭帮助文档，没什么看头。。。。
此问题已经解决，多谢。。。。

starhan1314

cumtmbqq1 发表于 2012-1-2 17:03
自重计算为什么不能考虑大变形，小变形也是一种特殊的大变形，道理就跟：一次函数是特殊的二次函数一样。 ...

楼主说说是为什么，让小弟也学习学习，初来咋到

cumtmbqq1

starhan1314 发表于 2012-1-3 11:45
楼主说说是为什么，让小弟也学习学习，初来咋到

大变形与小变形的不同即在于，位移与应变的关系不同。
小变形，即位移与应变成线性关系，即应变等于位移对坐标的一阶导。。
大变形，位移与应变呈非线性关系，即我们平时所说的几何非线性，应变等于位移的一阶导加上其他高阶项。当高阶项等于0时，不就成为小变形了。。。就和二次函数AX^2+BX+C=0一样，当A等于0时，就成为了一次函数。
所以，二次函数的球根公式使用与一次函数，令A等于零就是了。同样，大变形的算法，也适用于小变形。。。

starhan1314

cumtmbqq1 发表于 2012-1-3 11:53
大变形与小变形的不同即在于，位移与应变的关系不同。
小变形，即位移与应变成线性关系，即应变等于位移 ...

谢谢，受教了

lforever

学习了！

GodPrays

这么说大变形和小变形没什么区别了?大变形包括小变形？不明白...

yuanlin0916

太深奥了看不懂，不过还是谢谢分享

yuanlin0916

楼主最后是怎么解决这个问题的呢，能不能八一下

laiyonghui

楼主能不能把你的已经成功解决的txt程序流发下，我最近也在做这个，73915212@qq.com，谢谢你啊