Flac3D单元等效结点力的计算与应用研究.pdfVIP

第30卷第6期 江西理工大学学报 Vo1．30，N。．6 2009年 12月 JOURNALOFJIANGXIUNIVERSITYOFSCIENCEANDTECHNOLOGY Dec．2009 文章编号：1007—1229(2009)06—0001—03 Flac3D单元等效结点力的计算与应用研究 袁海平 1,2a， 王金安 ， 赵 奎 知， 黄 晖 (1．北京科技大学土木与环境工程学院，北京 100083；2．江西理工大学，a-应用科学学院．b．工程研究院，江西赣州341000) 摘 要：针对Flac3D无法直接输出单元等效结点力的不足，采用虚位移原理和虚功等效原则，在 建立单元形函数和应变刚度矩阵的基础上，推导了Flac3D单元等效结点力的计算公式，编写了 相应计算机程序，并通过单轴压缩模拟试验算例验证，得到了满意的结果．该方法能够为岩体开 挖能量释放率计算和深部岩体岩爆分析提供有效途径． 关键词：等效结点力；等参数单元；Flac3D；能量释放率 中图分类号 ：TD32 文献标识码 ：A StudyonCalculationandApplicationofEquivalent NodalForceinFlac3D Element YUAN Hai-ping一，WANG Jin-anI,ZHAO KIli西，HUANG HIlj曲 (1．CivilandEnvironmentEngineeringSchool，BeijingUniversityofScienceandTechnology，Beijing100083，China； 2．JinagxiUniversityofScienceandTechnology，a．FacultyofAppliedScience；b．EngineeringResearchInstitute，Gnazhou341000，China) Abstract：Inordertocovertheshortagethatelementequivalentnodalforcecan’tbeputinFlac3D，principlesof virtualdisplacementandvirtualworkwereadoptedtoderiveasetofcomputingformulaforequivalentnodalforcein Flac3D elementbasedonbuildingelementshapefunctionandstrainstiffnessmatrixofelement．Hereby，acorre— spondingcomputerprogram waswritten．Acheckingexampleaboutuniaxialcompression testshowedthattheresult wassatisfaetory．Themethodprovidedaneffectivewaytocalculateenergyreleaserateofexcavatingrockmassandto analyzerockburstofdeeprockmass． Keywords：equivalentnodalforce；isoparametricelement；Flac3D；energyreleaserate 平均能量释放率作为评价指标，用等效结点力与结 0 引 言 点位移计算回采过程的动态开挖释放能量，取得了 单元的等效结点力是指将作用在单元上的各 良好的研究效果． 种荷载，包括表面力荷载、孔隙水压力荷载、自重荷 等效结点力主要应用于有限单元法 中71-91， 载、初始应力荷载等，等效移置到单元结点上，再逐 Flac3D目前还无法直接输出单元的等效结点力 ，也 点加以合成而得到的结点力 l-【-2]．t单元等效结点力 未给出单元等效结点力的具体算法．在单元应力 的计算对开挖岩体的能量释放、深部岩体岩爆等问 已知的情况下，如何将单元应力转换为等效结点 题的