# Undesired cylindrical constraints behaviour - Hydrostatic Analysis

**URL:** <https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405>\
**Category:** Bug Reports\
**Created:** [January 22, 2025, 6:04pm UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405 "2025-01-22T18:04:21Z")\
**Posts on this page:** 9\
**Page:** 1

<div class="post-metadata">

**Author:** ![jeblais](https://avatars.discourse-cdn.com/v4/letter/j/ac8455/32.png) [@jeblais](https://prepomax.discourse.group/u/jeblais)\
**Post date:** [January 22, 2025, 6:04pm UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/1 "2025-01-22T18:04:21Z")

</div>

Hello!

First post here 🙂 I have to say that I’m fairly new to PrePoMax/CalculiX but I love it! Thanks to the developpers for this free piece of software!

I’m currently testing PrePoMax by comparing analysis’ results with the Nastran In-CAD module for Autodesk Inventor. More specifically, I’m analysing the shell of an open-top water tank when full. To check the robustness of PrePoMax (v2.2.9-dev), I compared the results for the following analysis:

- Full 3D solid analysis - Shell nodes of tank free
- Partial (1 deg revolution) 3D solid analysis - cylindrical constraints
- Full shell analysis - cylindrical constraints on shell
- Partial (1 deg revolution) shell analysis - cylindrical constraints on shell
- 2D axisymmetric of tank of half cross-section. (On a side note, the draftsman did not generate the shell on the positive side of the radial axis, causing nonpositive Jacobian errors, I simply reversed the normal of the plane to solve the error, ignoring the “negative Z” warning.)

For the full shell model analysis, I get non-expected displacements results. I restricted the shell nodes in the revolution direction (theta or T) based on a cylindrical coordinate system but free in the R direction, and it seem the model reacts as if it’s a rectangular coordinate system. Am I doing something wrong? Is it a bug?

Thanks!

 ![Capture d’écran 2025-01-21 144954](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/c/c89ee4467d53a85c7d8e43a21f18883089c17610.png)

---

<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:** [January 22, 2025, 6:58pm UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/2 "2025-01-22T18:58:28Z")

</div>

Did you select the local CSYS when defining the BC ? Can you share the .pmx file ?

---

<div class="post-metadata">

**Author:** ![jeblais](https://avatars.discourse-cdn.com/v4/letter/j/ac8455/32.png) [@jeblais](https://prepomax.discourse.group/u/jeblais)\
**Post date:** [January 23, 2025, 1:27pm UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/3 "2025-01-23T13:27:42Z")

</div>

Yes I selected the local CSYS for the BC. PMX file is too big to be inserted as attachement, here is a link to the pmx file :

> **[Filebin | hzbe80vr3k6odsvh](https://filebin.net/hzbe80vr3k6odsvh)**
>
> Convenient file sharing. Registration is not required. Large files are supported.

Edit: updated hosting site to be less restrictive in download limit.

Also, I should have been more specific, this is what is expected in an axisymmetric manner in the full model, in terms of radial displacement:

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

Thanks!

---

<div class="post-metadata">

**Author:** ![Paul18fr](https://avatars.discourse-cdn.com/v4/letter/p/35a633/32.png) [@Paul18fr](https://prepomax.discourse.group/u/Paul18fr)\
**Post date:** [January 25, 2025, 9:19am UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/5 "2025-01-25T09:19:19Z")

</div>

Hi,

I played with the model and i can confirm such behaviour when:

- cylindrical coordinates sytem are applied on _ **all nodes** _
- and the hydrostatic pressure is implemented  
(see modif3 model)

I’va also (obviously) noted that when the pressure is implemented in the gobal system, everything is going correctly.

As the hydrostatic pressure is implemented using CLOAD (pressure recalculated following the 3 components) by PrePoMax:

- it is naturally modified by \*TRANSFORM
- as well as by the BCs (component Theta is null here)

I can imagine the equilibrium is modified, so I’m wondering if it the origin of the issue?

Paul

Link [here](https://transfert.free.fr/ZWateUB)  
(available 30 days from01/25/2025)

---

<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, 9:58am UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/6 "2025-03-20T09:58:37Z")

</div>

I can confirm that the hydrostatic pressure is implemented using CLOAD, which is the root of all problems.

---

<div class="post-metadata">

**Author:** ![Paul18fr](https://avatars.discourse-cdn.com/v4/letter/p/35a633/32.png) [@Paul18fr](https://prepomax.discourse.group/u/Paul18fr)\
**Post date:** [April 4, 2025, 6:31am UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/7 "2025-04-04T06:31:18Z")

</div>

Hi,

I can confirm the latest release fix the previous issue using CLOAD (now it is DLOAD); i basically reran the model and displacements are consistent

Thanks Matej

Link [here](https://transfert.free.fr/xQiV6uZ)  
(validity of the link: up to 05/02/2025)

 ![image_2025-04-04_082637316](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/d/d7097caeb6451c1788c45a86594ca942a1b0c5af.png)

---

<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:** [April 4, 2025, 7:45am UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/8 "2025-04-04T07:45:02Z")

</div>

Great! Thanks for the feedback. May I close this topic now?

---

<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:** [April 4, 2025, 10:13am UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/9 "2025-04-04T10:13:36Z")

</div>

I have compute reaction force in one cone hydrostatically loaded with the new Dload. Compared to analitycal solution deviation is below 0.002%. Structured mesh and medium refinement level.

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

---

<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:** [April 7, 2025, 1:05pm UTC](https://prepomax.discourse.group/t/undesired-cylindrical-constraints-behaviour-hydrostatic-analysis/2405/10 "2025-04-07T13:05:48Z")

</div>


