# Multiple Shell Sections in one Model not working?

**URL:** <https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892>\
**Category:** General Questions\
**Tags:** material-definition\
**Created:** [July 22, 2024, 12:47pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892 "2024-07-22T12:47:54Z")\
**Posts on this page:** 18\
**Page:** 1

<div class="post-metadata">

**Author:** ![hnrwagner](https://avatars.discourse-cdn.com/v4/letter/h/e8c25b/32.png) [@hnrwagner](https://prepomax.discourse.group/u/hnrwagner)\
**Post date:** [July 22, 2024, 12:47pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/1 "2024-07-22T12:47:54Z")

</div>

Hello,

At the moment I am testing PrePoMax 2.1.0 for the buckling analysis of a wind turbine tower, the tower is a shell with multiple different sections (different wall thicknesses

 ![Tower](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/1/12601d521dc10f96508d32479ad64098fc1722fa.png)  
).

I start with a “Buckle” Analysis and it seems to me that PrePoMax can only handle 2 different sections, If I use only 2 Sections I get about the expected result (buckling factor = 3.4), If I try more than 2 different sections I get (buckling factor = 1).

“Update: During writing of the post I tried the same analysis with quadelement and second order and it works ! in the first try I used only linear tri elements. Still strange error though”

---

<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:** [July 22, 2024, 12:56pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/2 "2024-07-22T12:56:35Z")

</div>

> [@hnrwagner](#):
>
> in the first try I used only linear tri elements

I would avoid them in CalculiX - one such example: [Calculix vs Nastran big difference in results for shells - #22 by synt](https://prepomax.discourse.group/t/calculix-vs-nastran-big-difference-in-results-for-shells/1534/22)

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 4, 2025, 6:39pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/3 "2025-06-04T18:39:36Z")

</div>

Hello, I am encountering similar problems. When I try to solve the model using S4/S4R elements it works correctly for a cylinder of constant thickness, but when I include several thicknesses it does not give adequate results. Can this be attributed to the computational capacity of my PC (8Gb RAM), or what is the reason?

---

<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 4, 2025, 6:54pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/4 "2025-06-04T18:54:32Z")

</div>

Did you try refining the mesh or using different element types ? Is it a linear buckling analysis ?

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 4, 2025, 7:51pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/5 "2025-06-04T19:51:32Z")

</div>

I am modelling the cylinder with shell elements, so I where using S4/S4R. I tried different meshes finer and coarser, but similar results. Yes, I have specially noticed this in LBAs, but it is true that I tried in other computer (with less RAM) and it has similar problems for LA, that’s why I asking about computational capacity.

---

<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:** [June 4, 2025, 8:53pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/6 "2025-06-04T20:53:35Z")

</div>

I solved a similar problem (looks like the same) some time ago based on a well documented result of a FEM contest.

> **[8-MW wind turbine tower computational shell buckling benchmark. Part 1: An...](https://www.sciencedirect.com/science/article/pii/S135063072300078X?via%3Dihub)**
>
> An assessment of the elastic-plastic buckling limit state for multi-strake wind turbine support towers poses a particular challenge for the modern fin…

Mesh was dense + S8R. See picture.  
In my experience from diferent problems solved along the time, S4R is not reliable for Buckling analisys.  
Pay special attention to boundary conditions on the base. You have to constrain rotations on a curved surface. That’s not trivial to achieve.

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

Deviation from reference value was minimum.

 ![imagen](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/d/d0723cf42c348f4779bd4ac31c97d2d2b712db0d.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:** [June 4, 2025, 9:29pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/7 "2025-06-04T21:29:13Z")

</div>

What matrix solver are you using ? There might be some issues in those analyses with Spooles and PaStiX. Pardiso is the most reliable.

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 5, 2025, 6:18am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/9 "2025-06-05T06:18:27Z")

</div>

Thanks for the comment. I have been working a lot with shell buckling (it is true that it is in ABAQUS), and I have obtained very good results using S4R. In fact, the second part of the paper you have shared is the reference solution and they proposed the use of S4R. So, my question here is, are the problems with S4R limitation of Calculix?

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 5, 2025, 6:19am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/10 "2025-06-05T06:19:21Z")

