相似模拟试验和数值模拟在岩土工程中的应用及实际案例

相似模拟与模型试验在岩土工程中的应用

相似模拟与其它一样是社会生产发展的必然产物。由于社会生产的不断发展，岩土工程所提出的问题日益复杂和繁琐。用数学方法很难得到精确的解析解，只能作一些假设与简化再求解，因而带来一些误差。于是人们不得不通过实验的方法来探求那些靠数学方法无法研究的复杂现象的规律性。但是直接的实验的方法有很大的局限性，其实验的结果只能推广到与实验条件完全相同的实际问题中去，这种实验方法常常只能得出个别量的表面规律性关系，难以抓住现象的内在本质。《相似模拟》正是为解决这些问题而产生的，它不直接的研究自然现象或过程的本身，而是研究与这些自然现象或过程相似的模型， 它是理论与实际密切相结合的科学研究方法，是解决一些比较复杂的生产工程问题的一种有效方法。

一、相似模拟与模型试验的方要研究内容

它是研究自然界相似现象的一门科学。它提供了相似判断的方法。并用于指导模型试验， 整理试验结果，并把试验结果用于原型的理论基础。

二、相似常数

设c 表示相似常数，x 表示原型中的物理量，x ' 表示模型中的物理量，则：

i i

i x x c '

=

其中i c 表示第i 个物理量所对应的相似常数。

物理量包含于现象之中。而表示现象的物理量，一般都不是孤立的，互不关联的，而是 处在自然规律所决定的一定关系中，所以说各种相似常数之间也是相互关联的。在许多的情况下这种关联表现为数学方程的形式。下面举例说明：

设两个物体受力与运动相似

则它们的质点的运动方程和力学方程均可用同一方程描述，即： 原型的运动方程与物理方程

dt ds v =

dt dv m f = ① 模型的运动方程与物理方程

t d s d v ''=

' t d v d m f '

'

'=' ② 因为两个物体的现象相似，其对应物理量互成比例，即

s c s s

='

t c t t =' t c v v

='

m c m m =' f c f f ='

③ ①，②，③联合得到

1==c c c c s

t

v ④

1==c c c c c v

m t f

⑤

由④，⑤可以说明，各相似常数不是任意选择的，它们之间是相互关联的。

三、相似三定理

1. 相似第一定理

相似第一定理是指出两个相似物体之间物理量的关系，具体可以归纳为二点。一、 相似现象可以用完全相同的方程组来表示。二、用来表征这些现象的一切物理量在空间相对应的各点在时间上相对应的各瞬间各自互成一定比例关系。

2. 相似第二定理

相似第二定理描述了物理体系中各个物理量之间的关系，相似准则之间的函数关系。 π关系式（准则方程）

