# Error in constraint reactions

**URL:** <https://prepomax.discourse.group/t/error-in-constraint-reactions/2552>\
**Category:** General Questions\
**Created:** [March 19, 2025, 1:48pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552 "2025-03-19T13:48:29Z")\
**Posts on this page:** 10\
**Page:** 1

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [March 19, 2025, 1:48pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/1 "2025-03-19T13:48:29Z")

</div>

I analyzed a symmetrical shaped tank, with symmetrical loads and symmetrical constraints.  
The value of the constraint reactions RF1 and RF2 should be equal to zero, instead they are different from zero and different from each other. I don’t understand where I could have gone wrong.  
Thanks for the help.

> **[cono](https://wetransfer.com/downloads/5f7531f0ca4315b92483f46d22de880020250319134529/5dc17edc2466dae934a9b683b30d8fba20250319134529/1af874?t_exp=1742651129&t_lsid=2dba6f64-93e2-4491-a10c-6659bd3e5129&t_network=email&t_rid=ZW1haWx8NjdjMDg2OGMwMmI0Mzc4ZjAwNzBjM2Vi&t_s=download_link&t_ts=1742391929&utm_campaign=TRN_TDL_01&utm_source=sendgrid&utm_medium=email&trk=TRN_TDL_01)**
>
> 1 file sent via WeTransfer, the simplest way to send your files around the world

---

<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:** [March 19, 2025, 3:24pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/2 "2025-03-19T15:24:16Z")

</div>

Why are you using tetrahedrons instead of shell elements for such a highly thin-walled part ?

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [March 19, 2025, 5:27pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/3 "2025-03-19T17:27:05Z")

</div>

Reading the forum I thought I understood that calculix transforms shell elements into solid elements. I redid the model with Shell elements and the problem is not solved, the RF1 and RF2 constraint reactions are different from each other and always different from Zero and the deformed shape is not symmetric. I can’t understand the reason.  
[Cono\_Shell\_B.pmx](https://prepomax.discourse.group/uploads/short-url/gYcjbgbj3xxGweMoOECzk8fqWms.pmx) (3.5 MB)

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [March 20, 2025, 6:48am UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/4 "2025-03-20T06:48:51Z")

</div>

For the shell model your pressure definition is not symmetric. Use Right click → Preview to visualize it:

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/b/bb08bc0394ec28255aae7305830fbaad773b5aba.png)

---

<div class="post-metadata">

**Author:** ![Agazzotti](https://avatars.discourse-cdn.com/v4/letter/a/3da27b/32.png) [@Agazzotti](https://prepomax.discourse.group/u/Agazzotti)\
**Post date:** [March 20, 2025, 1:20pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/5 "2025-03-20T13:20:21Z")

</div>

Thanks for the reply. I understood my mistake: in the hydrostatic load field “Pressure change direction” by selecting two nodes of the loaded surface, the values ​​of N1 and N2 are also reported, which in this specific case must be equal to zero; therefore it is necessary to pay close attention to this step.  
I did not understand to what extent the values ​​of N1, N2 and N3 affect the load defined in the field: “Magnitude”

---

<div class="post-metadata">

**Author:** ![Gunnar](https://avatars.discourse-cdn.com/v4/letter/g/ed655f/32.png) [@Gunnar](https://prepomax.discourse.group/u/Gunnar)\
**Post date:** [March 20, 2025, 5:46pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/6 "2025-03-20T17:46:27Z")

</div>

I recently had a similar problem with a questionable reaction force when applying a pressure load (albeit using “uniform pressure”). I still don’t understand why I get a reaction force in the vertical direction with the press frame shown:

 ![frame](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/b/bdcbf840f6ae4319a022d8cf037e00b0e93b3903.png)

The uniform pressure multiplied by the cylinder base area results in exactly the same force as I applied to the top of the tabletop (surface traction). Nevertheless, I get a relatively large reaction force at the support edges.

On the other hand If I replace the pressure on the cylinder base with a equivalent force (the radial pressure on the cylinder wall remains unchanged), I get practically no reaction force in z-direction, as expected. The results of both calculations (stresses and displacements) are exactly the same… 🤔

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [March 20, 2025, 9:23pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/7 "2025-03-20T21:23:48Z")

</div>

Hi,

I assume you are using NLGEOM=ON. One reason could be that Traction (CLOAD) is not update while Pressure (DSLOAD) =F/S does. If the surface changes (due to Poisson effect for example), the overal reaction will end up unbalanced. This is very noticeble on hiperelastic materials for example.

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [March 24, 2025, 1:45pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/8 "2025-03-24T13:45:06Z")

</div>

Hi Gunnar,

May I ask you something ?, did those pretension arrows shown on your picture belong to unconnected opposite faces. ¿Did that work?. I have recently seen another case in which they are not able to pretension correctly oposite faces if there is not a continuous material inbetween. At least for me.

 ![imagen](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/1/19422895973f001e1ff90964c80db06ed04c4b59.png)

 ![imagen](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/7/7206c626fa89b887b2d750b8c0ae4f591cbf749e.jpeg)

---

<div class="post-metadata">

**Author:** ![Gunnar](https://avatars.discourse-cdn.com/v4/letter/g/ed655f/32.png) [@Gunnar](https://prepomax.discourse.group/u/Gunnar)\
**Post date:** [March 24, 2025, 7:21pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/9 "2025-03-24T19:21:32Z")

</div>

Hi,  
no, the preload is applied to the bolts, they were just not visible 🙂

 ![l](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/4/45a136ed8925d79312b2b60bb4eb1fa6847e6e98.png)

---

<div class="post-metadata">

**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\
**Post date:** [March 24, 2025, 7:30pm UTC](https://prepomax.discourse.group/t/error-in-constraint-reactions/2552/10 "2025-03-24T19:30:55Z")

</div>

aaaaaaa, Thank you 🙂
