# Suitability of tetrahedral & wedge element for contact analysis

**URL:** <https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642>\
**Category:** Uncategorized\
**Created:** [May 7, 2023, 11:04am UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642 "2023-05-07T11:04:14Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\
**Post date:** [May 7, 2023, 11:04am UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642/1 "2023-05-07T11:04:14Z")

</div>

Hi all,

I read some argument from HKS Abaqus document, it’s related to suitability using tetrahedral and wedge element for contact analysis. However, there’s no further test or verification cited.

> Second-order tetrahedra are not suitable for the analysis of contact problems: a constant pressure on an element face produces zero equivalent loads at the corner nodes. In contact problems this makes the contact condition at the corners indeterminate, with failure of the solution likely because of excessive gap chatter. The same argument holds true for contact on triangular faces of a wedge element.

Above limitation did not notify in CalculiX documentation, so is there any comment and suggestion?

Regards,

---

<div class="post-metadata">

**Author:** ![Calc\_em](https://avatars.discourse-cdn.com/v4/letter/c/43a26b/32.png) [@Calc\_em](https://calculix.discourse.group/u/Calc_em)\
**Post date:** [May 7, 2023, 11:26am UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642/2 "2023-05-07T11:26:17Z")

</div>

It’s more about the second order but it also depends on contact formulation. For node-to-surface contact, it’s advised to avoid second-order tetrahedrons since they converge poorly and cause contact pressure noise. Modified versions of those elements are recommended instead. Of course, first-order tetrahedrons and wedges are bad in general, they should be used only as filler elements.

---

<div class="post-metadata">

**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\
**Post date:** [May 7, 2023, 12:20pm UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642/3 "2023-05-07T12:20:17Z")

</div>

Thanks for commenting, as i was thinking and guess before it’s specific for node to surface contacts. Not applicable for surface to surface contact. Abaqus [documentation](https://classes.engineering.wustl.edu/2009/spring/mase5513/abaqus/docs/v6.6/books/stm/default.htm?startat=ch03s02ath64.html) in element sections may lead to a user doubt since it’s not clear and notify an exception.

---

<div class="post-metadata">

**Author:** ![Fabrice](https://avatars.discourse-cdn.com/v4/letter/f/5fc32e/32.png) [@Fabrice](https://calculix.discourse.group/u/Fabrice)\
**Post date:** [May 10, 2023, 4:14pm UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642/4 "2023-05-10T16:14:29Z")

</div>

Actually there are comments on element formulation page 218 in contact paragraph.

---

<div class="post-metadata">

**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\
**Post date:** [May 10, 2023, 7:11pm UTC](https://calculix.discourse.group/t/suitability-of-tetrahedral-wedge-element-for-contact-analysis/1642/5 "2023-05-10T19:11:03Z")

</div>

When subjected to uniform pressure, quadratic element have some amount of negative value (tension) at corner, but Abaqus notify as zero values. That’s lead me some doubt, also document placed at element section without except.

CalculiX seems done properly, it did not give any notify in element sections document since it’s contact type dependent. More clear at contact sections document by explanation, also improved the convergency by small amount of positive values (compression).
