巷道围岩初始应力场和弹性模量的区间反演

1巷道围岩初始应力场和弹性模量的区间反演方法摘要从有限元数值计算的角度对岩土工程参数的可辨识性进行了分析。在可辨识性分析的根底上，考虑观测数据的不确定性，引入不确定问题求解的区间分析思想，建立了联合采用围岩位移和应力观测数据同时识别地应力和围岩弹性模量的反演模型。并采用约束变尺度方法求解该反演模型，获得地应力和弹性模量所在区间范围。关键词区间分析，约束变尺度方法，初始应力场，弹性模量1．引言输入参数〔岩体力学参数、初始应力场〕和岩体本构关系给不准是岩土工程数值模拟的两个瓶颈问题[1]。采用位移等观测数据推断岩体弹性参数或初始地应力，这类弹性问题的位移反分析，由于公式推演和计算过程的实施较为简单，且所得结果通常能较好地近似反映实际工程问题的主要特征[2]，70年代以来受到人们的重视[2-6]。实际工程中，地应力和围岩弹性参数一般均是未知的，寻找能够同时确定它们的方法更具有实用性。文[3]曾就同时辨识围岩弹性参数和地应力可行性进行过研究，杨林德等率先将位移和扰动应力增量量测信息同时用作反演分析的依据，提出了可同时反求初始地应力参数和地层弹性参数的位移、应力反演分析法[2]；文[4]曾在引入所谓“预估约束〞的条件下同时识别弹性参数和地应力。但是，由于观测误差客观存在的不确定性，据此推演出的岩土工程参数也应为不确定的，所以，反演问题本质上是一个不确定问题，本文从该角度对此类问题进行研究。从有限元数值计算的角度，本文首先就同时辨识围岩弹性参数和地应力可行性进行了探讨。然后在可辨识性分析的根底上，采用区间分析方法的思想建立了联合利用位移和应力观测数据同时反演围岩弹性模量和地应力的区间反演计算模型。在该模型中，用有限元法建立系统的控制方程，把待求变量看作区间数，以量测信息作为设计变量的边界约束，运用约束变尺度方法求出地应力和弹性模量的最大值和最小值。2．同时辨识巷道围岩弹性模量和初始地应力的可行性分析当巷道埋深较大时，其受力情况可以简化为图1示。巷道开挖所引起的巷道围岩位移和扰动应力可以看作是，在水平压力P=σ0x和垂直压力Q=σ0y作用下与巷道形状相同的弹性孔的位移和应力减去没有开挖巷道前岩体已经完成的位移和应力。即相当于在巷道开挖边界反向作用与地应力大小相等的边界力时，所引起的围岩的位移和应力〔为作图方便，以开挖边界为矩形为例如图2示〕。采用有限元计算时，假设初始地应力作用下开挖边界上的作用力为f0，开挖引起的位移和应力增量相当于在开挖边界上反向施加等效的节点力-f0引起的围岩围岩和应力值[6]。(1)慧τ嗓0闭xy温τ劫0芝xy斧σ脸0馒x孝σ赚0跑y表P=疗σ幼0盏y牢Q=常σ尼0屋x掠y僵x

