复合顶板受力情况分析

复合顶板情况分析根据煤层综合柱状图：在UDEC中初步建立80*160的工作面，按照上图的煤层分布，设计成共计10层的工作面，具体如下：（在己15和己16之间有斜面的砂质泥岩，开采处厚度为0.2M，太薄，所以不考虑）为了模拟现实情况，开挖高度为1.4米，横向距离从30米到40米，40米到60米，60米到80米，80米到100米，100米到130米，130米到160米共计6个测点。

下面是程序运行2000步时复合顶板受力情况图：块体在X方向的位移图：块体在Y方向的位移图：应力分析图：应力云图：塑性区变形：具体程序如下：ro 0.1set ovtol 2bl 0 0 0 80 160 80 160 0;建立模型cr 0 20 160 20 ;1cr 0 30 160 30 ;2cr 0 32 160 32 ;3cr 0 37 160 37 ;4cr 0 38.4 160 38.4 ;5cr 0 40.4 160 40.4 ;6cr 0 43.4 160 43.4 ;7cr 0 53.4 160 53.4 ;8cr 0 56.4 160 56.4 ;9cr 0 80 160 80 ;10jreg id 1 0 0 0 20 160 20 160 0 deljset 0,0 160,0 0,0 5,0 0,0 range jreg 1 jset 90,0 5,0 5,0 8,0 0,0 range jreg 1 jset 90,0 5,0 5,0 8,0 3 5 range jreg 1jreg id 2 0 20 0 30 160 30 160 20 del jset 0,0 160,0 0,0 2,0 0 20 range jreg 2 jset 90,0 2,0 2,0 3,0 0,20 range jreg 2 jset 90,0 2,0 2,0 3,0 1 22 range jreg 2jreg id 3 0 30 0 32 160 32 160 30 deljset 0,0 160,0 0,0 1,0 0,30 range jreg 3jset 90,0 1 0 1,0 2,0 0 30 range jreg 3jset 90 0 1 0 1 0 2 0 0.5 31 range jreg 3jreg id 4 0 32 0 37 160 37 160 32 deljset 0,0 160,0 0,0 2,0 0,32 range jreg 4jset 90,0 2,0 2,0 3,0 0,32 range jreg 4jset 90,0 2,0 2,0 3,0 2 34 range jreg 4jreg id 5 0 37 0 38.4 160 38.4 160 37 deljset 0,0 160,0 0,0 0.75,0 0,37 range jreg 5jset 90,0 0.75,0 0.75,0 1,0 0,37 range jreg 5jset 90,0 0.75,0 0.75,0 1,0 0.5,37.75 range jreg 5jreg id 6 0 38.4 0 40.4 160 40.4 160 38.4 deljset 0,0 160,0 0,0 1,0 0,38.4 range jreg 6jset 90,0 1,0 1,0 2,0 0,38.4 range jreg 6jset 90,0 1,0 1,0 2,0 0.5 39.4 range jreg 6jreg id 7 0 40.4 0 43.4 160 43.4 160 40.4 deljset 0,0 160 0 0,0 1.5,0 0,40.4 range jreg 7jset 90,0 1.5,0 1.5,0 2,0 0 40.4 range jreg 7jset 90,0 1.5,0 1.5,0 2,0 0.6 41.9 range jreg 7jreg id 8 0 43.4 0 53.4 160 53.4 160 43.4 deljset 0,0 160,0 0,0 5,0 0,43.4 range jreg 8jset 90,0 5,0 5,0 8,0 0,43.4 range jreg 8jset 90 0 5 0 5 0 8 0 3 48.4 range jreg 8jreg id 9 0 53.4 0 56.4 160 56.4 160 53.4 deljset 0,0 160,0 0,0 1.5,0 0,53.4 range jreg 9jset 90,0 1.5,0 1.5,0 2,0 0,53.4 range jreg 9jset 90,0 1.5,0 1.5,0 2,0 0.5,54.9 range jreg 9jreg id 10 0 56.4 0 80 160 80 160 56.4 deljset 0,0 160,0 0,0 1.2,0 0,56.4 range jreg 10jset 90,0 1.2,0 1.2,0 2,0 0,56.4 range jreg 10jset 90,0 1.2,0 1.2,0 2,0 0.8 57.6 range jreg 10gen quad 6 range 0,160 0,80zone model mo range 0,160 0,80change jcons=2 range 0,160 0,80;给岩体岩块赋属性编号chang mat=1 range reg 0 0 0 20 160 20 160 0chang mat=2 range reg 0 20 0,30 160 30 160 20chang mat=3 range reg 0 30 0 32 160 32 160 30chang mat=4 range reg 0 32 0 37 160 37 160 32chang mat=5 range reg 0 37 0 38.4 160 37 160 38.4chang mat=6 range reg 0 38.4 0 40.4 160 40.4 160 38.4chang mat=7 range reg 0 40.4 0 43.4 160 43.4 160 40.4chang mat=8 range reg 0 43.4 0 53.4 160 43.4 160 53.4chang mat=9 range reg 0 53.4 0 56.4 160 56.4 160 53.4chang mat=10 range reg 0 53.4 0 80 160 80 160 53.4chang jmat=1 range reg 0 0 0 20 160 20 160 0chang jmat=2 range reg 0 20 0,30 160 30 160 20chang jmat=3 range reg 0 30 0 32 160 32 160 30chang jmat=4 range reg 0 32 0 37 160 37 160 32chang jmat=5 range reg 0 37 0 38.4 160 37 160 38.4chang jmat=6 range reg 0 38.4 0 40.4 160 40.4 160 38.4chang jmat=7 range reg 0 40.4 0 43.4 160 43.4 160 40.4chang jmat=8 range reg 0 43.4 0 53.4 160 43.4 160 53.4chang jmat=9 range reg 0 53.4 0 56.4 160 56.4 160 53.4chang jmat=10 range reg 0 53.4 0 80 160 80 160 53.4;赋值prop mat=1 dens=2683prop mat=2 dens=2673prop mat=3 dens=1400prop mat=4 dens=1456prop mat=5 dens=1773prop mat=6 dens=2500prop mat=7 dens=2683prop mat=8 dens=1673prop mat=9 dens=2400prop mat=10 dens=2456zone k=24e9,g=11e9,fric=29,coh=2.42e6,ten=3.1e6 range mat=1 zone k=8e9,g=6e9,fric=15,coh=1.87e6,ten=2.4e6 range mat=2 zone k=5e9,g=1.92e9,fric=50,coh=4.77e6,ten=0.70e7 range mat=3 zone k=6e9,g=4.5e9,fric=31,coh=3.6e6,ten=1.0e6 range mat=4 zone k=24e9,g=14e9,fric=10.00,coh=2.42e6,ten=4.2e6 range mat=5 zone k=3e9,g=2e9,fric=13.00,coh=1.2e6,ten=0.7e6 range mat=6 zone k=24e9,g=11e9,fric=21.00,coh=6.42e6,ten=3.1e6 range mat=7 zone k=8e9,g=6e9,fric=22.00,coh=1.87e6,ten=2.4e6 range mat=8 zone k=5e9,g=1.92e9,fric=20,coh=1.77e5,ten=7e7 range mat=9 zone k=6e9,g=4.5e9,fric=31,coh=1.6e6,ten=1.0e6 range mat=10prop jmat=1 jkn=12e9,jks=6e9,jcoh=1.87e5,jfric=20,jten=2e6 prop jmat=2 jkn=1.5e9,jks=1e9,jcoh=0.8e5,jfric=15,jten=1.6e5 prop jmat=3 jkn=1.4e9,jks=1e9,jcoh=1.5e5,jfric=11,jten=1e6 prop jmat=4 jkn=2e9,jks=1.3e9,jcoh=5.80e5 jfric=20,jten=1.2e6 prop jmat=5 jkn=7.2e9,jks=5e9,jcoh=2.2e5,jfric=15,jten=3e6 prop jmat=6 jkn=2e9,jks=1.5e9,jcoh=0.8e6,jfric=20,jten=0prop jmat=7 jkn=12e9,jks=6e9,jcoh=1.87e5,jfric=30,jten=2e6 prop jmat=8 jkn=1.5e9,jks=1e9,jcoh=0.8e6,jfric=6,jten=1.6e6 prop jmat=9 jkn=1.4e9,jks=1e9,jcoh=0.5e3,jfric=5,jten=1e5prop jmat=10 jkn=2e9,jks=1.3e9,jcoh=1.80e5 jfric=14,jten=1.2e5set gravity 0,-10bound stress 0,0,-5.75e6 range 0 160 78.8 80.2bound xvel=0 range -0.1,0.1 0,80bound xvel=0 range 159.9 160.1 0,80bound yvel=0 range 0,160 -0.1,0.1step 2000resetsave 1.savset pline 0 38.4 160 38.4 81set pline 0 37 160 37 81hist ydisp 120 37hist ydisp 120 38.4；多加几个测点del range 30 130 38.4 40.4step 5000sav 2.savdel range 30 40 37 38.4step 2000save 3.savdel range 40 60 37 38.4step 2000save 4.savdel range 60 80 37 38.4step 2000save 5.savdel range 80 100 37 38.4step 2000save 6.savdel range 100 130 37 38.4step 2000save 7.savdel range 130 160 37 38.4 step 2000save 8.sav。