网上flac3d学习笔记.docx

上传人:b****6 文档编号:6273327 上传时间:2023-01-05 格式:DOCX 页数:42 大小:206.24KB
下载 相关 举报
网上flac3d学习笔记.docx_第1页
第1页 / 共42页
网上flac3d学习笔记.docx_第2页
第2页 / 共42页
网上flac3d学习笔记.docx_第3页
第3页 / 共42页
网上flac3d学习笔记.docx_第4页
第4页 / 共42页
网上flac3d学习笔记.docx_第5页
第5页 / 共42页
点击查看更多>>
下载资源
资源描述

网上flac3d学习笔记.docx

《网上flac3d学习笔记.docx》由会员分享,可在线阅读,更多相关《网上flac3d学习笔记.docx(42页珍藏版)》请在冰豆网上搜索。

网上flac3d学习笔记.docx

网上flac3d学习笔记

网上FLAC3D学习笔记

FLAC2d命令的FAQ:

1、怎样开始入门?

答:

初步学习首先看说明书User'sgurid的Section2Gettingstarted和Section3ProblemslovingwithFLAC

2、哪些命令需要掌握?

答:

gridgenpropinterfaceinifixsethiststepsolvepausecontsav等,在学习的过程中遇到这些命令和不会的命令应该查看说明书的CommandReference,里面讲的很清楚,特别需要注意的是一定要仔细看一下CommandReference开头的OptionalRangePhrases,这是讲范围的,因为很多命令的施加都需要加到模型的一部分,这里就告诉你怎样来表示这一部分。

这一原则同样适合与其他5个ITASCA软件。

还有一点就是六个软件(FLACFALC3DUDEC3DECPFCPFC3D)的range大不相同,要注意区分。

3、FLAC和FLAC3D建模有什么区别?

答:

FLAC和FLAC3D基本的建模思想大致是一样的,不同的是FLAC是首先用GRID命令把整体的网格画出来,然后用MODELLNULL命令把网格分成许多部分,再用GEN命令把这些部分一块一块搭接起来。

然后再用GENLINE,GENARC等命令画出要开挖的边界,在FLAC中产生网格的命令只有一个GRID,而且只能用一次。

而FLAC3D是直接用GEN命令一块一块搭接,GEN可以用很多次。

4、怎样查看模型?

答:

plotgrid可以查看网格,plotgridnum可以查看节点号。

5、怎样查看边界条件?

答:

plotfix

6、怎样查看预先留设的开挖边界?

答:

plotgridmark。

7、怎么查看模型的云图?

答:

plsyy(sxxxdisydissig1..)fill,注意,带fill关键字就是查看云图,不带fill关键字就是查看等直线图。

8、怎样查看模型的矢量图?

答:

plgridvel(disp)red

9、怎样查看INTERFACE?

答:

plgridifaceidred

    printifaceid;id是interface的编号

10、FLAC怎么调整视图?

答:

ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小,双击屏幕。

1.请问在圆柱体四周如何施加约束条件?

可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其中r2

命令流如下:

new

genzonecylp1100p2040p3001size123012

genzonereflenorm100

genzonereflenorm001

moss

propden2500bulk2e8she1e8co2e6fric45ten1e6dil10

propftab1ctab2dtab3

table10450.05420.140140

table202e60.051e60.15e515e5

table30100.0530.10

fixxrangeend1100end2140rad1end10.500end20.540rad0.5

回答:

加上的约束不对!

fixxrangeend1100end2140rad1end10.500end20.540rad0.5

回答:

上面的rad是不是相差太大了?

是不是在1左右来限制了

2.关于计算不收敛问题

;初始网格

genzoneradcylp0150-15.75p1add1500p2add060p3add0015.75&size3485dim2.252.252.252.25rat1111.1fill

genzoneradtunsize3435dim0.752.250.752.25rat1111.1& p0150-15.75p1add00-14.25p2add060p3add1500fill

genzonereflectdip90dd270ori150-15.75

;摩尔库仑模型

modelmohr

macroniyan'bulk3.2825e3shear1.5150e3fric20coh0.07'

macromeiceng'bulk2.1735e3shear8.8917e2fric19coh0.05'

propniyanrangez=-13.50

propmeicengrangez=-15.5-13.5

propniyanrangez=-30-15.5

;为开挖定义整体部分

rangenameinz-15.750cylend1150-15.75end2156-15.75rad2.25

rangenamehangdaox12.7517.25y06z-16.50-15.75

;进行加载

fixxyzrangez-30.001-29.999

appsxx=-20rangex-0.0010.001

appsxx=-20rangex29.99930.001

