A phase-field fracture model for brittle materials based on COMSOL
Received:April 13, 2021  Revised:August 03, 2021
View Full Text  View/Add Comment  Download reader
DOI:10.7511/jslx20210413002
KeyWord:brittle materials  phase-field method  COMSOL  numerical simulation  sensitivity analysis
                 
AuthorInstitution
张文兵 上海海事大学 海洋科学与工程学院, 上海 ;河海大学 水利水电学院, 南京
沈振中 河海大学 水利水电学院, 南京
徐力群 河海大学 水利水电学院, 南京
张宏伟 河海大学 水利水电学院, 南京
甘磊 河海大学 水利水电学院, 南京
李舸航 上海市堤防泵闸建设运行中心, 上海
Hits: 970
Download times: 321
Abstract:
      The phase-field method (PFM) describes the material fracture process through a series of differential equations,which avoid tedious crack surface tracing and have advantages in simulating crack initiation,propagation,and bifurcation.This paper introduced a fracture model for brittle materials based on the PFM,gave the derivation process of the governing equations of the PFM for the brittle material fracture problem,and proposed a method to implement the phase-field fracture model for brittle materials in COMSOL based on the staggered algorithm.The cracking process of the brittle material one element model and the single-edge notched plate subjected to tensile and shear are reproduced,and the simulated crack propagation path is consistent with the results of the existing literature,which verifies the rationality of the computer program.For many parameters involved in the phase-field fracture model of brittle materials,a global sensitivity analysis based on the Morris method is carried out to identify the main factors of the phase-field fracture model pertinent to the load-displacement relationship.The results show that the Young’s modulus (E),critical energy release rate (Gc) and displacement increment (Δux) are the main parameters that affect the output results of the load-displacement relationship.The phase-field fracture model implemented in COMSOL can effectively simulate the crack initiation and propagation in brittle materials,and the model parameters E,Gc and Δux have an important impact on the improvement of the fracture performance of the material and the efficiency of the inversion of the model parameters.