Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA phase-field method for modeling cracks with frictional contact
77
0
0.0
(
0
)
Added by
Jinhyun Choo
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We introduce a phase-field method for continuous modeling of cracks with frictional contacts. Compared with standard discrete methods for frictional contacts, the phase-field method has two attractive features: (1) it can represent arbitrary crack geometry without an explicit function or basis enrichment, and (2) it does not require an algorithm for imposing contact constraints. The first feature, which is common in phase-field models of fracture, is attained by regularizing a sharp interface geometry using a surface density functional. The second feature, which is a unique advantage for contact problems, is achieved by a new approach that calculates the stress tensor in the regularized interface region depending on the contact condition of the interface. Particularly, under a slip condition, this approach updates stress components in the slip direction using a standard contact constitutive law, while making other stress components compatible with stress in the bulk region to ensure non-penetrating deformation in other directions. We verify the proposed phase-field method using stationary interface problems simulated by discrete methods in the literature. Subsequently, by allowing the phase field to evolve according to brittle fracture theory, we demonstrate the proposed methods capability for modeling crack growth with frictional contact.
rate research
Read More
We present a barrier method for treating frictional contact on interfaces embedded in finite elements. The barrier treatment has several attractive features, including: (i) it does not introduce any additional degrees of freedom or iterative steps, (ii) it is free of inter-penetration, (iii) it avoids an ill-conditioned matrix system, and (iv) it allows one to control the solution accuracy directly. We derive the contact pressure from a smooth barrier energy function that is designed to satisfy the non-penetration constraint. Likewise, we make use of a smoothed friction law in which the stick-slip transition is described by a continuous function of the slip displacement. We discretize the formulation using the extended finite element method to embed interfaces inside elements, and devise an averaged surface integration scheme that effectively provides stable solutions without traction oscillations. Subsequently, we develop a way to tailor the parameters of the barrier method to embedded interfaces, such that the method can be used without parameter tuning. We verify and investigate the proposed method through numerical examples with various levels of complexity. The numerical results demonstrate that the proposed method is remarkably robust for challenging frictional contact problems, while requiring low cost comparable to that of the penalty method.
Simulation of contact mechanics in fractured media is of paramount important in the scope of computational mechanics. In this work, a preconditioned mixed-finite element scheme with Lagrange multipliers is proposed in the framework of constrained variational principle, which has the capability to handle frictional contact mechanics of the multi-crossing fractures. The slippage, opening and contact traction on fractures are calculated by the resulted saddle-point algebraic system. A novel treatment is devised to guarantee physical solutions at the intersected position of crossing fractures. A preconditioning technique is introduced to re-scale the resulting saddle-point algebraic system, to preserve the robustness of the system. An iteration strategy, namely monolithic-updated contact algorithm, is then designed to update the two primary unknowns (displacement and Lagrange multiplier) in one algebraic block. A series of numerical tests is conducted to study the contact mechanics of single- and multi-crossing fractures. Benchmark study is presented to verify the presented numerical method. Two tests with crossing fractures are studied, in which the slippage and opening can be calculated. The effects of crossing fractures on the deformation field can be observed in the calculated results, in which the variation of slippage/opening is analyzed by different loading conditions.
In this paper, we present an efficient numerical algorithm for solving the time-dependent Cahn--Hilliard--Navier--Stokes equations that model the flow of two phases with different densities. The pressure-correction step in the projection method consists of a Poisson problem with a modified right-hand side. Spatial discretization is based on discontinuous Galerkin methods with piecewise linear or piecewise quadratic polynomials. Flux and slope limiting techniques successfully eliminate the bulk shift, overshoot and undershoot in the order parameter, which is shown to be bound-preserving. Several numerical results demonstrate that the proposed numerical algorithm is effective and robust for modeling two-component immiscible flows in porous structures and digital rocks.
In this paper, we focus on modeling and simulation of two-phase flow with moving contact lines and variable density. A thermodynamically consistent phase-field model with General Navier Boundary Condition is developed based on the concept of quasi-incompressibility and the energy variational method. Then a mass conserving and energy stable C0 finite element scheme is developed to solve the PDE system. Various numerical simulation results show that the proposed schemes are mass conservative, energy stable and the 2nd order for P1 element and 3rd order for P2 element convergence rate in the sense of L2 norm.
This paper studies a model of two-phase flow with an immersed material viscous interface and a finite element method for numerical solution of the resulting system of PDEs. The interaction between the bulk and surface media is characterized by no-penetration and slip with friction interface conditions. The system is shown to be dissipative and a model stationary problem is proved to be well-posed. The finite element method applied in this paper belongs to a family of unfitted discretizations. The performance of the method when model and discretization parameters vary is assessed. Moreover, an iterative procedure based on the splitting of the system into bulk and surface problems is introduced and studied numerically.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
