某煤矿放顶煤开采数值模拟报告

1、某煤矿10#煤层综放开采有限差分正交试验数值模拟研究报告通过对某煤矿 10 #煤层围岩物理力学特性研究、10 #煤层综合机械化放顶煤开采的相似模拟试验研究，得出了该矿10井煤层基本适合于综合机械化放顶煤开采的初步结论，以及围岩活动规律及参数，摸清了综放工作面煤岩组合结 构，从而获得了某煤矿综放开采的宝贵资料；为了进一步掌握工作面工艺参数 对顶煤回收效果的影响规律，优化综放工作面工艺参数，同时为采煤机、液压 支架等设备选型提供依据，在该课题前述研究基础上，利用相似模拟试验得到 的综放采场围岩物理结构，建立相应的数值模型，利用著名有限差分程序的 flaCD进行某煤矿有限差分数值模拟正交试验研究，研

2、究各工艺参数对工作面顶煤回收效果的影响规律及其最佳组合。i4、有限差分法基本原理及特点由于计算机的不断发展，数值模拟已成为研究岩土工程的强有力工具。起 初，岩石力学分析的进展在很大程度上有赖于其它工程领域的进展，特别是力 学在这些领域中的应用。自1962年以来岩石力学的一个主要进展是引入计算机的数值模拟。岩石力学中数值模拟的最早应用集中于有限单元法，该法最初是 为应用于结构材料变形性质的研究而开发的，材料的力学性质可以精确地加以 确定，结构的范围也较岩石结构的范围要小得多。而岩石工程设计的问题则大 为不同。岩石工程设计问题的特点主要有以下几点：(1) 可变性和不确定性；(2) 主要是压应力，包

3、括一些结构的预应力，如锚杆、锚索等支护；还有在 岩体开挖引起卸载或加载所产生的压应力；(3) 岩石的尺寸大小和时间的变化范围大；(4) 地质条件复杂；岩石工程一般是三维的、非均质、不连续的；(5) 岩体的强度和变形特性未知，“峰值后”的力学性质尤为重要；(6) 边干边设计一一即要根据开挖或建设时发生的不可预见的条件迅速改 正设计；(7) 解决问题的数据有限；(8) 分析的主要目的为“理解机理”而并非精确的定量计算。由于岩体的力学特性是很复杂的，以致当今最有力的分析和数值方法也不 能描述岩石变形的机理，因而只能从分析中得到对具有内在不确定性和多变性 的岩体力学性质的定性的结论，即得到对物理过程的

4、总的了解。然而，这种了 解的价值对于设计来讲是至关重要的。一旦内在的机理清楚了，工程师就能根 据预计条件的变化来进行设计。能否根据条件的变化来进行设计，是一个成功的和经济的设计同费用高甚至会发生事故的设计的区别所在岩石工程中早期的数值模型主要系采用原位结构分析所开发的数值方法，其特点是：（a）所用材料的力学性质可以由实验室试验来确定，（b）整个结构的性能可以直接根据原型试验来确定。在很多场合，变形可以限制在线弹性范围 内的小变形。在岩石结构中则往往必须考虑大变形、大应变和非线性，在很多 场合中还要研究“峰值后强度”的变化，这时结构可能破坏而呈应变软化的状 态。实际上即使岩石的强度为已知，也建议

5、去考虑强度高估的效应。由于认识到岩石和岩体的这种特点，许多专家认为在岩石变形的模型中用显式的有限差分法较在其它领域中广泛应用的隐式的有限单元法为好。有限差分法是给定初值和（或）边值求解微分方程组的一种数值方法。在用有限差分法时，控制方程组中的导数直接用空间离散点的场变量（例如应力或位移）的代数表达式来代替，这些变量在单元中是未定义的。相反，有限元法的一个根本的要求就是在每个单元中场变量（应力，位移）要用一定的参数函数来描述。在有限元法的列式中，这些参数要根据误差项或 能量项为最小的原理来确定。这两种方法都要解一个代数方程组。虽然这两个方程组系用不同的方法得到，但是容易证明（在特殊情况下）彼此的

