如何在衬砌和围岩之间加一层壳体结构？

Gracezhong_8

gen zon radcyl
p0 0 0 0 p1 0 0 -9
p2 0 1 0 p3 10.55 0 0 &
size 5 2 10 4 dim 4 5.55 4 5.55 rat 1 1 1 1.2 group outsiderockxy
gen zone cshell p0 0 0 0 p1 0 0 -4 p2 0 1 0 p3 5.55 0 0 &
size 1 2 10 4 dim 3.5 5.05 3.5 5.05 rat 1 1 1 1 group concretlinerxy fill group insiderockxy


我写的命令老不能把下部围岩上面也加上壳结构，究竟是哪里出了问题，还请高手指点，谢谢！

sel shell id 1 group outsiderockxy &

range cyl end1 0 0 0 end2 0 1 0 rad 5.6 &

cyl end1 0 0 0 end2 0 1 0 rad 3.9 not &

z -3.51 0.01 x -5.56 0.01

haiqsim

是你最后的z 和x 范围设置不正确！
首先，你肯定没有搞清楚模型的坐标原点在什么位置。显示网格组后，应添加显示坐标轴，这样才好设定范围。

你那个z的范围，应该比 —4 更靠下才对啊，下部围岩左侧最高点的z坐标才是-4，你得把这个包括进去；

x的范围就更离谱了，完全反了。模型的坐标系原点在左上角，你的x范围应为：(-0.01, 5.56)。

把这两点都改了，是可以生成壳的。显示一下看看，但是，生成的壳单元有点多，说明这个圆圈以及x、y范围还是限制的不够，你自己想想看是什么原因。从那个浅埋开挖支护的例子里可以找到类似的吧？

Gracezhong_8

谢谢你的指点，我把z坐标改过来了，不过x坐标应该是-5.56，0.01，因为那个坐标原点是在右上角，我只开挖了隧道的左下部分进行支护，这个例子就是你推荐的那本李围的书中第六章的例子。他里面有关这部分加壳结构的名利根本就不能用。
想请问一下，这两个圆圈的范围我试了很多次，总还是不对。我的想法是：外圈的半径必须比5.55大，而内圈的半径要比4.0小，不过怎么试外圈都会把另一部网格包括进去，还请您指点一下，谢谢！

haiqsim

1.  你这段命令流建立的模型，坐标原点明明在左上啊，难道我眼花？晕
2.  这个范围限制么，你目前只考虑半径大小和x方向、z方向限制，难道没考虑下y方向也给限制一下？？

Gracezhong_8

多谢指点，那我再考虑一下y方向试试看！

Gracezhong_8

谢谢你的指点，我又参考了一下那个浅层开挖的例子，发现我的圆心坐标取的有点问题，而且我当时给你的图像圆心确实是在左上角，不好意思，我自己图没给完全，其实我是在开挖隧道的左下部分，经过圆心的调整以及多次调整半径，终于形成了完整的壳体结构。
sel shell id 1 group outsiderockxy &
range cyl end1 0 0 1.55 end2 0 1 1.55 rad 6.0 &
cyl end1 0 0 1.55 end2 0 1 1.55 rad 5.2 not &
z -4.01 0.01 x -5.56 0.01

另外，我把你推荐的那本李围的书上关于隧道开挖的例子看完了，发现命令流存在一些错误，不过都是些小错误，自己能够看懂解决，就是关于这个壳结构的施加上他两个隧道例子里面都没有讲明白，我不知是否是他疏忽了。这几个隧道的例子在模拟过程中都做了很大的简化，看起来通俗易懂，而且三个例子大同小异，手册上的例子着实比他复杂了许多。看了这几个例子，对于模拟隧道问题信心大增。相信在实际工程方面做模拟也是比较常规的那种吧。

这本书上关于初始条件及初始应力条件时，并没有考虑到隧道埋深的作用，所以他并没有涉及隧道中地应力的施加，想请问一下，这会不会影响后面的计算。谢谢！

另外，我还有个问题想咨询一下，请问你有没有关于在隧道开挖中碰到地下水方面的资料，我看的那些隧道例子中好像都没有涉及到这方面的问题，我想知道在碰到这方面的问题时又是怎样的，应该怎样解决，谢谢了！

haiqsim

1. 很惭愧，我也只是个门外汉，您说的这个“隧道埋深”的作用是指的什么？我不大明白。还有“隧道中地应力的增加”是指什么？我也不大清楚。想必他的书中在计算的时候就是单纯考虑了自重应力场，这个自重应力场不就是随埋深增大而增大么？您说的莫非是指构造应力？
2. 很抱歉，我没有做过类似的工程。论坛上应该有不少流固耦合的资料与讨论吧？可以搜索看看。另外，上期刊网搜索些文章，应该会有不少相关的。

Gracezhong_8

哦，我所指的地应力就是围岩内部的应力，常常用initial命令施加的szz，syy和sxx，我在论坛上的一些帖子上也看到过有类似的内容，因为在那本隧道例子的书中，他并没有指明隧道的埋深，因此我们也不能用initial命令来施加自重应力以及侧向应力，所以我想问一下看看。不过还是要谢谢你，我再看看论坛上的其他帖子或是查查相关的文章吧。谢谢了！

haiqsim

我想有一点你可能误会了吧？？
那个自重应力、侧向应力等等，你都是指的初始应力吧？
这个初始应力并不是必须由initial命令来施加的吧？？
我的理解是：
用initial命令来施加初始应力的情况，一般是为了让模型尽早达到平衡状态，先给他施加一个梯度变化的初始应力，然后再进行计算，这个是认为施加的，为加速平衡用的，但好像也不是随便加的，我没有用过，版上好像有相关讨论，可以查看。
有不对的地方请指正。

hwaust

1# Gracezhong_8
李围书中好多错误，我刚刚开始学习FLAC，好多不会修改，兄弟能否把修改好的李围书中的命令流发给我一份，谢谢，hwaust@163.com