应变能释放率是粘弹性裂纹扩展判断的重要依据.为了解决低质量网格条件下的高精度粘弹性应变能释放率的数值计算问题,提出了裂纹尖端虚拟网格方法.在原始网格位移场基础上通过插值手段获得虚拟网格位移场、应力及应变场分布;基于虚拟网格信息并结合虚拟裂纹闭合方法,开展应变能释放率的数值计算.虚拟网格方法在应对裂纹扩展过程中任意复杂网格形式的同时,不再需要对原始裂纹尖端进行精细网格划分.两种典型断裂模式下的算例仿真结果表明,全积分形式的虚拟网格方法可以实现低泊松比条件下应变能释放率的高精度数值计算,相对误差均在 5%以内.为了应对高泊松比下的断裂问题,设计了虚拟网格方法的缩减积分方案;缩减积分方案下的应变能释放率相对误差在 1%左右,计算精度较全积分方案大幅提高.
The strain energy release rate is an important basis for the judgment of viscoelastic crack propagation.In order to solve the problem of numerical calculation of high-precision viscoelastic strain energy release rate under the condition of low-precision mesh,a virtual grid method for crack tip was proposed.Based on the original mesh displacement field,the virtual grid displacement field,stress and strain field distribution were obtained by interpolation.Based on the virtual grid information and the virtual crack closure method,the strain energy release rate was calculated.The virtual grid method can deal with any complex mesh form in the process of crack propagation,and at the same time,fine meshing of the original crack tip was no longer required.The simulation results show that the fully integrated virtual grid method can realize the high-precision numerical calculation of strain energy release rate under low Pois-son's ratio,and the relative error is less than 5%.To deal with the fracture problem at high Poisson's ratio,the reduced integration scheme of the virtual grid method was designed.Examples show that the relative error of strain energy release rate in the reduced inte-gration scheme is about 1%,the calculation accuracy is greatly improved compared with the full integration scheme.