6、结果是等同的。表1所示为有限差分方法和有限元隐式方法两者的比较。有限差分方法的缺点是需要小的时步，也就是时步数要很多。但有限差分方法更适用于病态的 系统，例如非线性，大变形和非稳定的情况，而对于线性，小变形的问题则效表1有限差分法和有限元法的比较有限差分有限元计算时步要取得比为稳定所需的临界 值大在用无条件稳定的格式时时步可任意大每个时步的计算开销小每个时步的计算开销大对于动态问题没有显著的数值阻尼在用无条件稳定的格式时数值阻尼和时步有关对于非线性本构方程无需迭代对于非线性本构方程需要迭代不用形成矩阵，要求内存小，无带宽 的限制必须存贮刚度矩阵，必须要解决随之而来 的例如带宽冋题，内存要求大由

7、于无需形成矩阵，大位移和大应变无需附加的机时为跟踪大位移和大应变需要附加的机时由于有限差分不需要形成刚度矩阵，所以在大变形模式时在每一个时步更 新座标是一件常规的事。位移的增量加在座标上，因此网格和它所代表的材料 一起移动和变形。这称之为“拉格朗日”表示法，与之相对的是材料参照固定 网格移动和变形的“欧拉”表示法。在每一个时步时的变形增量是一个小数， 但在许多时步后变形则是大变形。由于计算技术（特别是微机和工作站）的迅猛和持续的发展， 数值模拟在近十 几年中发展得非常快。这样，在计算程序中包括影响岩体力学响应的许多现实 的复杂性（非均质，节理，时间因素和三维效应等）就愈加成为可能。有限差分法数

8、值计算的建模原理：有限差分的计算方法，应用节点位移连续的条件，可以对连续介质进行大 变形分析。假定在有限的时间和空间间隔内，变量是线性变化的，变量关于空 间和时间的一阶导数均用有限差分来近似。利用动态逼近的求解方法，应用质 点运动方程求解，用惯性来衡量所考虑的系统是否达到平衡状态。利用有限差分法求解一个静力问题的解时，在计算式中要包括运动的动力 学方程。这样做的好处是即使所模拟的物理系统是不稳定的也要保证数值方法 的稳定。材料的非线性往往会导致系统的失稳，例如矿柱的突然破坏。在现实 生活中，系统中应变能的一部分转变为动能，然后辐射出去消散掉。有限差分 法直接模拟了这一过程，因其包括惯性项，动能

9、产生又消散。利用有限差分法进行数值模拟 ，首先，要分析模拟对象的几何特征，按照研 究的目的进行网格剖分，三维模型的网格一般剖分为六面体单元，对需要详细研究和对模型影响较大的部分要进行网格细化，使网格按模型的计算和分析要求合理分布。其次，要为模型的计算设定初始条件和边界条件，边界条件包括 应力边界条件（如面力、体力、集中载荷）和位移边界条件；初始条件为初始 应力条件。初始应力和集中载荷要加载在给定表面的节点上。最后，为模型的 计算过程，如图 3.1所示。先是用运动方程由应力和力求出新的速度和位移。然后由速度推导出应变率，并由此求得新的应力。一个时步为沿回线的一个循 环。平衡方程（运动方程）图1三

10、维有限差分法的基本计算循环二、正交试验设计正交试验设计法是研究与处理多因素试验的一种实用的试验设计法，它是在实际经验与理论认识的基础上，利用一种排列整齐的规格化表一“正交表”来安排试验，能够通过少量的试验次数，找到较好的生产条件，同时还能做许多进一步的分析。正交表有如下基本性质：1、正交性（1）在任何一列中各水平都出现，且出现的次数相等。（2） 任何两列之间各种不同水平的所有可能搭配都出现，且出现的次数相 由正交表的正交性可以看出，表中各列、各行、同一列的水平数均可相互置换，分别称为列间置换、行间置换、水平置换。2、代表性代表性的含义之一，在于正交表正交性中，（1）任一列的各水平都出现时的部分

