爆炸荷载下岩石破坏的数值流形方法模拟

为了更好地利用数值流形方法对动力学问题进行分析,在对原数值流形方法中的动力学问题求解思想进行分析的基础上,采用动力有限元方法中的Newmark法对该算法进行了改进。改进后的数值流形方法与原来相比具有三个明显的优势：（1）当选择合适的参数后,该方法能够保证解的无条件收敛;（2）可以采用比原算法大得多的时间步长;（3）充分考虑了动力学问题中的阻尼效应。最后通过一个算例说明了改进后的数值流形方法能够很好地模拟岩石在冲击载荷作用下破坏的全过程,克服了有限元法不能模拟岩石破坏后块体运动情况的不足。     

数值流形方法及其研究进展

数值流形方法是当前岩体力学与工程界十分关注的一种新的数值方法.本文简述了数值流形方法的基本概念及其研究与应用进展,指出了它与有限元法的主要区别,提出了作者的几点粗浅看法.      

冲击载荷作用下裂纹动态响应的数值模拟

对垂直、剪切以及斜向等各种冲击载荷作用下裂纹的动态响应进行了数值模拟,得到了一系列随时间变化的动态应力场以及应变场图;根据其定义,计算出了相应的动态应力强度因子,进而分析了斜向载荷作用下裂纹起裂情况,并对最优断裂问题进行了阐述。     

冲击载荷下导线覆冰脱落过程的数值模拟

采用ABAQUS有限元软件建立输电线路有限元模型,模拟研究单档覆冰导线在冲击载荷作用下的除冰过程。覆冰破坏采用拉伸破坏准则,利用ABAQUS的用户材料子程序VUMAT,定义覆冰的本构关系并实现覆冰破坏单元删除的判定。进一步模拟研究在不同冲击载荷幅值、线路档距和覆冰厚度等条件下的除冰率,结果表明:除冰率随冲击载荷幅值的增加而增大;覆冰厚度和档距越小,除冰率越大;连续档线路需要更大的冲击载荷才能达到较好的除冰效果。所得结果可以为机械式除冰装置的设计提供参考。     

冲击载荷下软钢梁早期响应的数值模拟和简化模型

冲击载荷作用下，梁的早期响应既有弹性变形也有塑性变形，两者相互耦合．有限元数值模拟的结果表明，弹性弯曲波的传播是梁早期变形的主要机制，刚塑性简化理论预言的初始阶段中梁的“移行塑性铰”实际上是不存在的．本文提出的弹性 理想塑性简化模型可以很好地模拟固支软钢梁的早期响应     

单孔台阶爆破中保留岩体损伤特征的数值模拟

为了描述露天矿单孔台阶爆破中预留岩体在爆炸冲击载荷下的损伤演化特征,建立了一个各向异性的损伤模型。模拟结果表明,该损伤模型能够模拟岩石材料在冲击载荷下的损伤特征,可以为岩土工程中经历过冲击扰动作用的岩体的稳定性分析提供理论依据。     

冲击载荷下孔隙塌缩过程的数值模拟

爆炸压实过程中多孔体的孔隙闭合程度对压实效果起着决定性作用。利用LS-DYNA有限元程序,对无氧铜中的圆形孔隙塌缩过程进行了数值模拟。根据模拟结果分析得出,在6 GPa的冲击压力下,孔隙闭合时不同边界区域会发生爆炸焊接和射流侵彻,这2种结合机理可以使材料结合更牢固,提高材料的致密度和机械强度,实现高质量的爆炸压实。     

爆炸冲击震动对砖墙破坏作用的数值模拟

砖墙在爆炸冲击震动作用下的动力反应非常复杂,本构关系很难精确建立。本文阐述了砖墙几种常用的有限元模拟方法,分析各种方法的优缺点,确定采用一种砖块和砂浆分开的精细化建模的三维砖墙有限元模型;通过LS-DYNA软件,得到砖墙在水平爆炸冲击震动荷载下的破坏过程,计算结果与实验现象很好。研究表明：该种分析模型综合考虑了砖块和砂浆之间复杂的相互作用,并且对砂浆层进行了单独建模,保证了砖墙在数值模拟上的真实性和正确性,因此可以准确地模拟出实验中砖块之间砂浆层的损伤积累破坏。

     