0),,,(211=n f πππ

π关系式的性质

① 对于彼此相似的现象，π关系式相同。

② π关系式中的π项在模型试验中有自变项与应变项之分。自变项是由单值条件的 物理量所组成的定性准则，应变项是包含非单值条件的物理量的非定性准则。

③ 若能做到原型与模型中的自变π项相等，由应变π项与自变π项之间的关系式可以得到应变π项，然后推广到原型中去，作为工程设计的各种参数。

3. 相似第三定理

相似第三定理是解决两个同类物理现象满足什么样的条件才能相似的问题。

第一条件：由于相似现象服从同一的自然规律，因此，可被完全相同的方程能所描述 第二条件：具有相同的文字方程式，其单值条件相似，并且从单值条件导出的相似准则 的数值相等。

所谓的单值条件是指从一群现象中，根据某一个现象的特性，把这个具体的现象从一群现象中区分出来的那些条件，单值条件中的物理量又称为单值量。单值条件包括几何条件、物理条件、边界条件和初始条件。

4. 相似三定理之间的关系

相似第一和第二定理是从现象已经相似这一基础上出发来考虑问题，第一定理说明了相似现象各物理量之间的关系，并以相似准则的形式表示出来。第二定理指出了各相似准则之间的关系，便于将一现象的实验结果推广到其它现象。相似第三定理直接同代表具体现象的单值条件相联系，并且强调了单值量相似，所以显于出了科学上的严密性，是构成现象相似的充要条件。是一切模型试验应遵守的理论指导原则。

但是在一些复杂的现象中，很难确定现象的单值条件，仅能借经验判断何为系统最主要的参量，或者虽然知道单值量，但是很难做到模型和原型由单值量组成的某些相似准则在数值上的一致，这使得相似第三定理真正的实行，并因而使模型试验结果带来近似的性质。

一、 同类相似与异类相似

同类相似是指相似的物体是同类物质，模型与原型的全部物理量相等，物理本质一致， 区别在于各物理量的大小比例不同。异类相似是指相似的物体不同类。仅因为对应量都遵循相同的方程式，具有数学上的相似性。

五、相似准则的导出方法

相似准则的导出方法有三种：定律分析法，方程分析法和因次分析法。从理论上说，三种方法可以得到同样的结果，只是用不同的方法对物理现象作数学上的描述。但是作为三种

不同的方法，又有各自的适用条件。

1. 三种方法的介绍

定律分析：这种方法是建立在全部现象的物理定律已知的基础上的，通过剔除次要因素，从而推算出数量足够的，反映现象实质的π项。这种方法的缺点上：

1)流于就事论事，看不出现象的变化过程和内在联系，故作为一种方法，缺乏典型意

义

2)由于必须找出所有的物理定理，所以对于未能掌握其全部机理的，较为复杂的物理

现象，运用这种方法是不可能的，甚至无法找到近似解

3)常常有一些物理定理，对于所讨论的问题表面上看去关系不密切，但又不宜于妄加

剔除，而必须通过实验找出各个定律间的制约关系，决定其重要因素，这实际问题

的解决带来不便。

优点：对于模型制作有指导性意义。

方程分析法：根据已知现象的微分或积分方程推出π项。此方法的的优点：

1)结构严密，能反映出现象的本质，故可望得到问题的可靠性结论

2)分析程序明确步骤易于检查

3)各种成份的地位一览无遗，有利于推断，比较和校验

缺点：对现象的机理不清楚，没有建立方程的问题，无法解决

因次分析法：是根据正确选定参量，通过因次分析法考察各参量的因次，求出和π定理一致的函数关系式，并据此进行相似现象的推广。因次分析法的优点，对于一切机理尚未彻底弄清，规律也未充分掌握的现象来说，尤其明显。它能帮助人们快速地通过相似性实验核定所选参量的正确性，并在此基础上不断加深人们对现象机理和规律性的认识。

以上各种方法，日前应用最广泛的是因次分析法，但是也不排除将各种方法结合使用的可能性。

六、相似准则导出方法的解题步骤

1. 三种方法的解题步骤

1）定律分析法的步骤

①分析现象，抓住主要矛盾，排除次要因素

②写出主要矛盾的物理表达式

③作等效变化，转化为具有相同因次的物理量

④两两作比值，求出相似准则π

2）方程分析法

通常的方程分析法有：相似转换法和积分类比法

相似转换法的步骤

①写出现象的基本微分方程

②写出全部的单值条件，并令其二现象相似

③将微分方程按不同现象写出

④进行相似转换

⑤求出相似准则π

积分类比法的步骤

①写出现象的基本微分方程和全部的单值条件

②用方程的任一项，除其它各项

③进行积分类比转换，求出相应的准则

3) 因次分析法

因次分析法一般分为两种：指数分析法和矩阵分析法。这两种方法的基本原理一样，运

算步骤稍有不同。指数分析法主要用于现象的物理量较少的情况，而矩阵分析法主要用于现象物理量较多的情况。

指数分析法

① 列出相似准则的表达式

② 根据方程两边因次相等列出物理量参数的方程K 个 ③ 设物理量有M 个，任选其中的M －K 个物理量为已知量

④ 将这M －K 个物理量，依次用M －K 个单位向量代入方程，得到M －K 组解 ⑤ 把这M －K 组解代入相似准则的表达式中，可以得出M －K 个独立的相似准则 矩阵分析法