11、实验中包含所有水平；（2）任意两列间所有搭配都出现，使得任意两列因素间都是全面试验。因此，部分实验中，所有因素的所有水平信息及两两因素间的所有搭配信息无 一遗漏。这样，虽然正交表安排的只是部分实验，但却能了解到全面的试验 情况，在此意义上，部分试验可以代表全面试验。3、综合可比性保证了任一因素各水平的试验条件都相同， 使得每列因素各水平的效果 中，最大限度地排除了其它因素的干扰， 从而可以综合比较各因素不同水平 对试验指标的影响，此性质即为综合可比性。正交表的三个基本性质中，正交性即均衡性是核心，是基础，代表性和综合性是正交性的必然结果，从而使正交性得以具体应用。正交表集三个性质于一体，成为正

12、交试验设计的有效工具。三、综放面工艺参数分析综放采场影响顶煤放出率的重要技术因素主要有：支架的支撑高度、控顶距、支架工作阻力、支撑合力作用点位置和放煤步距等。（1）支架支撑高度对于具体的煤层和工作面来说，支架的支撑高度一定， 底层机采高度就基本确定了。对于煤厚 8m以上的煤层，底层采高在0.5m内变化，其放采比变化不大，但对于小于8m煤层，其放采比就会有明显变化。例如，古书院矿3号煤层，平均煤厚6m ,f > 3.0，如底层采高2.5m，则顶煤厚度3.54.0m，放采比为1.4 1.6，观测表明，可出现 3个台阶的悬臂，其上位1m的顶煤放不出来，顶煤放出率不足70%，若将底层采高提高到3

13、.0m，放采比不大于 1,此时一般为 2个台阶悬臂，顶煤放出率可达85 %以上，因此，对强度较大的煤层，在一定范围内的底层采高，亦即支架的合理支撑高度，将是影响顶煤冒放性的重要因素， 放采比愈大，即顶煤愈滞后，顶煤松动愈困难，其冒放性就愈差。(2)控顶距控顶距主要是影响顶煤松动区的前后宽度，对于厚而硬的煤层， 控顶距应当大些，以便增加支架在控顶区内的反复支撑范围和次数，有利于上部顶煤的充 分松动；而对于软而薄的煤层，控顶距应当适当减小，以防架前漏顶，出现难 以控制的局面。在某些情况下，控顶距可能影响到支承压力的分布。因此，选 择合理的控顶距也是提高顶煤放出率的重要因素。(3 )放煤步距放煤步距

14、是指两次放煤之间工作面推进的距离。合理的放煤步距应是顶煤放出率最高，含矸率最低，因此它与煤矸的块度大小、质量、运动阻力、运动方 向、混矸程度、到放煤口的距离等有关，也就是与煤岩的强度、块度、顶煤厚 度、冒落角、矸堆高度、安息角等有关。实际上放煤步距是采煤机截深的整倍 数，一般为一刀一放、二刀一放和三刀一放三种采放配合方式，因此对于滚筒 式采煤机，放煤步距仅变化在 0.5 1.5m或0.6 1.8m之间。 煤的强度和块度影响分析一般硬煤冒落块度大，冒落角小，软煤冒落块度小，冒落角大。大块煤容易挤压成拱，运动阻力大，冒落滞后，小块煤运动阻力小，冒落超前，因此对于 冒落角小，冒落块度大的硬煤应增加冒

15、落空间，增加冒放次数，易采用一刀一 放，放煤步距小的方式。例如古书院矿3号煤，平均煤厚 6m，f > 3.5，冒块©0.5m以上者占60%一 70%，冒落角55。，采用两刀一放，放煤步距1.1m时，顶煤放出率仅70%左右，改为一刀一放，放煤步距0.55m时，顶煤放出率79.6 %。 对于冒落角大，冒落块度小的软煤应减少冒放空间，降低混歼率，易采用两刀一放，放煤步距大的方式。例如：唐山矿开采煤厚6.3m，煤质较软，采用一刀一放，放煤步距0.5m时煤岩界面不均，混矸率高，灰分达30.95%，顶煤放出率81.4 %，当改为两刀一放，放煤步距1.0m时，顶煤冒落充分及时，煤岩界面下降平

