cumtsmile 发表于 2004-8-15 09:48:00

【原创】一个煤矿巷道变形的例子

new
;100m*100m
gr 120 100
;ROCK
m m j=1,30
m m j=32,43
m m j=45,56
m m j=57,60
m m j=62,65
m m j=67,76
m m j=78,79
m m j=81,84
m m j=86,100
;INTERFACE
m n j=31
m n j=44
m n j=61
m n j=66
m n j=77
m n j=80
m n j=85
        
;40m            
gen   0.00, -32.500.0,30.00 40.00,17.00 40.00,-32.50 i=1,41 j=1,31
gen   0.00, 30.000.0,38.00 40.00,25.00 40.00,17.00 i=1,41 j=32,44
gen   0.00, 38.000.0,44.00 40.00,31.00 40.00,25.00 i=1,41 j=45,57
gen   0.00, 44.000.0,46.00 40.00,33.00 40.00,31.00 i=1,41 j=57,61
gen   0.00, 46.000.0,48.00 40.00,35.00 40.00,33.00 i=1,41 j=62,66
gen   0.00, 48.000.0,55.50 40.00,42.50 40.00,35.00 i=1,41 j=67,77
gen   0.00, 55.500.0,57.30 40.00,44.30 40.00,42.50 i=1,41 j=78,80
gen   0.00, 57.300.0,61.30 40.00,48.30 40.00,44.30 i=1,41 j=81,85
gen   0.00, 61.300.0,67.50 40.00,67.50 40.00,48.30 i=1,41 j=86,101
;20m            
gen   40.00, -32.5040.0,17.00 60.00,10.50 60.00,-32.50 i=41,81 j=1,31
gen   40.00,17.0040.0,25.00 60.00,18.50 60.00,10.50 i=41,81 j=32,44

gen   40.00, 25.0040.0,27.50 47.90,24.7675 47.90,22.4325i=41,57 j=45,50
gen   40.00, 27.5040.0,31.00 47.90,28.4325 47.90,24.7675i=41,57 j=50,57

gen   47.90,22.432547.90,24.7675 52.1,24.7675 52.1,21.0675 i=57,65 j=45,50
gen   47.90,24.767547.90,28.4325 52.1,27.0675 52.1,24.7675 i=57,65 j=50,57

gen   52.1,21.067552.1,24.7675 60,21 60,18.5 i=65,81 j=45,50
gen   52.1,24.767552.1,27.0675 60,24.5 60,21 i=65,81 j=50,57

gen   40.00, 31.0040.0,33.00 60.00,26.50 60.00,24.50 i=41,81 j=57,61
gen   40.00, 33.0040.0,35.00 60.00,28.50 60.00,26.50 i=41,81 j=62,66
gen   40.00, 35.0040.0,42.50 60.00,36.00 60.00,28.50 i=41,81 j=67,77
gen   40.00, 42.5040.0,44.30 60.00,37.80 60.00,36.00 i=41,81 j=78,80
gen   40.00, 44.3040.0,48.30 60.00,41.80 60.00,37.80 i=41,81 j=81,85
gen   40.00, 48.3040.0,67.50 60.00,67.50 60.00,41.80 i=41,81 j=86,101

;40m            
gen   60.00, -32.5060.0,10.50 100.00,-2.50 100.00,-32.50 i=81,121 j=1,31
gen   60.00, 10.5060.0,18.50 100.00,5.50 100.00,-2.50 i=81,121 j=32,44
gen   60.00, 18.5060.0,24.50 100.00,11.50 100.00,5.50 i=81,121 j=45,57
gen   60.00, 24.5060.0,26.50 100.00,13.50 100.00,11.50 i=81,121 j=57,61
gen   60.00, 26.5060.0,28.50 100.00,15.50 100.00,13.50 i=81,121 j=62,66
gen   60.00, 28.5060.0,36.00 100.00,23.00 100.00,15.50 i=81,121 j=67,77
gen   60.00, 36.0060.0,37.80 100.00,24.80 100.00,23.00 i=81,121 j=78,80
gen   60.00, 37.8060.0,41.80 100.00,28.80 100.00,24.80 i=81,121 j=81,85
gen   60.00, 41.8060.0,67.50 100.00,67.50 100.00,28.80 i=81,121 j=86,101

  ;INTERFACE
