![有限元分析](http://www.nataid.com/uploads/allimg/2407/1_1905294641.jpg)
爆破引起的围岩损伤问题,一直是爆破工程界所关心的问题之一。对于爆破岩体损伤的研究方法主要有现场监测、相似试验和数值有限元分析等。由于爆破开挖过程极其复杂,进行物理模拟试验比较困难,一方面是设备复杂昂贵,另一方面是动力相似律不易满足。随着计算机技术和数值仿真软件的快速发展及完善,数值模拟研究方法已成为爆炸力学问题研究中的一种主要手段。
目前发展较为成熟的数值分析方法是有限差分法和有限元法。与有限差分法相比,有限元法在网格划分和边界条件处理等方面具有优点,现在已成为爆炸动力学中的主要数值分析方法。在有限元分析过程中,作为数值计算的基础,网格划分对计算结果有直接的影响。Johansson通过比较静态和动态载荷研究了网格的依赖性,表明如果本构模型中考虑应变率效应,应变率效应将取决于数值网格。因此不合理的网格划分,不仅会影响到计算精度及计算效率,甚至会导致计算结果偏离实际。
利用ANSYS/LS-DYNA软件,分别采用四边形映射网格(工况一)、四边形自由网格(工况二)和四边形自由网格局部加密(工况三)三种方式划分模型并进行有限元分析。通过对结果进行比较,研究不同网格划分方法对模拟结果的影响。
采用文献中的实验模型,利用非线性动力有限元软件ANSYS/LS-DYNA建立二维平面应变模型,炸药和混凝土靶选用Lagrange算法,给出了轴对称坐标下该模型的几何示意图。TNT装药半径8mm、长90mm、埋深100mm,置于半径500mm、长1000mm的混凝土靶中,反向起爆。
自由网格划分是自动化程度最高的网格划分技术之一,它在面上可以自动生成三角形或四边形网格,在体上自动生成四面体网格。通常情况下,可利用智能尺寸控制技术(smartsize)来自动控制网格的大小和疏密分布,系统默认为第6级,共分为10级,级别越高网格划分越粗,级别越低网格划分越细。自由网格划分简便,不必考虑划分细节,只要将划分设定为一定级别即可,省掉了手工划分的麻烦。本文smartsize取值为1,划分得到的网格尺寸为cm级,见图。
由图可发现,智能网格划分时虽然采用了最高精度,但是在炸药附近区域网格划分仍然过于稀疏。为了弥补智能网格划分的不足,ANSYS/LS-DYNA提供了几种网格局部细化的方法。比较简单的是一种“相容”网格加密技术EREFINE,该方法是在原有网格基础上进一步细化单元,共分5级,级别越高单元尺寸细化越小。但是,这种细化方法必须考虑与原有节点连接等一系列问题,细化级别太高时,可能会导致网格划分的失败。细化深度设为2级,在原有网格的基础上,将炸药附近区域的网格尺寸细化到mm级,如图所示。
映射网格是对规则模型的一种网格划分方法,对于不规则模型,通常做法是利用布尔运算对其切割,然后对切割好的模型进行映射网格划分。因此,对于复杂模型,映射网格划分方式相较于自由网格,需要更多的时间和精力。但是该方法可以控制划分细节,即对应力集中区、爆破损伤区等感兴趣的区域进行网格尺寸的控制。图是对模型进行映射网格划分的结果,采用5mm网格尺寸均匀划分。
混凝土靶采用HJC本构,该模型能够反映混凝土等脆性材料在大变形、高应变率和高围压下材料损伤失效动态响应。HJC模型考虑了材料强度的压力-硬化效应、损伤弱化效应以及应变累积效应,通过引入损伤变量D反映材料的非弹性变形积累对岩石的影响。其等效屈服强度是压力、应变率和损伤的函数,损伤量则是塑性体应变、等效塑性应变和压力的函数。
专业从事机械产品设计│有限元分析│CAE分析│结构优化│技术服务与解决方案
杭州那泰科技有限公司
本文出自杭州那泰科技有限公司www.nataid.com,转载请注明出处和相关链接!