16、稳，顶煤放出率 88.7 %，灰分下降为 23.9 %。 顶煤厚度的影响分析众所周知，顶煤的冒块一般是下小上大，下部顶煤松动充分，上部顶煤松动滞后，因此顶煤愈厚应增加其冒落时间和冒落宽度，以便上部原煤有充分松 动冒落的机会而放出，易采用放煤步距大的方式。例如：大雁二矿开采低变质 褐煤，煤质不硬，但煤厚平均14.03m，一刀一放，放煤步距0.6m时，上部顶煤冒落缓慢，混矸率很高，当改用两刀一放，放煤步距1.2m时，顶煤冒落充分及时，冒落的块度明显减小，顶煤放出率达81.8 %。对于厚度较薄的顶煤，由于下部冒落空间大，极易混矸，应减少自然冒落次数，使之在混矸前就及时放出， 易采用一刀一放，放煤步距

17、小的方式。例如：谢桥矿1121综放工作面，平均煤厚4.56m，f = O.35 1.7，采用两刀一放，放煤步距1.2m时，工作面采出率平均74.7 %，原煤含矸率达 19.11 %，当采用一刀一放，放煤步距0.6m时，工作面采出率达88.16%，原煤含矸率仅11.25%。5)(6)图2放煤步距与放煤口水平投影长度关系 放煤口位置和纵向尺寸的影响分析图3放煤步距与截深的关系对于低位插板式放煤口均为架间连通式，无脊背损失，由于顶煤冒落空间大，相对来说松散程度高，其放煤步距主要取决于放煤结构和顶煤冒落角。当冒落 角大，半拱式放煤结构时，放煤步距应大些（图3），此时尾梁摆动放煤的范围较大，若一刀一放，

18、放煤步距小，则过早混矸，上部顶煤放不出来，若两刀一放， 放煤步距大，顶煤能充分放出，但采空区丢失三角煤较多。这种情况下，当截 深大，0.81.0m时，一般采用一刀一放方式，若截深小，0.50.6m时，一般采用两刀一放方式。如：潞安矿务局煤厚6.57.0m，顶煤冒落角70°80°，截深0.8m，刀一放的情况下工作面采出率最高，达86.9%，含矸率较低，为13.1 %，若采用0.6m截深，两刀一放时，虽然含矸率稍低些，为10.7 %，但工作面采出率仅 82.1 %，故选用一刀一放，放煤步距0.8m的方式。鲍店矿煤厚8.5 9.0m，冒落角70。以上，采用一刀一放，放煤步距0.6

19、m时，冒落空间狭窄，上部顶煤不能松动而放不出来，当采用两刀一放，放煤步距1.2m时，工作面采出率达81%，含矸率4.77 %，若采用三刀一放，放煤步距1.8m时，工作面采出率仅77.6%，含矸率4.89%。当冒落角小、漏斗式放煤结构时，亦即坚硬煤层条件下，由于顶煤呈倒台阶 悬臂，放煤步距应小些，主要应防止上部悬臂顶煤冒落在采空区的矸堆上而放 不出来（图4），此时若采用一刀一放（图4a），冒落顶煤正好位于后部输送机的上方，收回插板，摆动层梁，呈漏斗状放出；若采用两刀一放（图4b），第一刀冒落的顶煤未放出，移架后均落在采空区的底板和矸堆上，第二刀冒落的顶煤位置 在拉移后的后部输送机上，此时，第一刀

20、冒落的顶煤与第二刀冒落的矸石混合，8m，将有一大部分在输送机的后方而不能放出。例如：忻州窑矿煤厚平均图4放煤步距与顶煤悬臂的关系冒落角50°，顶煤呈3个台阶悬臂，采用两刀一放，放煤步距1.0m时，可见大块煤均落在后部输送机的后方，摆动尾梁根本不能扰动这部分煤放出，顶煤放出率仅60%左右，当改为一刀一放，放煤步距0.5m时，顶煤放出率可达 73%,大大提高了工作面的采出率。综合以上分析，放煤步距的确定应遵循如下规律： 低位放煤条件下，煤厚在8m以下，一般采用一刀一放的方式，放煤步距0.5 0.8m，软煤取大值，硬煤取小值，煤厚大于8m者，可采用小截深的两刀一放方式，放煤步距 1.0 1

