# Boundary conditions of a deformable strip contact

**URL:** <https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797>\
**Category:** General Questions\
**Created:** [June 8, 2024, 12:11pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797 "2024-06-08T12:11:02Z")\
**Posts on this page:** 14\
**Page:** 1

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 12:11pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/1 "2024-06-08T12:11:02Z")

</div>

![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/5/5065e6e666e46b4ab398a5a5a53927fcd9536f9f.jpeg)  
Hi guys, so as you can see, I’m trying to simulate a deformable strip contact. In this, I have kept the lower strip as rigid, and fixed, as well as I have applied a tie contact between the two strips. Additionally I am applying a concentrated force on the upper strip so that it deforms and basically sticks to the lower strip as long as the strip is in contact, however I am not been able simulate it without getting errors like

\*ERROR in cascade: the DOF corresponding to  
node 2924 in direction 1 is detected on the  
dependent side of a MPC and a SPC

in usermpc: mpc of type  
is unknown

Edit 1: I have also kept the boundary of the upper strip at the left corner ( the vertical one as fixed) and given a displacement - rotation constraint on the upper strip such that it can move in x y plane, but not rotate

Hence, I needed your help in understanding the problem behind this, what am I missing?  
[Strip\_Contact.zip](https://prepomax.discourse.group/uploads/short-url/6tG1kzASxj5YzqKE0PjClf4YvoS.zip) (6.5 MB)

---

<div class="post-metadata">

**Author:** ![a.baer](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/a.baer/32/3382_2.png) [@a.baer](https://prepomax.discourse.group/u/a.baer)\
**Post date:** [June 8, 2024, 2:02pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/2 "2024-06-08T14:02:00Z")

</div>

You have overconstraint in model (several nodes belong to the fix, the tie and to contact. Try the boundary condition and contact to split.

 ![grafik](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/d/d810457e73b9f6fbe9811a694e80976893bbfced.png)

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 2:36pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/3 "2024-06-08T14:36:54Z")

</div>

Does that mean no rigid body? since I dont want the lower strip to be deformed that much at all

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [June 8, 2024, 2:43pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/4 "2024-06-08T14:43:02Z")

</div>

1. You can keep the rigid body constraint but boundary conditions should be applied to its reference node, not to the nodes of the rigid part.
2. Don’t use concentrated force this way: [Don't apply concentrated force load to surfaces](https://prepomax.discourse.group/t/dont-apply-concentrated-force-load-to-surfaces/1271)
3. You don’t have contact in this model, only surface interaction (that could be used for contact pair definition) and tie constraint.
4. If you want to use both tie and contact, make a partition on that edge so one segment will have tie constraint and the other will have contact (to avoid overconstraint).
5. You can’t model sticky contact in CalculiX, there’s no cohesive behavior. So once the contact is established, it will be able to separate. Just saying.

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 2:44pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/5 "2024-06-08T14:44:12Z")

</div>

Wow thanks a lot this helps so much!

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 3:01pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/6 "2024-06-08T15:01:30Z")

</div>

So, I followed your advice, ditched contact, applied the correct constraints, but somehow a problem we talked about earlier rears its head:

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/e/ec6963363105076b05b395d198a4756bccf0d641.png)  
Any way, I can make it that the upper strip instead of going through the lower strip kind of ‘squishes’ on top of the lower strip?

I used two different kinds of surface interactions in two cases:

Case 1: Hard

Case 2: I set type to linear and  
K= 50 times of adjacent material (according to calculix handbook)  
Stick slope of friction = 1000N/mm^3

However I got the same result both times. I watched the spur gear and snap fit tutorial for reference.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [June 8, 2024, 3:03pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/7 "2024-06-08T15:03:05Z")

</div>

Can you share the modified file ?

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 3:08pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/8 "2024-06-08T15:08:30Z")

</div>

Here you go. This file contains the second case surface contact details  
[Strip\_Contact.pmx](https://prepomax.discourse.group/uploads/short-url/cwdqwZ1Jf42CI9PhiB8Eb52HX4s.pmx) (2.3 MB)

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [June 8, 2024, 3:14pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/9 "2024-06-08T15:14:13Z")

</div>

Now it looks better but there’s still no contact so the beams won’t see each other and there will be this penetration. As I said, surface interaction itself is not enough. You have to create a contact pair between those 2 surfaces and assign that surface interaction to it.

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 8, 2024, 3:17pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/10 "2024-06-08T15:17:07Z")

</div>

Done, will work on it, thanks a lot, I know I have been a bit troublesome these past few days but I appreciate your patience and guidance, means a lot

---

<div class="post-metadata">

**Author:** ![synt](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/synt/32/4658_2.png) [@synt](https://prepomax.discourse.group/u/synt)\
**Post date:** [June 8, 2024, 7:29pm UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/11 "2024-06-08T19:29:37Z")

</div>

only looking from the sketch of problem, maybe it can be tried with both penalty and tied contact type not a contraint.

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 9, 2024, 7:20am UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/12 "2024-06-09T07:20:51Z")

</div>

Thank you, it worked wonderfully. I had a small question: Let’s say I gave an amplitude to the entire system, such that the switch goes up and down. Whenever I try to do so by using Direct stepping, the system does move, but the simulation fails the moment it comes into contact ( bear in mind that the system does have contacts enabled). However, when I use automatic stepping, the system not only disregards the movements, but also the contacts. Am I missing something? I couldn’t find anything in the documentation that could help me with this.

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/3/30e841395a004c58593f1f488576dbf110b62280.jpeg)  
 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/e/ec6963363105076b05b395d198a4756bccf0d641.png)  
This is in automatic stepping, → Error: too many cutbacks  
 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/3/3e74dd59cc7d222161fa144ce9a83eb3ced14fd5.jpeg)  
This is in direct stepping → Error : Divergence  
I have changed nothing but the stepping method here

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [June 9, 2024, 7:24am UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/13 "2024-06-09T07:24:27Z")

</div>

Direct incrementation is not recommended in most cases, unless used with a very low increment size to account for potential convergence issues usually caused by contact. Automatic incrementation should be sufficient but to properly account for the amplitude, you may have to use low initial and maximum increment size (at least on the level of the amplitude’s resolution, possibly lower).

---

<div class="post-metadata">

**Author:** ![lastairbender](https://avatars.discourse-cdn.com/v4/letter/l/9d8465/32.png) [@lastairbender](https://prepomax.discourse.group/u/lastairbender)\
**Post date:** [June 9, 2024, 7:26am UTC](https://prepomax.discourse.group/t/boundary-conditions-of-a-deformable-strip-contact/1797/14 "2024-06-09T07:26:24Z")

</div>

I see, will work on this as well, thank you! Also anything that you can tell me about the contacts being disregarded by calculix? I went through the discourse history and tried some things that I found were relatable but to no avail, just a mesh thing then ? maybe because the sim is dynamic instead of just static? Idk because the contacts seemed to work fine in static idk why they won’t in dynamic
