simwe论坛关于flac+ansysWord文档下载推荐.docx

simwe论坛关于flac+ansysWord文档下载推荐.docx

《simwe论坛关于flac+ansysWord文档下载推荐.docx》由会员分享，可在线阅读，更多相关《simwe论坛关于flac+ansysWord文档下载推荐.docx（47页珍藏版）》请在冰点文库上搜索。

这是flac3d为了避免求解中单元出现“沙漏”现象而采取的一种单元离散技术

Q:

genzonebricksize666

modelelas

propbulk1e8shear7e7

fixxrangex-0.10.1

table10,0100,-1e5/我不知道table在这里什么意思,有什么作用

applysxx1.0histtable1rangex5.9,6.1y0,6z0,2

histzonesxx6,0,0

step100

我对这个table不是很清楚

genzonecylp0000p1100p2020p3001size454

genzonereflectnorm1,0,0

genzonereflectnorm0,0,1

modelss

propbulk1.19e10shear1.1e10

propcoh2.72e5fric44ten2e5

propttab1

table10,2e52e-5,0

请问这里的table和后面的数字到底是什么意思啊?

A:

table相当于数已经赋值的数组

请问selcableid=1begin=（0.0,0.5,3.5）end=（8,0.5,3.5）nseg=8中的nseg是什么意思?

把cable分为几段,比如nseg=8,就是分为八段.

genzonereflectdip90dd270origin0,0,0

应该是以90和270度还有原点组成的平面进行对称拷贝吧？

？

inistate＝0

初始化单元为弹性。

一般来说，mc材料在初始化应力的时候，由于是手动设置的k0，因此局部单元可能已经进入了塑性状态，这会影响后继的计算步骤，因此需要inistate＝0

我这里编写了一个ansystoflac3d的前处理小程序，用起来比较方便。

希望能和大家共享和交流。

（可以自动形成flac开挖信息）

ansys-flac2002转换程序使用说明:

本程序用于将ansys生成的数据格式转换为flac计算所需的数据格式。

程序需要输入5个文件名。

1。

单元信息文件，格式如下：

总单元数

单元信息。

。

，例如：

10（总单元数）

11110181238713881688128122882688318827

21110181228826883188278127868777918456

31110188128816881788138827883188328828

41110188278831883288288456779177908455

51110188138817881888148828883288338829

61110188288832883388298455779077898454

71110188148818881988158829883388348830

81110188298833883488308454778977888453

91110188158819762976328830883476997704

101110188308834769977048453778877007703

2。

节点信息

节点总数

节点信息。

例如：

20（节点总数）

1210.750000000228.398805249-327.458594084

2210.750000000227.440685854-299.440511342

3210.750000000226.118354064-279.017920414

4241.820265858240.000000000-331.304948353

5246.293378743240.000000000-331.239755597

6252.148198534240.000000000-331.154425346

7258.555372995240.000000000-331.061044881

8266.830116362240.000000000-330.940445794

9274.869142734240.000000000-330.823282132

10284.683736459240.000000000-330.680240712

11295.556691364240.000000000-330.521774355

12306.689058734240.000000000-330.359527228

13318.078795786240.000000000-330.193529103

14329.460877328240.000000000-330.027642551

15340.835169125240.000000000-329.861869531

16352.899895166240.000000000-329.686033873

17367.026374896240.000000000-329.480149475

18380.592277707240.000000000-329.282435127

19390.569611461240.000000000-329.137021875

20242.519404127240.000000000-307.394667740

3。

输入开挖信息

开挖步数

第一步开挖单元数

开挖单元信息

第二步开挖单元数

例如

2（两步开挖）

5（第一步开挖5个单元）

1252751101176391774084308221178911786984578689

1252851101178911786984578689178941787281278126

12529511011763817739177401763917890178681786917891

12530511011789017868178691789117893178711787217894

12531511011763717738177391763817889178671786817890

15（第二步开挖15个单元）

12532511011788917867178681789017892178701787117893

12533511011764117742177381763717887178731786717889

12534511011788717873178671788917888178741787017892

12535511011764317730177421764117885178751787317887

12536511011788517875178731788717886178761787417888

12537511011764517732177301764317883178771787517885