21、.2m。 高位和中位开天窗放煤条件下，放煤步距与放煤口纵向水平投影长度相等，一般均采用两刀一放方式，放煤步距1.21.6m，个别放煤口纵向水平投影长度过大，而循环作业又允许的条件下，可用三刀一放的方式，但放煤步距也不要大于1.8m。(4)放采比的确定放采比是指放煤厚度与采煤厚度之比。采煤厚度就是底层采煤机的采高或底层炮采采高，放煤厚度就是顶煤厚度。放采比的大小直接影响工作面的采出率和放顶煤的经济效益，下面从两个方面分析放采比的确定。分析最大最小放采比最大放采比是指放顶煤最大可能的放煤厚度与底层采高之比。众所周知，任何松散体的冒落都是自然成拱的，且拱高与拱宽是成一定比例的，不可能无限地冒落，煤也

22、不例外。顶煤是一种半拱式的冒放，前方为实体煤，后方为冒落松散体，拱宽向前受到限制，向后虽无限制，但冒落在采空区的煤是不能回 收的，是一种狭窄的有限拱宽，因此冒落拱高就是有限的。根据顶煤冒放性理论，放顶煤的最大开采厚度为13m , 13m以上的顶煤是放不出来的。其中底层采高变化在2.53.5m，所以放采比为10/3.5 10/2.5 ,即最大放采比为 34。 通常在最大采厚条件下，为了获得最高采出率，底层采高也应最大，故放顶煤 的最大放采比一般不超过 3。最小放采比是指放顶煤最小可能的放煤厚度与底层采高之比。放采比过小， 混矸率高。单一长壁综采的最大采高以4.5m为宜，那么放顶煤的最小采厚就是4

23、.5m，此时的放采比一般为 2/ 2.5 1.7/2.8，即0.6 0.8，也就是说，低于 0.6 的放采比是不合理的。因为支架对顶煤的破碎高度一般不超过2m，因此小于1.5m厚的顶煤，在冒落和放出时混矸层比例加大，不利于顶煤的回收。确定已知煤层厚度的放采比通常放顶煤开采设计之前，煤层厚度是已知的，此时的放采比主要取决于 底层采高的确定。原则上认为冒放性好的煤层，为了获得最佳的经济效益，即 提高出块率，降低能耗等，放采比可大些；冒放性差的煤层，为了提高采出率， 放采比可小些，因为这种情况下加大采煤厚度，减小放煤厚度，在总体上可以 提高工作面采出率。例如：对于冒放性差的坚硬煤层，底层采高一般应选

24、 3.0 3.5m，充分发挥底层采出率高的优势来提高采出率。对于冒放性好的中硬 以下煤层，底层采高应选2.5 2.8m，既有利于控制片帮和架前漏顶，又充分发挥了放顶煤能耗低、成本低的作用，我国放顶煤实践证明，大多数情况下放采 比为1 2是合理的。(5)液压支架工作阻力和支撑合力作用点位置大量现场实践和模型实验表明，支架阻力仅对下位 2m左右的顶煤破碎有影响，但是来压时，工作阻力对顶煤起压缩和水平膨胀的破碎作用，因此提高支 架工作阻力对顶煤冒放性，特别是下位顶煤的冒落角有一定影响。支架支撑合 力作用点位置是指整架合力作用点到煤壁的距离，它与无立柱空间的宽度有关， 也与前后柱的间距有关。一般来说，

25、支撑合力作用点靠近前柱，有利于保持机 道上方顶煤的稳定性，防止架前漏顶，但顶煤的冒落块度大，放出率不高。如 果支撑合力作用点靠近后柱，顶煤易处于受拉状态，有利于破碎，冒落块度小， 放出率就高，但是在某些情况下，容易发生架前漏顶。例如，王庄矿4309工作面的高位单输送机放煤支架，其支撑合力作用点距煤壁仅1.5米，无须担心架前漏顶问题，但是冒落块度大，直径© 0.50.8m的大块占70%以上，直径© 1.0m的大块5分钟内在运输机上见到7块。而该矿6102工作面的低位双输送机放煤支架，支撑合力作用点距煤壁3.5m，工作面煤块70%以上为© 0.30.5m，基本上见不到