惯闭榨婆驻哭穗参图豆1慌醋异地应崭力作霞用下训的巷冤道井淹服傲万摔瘦降燕俊钥拣图榴2哥地下壁巷道窄开挖疮释放差荷载稳靠秤绵姿访裙F洗ig璃1赛幕A辣tu妈nn类el炼a存ct搞ed枕b络y摄in景s价it惯u脉st扭re节ss筹苦钞些吓F唱ig拥2等演Re抖le妙as擦ed境l注oa吐d粮du拦e煎to斥e低xc勉av栏at恭io染n向of炉t蛮un僻ne形l余式中城：桨趟为单垄元形费函数奏矩阵淘；晶饮为几数何矩骑阵；橡应为岩叫石的仍自重目体积挑力矢抚量，单当巷轮道埋悲深较塑大时要可近既似认菌为钞醉=0政；纤骗为边肥界力乘，怖菠=(朵朴)损；隙穷为单外元的自厚度吐；桐毛是具垮有开都挖边串界的月所有筑单元但。有纠限元凉平衡店方程筐为秋炒电梅筋肾接壶腰锯骑甩侍忆希帐椒(刻2)战其中盖疑，挺栗，拣龄。节甘点应迹变和糠应力恳通过艺式〔庆3蝶〕求倒出单聪元应姜变和柄应力馋值并哭采用奋某种歉光滑亮措施检〔一泊般相哲当于适乘一编常数咸矩阵怕〕后另获得叫。津猾壶能料萄绪瘦滴填具垦茫(3熊)偷秘执由于馅弹性炭矩阵绒D呈与弹符性模就量成绍正比侧，僚即与初敌始地鸣应力须成正筹比，交应变匠矩阵压B理仅与乞单元鼓几何轨特性烤有关蚁，由矛式〔到2敬〕知拾对于承任意摆常数横法，参厌数组劲端、讽考、虽假设，庭愈，都袭可以拉与一寻个位队移分妙布训吵相对迷应，竭因此幻，仅油采用呜位移摄观测嗽数据引来确赚定弹贼性模长量和捉初始更地应第力场存其解袭是不考唯一固的，榨只能背确定密地应泽力与皮弹性收模量恼的比疑值。蒙再由联公式孙〔捏3兄〕可伟见，嗽一点之的应价变值绑与其姨位移厉值二替者是欠相关期的，基同样罢可以豆有无绣数组雹〔甚E惧,畏σ帝0搭x早,歪σ匹0住y掘,秋τ借0虽xy已〕与滚一个救应变淡分布炸窗相对敏应，赵即使掀同时鸦采用特位移坊和应催变观盾测信往息也保不可卫能改包变解铁的唯答一性马。但趟是，绵由于哑应力熊分布主箩与位暖移分败布戚脂通过企弹性联模量稻妥联系姓起来条，并悬且由眼应力谢分布秒愈与位夜移分怨布蕉耐可望碎唯一可确定杯弹性京模量稼抹，因宰此，饱如果挎有扰判动应假设力观薯测信茶息时杨，那么掠可望扶唯一英确定目所有糖待求绕变量控。实谎际上逝仅采犁用应触力观膨测信诚息，肠或联进合采树用应杆力、瓜位移垫观测夜信息杯，或秤联合只采用梯应变葡、应氧力观猜想信厘息都移可望裙到达劣这种行辨识推的目执的。边当观雪测位到移信拿息足贿够多店时，蓄哪怕观只增推加一诞个应们力观粘测信管息，备便可搞望唯迁一确芝定围哑岩弹成性模夺量和明地应蛛力。撇考虑掩到反晚演中索采用瞎的观误测量踪，应体该是而可观避测、袋易观奇测的宰，且末观测笔精度固愈高顽愈好档等要洪求，萍而与泽应力颤观测欣相比罚，一阵般而教言，勺位移屑观测闪更易状于实岁施、券精度皆高等葬特点替，本艘文考藏虑联负合使雕用位狗移和意应力乐观测怪信息点来同报时识扩别位占移弹碍性模痛量和克初始惨地应沟力。奋3．防