12538511011788317877178751788517884178781787617886

12539511011764717734177321764517881178791787717883

12540511011788117879178771788317882178801787817884

12541511011763317697177341764717832178361787917881

12542511011783217836178791788117831178351788017882

1254351101838382831734917325841086851755917537

1254451101841086851755917537841182071756217540

12545511011732517349173481732417537175591755817536

12546511011753717559175581753617540175621756117539

4。

输入计算结果文件名：

该文件存放各种几何信息，包括材料分组信息

5。

输入开挖信息文件名：

该文件存放开挖分组信息。

所有材料的分组和开挖分组（group）都是自动定义的在ansys里将要分步开挖的单元分别找出来，按使用说明输入即可对不同开挖步信息分别给出它的group名

本code最大单元和结点数量可达1000000个。

转换后的文件可能是大了点，不过很好用，转换后的文件直接call就行了。

将开挖信息列出来,是为了给它们一个groupname,以便再开挖计算时使用。

在ansys里列开挖单元是不难的。

第一步，在ansys中建立三维模型，并给每个单元赋予相应的材料编号

第二步，将ansys中的单元信息和结点信息导出，格式如上述例子

第三步，在ansys中找出每一步开挖的单元号信息，格式如上

第四步，执行转换程序即可

材料分组和开挖单元分组由转换程序自动完成。

该程序绝对能用，如果出现问题，检查输入文件是否正确

在ＦＬＡＣ现有的模型，包括摩尔库伦和应力应变软化模型，运用到隧道等开挖支护的时候，并不能考虑到由于岩层土体应力变化对其材料参数的影响，很多情况下，没有锚杆支护与使用锚杆支护的效果并不明显，与实际出入较大．这是模型的缺陷．如果中间调整土体岩层材料的参数后，效果要明显许多，但参数的重新选取只能是经验来决定，误差较大．另外，锚杆支护时有个初锚力，然后再逐渐加载到一定值．模拟时，可以按此分步实现，但发现效果也并不明显．如果用ANSYS模拟锚杆支护，只能对杆端起到作用，不能准确模拟整个杆的受力状况．具体可以查阅中国矿业大学出版社＜锚杆支护的数值模拟＞．

一个煤矿巷道的例子：

new

;

100m*100m

gr120100

ROCK

mmj=1,30

mmj=32,43

mmj=45,56

mmj=57,60

mmj=62,65

mmj=67,76

mmj=78,79

mmj=81,84

mmj=86,100

INTERFACE

mnj=31

mnj=44

mnj=61

mnj=66

mnj=77

mnj=80

mnj=85

40m

gen0.00,-32.500.0,30.0040.00,17.0040.00,-32.50i=1,41j=1,31

gen0.00,30.000.0,38.0040.00,25.0040.00,17.00i=1,41j=32,44

gen0.00,38.000.0,44.0040.00,31.0040.00,25.00i=1,41j=45,57

gen0.00,44.000.0,46.0040.00,33.0040.00,31.00i=1,41j=57,61

gen0.00,46.000.0,48.0040.00,35.0040.00,33.00i=1,41j=62,66

gen0.00,48.000.0,55.5040.00,42.5040.00,35.00i=1,41j=67,77

gen0.00,55.500.0,57.3040.00,44.3040.00,42.50i=1,41j=78,80

gen0.00,57.300.0,61.3040.00,48.3040.00,44.30i=1,41j=81,85

gen0.00,61.300.0,67.5040.00,67.5040.00,48.30i=1,41j=86,101

20m

gen40.00,-32.5040.0,17.0060.00,10.5060.00,-32.50i=41,81j=1,31

gen40.00,17.0040.0,25.0060.00,18.5060.00,10.50i=41,81j=32,44

gen40.00,25.0040.0,27.5047.90,24.767547.90,22.4325i=41,57j=45,50

gen40.00,27.5040.0,31.0047.90,28.432547.90,24.7675i=41,57j=50,57

gen47.90,22.432547.90,24.767552.1,24.767552.1,21.0675i=57,65j=45,50

gen47.90,24.767547.90,28.432552.1,27.067552.1,24.7675i=57,65j=50,57

gen52.1,21.067552.1,24.767560,2160,18.5i=65,81j=45,50