26、© 1.0m以上的大块。以上架型参数的影响是定性分析，对于某具体煤层、采区或工作面，应通过数值分析，进一步优化选择对冒放有利的架型参数。四、数值模拟方案设计影响顶煤冒放性的技术因素除工作面长度外，主要是一些与架型有关的参数，合理选择这些参数，将会对顶煤冒放性达到最佳的控制效果，即可防止架前漏顶，又可达到最高顶煤放出率的要求。实现这一目的的方法是数值分析法，它是根据不同的顶煤顶板组合移动结构，建立相应的数值模型，运用正交试验法，找出各种架型参数的影响规律，进行优化选择的。数值计算用 5个架型参数的正交法确定，即支架的支撑高度、控顶距、支架工作阻力、支撑合力作用点位置和放煤步距。每个参数

27、根据实际需要取一定 范围的变化值，运用正交法得出若干个实验方案。各因素在正交表中的列位置随机确定。因素及其水平：放煤步距：0.6、0.9、1.2、1.5、1.8m。液压支架合力作用点到支架前方距离：2.1、2.7、3、3.6、4.5m。控顶距：3.3、3.9、4.5、5.1、5.7m。采煤机采高：2、2.4、2.8、3.2、3.6m。液压支架工作阻力：3800、4200、4600、5000、5400KN/架。支架阻力的模拟：阻力的一半为均布载荷，一半为合力作用点载荷。所用的正交表见表2。12表2 所用正交表13表3各方案参数及液压支架工作阻力换算列号试验放煤步距作用点距离控顶距采高工作阻力均布

28、载荷集中力mmmmNPaN10.62.13.323.8e63.84E+051.27E+0620.62.73.92.44.2e63.59E+051.40E+0630.63.34.52.84.6e63.41E+051.53E+0640.63.95.13.25.0e63.27E+051.67E+0650.64.55.73.65.4e63.16E+051.80E+0660.92.13.92.85.0e64.27E+051.67E+0670.92.74.53.25.4e64.00E+051.80E+0680.93.35.13.63.8e62.48E+051.27E+0690.93.95.724.2e62

29、.46E+051.40E+06100.94.53.32.44.6e64.65E+051.53E+06111.22.14.53.64.2e63.11E+051.40E+06121.22.75.124.6e63.01E+051.53E+06131.23.35.72.45.0e62.92E+051.67E+06141.23.93.32.85.4e65.45E+051.80E+06151.24.53.93.23.8e63.25E+051.27E+06161.52.15.12.45.4e63.53E+051.80E+06171.52.75.72.83.8e62.22E+051.27E+06181.53.

30、33.33.24.2e64.24E+051.40E+06191.53.93.93.64.6e63.93E+051.53E+06201.54.54.525.0e63.70E+051.67E+06211.82.15.73.24.6e62.69E+051.53E+06221.82.73.33.65.0e65.05E+051.67E+06231.83.33.92.85.4e64.62E+051.80E+06241.83.94.53.23.8e62.81E+051.27E+06251.84.55.13.64.2e62.75E+051.40E+06上表中方案10、15液压支架合力作用点距离超出支架控顶距，

31、为不可能方 案，为方便试验结果分析对数据数目的要求，均给予计算结果中最小破碎系数。五、数值模型的建立根据某煤矿10 #煤层综放开采相似模拟试验的煤岩结构建立数值模型。 I | 丁 T T 1图5某煤矿综放采场煤岩结构模型尺寸：水平方向工作面前方20m,工作面后方 20m,总共52m垂直方向 56m。15JOB TITLE :.(*10A1)16JOB TITLE :.(*10A1)FLAC (Version 4.00)LEGEND7.00019-May-03 11:13step 19825-1.600E+01 <x< 9.600E+01-2.791E+01 <y< 8.

32、410E+01Grid plotiI U ! i L u I i.u i i.iJ02E 1Boun dary plotI ! i L iL L1 I 1.J ! J02E 1User-defi ned Groupspeat coal10 #1 packsand main roof coarsegra ined coal9 coal10 #2Itasca Con sult ing Group, I nc. Minn eapolis, Minn esota USAJOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 11:13step 19825-1.09