矩阵分析法与指数分析法的基本原理一样，矩阵分析法把线性方程组的求解用矩阵的求 解来代替。其运算步骤不再此重复。

2. 证明指数分析法解出的独立π项的广泛代表意义

例设某现象由5个物理量A1，A2，A3，A4，A5组成，这5个基本物理的独立因次为L ，M ，N

物理量的表达式

i i i T M L Ai γβα=

5,4,3,2,1=i

相似准则的表达式

v u z y x A A A A A 54321=π

因为π项为零，故有 对于L 054321=++++V U Z Y X ααααα 对于M 054321=++++V U Z Y X βββββ

对于T

054321=++++V U Z Y X γγγγγ

固定U ，V 这两个参数，设U ＝0，V ＝1则可以得出一组解，设为X ＝X1，Y ＝Y1，Z ＝Z1，但若设U ＝0，V ＝N 则方程得出另一组解，设为X ＝X2，Y ＝Y2，Z ＝Z2这两组解之间存在着如下关系，即：

211X N X =

2N 1Y1Y = 21

1Z N

Z = 由上式可知，这个相似准则和前一个相似准则只差方次关系，又因为相似准则可以通过

加、减、乘、除、幂运算等进行相互变换，故这两个相似准则实为同一个无因次量群。

设U ＝1，V ＝0则可以得出一组解，设为X ＝X3，Y ＝Y3，Z ＝Z3，但若设U ＝1，V ＝1则方程得出另一组解，设为X ＝X4，Y ＝Y4，Z ＝Z4

这两组解之间存在着如下关系，即：

2314X X X +=

2Y3Y1Y4+= 2

3

14Z Z Z += 故U ＝1，V ＝1的相似准则可以用U ＝0，V ＝1和U ＝1，V ＝0的相似准则表示，所以说

U ＝0，V ＝1和U ＝1，V ＝0的相似准则可以表示U ，V 为任何实数的相似准则。

3. 三种方法解题

1）定律分析法

已知一个简支梁受有大小为4KN/M 均布荷载，简支梁的跨度为4M ，截面的高为0.5M ，

宽为0.4M ，跨中截面的最大正应力为4802

/M KN ，求当梁的跨度为2M ，截面尺寸相同受均布荷载为2KN/M 时的跨中截面的最大正应力。

跨中弯矩的公式

M ＝8

2ql

最大正应力公式

2

6bh

M

=

δ 解：由最大正应力公式可以推出6

2δ

bh M =

又因为2

8ql M =

所以2

2ql bh πδ

=

由m ππ=得

22

22

m m m m m

q l ql bh b h δδ= 又因为截面的尺寸相同所以可以简化为

2

2

m m m

q l ql δ

δ=

所以2

2

m m m q l ql

δδ==602/M KN 2）方程分析法

以弹性力学中的极坐标的平面应力问题为例说明 1.写出现象的基本微分方程 1）静力学平衡方程

10

210f f ρ

ρ?ρ?

ρ?

ρ?ρ??δτδδρρ?ρ

δττρ?ρρ

??-+

++=????+++=??

2）几何方程

()

11

111

111

111

2

3451111 ( ) 1

U l E E l

U E E

c u c c c c c c c c c c c c c c c c E

c c c c c c l m c c l l lU

l U E E ρερμδ

δεεερδ

ρμδ?μεδ

δρρ?ρρδερ

εδμδρδτδδπρρεεπεεπδδμμπρρπδδ----------?=

?==

==-===+='=='''==''

'=='''=='=='

1u u ρ

?

?ερρ??=+?

1u u u ρ??

ρ?γρ?ρρ

??=

+-?? 3）物理方程

()()1

1

2(1)E E

E

ρρ???ρρ?ρ?

εδμδεδμδμγτ=

-=-+=

4）边界条件（2个）

l m ρρ?ρδτ=+

另外一个类似

2.写出全部的单值条件，并令其二现象相似 1）几何单值条件相似

'l c ρρ= 'c εεε= 'U U c U