gen52.1,24.767552.1,27.067560,24.560,21i=65,81j=50,57

gen40.00,31.0040.0,33.0060.00,26.5060.00,24.50i=41,81j=57,61

gen40.00,33.0040.0,35.0060.00,28.5060.00,26.50i=41,81j=62,66

gen40.00,35.0040.0,42.5060.00,36.0060.00,28.50i=41,81j=67,77

gen40.00,42.5040.0,44.3060.00,37.8060.00,36.00i=41,81j=78,80

gen40.00,44.3040.0,48.3060.00,41.8060.00,37.80i=41,81j=81,85

gen40.00,48.3040.0,67.5060.00,67.5060.00,41.80i=41,81j=86,101

gen60.00,-32.5060.0,10.50100.00,-2.50100.00,-32.50i=81,121j=1,31

gen60.00,10.5060.0,18.50100.00,5.50100.00,-2.50i=81,121j=32,44

gen60.00,18.5060.0,24.50100.00,11.50100.00,5.50i=81,121j=45,57

gen60.00,24.5060.0,26.50100.00,13.50100.00,11.50i=81,121j=57,61

gen60.00,26.5060.0,28.50100.00,15.50100.00,13.50i=81,121j=62,66

gen60.00,28.5060.0,36.00100.00,23.00100.00,15.50i=81,121j=67,77

gen60.00,36.0060.0,37.80100.00,24.80100.00,23.00i=81,121j=78,80

gen60.00,37.8060.0,41.80100.00,28.80100.00,24.80i=81,121j=81,85

gen60.00,41.8060.0,67.50100.00,67.50100.00,28.80i=81,121j=86,101

interface1asidefrom1,31to121,31bsidefrom1,32to121,32

interface1c=0.2e+6f=28kn=588.00e+6ks=14.70e+6

interface2asidefrom1,44to121,44bsidefrom1,45to121,45

interface2c=0.2e+6f=28kn=588.00e+6ks=14.70e+6

interface3asidefrom1,57to121,57bsidefrom1,58to121,58

interface3c=0.05e+6f=20kn=480.00e+6ks=9.00e+6

interface4asidefrom1,61to121,61bsidefrom1,62to121,62

interface4c=0.1e+6f=25kn=588.00e+6ks=14.70e+6

interface5asidefrom1,66to121,66bsidefrom1,67to121,67

interface5c=0.1e+6f=25kn=588.00e+6ks=14.70e+6

interface6asidefrom1,77to121,77bsidefrom1,78to121,78

interface6c=0.1e+6f=25kn=588.00e+6ks=14.70e+6

interface7asidefrom1,80to121,80bsidefrom1,81to121,81

interface7c=0.1e+6f=25kn=588.00e+6ks=14.70e+6

interface8asidefrom1,85to121,85bsidefrom1,86to121,86

interface8c=0.1e+6f=25kn=588.00e+6ks=14.70e+6

BOTTOM

propD=2500C=2.50e+6F=24S=3.6e+9B=9.6e+9j=32,43

COAL

propD=1400C=0.80e+6F=18S=8.06e+8B=1.46e+9j=45,56

ROOF

propD=2500C=3.00e+6F=28S=4.2e+9B=13.1e+9j=57,60

propD=2600C=6.8e+6F=35S=13.2e+9B=19.4e+9j=62,65

propD=2500C=2.50e+6F=24S=3.6e+9B=9.6e+9j=67,76

propD=2600C=5.00e+6F=30S=7.8e+9B=16.5e+9j=78,79

propD=2500C=3.50e+6F=30S=4.2e+9B=13.1e+9j=81,84

propD=2550C=5.00e+6F=32S=7.8e+9B=16.5e+9j=86,100

setg=9.81

inisyy=-25.00e6var0,2.5e6

inisxx=-30.00e6var0,3.0e6

iniszz=-30.00e6var0,3.0e6

appsyy=-22.50e+6j=101

fixxi=1

fixxi=121

fixyj=1

setlarge

hisnstep=100

hisunbal

pa

solveforce1e3

savepdsyy.sav

calpdsold.txt

title

oldsupporttechnology

iniyd=0

inixd=0

iniyv=0

inixv=0

inix