appsyy=-20rangey-0.0010.001

appsyy=-20rangey5.9996.001

appszz=-20rangez-0.0010.001

solve          

(1)

modelnullrangein

modelnullrangehangdao

rangenameinz-15.750cylend1150-15.75end2156-15.75rad2.25

rangenamehangdaoliangce1x12.74912.751y06z-16.50-15.75

rangenamehangdaoliangce2x17.24917.251y06z-16.50-15.75

rangenamehangdaox12.7517.25y06z-16.501-16.499

applynstress-0.05planedip90dd0rangein

;给巷道底部加压力

appnstress-0.05planedip90dd0rangehangdao

;给巷道直墙部分(左侧)加压力

applynstress-0.05planedip90dd0rangehangdaoliangce2

;给巷道直墙部分(右侧)加压力

applynstress-0.05planedip90dd0rangehangdaoliangce1;对点进行监测

histunbal

histgpxdisp12.750-15.75 :

监测左帮

histgpzdisp150-13.5;监测顶板

histgpzdisp150-16.5;监测底板

solve

问:

第一个solve平衡了,是第二个solve不能平衡,这是为什么?

如果我把支护力加大,也不能平衡?

回答:

我遇到的不收敛的原因,有两个:

1)检查你的模型是不是没有建好,相互之间是不是没有连接好;

2)仔细检查你的参数,是不是已经发生了破坏。

回答:

力施加有问题,隧道底部隆起破坏,所以无法计算

问:

我这个模型就是该模型就是要出现巷道底臌现象,如果按照你的意思,力应该怎样施加才对?

回答:

如果无法收敛,应该用step命令,关键问题是要运行多少步。

3.问:

命令流如下:

g1010

me

defdistort

loopi(2,jzones)

loopj(2,jzones)

x(i,j)=x(i,j)+mag*(urand-0.5)

y(i,j)=y(i,j)+mag*(urand-0.5)

end_loop

end_loop

end

propd1000sh1e8bu2e8

其中:

x(i,j)=x(i,j)+mag*(urand-0.5)中的urand是否是随机生成数的命令?

回答:

urand:

randomnumberdrawnfromuniformdistributionbetween 0.0and1.0;

grand:

randomnumberdrawnfromnormaldistribution:

mean=0.0;standarddeviation=1.0

4.问:

怎么能把一个PLOT的图像数据导出来以便用其他软件绘图?

回答:

用setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。

5.问:

用命令建立模型后,如何显示点的坐标?

回答:

使用ploblogrogpnumon命令

6.interface问题

问:

我加了一个层面,初始平衡就出现了这么大的位移?

是怎么原因?

;No1

genzonebrickp0000p11000p20100p3008&

size886rat1.021.021.02

genzonebrickp0009p11009p20109p30011&

size882rat1.021.021.02

genzonereflectdip90dd90origin(0,0,6)rangez011

grouplayer1rangez911

;interface1facerangeplanenorm(0,0,1)origin(0,0,8)dist0.01

inizadd-1.0rangegrouplayer1

;No2

genzonebrickp00011p110011p201011p30013&

size882rat1.021.021.02

genzonereflectdip90dd90origin(0,0,1)rangez1113

grouplayer2rangez1113

;interface2facerangeplanenorm(0,0,1)origin(0,0,10)dist0.01

inizadd-1.0rangegrouplayer2

;No3

genzonebrickp00013p110013p201013p30014&

size881rat1.021.021.02

genzonereflectdip90dd90origin(0,0,1)rangez1314

grouplayer3rangez1314

;interface3facerangeplanenorm(0,0,1)origin(0,0,12)dist0.01

inizadd-1.0rangegrouplayer3

;No4

genzonebrickp00014p110014p201014p30015&

size881rat1.021.021.02

genzonereflectdip90dd90origin(0,0,1)rangez1415

grouplayer4rangez1415

interface4facerangeplanenorm(0,0,1)origin(0,0,13)dist0.01

inizadd-1.0rangegrouplayer4

;No5

genzonebrickp00015p110015p201015p30017&

size882rat1.021.021.02

genzonereflectdip90dd90origin(0,0,15)rangez1517

grouplayer5rangez1517

interface5facerangeplanenorm(0,0,1)origin(0,0,14)dist0.01

inizadd-1.0rangegrouplayer5

;No6

genzonebrickp00017p110017p201017p30018&

size881rat1.021.021.02

genzonereflectdip90dd90origin(0,0,1)rangez1718

grouplayer6rangez1718

interface6facerangeplanenorm(0,0,1)origin(0,0,16)dist0.01