= 式中：l c c ε U c 分别表示长度相似常数，应变相似常数和位移相似常数 说明?不为单值条件，且为无因次量 2）物理单值条件相似

'c μμ

μ

=

'E E

c E

=

'c ρρ

ρ

= 式中：c μ E c c ρ分别表示泊松比相似常数，弹性模量相似常数和容重相似常数 3）位移边界条件相似

'c δδδ= 'c ρρ

ρ

= 式中：c δ表示应力的相似常数，q c 为面力的相似常数 3.将微分方程按不同现象写出

第二现象的静力平衡方程（只写一个，另一个类似）

10f ρρ?ρ?ρδτδδρρ?ρ''''??-'+++='''

?? 几何方程（只写一个，其它类似）

u ρ

ρ

ερ'?'='

? 物理方程（只写一个，其它类似）

()1

E ρ

ρ?εδμδ''''=-'

边界条件（只写一个）

l m ρ

ρ?ρδτ'''=+ 4）进行相似转化

将有关的相似系数代入得 对平衡方程

111

1111

()10l l l c c c c f c c c δρ

δρ?δρ?ρρδτδδρρ?ρ

-------??-+++=?? 这了保证与原型方程的一致，必须使得

l l l

c c c c c c c δδδ

ρ===

即

1l c c c δ

ρ

= （从另外的一个方程也可以得到这个结果） 对几何方程

111U l c u c c ρερερ

---?=

?

为了保持与原方程的一致，可得

U

l

c c c ε=

即1l

U

c c c ε= （从另外的二个方程也可以得到这个结果） 对物理方程

()11

111

1E c c c c c E

ερδ

ρμδ?εδμδ-----=

-

为了保持与原方程的一致，可得

E E

c c c c c c μδ

δε==

即1c μ

=

1E

c c c εδ

= 从另外的二个方程也可以得到这个结果 对边界条件

11 ( ) c c l m δρρ?ρρδτ--=+

为了保持与原方程的一致，可得

1c c ρδ

=

5）求出相似准则

12

l l lU

l U δδπρρεεπ'==''

'==''

3

4

E E εεπδδ

μμπ'=='''==

5ρρπδδ'=='

与弹性力学的直角坐标系下的相似准则的比较可知是一样的，这同时也说明了相似准则与坐标系的选取没有任何关系。

3）因次分析法

设有半平面体，在其直边界上受有集中力，取单位厚度的部分来考虑，影响力F 作用点正下方h 深处的正应力δ的参量有F 与边界法线成角β，设单位厚度上所受有力为F ，埋深h

解：取基本因次M 、L 、T 此问题独立的相似准则m=4－3=1

A1 A2 A3 A4 δ

F h β

M 1 1 0 0 L -1 0 1 0 T -2 -2 0 0 根据因次和谐原则得 A1+A2＝0 -A1+A3＝0 -2A1－A2＝0 固定A1＝1

则可以得到A2＝-1 A3＝1 所以F h

πδ=

此问题的弹性力学解析解为2cos F

h

δβπ=-

相似准则与原问题的解保持了一致性

中盛上行开采相似模拟实验报告

1.相似模拟实验的目的

根据横向课题《近距离煤层群上行复采薄煤层可行性研究》的研究需要，并结合现场的实际要求，针对山西汾西矿业中盛煤矿的煤层地质赋存条件，进行近距离煤层群上行开采的三维相似材料模拟实验。

主要研究目的是：

●分析10#煤开采对上部9#、7#煤层的影响，判断上部煤层的整体性及可采性。

●测定10#煤开采过程中上部煤层的应力变化规律和位移变化，分析极近距离煤层群

上行开采的卸压特征。

实验时间：3D相似模型于2011年11月日开始铺设，历时日完成，2011年11月日进行模拟开挖实验。

主要实验设备：

3D模拟实验台

7V14数据自动采集系统及数据处理系统

位移计和应力传感器（压力盒和直角应变花）

照相机

实验原型基本条件：实验原型为山西汾西矿业中盛煤矿7#、9#、10#煤层。

1、7#煤层

