# Calculix\_2.20\_4win default solver spooles

**URL:** <https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257>\
**Category:** Uncategorized\
**Created:** [September 10, 2022, 5:38pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257 "2022-09-10T17:38:10Z")\
**Posts on this page:** 20\
**Page:** 2

<div class="post-metadata">

**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\
**Post date:** [September 12, 2022, 4:23pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/22 "2022-09-12T16:23:20Z")

</div>

> [@JuanP74](#):
>
> I think your boundary conditions have to be as detailed as your structure’s model.

right, someone need to consider warping support conditions were classical formulation ignoring this.

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 13, 2022, 3:49am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/23 "2022-09-13T03:49:49Z")

</div>

using \*RIGID BODY with the entire cross-section and translation at the supports  
in length direction is fixed is like fixed warping displacement.  
why is there no change of the result if i fix only one support in length direction?  
in LTBEAM the result is different, but ccx it has no influence!?  
also with \*DISTRIBUTING!?  
because both supports are fixed in cross / traverse direction.  
so with RIGID BODY and entire cross-section warping displacement is fixed,  
with \*DISTRIBUTING warping displacement is free.  
solution is with \*RIGIG BODY to create two support over the cross-section,  
one for the upper part one for the lower part,  
or fixing only the web, or better, create the web with thickness 2 elements and  
fix only the nodes in the middle and in the web.  
to be continued

---

<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:** [September 13, 2022, 9:13am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/24 "2022-09-13T09:13:24Z")

</div>

Hi dichtstoff,

I have read in the ReleaseInfo.txt file contained in the LTB folder this:

Version N°1.0.6 Release date: 07/2005

* * *

Fixed errors:

. correction in units of restraining springs for teta (torsion) and teta’ (warping) (“/rd” added)

That would suggest some springs are probably used to properly constrain the beam and help to make warping effect to emerge.

It should also be more noticeable in short beams.

I’m looking for the proper way to apply BC for warping to show up.

EDITED: Found this in the help.

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

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 13, 2022, 10:57am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/25 "2022-09-13T10:57:25Z")

</div>

Hello Disla,

i have fixed for \*RIGID BODY only the nodes from web,  
and the results are adorable:

\*\*\*\* centric pressure: 1.000,0  
\*\*\*\* lateral buckling: 4.385,0 ccx: 4.405,0 0,005 %  
\*\*\*\* lateral buckling: 15.877,0 ccx: 15.242,0 0.040 %  
\*\*\*\* torsional buckling: 10.427,0 ccx: 11.233,0 0.077 %  
\*\*\*\* bending moment: 1.000.000,0  
\*\*\*\* lateral torsional buckling: 1.502,0 ccx: 1.534,0 0,021 %

if fix the entire cross-section at support,  
what triggers the fixed warp displacement?

wbr

![image](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/4/4800c0cddaad8f306b6e79674eec130d00780499.png)

 ![image](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/8/85891d3b12dc336fc354f8b8af9eaa49f55fde97.png)

---

<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:** [September 13, 2022, 12:56pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/26 "2022-09-13T12:56:58Z")

</div>

Hi,

I have downloaded your files. Unfortunately, there was no msh file to identify where those BC were applied.  
I have also seen you use solids (_SOLID SECTION) and nonlinear geometric analysis (_\*\*\* proof of lateral buckling by nonlinear geometric calculation).

¿I’m I right.? ¿How many nodes are involved to get that accuracy?

None of the coupling cards works for me with shells.

Regards

---

<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:** [September 13, 2022, 1:13pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/27 "2022-09-13T13:13:17Z")

</div>

¿How do you compute the deviation?¿Is doesn’t seems percentage %?

\*\*\* lateral buckling: 4.385,0 ccx: 4.405,0 (4.385,0-4.405,0)/4.385,0\*100 = 0.45%

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 13, 2022, 2:17pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/28 "2022-09-13T14:17:29Z")

</div>

Hi,

i have upload all files.  
i have forgotten the 100 in calc.  
i have also upload the shell-beam example