岩石双孔爆破过程的流形元法模拟

在流形元基本理论基础上,引入断裂力学的裂纹产生及扩展判据,应用二阶流形元对岩石双孔爆破过程进行了研究,模拟了双孔同时起爆和毫秒延时起爆两种条件下裂纹的产生和扩展、块体的形成以及漏斗形成过程,分析了爆炸载荷作用下岩石的破坏规律以及起爆条件对该过程的影响,验证了流形元法在研究冲击载荷作用下岩石从连续体到不连续体破坏过程的准确性和有效性,为研究类似问题提供一个新的思路和方法。     

采动诱发岩体移动破坏过程数值模拟研究

本文基于岩石介质的宏观非线性主要是由非均质性和各向异性造成的, 应用新的数值计算软件RFPA(2D), 对采动引起岩体失稳破坏的全过程进了数值模拟研究。     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

THEORETICAL STUDY OF THREE-DIMENSIONAL NUMERICAL MANIFOLD METHOD

IntroductionNumerical manifold method is a new numerical method established on the basis of finitecover of manifold[1,2].By using continuous and non-continuous finite cover system,thenumerical manifold method includes the continuous and non-continuous as …     

数值流形单元法数学网格自适应

基于数值流形方法和有限覆盖技术,将有限元法的后验误差估计理论及h型网格自适应技术推广应用到数值流形单元法中,提出了数值流形单元法的后验误差估计方法和数学网格自适应技术,并编制了相应的程序。数值算例表明,经过网格自适应,可以在粗糙的初始网格基础上得到质量比较理想的网格,计算结果可达到用户要求的精度。     

非均匀泡沫金属材料在冲击载荷下的变形模拟

提出一种二维非线性弹塑性质量-弹簧-连杆模型,该模型将泡沫金属材料离散成许多质量块,质量块在受载方向由非线性弹塑性弹簧连接,垂直于受载方向由可延伸的弹性连杆铰接。采用该模型模拟并分析了层非均匀泡沫金属材料及局部不均匀泡沫金属材料在冲击载荷下的变形特性,说明了非均匀性对泡沫金属材料冲击变形的影响。     

数值流形方法的对象设计

研究了数值流形方法的对象设计方法和组织方式,为流形方法的理论研究和向三维问题的扩展打下良好的基础,研究发现数值流形方法具有编程和前后处理简单的特点,且仅用三角形流形单元和一阶近似的覆盖位移函数就可达到有限元多结点等参元的求解精度,具有深远的工程意义,计算结果表明,数值解与理论解吻合。     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

岩石巴西圆盘动态劈裂的流形元法模拟

利用二阶流形元法,通过引入裂纹产生与扩展判据,对冲击载荷作用下岩石平台巴西圆盘的动态拉伸劈裂过程进行了数值模拟,再现拉伸波作用下圆盘被劈裂的过程。模拟现象与实验结果相符,动态平衡时的应力分布与有限元结果基本一致。从而验证了流形元在模拟冲击载荷作用下材料动态破坏过程的有效性和可行性。     

冲击载荷作用下弹塑性板的反常动力响应研究

对弹塑性方板在横向脉冲载荷作用下动力响应的反直观行为进行了数值模拟及分析,给出了产生反直观行为的载荷范围和对相关问题的探讨。详细的分析发现,随着脉冲强度的增加,在几个窄的载荷区域,板的响应是反直观的,而且在此附近,结构参数、载荷等因素的微小改变将导致响应模式的很大差异,表明反直观行为对这些参数的敏感性。进一步的计算表明,这一特殊的动力行为主要与板的内力间的相互耦合作用密切相关,同时,卸载后的结构反弹到另一侧时发生较大的反向塑性变形,导致能量的进一步耗散,使板呈现反常的动力响应,这一现象是几何与材料两种非线性相互作用的结果。