采空區(qū)開挖分析第一頁，共十五頁，編輯于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_建立采空區(qū)模型.in*read03_定義邊界條件及荷載.in*read04_定義材料單元組及劃分網(wǎng)格.in*read05_定義采空區(qū)單元生死.in**--------------------------------------------------------------------------*ADINAOPTIMIZE=YESFILE='01_Initial-Stress.dat'*databasesave01_Initial-Stress.idbprompt=no**-----------------------------------------------------------------------------00_采空區(qū)開挖模型.in命令流文件求解設(shè)置分別讀入in命令流文件生成初始地應(yīng)力的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'**--------------------------------------------------------------------------*--------------------------------------------------------------------------求解設(shè)置生成開挖采空區(qū)的dat求解文件及idb模型文件對于采空區(qū)開挖,由于煤層自穩(wěn)性差且無支護,如果按照正常計算的話則無法收斂,此時需打開靜力計算的低速動力算法選項。00_采空區(qū)開挖模型.in命令流文件（續(xù)）第三頁，共十五頁，編輯于2023年，星期五打開低速動力算法

對于采空區(qū)開挖,由于煤層自穩(wěn)性差且無支護,如果按照正常計算的話則無法收斂,此時需打開靜力計算的低速動力算法選項。

對于一些靜力計算難以收斂或者是材料破壞的問題，此時可采用Low-SpeedDynamicDamping方式控制收斂； Dampingfactor會直接控制收斂過程，damping帶來的振蕩過濾后才是需要的平衡狀態(tài)；第四頁，共十五頁，編輯于2023年，星期五有限元模型采空區(qū)第五頁，共十五頁，編輯于2023年，星期五有限元模型采空區(qū)第六頁，共十五頁，編輯于2023年，星期五網(wǎng)格連續(xù)性的控制對于采空區(qū)上下土層而言，其可以進行映射網(wǎng)格劃分；對于采空區(qū)而言，其必須進行多次切割才能劃分映射網(wǎng)格，故在此僅進行四面體網(wǎng)格劃分，網(wǎng)格與上下土層不進行連續(xù)性控制，采空GlueMesh的方式與上下土層的網(wǎng)格進行網(wǎng)格連接。第七頁，共十五頁，編輯于2023年，星期五GlueMesh網(wǎng)格此處兩個面不能透視，也即網(wǎng)格不連續(xù)，此滿足GlueMesh對網(wǎng)格的要求第八頁，共十五頁，編輯于2023年，星期五網(wǎng)格連續(xù)性控制的技巧對煤層及采空區(qū)進行網(wǎng)格劃分時，可把煤層和采空區(qū)的Body體做成一個Domain，這樣的話在對煤層及采空區(qū)的Body進行網(wǎng)格劃分時選擇僅對Domain進行節(jié)點檢查，這樣就可以控制網(wǎng)格在煤層及采空區(qū)是連續(xù)的，在與其它土層相鄰處網(wǎng)格是不連續(xù)的。第九頁，共十五頁，編輯于2023年，星期五01_Initial-Stress模型的Z方向位移云圖第十頁，共十五頁，編輯于2023年，星期五01_Initial-Stress模型的最大主壓應(yīng)力云圖第十一頁，共十五頁，編輯于2023年，星期五02_Construction模型的初始地應(yīng)力場階段的Z方向位移云圖初位移得到了清零第十二頁，共十五頁，編輯于2023年，星期五02_Construction模型的初始地應(yīng)力場階段的最大主壓應(yīng)力云圖初應(yīng)力得到了繼承第十三頁，