> **[Ihr persönlicher Freigabelink | GMX Cloud](https://c.gmx.net/@324820432001700943/c_i7kTnkR2uFf15P9r2CMQ)**
>
> Mit diesem Link erhalten Sie Zugriff auf die mit Ihnen geteilten Daten.

wbr dichtstoff

---

<div class="post-metadata">

**Author:** ![JuanP74](https://avatars.discourse-cdn.com/v4/letter/j/b5a626/32.png) [@JuanP74](https://calculix.discourse.group/u/JuanP74)\
**Post date:** [September 13, 2022, 3:14pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/29 "2022-09-13T15:14:42Z")

</div>

> [@dichtstoff](#):
>
> using \*RIGID BODY with the entire cross-section and translation at the supports  
> in length direction is fixed is like fixed warping displacement.

In fact for thick profiles warping isn’t relevant most of the time.

> [@dichtstoff](#):
>
> why is there no change of the result if i fix only one support in length direction?

Is NLGEOM activated? buckling is a linear procedure so it doesn’t activates non linear effects by default.

---

<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:** [September 14, 2022, 7:42am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/30 "2022-09-14T07:42:01Z")

</div>

Hi,

I’m working on the IP 448x400x24x12.

The beam technique as support on the extreme wedges has shown to perform pretty well on shells without the need of Distributing.  
I would like to finish the analysis for comparison, but I don’t have clear what does it mean the BC called

\*\*\*\* bending moment 1.000.000 / \*\*\*\* lateral torsional buckling: 1.502,0  
\*\*\*\*BENDING MOMENT  
Nrotnode\_01,2, 1000000  
Nrotnode\_02,2,-1000000  
\*\*\*\*BENDING MOMENT

¿Are these opposite moments of +/- 1 KN/m on each extreme of the beam?  
Thanks

 ![screenshot.14](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/b/b638087be4b0068cfa7015250c14e1d67aea27bf.png)

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 14, 2022, 10:48am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/31 "2022-09-14T10:48:50Z")

</div>

yes, unit load with 1.0 kN x 1000 N 1.0 m x 1000 mm  
do you know how you can place here a text with scrollbar?  
wbr

---

<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:** [September 14, 2022, 11:08am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/32 "2022-09-14T11:08:48Z")

</div>

> [@dichtstoff](#):
>
> do you know how you can place here a text with scrollbar?

Click the “Preformatted text” option (or press Ctrl+E).

---

<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:** [September 14, 2022, 2:10pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/33 "2022-09-14T14:10:16Z")

</div>

This is my results.

Buckling factor in the strong axis differs the most in a 4% .

Apart from the FEA result I have serious doubts the Euler’s formula used in STALBAU III as a reference has applicability for this beam.

This section (l/k) ratio for the strong axis is too low for the Euler formula applicability.

I’m using a beam element as said in the LTBeam Validation on each side (only web) and it works well.

There is no need to use \*DISTRIBUTING or Rigid Body. Loads and BC are applied directly to that beams. That fact allows me to use ccx v2.19 and shells.

 ![imagen](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/5/59ca5c270e87f44aee28e81c28f643ddec5b0bf2.png)  
 ![screenshot.19](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/c/cb9d9e259db4f8847a8809bc7879b097ea0c65fc.png)

---

<div class="post-metadata">

**Author:** ![xyont](https://yyz2.discourse-cdn.com/free1/user_avatar/calculix.discourse.group/xyont/32/3553_2.png) [@xyont](https://calculix.discourse.group/u/xyont)\
**Post date:** [September 14, 2022, 2:37pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/34 "2022-09-14T14:37:41Z")

</div>

> [@Disla](#):
>
> Buckling factor in the strong axis differs the most in a 4% .

this may caused by offset in the model, another possibility are support modeling.

below figures & tables i have got in 2015 for case of distributed loads,

 ![2015-09-25-13_12_14-calculix-graphix](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/2/2a3085366636a197a331ad83fdaf874e0b54ca9e.png)

 ![2015-09-25-14_50_28-tabelmaxlengthcolumn-ods-openoffice-calc](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/14f89b8260d63717fbeb5ba9a3421c4195c61eda.png)

> [@Disla](#):
>
> There is no need to use \*DISTRIBUTING or Rigid Body.

stiff beam approach is similar to MPC or Rigid Body itself, this is not recommended since it may lead to numerical problems or conflict with another features (e.g contact, plasticity) but still can be use with cautions.

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 14, 2022, 4:39pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/35 "2022-09-14T16:39:27Z")

</div>

I have changed my system with solid elements,  
i fixed one end and applied 1000 N on the other end.  
fixed and load on nodes directly:

\*BOUNDARY  
Nsbound\_01,1,3,0  
\*CLOAD  
Nsbound\_02,3,-18.51852

total force (fx,fy,fz) for set NSBOUND\_01 and time 0.0000000E+00  
2.684040E-10 2.313234E-10 1.000000E+03  
total force (fx,fy,fz) for set NSBOUND\_02 and time 0.0000000E+00  
-5.163383E-12 -3.188363E-12 -1.000000E+03

lateral buckling weak axis: (1096-1107)/1096_100 = 1,00 %  
lateral buckling strong axis: (3969-3944)/3969_100 = 0.63 %

 ![image](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/2/20c544afe1cdc0bb2367d6731c36bcbf37fd0970.png)

 ![image](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/1515324d71eb5e1e430050e44d7c7610dd53676b.png)

---

<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:** [September 14, 2022, 5:43pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/36 "2022-09-14T17:43:30Z")

</div>

Hi dichtstoff,

I have scaled my beam to **measure twice** the actual length .¿ Same as you ?.

Lateral buckling weak axis: (1096-1096)/1096 x 100 = 0,00 % deviation

Lateral buckling strong axis: (3969-3932)/3969 x 100 = 0.93 % deviation.

The mixture of shells with small beams to help support works nice. I will use it with caution untill \*RIGID BODY and DISTRIBUTING are available for shells.

I would say that the deviations seen with the Stahlbau III pag 29 is because Euler’s formula is not appropriate for the beam used in the example. Slender ratio is too small.

Johnson’s parabolic would be a safety option if one wants a theoretical reference.

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 14, 2022, 6:24pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/37 "2022-09-14T18:24:11Z")

</div>

Hi Disla,  
thnx a lot. I’ll check these Johnson’s parabolic.  
Can you calc. by hand the values for buckling strong axis with johnson’s parabolic?  
or what is the criteria between euler and johnson?  
I take so much out with these discussion.  
I know how to handle now the support for lateral torsional buckling.  
Anyway, I still prefer to work with solid elements,  
I’ll try to update my example.  
thnx to all, I’m so happy. 🙂  
wbr dichtstoff

---

<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:** [September 14, 2022, 8:47pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/38 "2022-09-14T20:47:08Z")

</div>

You have some notes of Johnson’s parabolic in Wikipedia with a reference for the validity range of Euler’s formula.  
Strong axis (l/k) is about half below the critical value. If a column has (l/k)\<(l/k)cr you must pay attention to other failure mechanisms as it may reach the plastic range before buckling. Nonlinear buckling analisys with plasticity would be required.

| | | Strong Axis | | Weak Axis | |
| --- | --- | --- | --- | --- | --- |
| Pcr | | 15.877,27 | KN | 4.385,04 | KN |
| (l/k) | | 56,0 | | 106,5 | |
| (l/k)cr | | 108,1 | | 108,1 | |

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 14, 2022, 10:41pm UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/39 "2022-09-14T22:41:37Z")

</div>

Hello Disla,  
I’m not convinced with your thesis about Johnson’s parabolic,  
because it is empirically based equation.  
what have you in mind, when you say nonlinear buckling analysis?  
wbr

---

<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:** [September 15, 2022, 6:58am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/40 "2022-09-15T06:58:59Z")

</div>

Johnson’s parabolic considers that some beam geometries overpass the yield point way before the Ncr is reached.

If your first buckling mode would have been Ncr=15.877 KN , by that time, you would have 661 Mpa in the section.

Nonlinear buckling, I mean to introduce some geometric imperfections ( better with plasticity here ) and run an increasing load/displacement up to the point where the buckling shape starts to show.

If you want to see the aftermath of reaching 15.877KN you would need to adjust BC to remove the first two modes and make Ncr3 to become Ncr1.See pict.

 ![screenshot.26](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/3/30c3a9948d89ec10a461d73a413e3f1205fa95b4.png)

---

<div class="post-metadata">

**Author:** ![dichtstoff](https://avatars.discourse-cdn.com/v4/letter/d/dec6dc/32.png) [@dichtstoff](https://calculix.discourse.group/u/dichtstoff)\
**Post date:** [September 15, 2022, 10:56am UTC](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257/41 "2022-09-15T10:56:42Z")

</div>

Hello Disla,  
maybe we should open a new topics.

the reason for the discrepancy is not yet defined or spotted for me.  
you connect these with Johnson’s parabolic where you have a yield point,  
but we don’t use or have these in eigenvalue Buckling Analysis.

If you show me how to apply residual stress for these profile,  
i start immediately with nonlinear buckling 😉

![3-s2.0-B9780081029312000073-f07-24-9780081029312](https://global.discourse-cdn.com/free1/uploads/calculix/original/2X/1/11303b6cc2beef19ed35154dd81cc32f0c88b1b9.jpeg)

> **[ab\_stabilitaet.pdf](https://www.ifme.ovgu.de/ifme_media/TM/pdf/zusatzmaterial/ab_stabilitaet.pdf)**
>
> 1071.98 KB

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

[Previous page](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257.md?page=1)

[Next page](https://calculix.discourse.group/t/calculix-2-20-4win-default-solver-spooles/1257.md?page=3)
