# Overconstraint conflict between TIE and BOUNDARY: Matching Displacements but Divergent Eigenvalues

**URL:** <https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873>\
**Category:** Analysis issues\
**Created:** [April 10, 2026, 11:06am UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873 "2026-04-10T11:06:45Z")\
**Posts on this page:** 6\
**Page:** 1

<div class="post-metadata">

**Author:** ![aj\_warrier](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/aj_warrier/32/4213_2.png) [@aj\_warrier](https://calculix.discourse.group/u/aj_warrier)\
**Post date:** [April 10, 2026, 11:06am UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/1 "2026-04-10T11:06:45Z")

</div>

**Hi everyone,**

I am researching a composite beam modeling approach by comparing two cantilever setups:

1. A baseline beam modeled entirely with **C3D20** (20-node hexahedral) elements.

2. A composite beam where **S8** shell elements are “jacketed” onto the top and bottom surfaces of a C3D20 solid core.

**The Setup:**

- **Loading:** Surface traction on the free end (right face).

- **Boundary Condition:** The left face of the solid core is fixed using `*BOUNDARY`.

- **Interaction:** I used `*TIE` to bind the shell-solid interfaces, with the **solid face as the Master** and the **shell surface as the Slave**.

**The Issue:**

In the initial run, the shells “peeled” away from the solid core specifically at the fixed end. Even though the Master (solid) nodes were fixed, the Slave (shell) nodes at that specific boundary seemed to ignore the constraint.

**Attempted Fix:**

I applied an additional `*BOUNDARY` to the fixed-end edges of the shells. This stopped the peeling, but the solver issued warnings indicating that **tied constraints were not generated for those nodes** because they were already subject to an SPC. However-

**The “Paradox” in the Results:**

1. **Static Analysis:** Despite the deactivation of the tie at the fixed boundary, the **displacement results** for both the pure solid and composite beams match almost perfectly.

2. **Modal/Matrix Analysis:** When I extract the final global stiffness matrix (I wrote a custom code to export the stiffness vals at the end in mafillsm.f) and solved for the lowest positive **eigenvalues** , the results revealed that the composite beams were much softer than the pure solid one. The eigenvalues were consistently ~45% lower than those of the pure solid cantilever beam.

**The Question:** It seems the “overconstraint” is creating a localized region at the support where the shell and solid are no longer mathematically coupled in the stiffness matrix, even if the static displacements look “correct” due to the shared boundary condition.

1. How can I ensure the `*TIE` remains active across the entire interface, including the fixed boundary?

2. Is there a recommended way to “fix” the end of a tied shell-solid assembly without triggering the overconstraint/redundancy errors that deactivate the `*TIE` at the support?

I appreciate any insights on this. Have a nice day!

---

<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:** [April 10, 2026, 11:17am UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/2 "2026-04-10T11:17:58Z")

</div>

Did you try swapping the master/slave surfaces ? Or using tied contact instead of tie constraints ?

---

<div class="post-metadata">

**Author:** ![aj\_warrier](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/aj_warrier/32/4213_2.png) [@aj\_warrier](https://calculix.discourse.group/u/aj_warrier)\
**Post date:** [April 10, 2026, 2:35pm UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/3 "2026-04-10T14:35:28Z")

</div>

Hey, yeah I did.

-Swapping the master/slave surfaces solved the MPC-SPC conflict issue. But the stiffness matrix analysis result remained the same.

-I had also tried using tied contact instead of tie constraints, but the MPC/SPC conflict issue remained.

---

<div class="post-metadata">

**Author:** ![Disla](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/disla/32/3548_2.png) [@Disla](https://calculix.discourse.group/u/Disla)\
**Post date:** [April 10, 2026, 4:09pm UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/4 "2026-04-10T16:09:49Z")

</div>

I have tried modeling a wall. The beam nodes are merged with the wall. No need to fix. The fix is in the backside of the wall far from the contact areas. The shells are \*tie agains the wall and to the beam.

Beam and wall are master (Blue) and shell surfaces are slaves.

 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/16498af260c56e04037e903e90d78a3e7d4f7f99.png)

 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/5/5e861dd187bb8b7755fe04a6879a99f801966225.png)

---

<div class="post-metadata">

**Author:** ![aj\_warrier](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/aj_warrier/32/4213_2.png) [@aj\_warrier](https://calculix.discourse.group/u/aj_warrier)\
**Post date:** [April 10, 2026, 4:27pm UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/5 "2026-04-10T16:27:11Z")

</div>

That makes a lot of sense to move the SPCs away from the MPC zone. In my case, I’m trying to validate the ‘subtractive’ effect of shells, so I’m keeping the geometry as lean as possible.

**Follwo-up question:** Since displacements match between cases but the extracted K is ‘softer’ in Python, is there a specific way in CCX to export the **condensed/reduced** stiffness matrix where the slave DOFs are already eliminated? Or is `mafillsm.f` always going to provide the expanded system?

---

<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:** [April 14, 2026, 9:27am UTC](https://calculix.discourse.group/t/overconstraint-conflict-between-tie-and-boundary-matching-displacements-but-divergent-eigenvalues/3873/6 "2026-04-14T09:27:22Z")

</div>

> [@aj\_warrier](#):
>
> the Slave (shell) nodes at that specific boundary seemed to ignore the constraint

Seems related to this bug/limitation: [Tied constraint fails if a forces are prescribed to slave nodes](https://calculix.discourse.group/t/tied-constraint-fails-if-a-forces-are-prescribed-to-slave-nodes/1660)