位于太原组中部，K4与K3石灰岩之间，下距9#煤层19.50～26.30m，平均22.60m。煤层厚0.10～1.95m，平均1.11m。一般含0～1层夹石，局部达2层，结构简单，属不稳定的局部可采煤层。顶板一般为泥岩，底板多为中粒砂岩。

2、9#煤层

位于太原组下部，下距10#煤层4.80～7.83m，平均6.45m。煤层厚0～1.10m，平均0.90m，不含夹石，结构简单，该煤层可采范围分布于井田中东部，可采面积 2.139km2，属较稳定的局部可采煤层。顶板为石灰岩，底板一般为泥岩。

3、10#煤层

位于太原组下部，下距11#煤层6.76～13.82m，平均11.20m。煤层厚0.85～2.62m，平均1.86m。一般含0～1层夹石，局部达两层，结构简单。顶板一般为泥岩，底板为泥岩或细粒砂岩，属全区稳定可采煤层。该煤层于原矿区内已大部采空，仅西南部有小部分未开采地段。

复采煤层特征见表1所示。

表1复采煤层特征表

煤层号煤层厚度（m）煤层间距（m）

夹石

层数

顶板岩性

底板岩性

稳定程度

可采程度最小～最大

平均

最小～最大

平均

7 0.10～1.95

1.11 19.50～26.30

22.60

0～1层

局部2层

泥岩

中粒砂岩

不稳定

局部可采

9 0～1.10

0.90

石灰岩

泥岩

较稳定

局部可采

4.80～7.83

6.45

10

0.85～2.62

1.86

0～1层

局部2层

泥岩

泥岩或细粒砂岩

稳定

全区可采煤层综合柱状图如图1所示。

图1 煤层综合柱状图

2 相似模拟实验

2.1 模型相似比

本实验采用中国矿业大学（北京）的三维试验台，实验台尺寸为：长×宽×高为3000mm ×200mm ×1800mm 。设几何相似比为αL =50：1，设容重比为αγ=1.6：1，要求模拟与实体所有各对应点的运动情况相似，即要求各对应点的速度、加速度、运动时间等都成一定比例。所以，要求时间比为常数，即：L t a =

α= 7.07

式中： t α—时间相似比。

试验采用平面应变条件，各岩层在相似模型中的厚度为：

L H M a L L /11= =2000/50=40.0mm ； L H M a L L /22==6400/50=128.0mm ；

L H M a L L /33==900/50=18.0mm ； L H M a L L /44==8000/50=160.0mm ；

L H M a L L /55= =8000/50=160.0mm ； L H M a L L /66==4000/50=80.0mm ；

L H M a L L /77==3000/50=60.0mm ； 8M L =L H a L /8=1110/50=22.0mm ； L

H M a L L /99==4000/50=80.0mm

；

L

H M a L L /1010==8000/50=160.0mm ；

L H M a L L /1111==10000/50=200.0mm 1212/M H L L L α==10000/50=200.0mm

1205.2mm 其中M12为最上部覆岩。M1为最下部10#煤。

2.2模拟岩石的力学性质

根据试验目的，选择试验方案。本实验主要研究采场煤体变形及其上覆岩层的运移和应力的变化规律以及无煤柱条件下巷道变形规律以及顶板受力分析。研究现场实测的地质资料，根据煤岩层的分布，详细地整理了中盛煤矿各煤层顶、底板岩层物理力学性质的实测数据，详细情况见下表。

表2 模拟岩层物理力学性质表

2.３模型岩石的强度指标计算

逐层计算模型岩石的强度指标，由L a =50，γa =1.6得γσ

a a a L ?==８0。

由主导相似准则可推导出原型与模型之间强度参数的转化关系式，即：

[][][][]σ

γσσσγγσa a a L L c L c H c H M

H M M

c =?==· 式中： c σ—单轴抗压强度

根据上面的资料，可以求出煤质以及不同顶板岩层模型的单轴抗压强度c σ及容重M γ为：

第一层１０＃煤模型的抗压强度及容重为：

[]1M c σ=１４/８0=0.175MPa 1M γ =１.４/1.6=0.875g/cm 3

