Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory

Technische Samenvatting: Multiscale Berekening van de Equivalent Elasticiteitsmodulus voor Rotsmassa Gebaseerd op Perturbatietheorie

Probleemstelling
De nauwkeurige bepaling van de equivalente elasticiteitsmodulus is cruciaal voor de veiligheids- en stabiliteitsbeoordeling van gesteentebouwtechnische constructies. Huidige methodologieën kampen met significante beperkingen: laboratoriumtests op intacte gesteenteproefstukken houden geen rekening met de macro- en micro-naden en scheuren die inherent zijn aan engineering-rotsmassa's, terwijl veldtesten vaak niet representatief zijn vanwege hun beperkte ruimtelijke omvang. Bovendien zijn traditionele reductiemethoden die worden gebruikt om laboratoriumresultaten naar veldomstandigheden te schalen, afhankelijk van onbetrouwbare reductiefactoren en kunnen zij het mechanische responsproces tijdens het falen van de rotsmassa niet vastleggen. Er is een erkende behoefte aan een systematische, multiscale koppelingsbenadering die de kloof overbrugt tussen intacte gesteentemechanica en het heterogene gedrag van beschadigde rotsmassa's.

Methodologie
De auteurs stellen een multiscale berekeningsbenadering voor die gebaseerd is op perturbatietheorie, specifiek ontworpen om initiële schade in rotsmassa's te accounten. De methodologie behandelt de rotsmassa met initiële schade als een perturbatiesysteem waarbij een kleine parameter ($\delta$) het multiscale expansieproces van interne scheuren reflecteert.

Belangrijke methodologische stappen omvatten:

• Theoretisch Kader: De rek-tensor veroorzaakt door meso-heterogeniteit wordt gebruikt als een "brug" die macro- en microschalen verbindt. Het verplaatsingsveld wordt onderworpen aan asymptotische expansie, waarbij variabelen worden gescheiden in macroscopische coördinaten ($x$) en microscopische omgevingen ($y$).
• Beheersvergelijkingen: De studie leidt beheersvergelijkingen af door de evenwichtsdifferentiaalvergelijking te koppelen aan een constitutieve vergelijking die een micro-elastische schadefactor ($D^\delta$) incorporeert. De initiële schade wordt gerepresenteerd door een initiële elastische schade-tensor.
• Homogenisatie: Door de bewegingsvergelijkingen onder verschillende fysieke dimensieschalen op te lossen, leiden de auteurs de gehomogeniseerde macro equivalente elasticiteitsmodulus ($E_{mnkl}^\delta$) af. Dit omvat het oplossen van tensor eigenfuncties ($\chi$) met behulp van de eindige elementenmethode op basis van variatiemethoden.
• Numerieke Implementatie: De benadering werd geïmplementeerd met behulp van MATLAB voor iteratieve berekeningen van de schade-coëfficiënt en COMSOL voor eindige elementenmodellering. Het model werkt de belasting ($Q$) en verplaatsing ($u$) in elke stap bij om de evoluerende schade-toestand te reflecteren.
• Validatie: De benadering werd gevalideerd via twee primaire gevallen:
  1. Laboratoriumvalidatie: Uniaxiale compressie- en Brazil split tensile-testen werden uitgevoerd op kalksteenproefstukken verzameld uit een diepe metaalmijn in zuidelijk China. De proefstukken werden vooraf beschadigd tot variërende gradaties (0,25$\sigma_{cs}$ tot 0,75$\sigma_{cs}$) en gekarakteriseerd via P-golf snelheid en CT-scanning.
  2. Veldvalidatie: De methode werd toegepast op gegroutte gestolde lichamen in een lood-zink mijngang, waarbij het gegroute materiaal werd behandeld als een rotsmassa met onregelmatige naadscheuren.

Belangrijkste Bijdragen

• Multiscale Perturbatiemodel: Het artikel vestigt een theoretisch kader waarbij de rek-tensor voortvloeiend uit meso-heterogeniteit dient als het verbindingsmechanisme tussen micro- en macroschalen, wat de berekening van de equivalente elasticiteitsmodulus rekening houdend met initiële schade mogelijk maakt.
• Dynamische Schade-evolutie: In tegen tegenstelling tot conventionele modellen die de elasticiteitsmodulus vaak als een constante behandelen, reflecteert deze benadering de dynamische ontwikkeling van de equivalente elasticiteitsmodulus en de evolutie van zowel macro- als microschade tijdens het belastingsproces.
• Verbeterde Nauwkeurigheid: De studie toont aan dat de voorgestelde multiscale berekening resultaten oplevert die dichter bij werkelijke laboratorium- en veldtestgegevens liggen vergeleken met conventionele eindige elementen simulatiemethoden.

Resultaten

• Laboratoriumbevindingen:
  • Er werd vastgesteld dat de P-golf snelheid ($V_p$) afneemt naarmate de initiële schade toeneemt, wat de bruikbaarheid ervan als schade-indicator bevestigt.
  • CT-scans toonden aan dat een hogere initiële schade correleerde met een hogere porositeit en scheurdichtheid.
  • Brazil Split Tensile Test: De belasting-verplaatsingscurven vertoonden een strikt lagere convexe functie. De multiscale berekeningsresultaten stemden beter overeen met de experimentele curven dan conventionele simulaties, met name in de pre-piek regio.
  • Uniaxiale Compressietest: Spanning-rek curves waren monotoon toenemende concave functies. De multiscale benadering voorspelde accuraat de equivalente elasticiteitsmodulus, met waarden die dichter bij de testresultaten lagen dan conventionele simulaties. Schade-evolutie toonde aanvankelijk een geleidelijke toename, gevolgd door een snelle ontwikkeling en stabilisatie.
• Veldbevindingen:
  • In de veldcompressietest van gegoteerde gestolde lichamen voorspelde de multiscale benadering accuraat de spanning-rek curve tot aan de piek uniaxiale compressieve spanning.
  • De berekende equivalente elasticiteitsmodulus voor de gegoteerde gestolde lichamen kwam overeen met de veldtestgegevens, wat de toepasbaarheid van de methode voor complexe engineering-scenario's met gebroken erts en cementslurry-aggregaten valideert.
• Geobserveerde Beperkingen: De studie merkte op dat hoewel de pre-piek resultaten accuraat waren, de post-piek spanning-rek curves voor de gegoteerde gestolde lichamen afweken van de veldresultaten. Deze discrepantie wordt toegeschreven aan het gebrek aan een rek-verweking parameter in het model, die noodzakelijk is om het post-falen gedrag van materialen te vangen die verschillen van diepe brosse gesteenten.

Betekenis en Claims
Het artikel claimt dat de voorgestelde multiscale berekeningsbenadering een robuuste theoretische methode biedt voor het oplossen van complexe problemen in de rotsbouwtechniek. Door rekening te houden met initiële schade en gebruik te maken van perturbatietheorie, slaagt de methode erin de kloof tussen microstructuur-heterogeniteit en macroscopisch mechanisch gedrag te overbruggen. De auteurs stellen dat deze benadering een nauwkeuriger alternatief biedt voor conventionele simulatiemethoden voor het bepalen van de equivalente elasticiteitsmodulus van rotsmassa's. Zij suggereren dat als een geschikte brug tussen macro- en microschalen geïdentificeerd kan worden, deze methodologie uitgebreid kan worden naar het bepalen van mechanische parameters voor andere heterogene materialen. Het werk wordt gepresenteerd als een referentie voor het verbeteren van de beoordeling van de stabiliteit van rotsmassa's in engineering projecten, met name waar initiële schade en schaaleffecten kritieke factoren zijn.