愉基于锡区间强分析江思想神的反获演模宇型的隆建立点及其号求解勉方法进的选栋择屡甚众所渐周知穗，计吵算误专差一都直是闯数值误分析扁中一驴个比拢较麻询烦的饭问题遥，它利来源努于数饲据误琴差、舟截断执误差抛和舍闲入误静差。俊人们中努力嫂使计达算结而果能沈在所样要求抗的精晴度内纠，然春而，悲在许肯多问市题中拥，往即往是酿推测邻计算喷结果疮的某鄙种精贤度或激者使卫用高惯精度压的运辈算以槽保证胸计算构结果岁的精阴度。旱但是甩，由付于计捎算误拥差的香积累茶，可卖能使末计算州结果恳失去卷意义究。巷站道施股工期犯间的奇监控变量测律，由角于仪衫器精差度和丈其它套偶然贼误差创，所胸测变样形是狸一随涛机值泄，因耐而反疮推的脱围岩驼参数率和初奸始地课应力萍也是刷随机宪变量址。当仅根据跟带误趣差观逗测数忙据采金用确悠定性居方法蚊进行内反演帽时，谅反演箩结果娱必然渴不稳沈定，来而且己观测榴数据蹲的少勇量误礼差常幻可以借导致臂反演虽结果元有大讨幅度姿的波撑动参[2鸡]嫩。广随机梢理论乞、模驴糊集厌理论糠和区翅间分察析是巴解决司不确说定性愿问题税的三观种方视法驱[6萍-7辩]魂。用殿随机既理论有或模驳糊集绳理论抖求解技问题仅时需野要知鼓道不娇确定谋参数摆的概汽率密弯度函顶数或殃隶属析函数粘。而蔽实际揪中这麻些函结数往川往难报以确为定，爹常常溉通过燥人为滔方式衣取舍机。为认了真刑实地单反映跑客观林实际份，减奔少人文为因凶素的绪影响拥，提信高计功算结杯果的蚊可靠误性，烧运用设区间坦分析钥求解番不确见定问燕题，纯粹是轮这一鸦思想是的体讨现，迎它提估供了赔一种尤简便果的方射法。代考虑栋到各湿种计血算误弯差，惠作为腔计算绝结果槽，区剪间分化析方缝法得趴到的救是包任含其滚精确甩解的暑一个膊区间染和误惯差界团，这挽就为两人们躺能更编准确撞地估兆计和透评价嚷所得侵的结败果提泡供了客一定瓦的依绵据。缘近年朗来，截采用鸽区间练分析绿求解远不确障定问漠题的末方法友已经禽有所瘦研究板。区蜜间数族与普欲通数池不同姑，四廉那么运敏算规乡那么也资完全鸣不同曲，对捕于两仇个区马间数脆仁和掏芽，四遵那么运搂算规阅定为锣：懒贤，澡奉，甚匀，瑞奖楼，区脾间数今运算魄容易秃导致范区间午扩张渠。在谜用区尾间分恢析求复解不祖确定新问题帅时，怖有许占多新悟问题蝶需要悼解决荷，如痰如何益保证辩白区近间不阶被放本大或绵缩小馒。文闯[9迅]利提出弯采用江直接乱优化块法进澡行非刑确定绣结构根系统雄区间役分析扛，算透例表窑明采碧用直旦接优六化方台法与摘穷举俩组合旅法所忍得结钓果相常同，刑而直仓接采区用区乏间数瘦运算酸规那么蛮所得勺结果跨区间暮有较究大扩沫张，般这里播仍然那采用拌这种损方法张处理佛本文对问题闷。络当把醋待求企变量赞〔某E建,夹σ况0稍x得,此σ俯0阳y吓,割τ鸡0漏xy园〕作颂为一最个区简间数姻时，茎相应尊地采暗用有俊限元绣法要胃区间毁的区两间线匙性结躁构静巴力学炼方程范形式替为穿鸽被匆席遣涂捆扔守送吐栋巷悬少够朽(满4)石与公细式恳(2三)爸不同蒜的是杨式抚(4漠)袍中蹈镰、金宾、满为都是秀区间述矩阵铁和区孤间向燃量，搬其元挂素用苗区间池数表舞示为别，棋抱，侵闸，传吗。宵在已论知位河移和完应力是观测傻精度丝的情遍况下冷，由文位移鄙和扰黄动应顽力观凭测信不息可途得含港植口宵江突掠矩柄收懂腰勉段踩袍至防.奶皂弓瞒伍支赞该倾旺冈粉犹豪滋(5形)愿式中骗：绸悔、押敢分别御是第秆歪个观歌测点识位移沉值和烈第呼豆个应聋力观踢测点敌扰动潜应力冻值；婶别、岁膨分别精是相否应观奇测点蹄的位集移观拾测值加和扰驰动应肚力观把测值不；暴±举升和荡±刑补分别便相应露观测传点的激位移金观测翠精度漂和扰旗动应灾力观蚀测精线度；雹令、愁撑分别链是位陆移和票扰动搁应力照观测对数量感。假设歪再考平虑对床于实谊际工洲程问君题，肌一般蹦都能竖给出币如式惰(6佳)阅所示奸的待瞎求变弃量较恋为宽及松的投上下仙界限推，以殿保证忙问题炭的解赛与实骑际物拉理意共义相洞符。大,指,眯,烈稿星古拢请脾(6害)丢从而挪，可激以建尽立同静时确燃定围狸岩参袍数和偶原始窜地应由力的窝区间回直接刃优化章计算婆模型口：究谢躲越岛声毁启阻茫务M英in层/抄M冻ax酷速筋谷判借死清滤指债s.孙t.误粥珠座宝啄雁合惑与锄不四托谎裕区财再嗽换虾圆喘留犁(7嘉)阳跑杂屋牺可归兆穷雪挎懒虎爬裹婶丽壤灵纤.俊通兰骑鸦天南技奥惰死眯弟尚,掠,妖,酷漏陷模型顾(7过)换是一肥个复案杂的毕约束勉非线揉性优包化问枕题，橡本文婶采用温约束胁变尺山度方怕法求写解。告约束辞变尺盲度方德法的终的基诚本思什想是辞迭代氏和逼葱近。敌 忆约束脉变尺曲度方慌法具遗有收枕敛快狗、效编率高隆、可云靠性时好与瓶整体坝收敛略性好舅、适激应能贺力强厚等一辛系列线优点箱。求袖解模巾型形(7逮)着的约广束变霜尺度辨方法茅程序泽是在殿Po棕we悉ll哈型约瞧束变拍尺度胞方法壳的基圣础上奥，懂参考目文献召[1研1]离采取成的监邀控搜话索技剂术、凭尺度识矩阵暴保持珠正定松等改恒进措便施编鸦制而垄成。非4．造数值修算例劝乌由疏某地丸下巷男道，假设围岩丘泊桑告比瓦u滤=0兆.2础，假泛设初杠始地竖应力通场均爬匀分萍布，牙且剪粉应力握τ春0产xy晕=0族〔剪烫应力帽τ铃0讽xy膛不参柿加反完演〕愉，由泛观测面位移段值和肺扰动疼应力闪观测冷值，泡确定锄围岩签弹性仁模量很E才和初裕始地奔应力歌分量穗σ续0傅x淹荣，艳σ奖0恭y材。浙公1促2臭①雕②果③攻数值谅模拟刻时首救先选颤取一辽组参授数值手(章E纹,锄σ虎0堂x俭,吉σ克0瞧y类)=侍(6丰0架00逗0,坟1吉0,翠4葛0)锦MP驰a何，用裁有限折元进窜行正究演计强算，驴得到煎图示铺测线糖①眠、尤②粪、坚③望的收页敛位捏移值蜻，以怎及观锁测点瞎1染和点刷2扔的水责平方修向应得力扰炼动分素量〔约见表绣1轿〕。此然后饺把这库这映4挤个收往敛位岛移值访和扰妨动应北力值款作为泥观测顾值，控并考庭虑观床测仪母器的惑精度拌确定俯围岩例弹性元模量您E恒0邮和原案始应投力分蚕量长σ倘0融x起,狐σ寻0笼y轮所在疮区间摇范围裳。用煌有限窃元法榆计算普时，能由于石对称恳，只馆计算染半个未区域绳，共因分为恼26欺5甘个三雁角形塘单元槐。各景参数茂上下垂限取状集=1尺0灰00寇0M愉pa叮，拌护=9桐0什00役0M拌pa立，勒价=罩桐=2柜M耻Pa牙，贷,念册=祸欲=6注0M盼pa疫。计漂算结豪果如绒表考2邮所示三，对极于某木一给唯定初伸值在泡奔腾员13绵3互微机革上运墓行计滩算时近间为州40室~7纱0s最。喷明茫透状铲

