空间轴对称极限平衡应力数值分析.docxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
空间轴对称极限平衡应力数值分析

1998年8月JOURNALOFCHINACOALSOCIETYAug.1998空间轴对称极限平衡应力数值分析李国强袁辉韩光立苏立凡(煤炭科学研究总院北京建井研究所)摘要以莫尔-库仑准则和全塑性假设为基础的空间轴对称问题的基本方程化为无量纲形式,利用比较精确的应力方程数值积分方法对空间轴对称问题的基本方程作了数值解,编制了程序CALC,在此基础上,将立井围岩视为松散介质,进行模拟计算,得出了井壁压力不是直线型的,而是曲线型的;别氏公式在浅部适用性较好,而在深部不适用等结论,并对滑移线形状进行了描述.关键词松散介质全塑性假设应力方程数值积分方法极限平衡中图分类号TD321O24114松散介质空间轴对称极限平衡问题的基本方程是由前苏联著名学者B1Γ1别列赞采夫首先提出的,由于直接求解该方程的解析解的难度较大,因此,B1Γ1别列赞采夫将该方程的滑移线假设为直线,在此假设的基础上对空间轴对称问题的基本方程做了求解,所得结果即为著名的别氏公式.由于别氏公式将滑移线视为直线,因此限制了解的应用范围1.80年代后,由于计算机的飞速发展,对于空间轴对称问题的基本方程,用数值方法进行模拟计算成为可能,梁力群2曾用有限差分方法对空间轴对称问题的基本方程做了数值求解,但由于差分方法累计误差较大,因此计算结果不够精确,同时,差分方法无法对滑移线进行详细描述.本文采用了一种比较新的数值方法———应力方程数值积分方法———来求解空间轴对称问题的基本方程.与有限差分方法相比较,该方法从理论上讲精确度比差分方法要好的多,不仅求出应力分量,而且能够对滑移线进行详细描述,作出滑移线图.在此基础上,将立井围岩视为松散介质,利用该方法做了数值求解,得出了一些对立井设计、施工有一定指导作用的结论.立井空间轴对称问题基本方程的无量纲化立井空间轴对称问题在岩体子午面内的基本方程3为15σφ5δσtgφ(1+sinφ)[sin(δ+φ)-cosδ]ρsin(δ+φ)σ5S1+2tg5S1+=,Rcosφσtgφ(1+sinφ)[sin(δ+φ)cosφ(1)5σφ5δ-cosδ]ρcosδσ5S2-2tg5S2-=-cosφ,Rcosφ式中,σ=(σ1-σ3)2sinφ=(σR+σZ)2+Cctgφ;φ为岩体的内摩擦角;C为岩体内聚力;ρ为岩体密度;S1,S2为岩体子午面内的特征线.式(1)中未知数为σ,δ,且为双曲线型的一阶拟线性偏微分方程组(图1),特征线的物理意义即为滑移线,两组特征线分别为S1和S2线,这一特性为数值积分方法提供了很大的方便.求解之前收稿日期:1996-12-17煤炭学报1998年第23卷358先将式(1)化为无量纲形式.由图1可以看出,滑移线方程dZdR=tgδ,dZdR=tg(δ-π2+φ).沿S1方向沿S2方向将式(1)改写为(2)5σ5S1+2σtgφ5δ5S1+σtgφ[sinδ-ρsin(δ+φ)cosφ,cos(δ+φ)]R=(3)5σ5S2-2σtgφ5δ5S2+σtgφ[cos(δ+φ)-ρcosδcosφ.由图2可以得到sinδ]-R=图1岩体子午面内点的应力和滑动面Fig11StressandslidingplaneofapointinthemeridianplaneofrockmassdRdS1=cosδ,dZdS1=sinδ,(4)dRdS2=cos(δ-2μ),dZ式中,μ=π4-φ2.令r=RρC,z=ZρC,Ω=σC,dS2=sin(δ-2μ),通过量纲分析可以看出,r,z,Ω均为无量纲形式,则R=rCρ,Z=zCρ,σ=ΩC.(3),并令dZ=dz,dR=dr,R=r,dσ=dω,σ=ω.将其代入式可以得到(2),dZ=dRtgδ,dZ=dRtg(δ-2μ).沿S1方向沿S2方向(5)则沿S1方向dσ+2σtgφdδ=(tgφ+σsinφR)1-σtgφ(1+sinφ)RdZ,沿S2方向dσ-2σtgφdδ=(-tgφ+σsinφR)1+σtgφ(1+sinφ)RdZ.dR+(6)图2计算转化Fig12ConversionincalculationdR+式(5),(6)即为立井空间轴对称问题的基本方程.立井空间轴对称问题的数值方法求解Masseau4提出了对平面问题双曲线方程的最简单解法.2将Masseau提出的方法通过下列运算应用到立井空间轴对称问题的数值方法求解中.211域内求解如图3所示,已知在子午面内A,B两点的σ,δ值,欲求P点的σ,图3域内求解的特征线Fig13Characteristiclinewhensolvingtheequationwithinthefieldδ值.P点是通过A点的α线(S线)和B点的β线(S线)的交点.可12以看出,P点的位置未知,需在运算中确定.很明显,若不知σ,δ的精确变化,

文档评论(0)

pangzilva + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档