</div>

Yes, I am using Pardiso. It was the fastest one too

---

<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 5, 2025, 6:49am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/11 "2025-06-05T06:49:47Z")

</div>

Check this benchmark study: [Buckling of composite panels - #34 by lucas\_bueno - CalculiX (official versions are on www.calculix.de, the official GitHub repository is at https://github.com/Dhondtguido/CalculiX).](https://calculix.discourse.group/t/buckling-of-composite-panels/1624/34)

That whole thread can be interesting too.

---

<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:** [June 5, 2025, 8:33am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/12 "2025-06-05T08:33:06Z")

</div>

In CalculiX, S4 and S4R four-node shell elements are expanded into three-dimensional C3D8I and C3D8R elements, respectively.

By other hand a node leads to a knot if

• the direction of the local normals in the elements participating in the node

diﬀer beyond a given amount. Notice that this also applies to neighbor-

ing elements having the inverse normal. Care should be taken that the

elements in plates and similar structures are oriented in a consistent way

to avoid the generation of knots and the induced nonlinearity.

_• the thickness is not the same in all participating elements._\<------

• the oﬀset is not the same in all participating elements.

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 5, 2025, 9:19am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/13 "2025-06-05T09:19:42Z")

</div>

Thank, it is a very interesting post and quite helpful

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 5, 2025, 9:29am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/14 "2025-06-05T09:29:01Z")

</div>

Ok, thank you. So, now I understand why this problem occurs when considering different thicknesses. So, how do you solve this by considering S8R since in this case it is not possible to avoid the generation of such knots because of the variable thickness?

---

<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 5, 2025, 9:38am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/15 "2025-06-05T09:38:08Z")

</div>

Sometimes it makes sense to try with corresponding solid elements in the case of such issues with expanded shells in CalculiX. PrePoMax has a tool for that - Thicken Shell Mesh.

---

<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:** [June 5, 2025, 1:15pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/16 "2025-06-05T13:15:40Z")

</div>

I have translated my old file with S8R from Mecway to most recent Prepomax version and it has perform excellent too.

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

S4 and S4R results are empty. It reports Buckling factor =1 but mode shapes are almost zero.

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

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 5, 2025, 1:49pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/17 "2025-06-05T13:49:36Z")

</div>

Yes, that’s exactly what I am obtaining for S4R but I am not able either to obtain the results for S8R. This is how they are looking

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/4/46b96fee0ed57a596bdfc3a6a26d8bd53b3fe5be.jpeg)

is it possible to know how you made the model??

---

<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:** [June 5, 2025, 2:36pm UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/18 "2025-06-05T14:36:44Z")

</div>

> [@jacabrera96](#):
>
> is it possible to know how you made the model??

You are very close.  
Try to increase the accuracy in the Buckling step window. That should improve your mode shape.  
Regarding your buckling factor (slightly low) , read the paper carefully because you are making a mistake that is described in it.

---

<div class="post-metadata">

**Author:** ![jacabrera96](https://avatars.discourse-cdn.com/v4/letter/j/67e7ee/32.png) [@jacabrera96](https://prepomax.discourse.group/u/jacabrera96)\
**Post date:** [June 6, 2025, 9:32am UTC](https://prepomax.discourse.group/t/multiple-shell-sections-in-one-model-not-working/1892/19 "2025-06-06T09:32:59Z")

</div>

Thank you so much. I think that more or less I got it.

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/0/0f842877dc7e79f8dd76d8faabc92352662a48f8.jpeg)

The only question I still have is why with the S8R elements it works but with the S4R didn’t? the knots are also appearing
