FLAC新本构模型开发成功：层状介质Cosserat本构

sunnyRSM

报告各位同仁，经过半年多的努力，将层状介质Cosserat本构嵌入Mohr-Coulomb模型，我们开发成功了一种FLAC新本构，可以分析软或者硬夹层对岩体破损的影响，以及对地下硐室稳定性的影响。
如图所示为Cosseart介质微元，首先算出平均宏观应力，然后利用微元内位移协调，分别求出各层应力，利用各自破坏准则判断是否破坏，再综合回复合体。
现模型刚刚调试成功，各位若有兴趣，欢迎一起探讨。

flacfighter

不错！

ying111

he

lakewater

利用FISH还是C++？

sunnyRSM

回版主，我们的模型是用C＋＋开发的，建立了动态链接库文件。

和内嵌MC或DP模块比，我们模块能考虑夹层破坏对材料进一步破损的影响，运行速度略慢一些。

但由于我们将薄层状体概化为均匀介质，可以不必细究各层厚度，而只考虑夹层以及含量便可，这样单元可以取得大一些，比如和层厚一般。综合起来速度快多了，建模过程也方便多了。

lakewater

可以将简单算例的结果公布开来，让大家看看计算效果。

sunnyRSM

好的，正在准备。
看来版主把关很严哦 :-)

要获得积分不易阿:-P

hja2000415

我的论文就是层状岩体的，希望楼主快点公布算例啊，我的qq：4776689希望和楼主交流。

dynamax

不错不错，这是我看到的国内第1个用C++开发的本构，赞一个

starsmoon

有机会象楼主请教啊。

geostability

lakewater wrote:
可以将简单算例的结果公布开来，让大家看看计算效果。

同意!!!!

lakewater

dynamax wrote:
不错不错，这是我看到的国内第1个用C++开发的本构，赞一个

岩土力学马上要看出一个用C++开发的流变本构，另外本人实现的也在审稿中。

sunnyRSM

现献上一计算结果，请各位指正。

红色区域为含硬夹层岩体（Cosserat 复合介质），夹层体积含量30％，夹层初始弹模（einter）21.78GPa，层面法向为铅垂方向，岩体弹模为5.15GPa。洞室内压为5.6MPa，可见洞周夹层弹模部分（蓝色部分）退化为初始弹模的1/2，即10.89GPa。

绿色区域为大夹层，和复合介质内夹层性质一样，采用内嵌MC模型，蓝色域为破坏单元。

我们的模型有待进一步完善，请各位多提意见！

yycgy

请教楼主，你的模型和复合材料中的层板模型有哪些异同之处？

sunnyRSM

谢谢版主的加分鼓励！我总算成为仿真论坛的一员了。

Cosserat介质本构的特点是考虑了局部应力场，局部的弯曲效应，引入了偶应力，这样导致应力张量和应变张量不一定对称。该本构能考虑材料的内尺度（internal scale length）的影响，如岩层厚度。

当然目前的FLAC是不能处理非对称应力张量的，故而作为第一步，我们先不考虑偶应力的影响，仅考虑各层细观位移协调导致各层应力状况不同，分别采用塑性流动法则来分析，最后综合各层应力得到平均（cosseart）应力，然后依次迭代。这个模块对于含硬脆性夹层的岩体变形和破损分析是很有用的，因为刚度较大而变形能力差的夹层往往会先行破坏，从而对岩体进一步破损产生影响。

复合材料中层板模型，应该是一个纯复合材料模型，如宏观上为横观各向同性，和我们的模型是有区别的。

yun

可否将编译dll的方法介绍一下呢？

wishes188

好东西，厉害

sy1884

恭喜楼主！

lsxlsxlsx

看不到图片啊

zhangyewei720

楼主半年的努力没有白费！再接再厉啊！