interface 1 aside from 1,31 to 121,31 bside from 1,32 to 121,32
interface 1 c= 0.2e+6 f=28 kn=588.00e+6 ks=14.70e+6
interface 2 aside from 1,44 to 121,44 bside from 1,45 to 121,45
interface 2 c= 0.2e+6 f=28 kn=588.00e+6 ks= 14.70e+6
;interface 3 aside from 1,57 to 121,57 bside from 1,58 to 121,58
;interface 3 c= 0.05e+6 f=20 kn=480.00e+6 ks= 9.00e+6
interface 4 aside from 1,61 to 121,61 bside from 1,62 to 121,62
interface 4 c= 0.1e+6 f=25 kn=588.00e+6 ks= 14.70e+6
interface 5 aside from 1,66 to 121,66 bside from 1,67 to 121,67
interface 5 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6
interface 6 aside from 1,77 to 121,77 bside from 1,78 to 121,78
interface 6 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6
interface 7 aside from 1,80 to 121,80 bside from 1,81 to 121,81
interface 7 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6
interface 8 aside from 1,85 to 121,85 bside from 1,86 to 121,86
interface 8 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6
  ;BOTTOM
prop D=2500C= 2.50e+6F=24S=3.6e+9B=9.6e+9 j=32,43

;COAL
prop D=1400C= 0.80e+6F=18S=8.06e+8B=1.46e+9 j=45,56

;ROOF
prop D=2500C= 3.00e+6F=28S=4.2e+9B=13.1e+9 j=57,60
prop D=2600C= 6.8e+6F=35S=13.2e+9B=19.4e+9 j=62,65
prop D=2500C= 2.50e+6F=24S=3.6e+9B=9.6e+9 j=67,76
prop D=2600C= 5.00e+6F=30S=7.8e+9B=16.5e+9 j=78,79
prop D=2500C= 3.50e+6F=30S=4.2e+9B=13.1e+9 j=81,84
prop D=2550C= 5.00e+6F=32S=7.8e+9B=16.5e+9 j=86,100
set g=9.81

ini syy=-25.00e6var 0,2.5e6
ini sxx=-30.00e6var 0,3.0e6
ini szz=-30.00e6var 0,3.0e6
app syy=-22.50e+6 j=101
fix x i=1
fix x i=121
fix y j=1
set large
his nstep=100
his unbal
pa

solveforce 1e3
save pdsyy.sav

cal pdsold.txt
title
old support technology
ini yd=0
ini xd=0
ini yv=0
ini xv=0
ini x=47.9 i=57
ini x=52.1 i=65
his xd i=57 j=53
his xd i=65 j=53
his yd i=61 j=57
his yd i=61 j=50

M N i=57,64 j=50,56

INTERFACE 9 ASIDE FROM 57,50 TO 65,50 BSIDE FROM 65,50 TO 65,57
INTERFACE 10 ASIDE FROM 65,50 TO 65,57 BSIDE FROM 65,57 TO 57,57
INTERFACE 11 ASIDE FROM 65,57 TO 57,57 BSIDE FROM 57,57 TO 57,50
INTERFACE 12 ASIDE FROM 57,57 TO 57,50 BSIDE FROM 57,50 TO 65,50
interface 9 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 10 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 11 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 12 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6

;maogan

;left
;beam (tuo pan)
STR BEAM BEG 47.90,25.15 END 47.90,25.25 SEG=2 PROP=1
INTERFACE 21 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 1 TO NODE 3
INTERFACE 21 c=0.7e6 f=12 ks=1e6 kn=3e7glued

STR BEAM BEG 47.90,25.9 END 47.90,26.0 SEG=2 PROP=1
INTERFACE 22 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 4 TO NODE 6
INTERFACE 22 c=0.7e6 f=12 ks=1e6 kn=3e7 glued

STR BEAM BEG 47.90,26.65 END 47.90,26.75 SEG=2 PROP=1
INTERFACE 23 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 7 TO NODE 9
INTERFACE 23 c=0.7e6 f=12 ks=1e6 kn=3e7 glued

STR BEAM BEG 47.90,27.4 END 47.90,27.5 SEG=2 PROP=1
INTERFACE 24 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 10 TO NODE 12
INTERFACE 24 c=0.7e6 f=12 ks=1e6 kn=3e7glued

STR BEAM BEG 47.90,28.15 END 47.90,28.25 SEG=2 PROP=1
INTERFACE 25 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 13 TO NODE 15
INTERFACE 25 c=0.7e6 f=12 ks=1e6 kn=3e7glued

;cable

STR CAB BEG node 2 END 47.8034,25.1741SEG=1 PROP=2
STR CAB BEG NODE 16 END 46.9341,24.9412SEG=4PROP=3
STR CAB BEG NODE 20 END 45.9681,24.6824 SEG=5 PROP=4 ten 10000