33、1E+01 <x< 9.091E+01-2.282E+01 <y< 7.901E+01Grid plot02E 1Boun dary plot02E 1Itasca Con sult ing Group, I nc. Minn eapolis, Minn esota USA0.0002.0004.0006.0008.0005.0003.0001.000-1.000(*10A1)图6采场围岩(*10A1)7.0005.0003.0001.000_ -1.0000.0002.0004.000(*10A1)6.0008.000图7数值模型网格剖分17X S T TLE12 W

34、5ieQ 01 埔锂 旳1<?« £ «JJE+O1-$«»e+<o 呵12 ii7e*oi繼tR IE Fotsas|i*vectoi- Z *1OE-+CI6Q3E S桂t邛P IE FOTJESVecfcT - Z 貂 OEM$QSE SB: Jtiary plotzm» Ccni.uirgnwp me im wsii tmwR* ar = W18#图8液压支架支护作用模拟工作面采深175m 数值模型加载载荷：(175 56.36 )X 2.5 = 2.96MPa。六、有限差分数值模拟结果分析根据综放开采支承压力破碎

35、顶煤的现象，有限差分数值模拟注重考虑工作面前方顶煤在不同工艺参数影响下的破碎状况，即数值模拟时综放采场数值模 型煤壁前上方 24卅顶煤区域内破坏单元面积所占比率，称之为顶煤宏观破碎系 数Y;显而易见，顶煤宏观系数Y值越大，说明顶煤破碎效果越明显，顶煤回收效果越好。图931为有限差分程序 FLACT (Version 4.00 )计算的各方案顶煤 破坏单元分布图，表4为各方案破碎系数计算结果及其极差分析法、方差分析法处理结果，明确了各工艺参数对顶煤破碎状况的影响程度。#JOB TITLE 一FLAC (Version 4.00)LEGEND19-May-03 11:13step 198251.9

36、74E+01 <x< 3.698E+013.050E+00 <y< 2.029E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in pastGrid plotI I I I I I05E 0Boundary plot05E 0Itasca Consulting Group, Inc. Minneapolis, Minnesota USA图9方案1顶煤破碎状况JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 11:52step 1919

37、92.003E+01 <x< 3.670E+013.337E+00 <y< 2.000E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in past Grid plotBoundary plotItasca Consulting Group, Inc. Minneapolis, Minnesota USAJOB TITLE :.FLAC (Version 4.00)LEGEND22-May-03 21:38step 207012.003E+01 <x< 3.670E

38、+013.337E+00 <y< 2.000E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in pastGrid plot2.0001.8001.6001.4001.2001.0000.8000.6000.400Il III III Hl III III IIIs I1： . r-Mil I-mi- i I u u s_ nJub Hmu M u IIB! 空测日肚划吐日 阳"J1.8001.6001.4001.2001.0000.8000.6000.4001詢1:图10 方案

39、2顶煤破碎状况05E 0Boundary plot05E 0Itasca Consulting Group, Inc. Minneapolis, Minnesota USA1.8001.6001.4001.2001.0000.8000.6000.400图11方案3顶煤破碎状况20JOB TITLE 一(*10n1)FLAC (Version 4.00)LEGEND19-May-03 15:48step 193951.864E+01 <x< 3.530E+013.364E+00 <y< 2.003E+01Plasticity Indicator* at yield in

40、shear or vol.X elastic, at yield in pasto at yield in tensionGrid plot05E 0Boundary plot05E 0Itasca Consulting Group, Inc.Minneapolis, Minnesota USAEIHH n M H m M HI HHM UU&- UUU BJMU k H MlM J bJuiu uui i 竺凹仝亞 En_Egiiii risswrI1.8001.6001.4001.2001.0000.8000.6000.400JOB TITLE :.FLAC (Version 4.

41、00)LEGEND图12 方案4顶煤破碎状况(*10U)19-May-03 16:50step 225651.864E+01 <x< 3.530E+013.364E+00 <y< 2.003E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in pasto at yield in tensionGrid plot05E 0Boundary plot05E 0Itasca Consulting Group, Inc. Minneapolis, Minnesota USA1.8001

