Partition of unity finite element for analysis of contact interface stress singularity
Singular mode shapes around the sharp frictional contact interfaces of different materials with linear elastic behavior are derived semi-analytically using the asymptotic eigen-function expansion method. The characteristic equation which is not proposed in a short explicit expression in the literatu...
Gespeichert in:
Veröffentlicht in: | European journal of mechanics, A, Solids A, Solids, 2020-05, Vol.81, p.103957, Article 103957 |
---|---|
1. Verfasser: | |
Format: | Artikel |
Sprache: | eng |
Schlagworte: | |
Online-Zugang: | Volltext |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Zusammenfassung: | Singular mode shapes around the sharp frictional contact interfaces of different materials with linear elastic behavior are derived semi-analytically using the asymptotic eigen-function expansion method. The characteristic equation which is not proposed in a short explicit expression in the literature is extracted in symbolic terms; then solved numerically and the singularity order is derived using Newton's iterative method for the first time in the literature.
The extracted eigen-shapes of the high gradient stress fields and the corresponding displacement fields satisfy the compatibility of deformations and tractions on the contact interface with arbitrary frictional interface equations. These asymptotic mode shapes are then implemented into the finite element interpolation using Partition of Unity Finite Element Method (PUFEM). Since the additional enrichment degrees of freedom on the singular contact nodes have the physical meaning of generalized stress intensity factor (GSIF), the direct GSIFs obtained are more accurate and the contour integral manipulations need not be maintained in the post processing phase of the contact analysis.
Finally, some engineering problems are solved which demonstrate the capacities of the approach concerning the accuracy of GSIF obtained from the direct-PUFEM proposed compared with that derived by reciprocal work contour integral method used in traditional finite element analyses.
•The order of singularity and shape functions for frictional sliding kinked contact interface are calculated numerically.•The number of additional DOFs is decreased substantially which are the same as Stress Intensity Factors (direct GSIF).•The sensitivity of the contour integral results out of the enrichment area is assessed.•Using direct-GSIF, the need for contour Integral formulation for finding the GSIFs in contact analysis is omitted.•It is shown that, the superposition of singular modes is acceptable for early steps of unidirectional contact. |
---|---|
ISSN: | 0997-7538 1873-7285 |
DOI: | 10.1016/j.euromechsol.2020.103957 |