inizadd-1.0rangegrouplayer6

;No7

genzonebrickp00018p110018p201018p30020&

size882rat1.021.021.02

genzonereflectdip90dd90origin(0,0,1)rangez1820

grouplayer7rangez1820

;interface7facerangeplanenorm(0,0,1)origin(0,0,17)dist0.01

inizadd-1.0rangegrouplayer7

;No8

genzonebrickp00020p110020p201020p30024&

size883rat1.021.021.02

genzonereflectdip90dd90origin(0,0,17)rangez2024

grouplayer8rangez2024

;interface8facerangeplanenorm(0,0,1)origin(0,0,19)dist0.01

inizadd-1.0rangegrouplayer8

;No9

genzonebrickp00024p110024p201024p30029&

size884rat1.021.021.02

genzonereflectdip90dd90origin(0,0,0)rangez2429

grouplayer9rangez2429

;interface9facerangeplanenorm(0,0,1)origin(0,0,23)dist0.01

inizadd-1.0rangegrouplayer9

;----------边界条件------------------------

fixxrangex-20.1-19.9

fixxrangex19.920.1

fixyrangey-0.10.1

fixyrangey9.910.1

fixzrangez-0.10.1

;-----------初始条件-------------------------

applyszz-26e6rangez27.928.1

initialszz-1.08e6grad002.6e4rangez028

setgravity00-10

initialdensity2.6e3

initialsxx-45.8e6grad002.6e4rangez028

initialsyy-45.8e6grad002.6e4rangez028

;---------------------------------------------

modelmohr

;-----------材料属性--------------------------

;-----------No1-----------------------

probulk8e9shear7e9fric25coh7.5e6ten2.6e6rangez08

;-----------No2-------------------------

probulk7.8e9shear3.5e9fric30coh6.5e6ten2.4e6rangez810

;------------No3---------------------

probulk7e9shear3.0e9fric22coh5.5e6ten2.2e6rangez1012

;------------No4-----------------------------

probulk6e9shear2.3e9fric26coh5.5e6ten2.0e6rangez1213

;-----------No5-----------------------------------

probulk6.2e9shear2.0e9fric28coh5.2e6ten1.2e6rangez1314

;-------------------No6------------------------

probulk5.8e9shear2.5e9fric35coh6.0e6ten1.3e6rangez1416

;-------------------No7------------------------

probulk9e9shear3.2e9fric30coh10.3e6ten1.2e6rangez1617

;-------------------No8------------------------

probulk8.1e9shear3.1e9fric26coh8.0e6ten1.5e6rangez1719

;-------------------No9------------------------

probulk7.6e9shear2.6e9fric23coh8.5e6ten1.8e6rangez1923

;-------------------No10------------------------

probulk7.5e9shear2.1e9fric26coh8.0e6ten1.9e6rangez2328

;int1propks8e10kn8e10fric30dil6tension2e7bslip=on

;int2propks8e10kn8e10fric30dil6tension2e7bslip=on

;int3propks8e10kn8e10fric35dil8tension2e7bslip=on

int4propks9e11kn9e11fric32dil5tension2e7bslip=on

int5propks9e11kn9e11fric32dil5tension2e7bslip=on

int6propks8e11kn8e11fric33dil5tension2e7bslip=on

;int7propks8e10kn8e10fric33dil5tension2e7bslip=on

;int8propks8e10kn8e10fric33dil5tension2e7bslip=on

;int9propks8e10kn8e10fric33dil5tension2e7bslip=on

histunbal

step5000

applyszz-26e6rangez27.928.1

initialszz-1.08e6grad002.6e4rangez028

setgravity00-10

initialdensity2.6e3

initialsxx-45.8e6grad002.6e4rangez028

initialsyy-45.8e6grad002.6e4rangez028

initialszz-1.08e6grad002.6e4rangez028

initialdensity2.6e3

initialsxx-45.8e6grad002.6e4rangez028

initialsyy-45.8e6grad002.6e4rangez028

这是给的模型加的初始地应力,难道这样加地应力,模型就自动平衡了吗?

回答:

你用ini平衡,为什么还要计算5000步?

问:

我的INI是用来施加初始地应力的,难道用INI那样加应力边界条件就不用step或solve算平衡了吗?

模型就能自动平衡吗?

;-----------初始条件-------------------------

applyszz-26e6rangez27.928.1

initialszz-1.08e6grad002.6e4rangez028

setgravity00-10

initialdensity2.6e3

initialsxx-45.8e6grad002.6e4rangez028

initialsyy-45

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 表格模板 > 合同协议

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1