第二层泥岩模型的抗压强度及容重为：

[]2M c σ

=27/８0=0.34MPa

2M γ=2.6/1.6=1.63g/cm 3

第三层９＃煤的模型抗压强度及容重为：

[]3M c σ=１４/８0=2.17kg/cm 2

=0.175MPa

3M γ =１.４/1.6=0.875g/cm 3

第四层石灰岩的模型的抗压强度及容重为：

[]4M c σ=56/８０=0.7MPa

4M γ =2.７3/1.6=1.71g/cm 3

第五层砂质泥岩的模型抗压强度及容重为：

[]5M c σ=37/８0=0.46MPa

5M γ =2.6/1.6=1.63g/cm 3

第六层石灰岩的模型的抗压强度及容重为：

[]6M c σ=63.2/８０=0.79MPa

6M γ=2.７3/1.6=1.71g/cm 3

第七层泥岩的模型抗压强度及容重为：

[]7M c σ=27/８0=0.34MPa

7M γ=2.6/1.6=1.63g/cm 3

第八层７＃煤模型的抗压强度及容重为：

[]8M c σ=１４/８0=0.１75MPa

8M γ =１.４/1.6=0.875g/cm 3

第九层石灰岩模型的抗压强度及容重为：

[]9M c σ=68.74/８０= 0.859 MPa

9M γ =2.70/1.6=1.69g/cm 3

第十层砂质泥岩模型的抗压强度及容重为：

[]10M c σ=47/８０=0.59MPa

10M γ =2.6/1.6=1.63g/cm 3

第十一层中粒砂岩模型的抗压强度及容重为：

[]11M c σ=61.04/８0 =0.76MPa

11M γ =2.67/1.6=1.67g/cm 3

表3 相似模拟试验材料配比表

2.４相似实验材料的制备

根据汾西矿业中盛煤矿煤岩层的实际地质资料，选择组成相似模拟材料的成分，相似模拟材料主要由两种成分组成－骨料和胶结料。骨料在材料中所占的比重较大，是胶结料胶结的对象，其物理力学性质对相似材料的性质有重要的影响。骨料主要有细砂、石英砂、岩粉等，本试验骨料采用细砂。

胶结料是决定相似材料性质的主导成分，其力学性质在很大程度上决定了相似材料的力学性质，常用的胶结材料主要有石膏、水泥、碳酸钙、石灰、高岭土、石蜡、锯末等。根据试验及地质成分，本试验胶结料采用石灰和石膏。

不同胶结料与骨料混合组成不同种类的相似材料，其力学性能不同。根据已计算出的模型的力学参数，选定骨料及胶结料进行配比试验，为了精确选定与计算参数一致的配比，经过了多次配比试验，做出了各种配比表，最后选择出满足试验要求的一种，见表3。

表4 模型铺设分层材料用量表

上覆岩层厚度约为640m ，应力为16Mpa ：

根据模型的尺寸，以及预定比例，实际加载压力为： σ/αL αγ=16/（7.47*1.6）=1.339Mpa 2.5制作模型的步骤

本模型的架子的主体，是由,以便固定模板，模板用厚1厘米的钢模板制成，为了防止装填材料时模板向外凸起，架子中部可用竖向小槽钢加固，模型架的长度保证围岩尺寸为巷道尺寸的10倍以上。

（1）上模板：将模型后面模板全部上好，前面边砌模型边上模板。

（2）配料：按已计算好的各分层材料所需量，把水泥、石膏、沙子、缓凝剂及水用天平、量杯称好；其中沙子、石膏、水泥可装在一个搅拌容器内，但需将水泥倒在石灰上面，以免与沙子中所含的水分化合而凝固，缓凝剂融解后可放入已经称量好的水中，并搅拌均匀使无沉淀。其中，模拟煤层加墨汁，使其变黑，以便区别。

（3）搅拌：先将干料拌匀，再加入含缓凝剂的水中，并迅速搅拌均匀，防止凝块。 （4）装模：将搅拌均匀的材料倒入模子内，然后夯实，以保持所要求的容重，压紧后

