simwe论坛关于flac+ansys.docx
- 文档编号:11351469
- 上传时间:2023-02-28
- 格式:DOCX
- 页数:47
- 大小:75.37KB
simwe论坛关于flac+ansys.docx
《simwe论坛关于flac+ansys.docx》由会员分享,可在线阅读,更多相关《simwe论坛关于flac+ansys.docx(47页珍藏版)》请在冰豆网上搜索。
simwe论坛关于flac+ansys
关于surfer绘制等直线,由于需要结合自己实际需要绘制,sufer只能绘制2D等直线,所以首先确定自己绘制等直线的面,然后flac3d导出的单元信息(坐标及应力,应变),通过简单的fortran程序找出所要求的面上的单元的,并记录其坐标位置,即应力应变值。
输出格式为surfer认可的格式即可以了。
授人以鱼,不如授人以渔。
我编的程序只能针对自己的问题,所以大家不要关心程序,如果,确实需要,按照上述思路实现应该很容易。
建立模型以后直接solve让程序产生应力分布!
虽然时间花费上多一点,可是不用那么伤脑筋了。
我得懒人做法是:
1)现把模量×10,c,Fi取很大,然后solveforbalance;2)inidisp=0,inistate=0,;3)回复原来得E,G,C,Fi等,
subzone和zone哪个是“最小几何单位”?
二者有和区别呢?
我们通常说的“单元element”指的是哪一个?
用户接触到的最小几何单位应该是zone
通常说的element也是指zone
subzone和FLAC3d的MixedDiscretization算法有关
比如说一个brick
计算时又自动划分为2套overlay
每套overlay中各有5个tetra
即subzone
据说这样可以更好地模拟材料的塑性变形
这是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中找出每一步开挖的单元号信息,格式如上
第四步,执行转换程序即可
材料分组和开挖单元分组由转换程序自动完成。
该程序绝对能用,如果出现问题,检查输入文件是否正确
在FLAC现有的模型,包括摩尔库伦和应力应变软化模型,运用到隧道等开挖支护的时候,并不能考虑到由于岩层土体应力变化对其材料参数的影响,很多情况下,没有锚杆支护与使用锚杆支护的效果并不明显,与实际出入较大.这是模型的缺陷.如果中间调整土体岩层材料的参数后,效果要明显许多,但参数的重新选取只能是经验来决定,误差较大.另外,锚杆支护时有个初锚力,然后再逐渐加载到一定值.模拟时,可以按此分步实现,但发现效果也并不明显.如果用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
;40m
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
;INTERFACE
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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- simwe 论坛 关于 flac ansys