STR CAB BEG node 5 END 47.80,25.95SEG=1 PROP=2
STR CAB BEG NODE 26 END 46.90,25.95SEG=4PROP=3
STR CAB BEG NODE 30 END 45.90,25.95 SEG=5 PROP=4 ten 10000

STR CAB BEG node 8 END 47.80,26.7SEG=1 PROP=2
STR CAB BEG NODE 36 END 46.90,26.7SEG=4PROP=3
STR CAB BEG NODE 40 END 45.90,26.7 SEG=5 PROP=4 ten 10000

STR CAB BEG node 11 END 47.80,27.45SEG=1 PROP=2
STR CAB BEG NODE 46 END 46.90,27.45SEG=4PROP=3
STR CAB BEG NODE 50 END 45.90,27.45 SEG=5 PROP=4 ten 10000

STR CAB BEG node 14 END 47.8034,28.2259SEG=1 PROP=2
STR CAB BEG NODE 56 END 46.9341,28.4588SEG=4PROP=3
STR CAB BEG NODE 60 END 45.9681,28.7176SEG=5 PROP=4 ten 10000

;ROOF
;(gangjinti)
STR BEAM BEG 47.90,28.4325 END 52.10,27.0675 SEG=10PROP=5
INTERFACE 26 ASIDE FROM 57,57 to 65,57 BSIDE FROM NODE 66 TO NODE 76
INTERFACE 26 c=2e6 f=22 ks=1e7 kn=3e8 glued

;beam (tuo pan)
STR BEAM BEG 48.00,28.4 END 48.10,28.3675 SEG=2 PROP=1
INTERFACE 27 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 77 TO NODE 79
INTERFACE 27 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 48.65,28.1887 END 48.75,28.1563 SEG=2 PROP=1
INTERFACE 28 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 80 TO NODE 82
INTERFACE 28 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 49.30,27.9775 END 49.40,27.9450 SEG=2 PROP=1
INTERFACE 29 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 83 TO NODE 85
INTERFACE 29 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 49.95,27.7662 END 50.05,27.7337 SEG=2 PROP=1
INTERFACE 30 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 86 TO NODE 88
INTERFACE 30 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 50.60,27.5550 END 50.70,27.5225 SEG=2 PROP=1
INTERFACE 31 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 89 TO NODE 91
INTERFACE 31 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 51.25,27.3437 END 51.35,27.3113 SEG=2 PROP=1
INTERFACE 32 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 92 TO NODE 94
INTERFACE 32 c=2e6 f=22 ks=1e7 kn=3e8 glued

STR BEAM BEG 51.90,27.1325 END 52.00,27.10 SEG=2 PROP=1
INTERFACE 33 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 95 TO NODE 97
INTERFACE 33 c=2e6 f=22 ks=1e7 kn=3e8 glued

  ;cable
STR CAB BEG node 78 END 48.05,28.4889 SEG=1 PROP=6
STR CAB BEG NODE 98 END 48.05,29.0889 SEG=2PROP=7
STR CAB BEG NODE 100 END 48.05,30.5889 SEG=7 PROP=8 ten 20000

STR CAB BEG node 81 END 48.7309,28.2677SEG=1 PROP=6
STR CAB BEG NODE 108 END 48.9163,28.8383 SEG=2   PROP=7
STR CAB BEG NoDE 110 end 49.3798,30.2648 SEG=7 PROP=8 ten 20000

STR CAB BEG node 84 END 49.3809,28.0565SEG=1 PROP=6
STR CAB BEG NODE 118 end 49.5663,28.6271 SEG=2   PROP=7
STR CAB BEG NoDE 120 END 50.0298,30.0536 SEG=7 PROP=8 ten 20000

STR CAB BEG node 87 END 50.0310,27.8453SEG=1 PROP=6
STR CAB BEG NODE 128 END 50.2164,28.4159 SEG=2 PROP=7
STR CAB BEG NoDE 130 END 50.6798,29.8423 SEG=7 PROP=8 ten 20000

STR CAB BEG node 90 END 50.6810,27.6341SEG=1 PROP=6
STR CAB BEG NODE 138 END 50.8664,28.2047 SEG=2   PROP=7
STR CAB BEG NoDE 140 end 51.3298,29.6311 SEG=7 PROP=8 ten 20000

STR CAB BEG node 93 END 51.3310,27.4229SEG=1 PROP=6
STR CAB BEG NODE 148 end 51.5164,27.9935 SEG=2   PROP=7
STR CAB BEG NoDE 150 END 51.9798,29.4198 SEG=7 PROP=8 ten 20000

STR CAB BEG node 96 END 52.0066,27.20333SEG=1 PROP=6
STR CAB BEG NODE 158 END 52.3449,27.7243 SEG=2 PROP=7
STR CAB BEG NoDE 160 END 53.1482,28.9613 SEG=7 PROP=8 ten 20000