42、.6001.4001.2001.0000.8000.6000.400图13 方案5顶煤破碎状况JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 17:04step 193091.962E+01 <x< 3.629E+012.571E+00 <y< 1.924E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in pastGrid plotlii11I105E 0Boundary plotLkiiii05E 0Itasca Consult

43、ing Group, Inc. Minneapolis, Minnesota USA1.7001.5001.3001 100(*10n1):J-沖"(*10n1)1.900335IS3H3d=BirUZLLEIIn咖h-jtli图14 方案6顶煤破碎状况22JOB TITLE ：.FLAC (Version 4.00)LEGEND19-May-03 18:17step 180991.962E+01 <x< 3.629E+012.571E+00 <y< 1.924E+01Plasticity Indicator* at yield in shear or vol

44、.X elastic, at yield in past Grid plot_I5E 0Boundary plot5E 0Itasca Consulting Group, Inc. Minneapolis, Minnesota USA1.7001.5001.3001.1000.9000.7000.5000.300("10n1)("10n1)1.900图17 方案9顶煤破碎状况23JOB TITLE ：.图17 方案9顶煤破碎状况#JOB TITLE ：.图15 方案7顶煤破碎状况图17 方案9顶煤破碎状况#JOB TITLE ：.图17 方案9顶煤破碎状况#JOB TITL

45、E ：.JOB TITLE :.("10n1)FLAC (Version 4.00)LEGEND19-May-03 18:38 step 212801.864E+01 <x< 3.530E+013.255E+00 <y< 1.992E+01Plasticity Indicator* at yield in shear or vol. X elastic, at yield in past o at yield in tension Grid ploti i i n 05E 0Boundary plot05E 0Itasca Consulting Group,

46、Inc.Minneapolis, Minnesota USA1.8001.6001.4001.2001.0000.8000.6000.400图16 方案8顶煤破碎状况JOB TITLE :.("10n1)FLAC (Version 4.00)LEGEND19-May-03 18:54step 202841.858E+01 <x< 3.525E+013.282E+00 <y< 1.995E+01Plasticity Indicator* at yield in shear or vol.X elastic, at yield in pasto at yield

47、in tensionGrid ploti i u n 05E 0Boundary plotI. L li ii ii 05E 01.800Itasca Consulting Group, Inc.Minneapolis, Minnesota USA2.0002.2001.6001.4001.2001.0000.8000.600(*10U)图17 方案9顶煤破碎状况24JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 19:07step 162721.858E+01 <x< 3.525E+013.282E+00 <y< 1.99

48、5E+01Plasticity Indicator* at yield in shear or vol. X elastic, at yield in past o at yield in tensionGrid plotBoundary plotItasca Consulting Group, Inc.Minneapolis, Minnesota USA1.8001.6001.4001.2001.0000.8000.6000.400图20方案13顶煤破碎状况25JOB TITLE :.图20方案13顶煤破碎状况#JOB TITLE :.图18方案11顶煤破碎状况图20方案13顶煤破碎状况#J

49、OB TITLE :.图20方案13顶煤破碎状况#JOB TITLE :.JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 19:23step 201011.799E+01 <x< 3.509E+013.070E+00 <y< 2.017E+01Grid plot0 I I I 5E 0Boundary plot1IhhI105E 0Plasticity Indicator* at yield in shear or vol.X elastic, at yield in past o at yield in tensionIt

50、asca Consulting Group, Inc. Minneapolis, Minnesota USA1.8001.6001.4001.2001.0000.8000.6000.400图19 方案12顶煤破碎状况JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 21:53step 226481.799E+01 <x< 3.509E+013.070E+00 <y< 2.017E+011 1 11iI05E 0Grid plot1.8001.6001.4001.2001tthii05E 0Boundary plotPlasti

51、city Indicator* at yield in shear or vol.X elastic, at yield in past o at yield in tensionItasca Consulting Group, Inc.Minneapolis, Minnesota USA.IHk+lH-Sh-l1.0000.8000.6000.400("10n1)图20方案13顶煤破碎状况26JOB TITLE :.FLAC (Version 4.00)LEGEND19-May-03 20:00step 180642.077E+01 <x< 3.632E+013.681E+00 <