的高度应基本符合计算时的分层高度。分层间撒一层云母以模拟层面，每一分层的制作工作应在20min内完成。

（5）风干：通常在制模后一天开始风干，风干后，开始采煤。

（6）加重：由于考虑到地层很深，而所要研究的问题仅涉及煤巷附近一部分围岩，可施加面力的方法来代替研究范围以外的岩石自重。

2.6 监测点布置

模型铺设总高度120.5cm，在老顶上方依次铺设干砂、铁板、槽钢，通过液压油缸对槽钢加压，并通过铁板、干砂均匀传递载荷。

材料水灰比：材料：水=100：8

铺设模型时，每铺设4cm厚度材料后，进行夯实。

图9 3D相似模拟实验台

各岩层之间加铺云母粉。

图10为3D相似模拟实验台，标准尺寸为长×宽×高=300×200×180cm，模型四周采用槽钢封闭，作为位移边界条件；模型上部采用油缸施加压力模拟上覆岩层的载荷，压力的大小由操作控制台（图9中位于左侧的为控制台）调节；模型底部由两部分组成（见图10）：图中阴影部分组成的活动区和活动区以外的固定区，固定区由钢板制成，而活动区尺寸为180×100cm，由90条长100cm宽2cm的活动钢板组成，活动钢板两侧有螺栓可以拆卸，铺设模型时将全部活动钢板装上，模拟开挖时依次卸下活动钢板，就可以从模型底部挖出煤体，而且同时保留固定区的煤体，因此该模型更加真实地模拟实际的开采过程。该模型的唯一缺点是无法直接目视观测。安装各种传感器时，均以活动区为参照系。

图10 3D模拟实验台底部结构示意图

3.2 传感器位置

(1). 上覆煤层压力变化观测：压力盒布置2个层位，第一层为底板以上17.7cm处，即9#煤中。第二层为底板以上65.7cm处即7#煤中每个煤层内沿走向方向布置3测线，每条测线布置3测点。每层设置12个压力盒。共布置压力盒24个。具体位置见图11（图中虚线交点位置即为压力盒布设位置）。

图11 煤体压力力监测布置平面图

省去(2). 位移观测：位移观测布置3个层位；第一层为煤层底板以上17.42cm处即9#煤中，沿走向方向布置3条测线，每条测线布置5个测点，其中中央3个点用于测量煤层沿走向方向的位移，两侧不开挖煤体中各布置1个测点用于测量巷道两侧固定帮煤体向采空区的移动情况，9#煤煤层测点的引线方向为垂直方向（见图13）；第二层为煤层底板以上42.32cm处即9#煤顶板砂质泥岩中，设置三条测线，每条三个测点；第三层为煤层底板以上65.42cm处即7#煤中，沿走向方向布置3条测线，每条测线布置3个测点，7#煤的测点引线方向为垂直方向，具体位置见图14。

位移测点共布置33个。

图12第一层煤层位移测点布置平面图

图13第二、三层煤岩层位移测点布置平面图

(3). 上覆煤层应力变化观测：应变片3个层位，第一层为底板以上17.7cm处即9#煤中，第二层为底板以上42.6cm即9#煤层上部泥质砂岩中，第三层为底板以上65.7cm处即7#煤中；每个煤层内沿走向方向布置4条测线，每条测线布置4个测点，每个点布置两组应变花，中间三条侧线每个测点的2个应变花采用如下布置：1个为垂直于工作面走向并平行

于工作面倾向方向（即垂直放置）；另一个为平行于工作面走向并平行于煤层的方向（即水平放置）。外侧测线每个测点的2个应变花采用如下布置：1个为平行于于工作面走向并垂直于工作面倾向方向（即垂直放置）；另一个为平行于工作面走向并平行于煤层的方向（即水平放置）。

编号为7-1-1-1x,表示7#煤一测线，一测点，一测组x方向应变片。共设置96个应变花。具体位置见图13（图中虚线交点位置即为应变花布设位置）。

三维相似模拟模型共布置各类传感器个。