蜘图悼3系浙瑞测点映布置驰跳内Fi福g少3慢A劣ll等oc欺at逝io警n石of烫o毯bs授er帆va蔽ti草on辨s使表费1停位移旷和扰森动应喜力观音测数抗据廊Ta膀bl先e苍1忆Me透as撒ur挠ed役d麻at炊a痰of鲁d晨is裳pl南ac歇em尚en样ts剂a昆nd使d糟is捡tu佳rb瞒ed柜s娘tr煎es筒se理s械置轮毯止观航测症位队置录尖源售猜标垦趟①懂涨煎室献②责益碰甚胆③绣匪朱辰沫俯悟1净绪颤与惹岭2猎宏萍缘知收赌合敛菊竿位杰倦移趁/庭10效-3冒m勾膀丛筋低偏熔冤附挤丑6.锋00晚92湿讽去惰3.晋06汤38似6型绵副莲-2后.3钩52红80从阀胃葬愤晒扰动躁应力硬σ涝x婆/伯MP牧a渗仍尤拨激田锣茫蜜扶朋哪材茅芳疾吃伸罚半杆盛勺爸质敌毙萝丈恢至救猫攻策调送躺唯趋罚肯芽响必鼓胁终8.快64亏51招00辉挑臂1株1.络02崇07姻0忘表爸2动类反演盼识别彩结果忠Ta铃bl族e裕2槐E飘st颜im贤at勿ed胸r怨es毯ul孩ts柴o疑f驳in呀ve缴rs思e挎an鼠al傲ys带is秀损须使妹/1吃0抱3裂MP指a军售毒之乳σ含0哨x贞/狠MP洪a蛙侮直宝族宴σ山0南y贷/M脂Pa凉乘精确思值劫勺逆雁欺摆汽辈粮6订0.基00我0墨卡窝村沈荷忍危割孙杀催尾痛思1堆0.败00坏0浓外铲驳孝充讽烦还炸激程呀柄40裕.0扁00举留胃街强识别越值痕(1久)足即横白朽萌5每9.给74呼81慰2~伍62忠.6狡74勾81舰底捆表咱9.注33查40齿1~尘10画.9稻31咐87悄抖毯薯欠沈38暴.4园36圈05惹~4腰2.至40任22愈8欠

菌语晨归识别浑值爷(2亮)梳豪奋撤授顷5窜7.征05设52剖3~赌63涛.1悉06件62损筒侧菊8.膝99啦39元4~步11这.0戚50阳59帆谁