- 积分
- 0
- 注册时间
- 2012-2-22
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2012-5-29 22:42:57
|
显示全部楼层
来自 湖南湘潭
本帖最后由 struggle_ycc 于 2012-5-29 22:46 编辑
最近研究大神的这个帖子,终于搞出来了,看了很多次楼主附件和跟帖,搞到最后去看陈育民老师那本书就有这个例子,大神就是大神,膜拜。在大神面前不敢造次,下面我来说下学习体会,说的不对大家见谅。
首先说关于z1/z2...多层的问题,其实这个只要控制不同层的x、y的范围,详细点就是把每层的x范围和y的范围控制成一样,如果是CAD地形图就直接把地形图范围控制成一样就可以了,然后在surfer插值的时候把x和y方向坐标的间隔都调成一样,这样就OK了。
再说模型建立出现的问题。第一次把surfer_Data换成我建模的坐标并把command文件改完之后,程序运行就会提示单元体不符合节点排列规则。一直从01-03和生产的table,Flac3D_Model都认真研究了一遍,最后发现每个table的x值都只到70,这里明显的少了一个值(72),原因是02_Creat_Table.fis里面table的第二层循环是从1-col_Num,只有col_Num个值,应该是col_Num+1才对,所以这个地方是有问题的,这里的问题会直接影响到03_Gen_FlacModel.fis这个文件在引用这些数据建立三菱柱。
再说03_Gen_FlacModel.fis这个文件,其思路就是通过两个循环不断累加直立的wedge单元体建两个不同的group。其中group soil的zr1-zr4的table对应的表号为j+n_zon_row(+1),而这里的n_zon_row = Row_Num - 1,这里n_zon_row其实是x方向上坐标个数,要比x方向上的单元体多1个,
zr1=table(j+n_zon_row,x1)
zr2=table(j+n_zon_row+1,x1)
zr3=table(j+n_zon_row+1,x2)
zr4=table(j+n_zon_row,x2)
而循环式从1开始,这样造成table的表好小1,就是说zr1=table(j+n_zon_row,x1)=table(26,x1),应该是table(27,x1)这样,所有单元体就建立不起来。group rock也是同样的道理,楼主之所有能够建立起来就是因为在每个table本来就少一个数值(72),因此是OK的。到这里模型没问题了,因为我只要flac的,其他两个文件没有去研究。
这个帖子虽然很久之前就发了,都被陈老师收录到书里面去了,之所以花那么多时间回这个帖子是想希望和我一样想学习这个方法建模的童鞋少走弯路,少花费时间。写得很凌乱,大家将就看下,祝大家和我一样好运,谢谢大神。
|
|