采空区开挖分析

采空区开挖分析第一页，共十五页，编辑于2023年，星期五*--------------------------------------------------------------------------*MASTERANALYSIS=STATICAUTOMATI=ATSINSITU-D=YES*KINEMATICSDISPLACE=largeINCOMPAT=no*ITERATIONMETHOD=FULL-NEWTONLINE-SEA=YESMAX-ITER=50*PORTHOLESAVEDEFA=YESFILEUNIT=60FORMATTE=NOINPUT-DA=1DISPLACE=YESVELOCITI=NOACCELERA=NOTEMPERAT=NO*TIMESTEPNAME=DEFAULT@CLEAR10.1**--------------------------------------------------------------------------*read01_建立地层模型.in*read02_建立采空区模型.in*read03_定义边界条件及荷载.in*read04_定义材料单元组及划分网格.in*read05_定义采空区单元生死.in**--------------------------------------------------------------------------*ADINAOPTIMIZE=YESFILE='01_Initial-Stress.dat'*databasesave01_Initial-Stress.idbprompt=no**-----------------------------------------------------------------------------00_采空区开挖模型.in命令流文件求解设置分别读入in命令流文件生成初始地应力的dat求解文件及idb模型文件第二页，共十五页，编辑于2023年，星期五*--------------------------------------------------------------------------*--------------------------------------------------------------------------** 不使用低速动力算法则无法计算收敛*MASTERMODEX=RESTARTTSTART=0.1*AUTOMATICTIME-STEPPINGMAXSUBD=100RESPS=YESRESFAC=0.005*TIMESTEPNAME=DEFAULT@CLEAR1 0.11 1.0**--------------------------------------------------------------------------*--------------------------------------------------------------------------*databasesave02_Construction.idbprompt=no*ADINAOPTIMIZE=YESFILE='02_Construction.dat'**--------------------------------------------------------------------------*--------------------------------------------------------------------------求解设置生成开挖采空区的dat求解文件及idb模型文件对于采空区开挖,由于煤层自稳性差且无支护,如果按照正常计算的话则无法收敛,此时需打开静力计算的低速动力算法选项。00_采空区开挖模型.in命令流文件（续）第三页，共十五页，编辑于2023年，星期五打开低速动力算法

对于采空区开挖,由于煤层自稳性差且无支护,如果按照正常计算的话则无法收敛,此时需打开静力计算的低速动力算法选项。

对于一些静力计算难以收敛或者是材料破坏的问题，此时可采用Low-SpeedDynamicDamping方式控制收敛； Dampingfactor会直接控制收敛过程，damping带来的振荡过滤后才是需要的平衡状态；第四页，共十五页，编辑于2023年，星期五有限元模型采空区第五页，共十五页，编辑于2023年，星期五有限元模型采空区第六页，共十五页，编辑于2023年，星期五网格连续性的控制对于采空区上下土层而言，其可以进行映射网格划分；对于采空区而言，其必须进行多次切割才能划分映射网格，故在此仅进行四面体网格划分，网格与上下土层不进行连续性控制，采空GlueMesh的方式与上下土层的网格进行网格连接。第七页，共十五页，编辑于2023年，星期五GlueMesh网格此处两个面不能透视，也即网格不连续，此满足GlueMesh对网格的要求第八页，共十五页，编辑于2023年，星期五网格连续性控制的技巧对煤层及采空区进行网格划分时，可把煤层和采空区的Body体做成一个Domain，这样的话在对煤层及采空区的Body进行网格划分时选择仅对Domain进行节点检查，这样就可以控制网格在煤层及采空区是连续的，在与其它土层相邻处网格是不连续的。第九页，共十五页，编辑于2023年，星期五01_Initial-Stress模型的Z方向位移云图第十页，共十五页，编辑于2023年，星期五01_Initial-Stress模型的最大主压应力云图第十一页，共十五页，编辑于2023年，星期五02_Construction模型的初始地应力场阶段的Z方向位移云图初位移得到了清零第十二页，共十五页，编辑于2023年，星期五02_Construction模型的初始地应力场阶段的最大主压应力云图初应力得到了继承第十三页，