;Right
;beam (tuo pan)
STR BEAM BEG 52.10,25.3175 END 52.10,25.4175 SEG=2 PROP=1
INTERFACE 34 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 168 TO NODE 170
INTERFACE 34 c=0.7e6 f=12 ks=1e6 kn=3e7glued

STR BEAM BEG 52.10,26.0675 END 52.10,26.1675 SEG=2 PROP=1
INTERFACE 35 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 171 TO NODE 173
INTERFACE 35 c=0.7e6 f=12 ks=1e6 kn=3e7 glued

STR BEAM BEG 52.10,26.8175 END 52.10,26.9175 SEG=2 PROP=1
INTERFACE 36 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 174 TO NODE 176
INTERFACE 36 c=0.7e6 f=12 ks=1e6 kn=3e7 glued

  ;cable
STR CAB BEG node 169 END 52.20,25.3407 SEG=1 PROP=2
STR CAB BEG NODE 177 END 53.10,25.0996 SEG=4   PROP=3
STR CAB BEG NODE 181 END 54.0319,24.8499 SEG=5 PROP=4 ten 10000

STR CAB BEG node 172 END 52.20,26.1175SEG=1 PROP=2
STR CAB BEG NODE 187 END 53.10,26.1175 SEG=4PROP=3
STR CAB BEG NoDE 191 end 54.10,26.1175 SEG=5PROP=4 ten 10000

STR CAB BEG node 175 END 52.20,26.8943SEG=1 PROP=2
STR CAB BEG NODE 197 end 53.10,27.1354 SEG=4   PROP=3
STR CAB BEG NoDE 201 END 54.0319,27.3851 SEG=3 PROP=4 ten 10000

;tuo pan prop=1(steel size:100*100*10mm)
STR PROP=1 E=2.0E11 I=1.25e-9 area=5.5e-4d=7800.00 wid=0.125
;cable prop (20mm*2m)
STR PROP=2 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=2 d=7800.0 kb=10e10 ks=9e7 yield=0.60e5 ycom=3e8
STR PROP=3 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=3 d=7800.0 kb=6e2 ks=7e3yield=0.6e5 ycom=3e8
STR PROP=4 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=4 d=7800.0 kb=10.0e9 ks=9.0e6 yield=0.6e5 ycom=3e8
;gangjinti prop
STR PROP=5 E=2.0E11 i=5.2e-9 area=2.0e-4 d=7800.00wid=0.25

STR PROP=6 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=6 d=7800.0 kb=10e10 ks=9e7 yield=0.60e5 ycom=3e8
STR PROP=7 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=7 d=7800.0 kb=6e2 ks=7e3yield=0.6e5 ycom=3e8
STR PROP=8 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=8 d=7800.0 kb=10.0e9 ks=9.0e6 yield=0.6e5 ycom=3e8

step 11756
save pdsold.sav

qianghuicm 发表于 2004-8-25 13:03:03

landslide 发表于 2004-8-26 09:51:09

回复: 【原创】一个煤矿巷道变形的例子

不错,支持原创

starsmoon 发表于 2004-8-26 21:12:28

回复: 【原创】一个煤矿巷道变形的例子

是两维的东西吧,大家用的是什么版本啊?

cjiao2000 发表于 2004-8-28 11:06:54

回复: 回复: 【原创】一个煤矿巷道变形的例子

starsmoon wrote:
是两维的东西吧,大家用的是什么版本啊?
I have 3.3 version,do you need?!

rocksoul 发表于 2004-8-28 19:51:01

starsmoon 发表于 2004-8-28 23:44:42

回复: 【原创】一个煤矿巷道变形的例子

cjiao2000 兄,is your version "education version"? if not ,how can i connect with you?
My mail: starsmoon@citiz.net
msn:YingChee982597@hotmail.com
popo:starsmoon
在隧道方面你可是前辈啊,呵呵。

q1 发表于 2004-9-6 08:07:56

liujunfeng 发表于 2005-6-6 18:01:28

tg0215 发表于 2005-7-5 22:29:13

Re:【原创】一个煤矿巷道变形的例子

我有4.0版的。

lc_rw 发表于 2005-7-19 10:44:55

Re:【原创】一个煤矿巷道变形的例子

碰到一个同行前辈!

yuanyuekafu 发表于 2012-4-12 18:22:56

强大啊,有没有模拟的效果图呢!

快乐 发表于 2018-3-9 20:18:35

厉害呀!
页: [1]
查看完整版本: 【原创】一个煤矿巷道变形的例子