i have maybe a possible bug in the load

About the development of the FEM module/workbench.

Moderator: bernd

Jee-Bee
Veteran
Posts: 2566
Joined: Tue Jun 16, 2015 10:32 am
Location: Netherlands

i have maybe a possible bug in the load

Post by Jee-Bee »

I created a plate created of shells with in the center an (point)load
after creating the inp i checked the inp file and there was no load at all
see attachments
Attachments
Compound_Mesh.inp.txt
(51.83 KiB) Downloaded 47 times
quadplate.FCStd
(43.11 KiB) Downloaded 50 times
User avatar
bernd
Veteran
Posts: 12851
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland
Contact:

Re: i have maybe a possible bug in the load

Post by bernd »

your analysis type is frequency. For frequency analysis type neither force nor pressure nor selfweight constraint are written to input file. FEM and CalcluiX experts around. Should load constraints be written on analysis type frequency?

frequency:
https://github.com/berndhahnebach/FreeC ... cx.py#L129
Jee-Bee
Veteran
Posts: 2566
Joined: Tue Jun 16, 2015 10:32 am
Location: Netherlands

Re: i have maybe a possible bug in the load

Post by Jee-Bee »

beamf2 one of the standard calculix examples have a force applied.
Attachments
Archive.zip
(4.6 KiB) Downloaded 53 times
fandaL
Posts: 440
Joined: Thu Jul 24, 2014 8:29 am

Re: i have maybe a possible bug in the load

Post by fandaL »

Seems that CalculiX enables frequency analysis of loaded structure in different way: first compute step with loading and then compute frequencies in the next step where parameter PERTURBATION means taking loaded state from the previous step.
Details from the manual but rather look at example code
http://feacluster.com/CalculiX/ccx_2.11 ... de245.html
... If the PERTURBATION parameter is used in the *STEP card, the load active in the last *STATIC step, if any, will be taken as preload. Otherwise, no preload will be active.
and http://feacluster.com/CalculiX/ccx_2.11 ... de147.html
If the perturbation parameter is not activated on the *STEP card, the frequency analysis is performed on the unloaded structure, constrained by the homogeneous SPC's and MPC's. Any steps preceding the frequency step do not have any influence on the results.
If the perturbation parameter is activated, the stiffness matrix is augmented by contributions resulting from the displacements and stresses at the end of the last non-perturbative static step, if any, and the material parameters are based on the temperature at the end of that step. Thus, the effect of the centrifugal force on the frequencies in a turbine blade can be analyzed by first performing a static calculation with these loads, and selecting the perturbation parameter on the *STEP card in the subsequent frequency step. The loading at the end of a perturbation step is reset to zero.
In the beamf2.inp there are 2 steps

Code: Select all

...
*STEP
*STATIC
*CLOAD
LAST,3,-48.155
*END STEP
*STEP,PERTURBATION
*FREQUENCY
10
*END STEP
btw a formal detail
Is it needed to print into .inp file line *STATIC before *FREQUENCY, as FreeCAD currentli does? I thought these are 2 different analysis types with own keyword *STATIC or *FREQUENCY
Jee-Bee
Veteran
Posts: 2566
Joined: Tue Jun 16, 2015 10:32 am
Location: Netherlands

Re: i have maybe a possible bug in the load

Post by Jee-Bee »

fandaL wrote:btw a formal detail
Is it needed to print into .inp file line *STATIC before *FREQUENCY, as FreeCAD currentli does? I thought these are 2 different analysis types with own keyword *STATIC or *FREQUENCY
Is not needed. in the frequency examples they use just *FREQUENCY. I think only on case like described above it is needed.
however i have seen some example dat files with just the eigen frequencies. so it could be that static is needed just for displacement and stresses etc
User avatar
bernd
Veteran
Posts: 12851
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland
Contact:

Re: i have maybe a possible bug in the load

Post by bernd »

for reference:

http://feacluster.com/CalculiX/ccx_2.11 ... de289.html

http://feacluster.com/CalculiX/ccx_2.11 ... de245.html
http://feacluster.com/CalculiX/ccx_2.11 ... de287.html
http://feacluster.com/CalculiX/ccx_2.11 ... de215.html

It will take some changes in the input file writer to change this. But we will need them anyway if we would like to support multiple steps at one point in the future. I will come up with some changes ...
User avatar
bernd
Veteran
Posts: 12851
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland
Contact:

Re: i have maybe a possible bug in the load

Post by bernd »

I reimplemented the analysis type handling. It is all in one def now. This makes it easier once we would like support multiple steps. Any testing is highly appreciated.

branch:
https://github.com/berndhahnebach/FreeC ... femccxstep

commit:
https://github.com/berndhahnebach/FreeC ... dfc9b681c6

cheers bernd
Jee-Bee
Veteran
Posts: 2566
Joined: Tue Jun 16, 2015 10:32 am
Location: Netherlands

Re: i have maybe a possible bug in the load

Post by Jee-Bee »

Hi was by accident looking true CCX writer.
I saw by accident that (still?) no constraints are applied. see https://github.com/FreeCAD/FreeCAD/blob ... er.py#L137
I would asume this is (still) a bug that never made it to master...
thschrader
Veteran
Posts: 3157
Joined: Sat May 20, 2017 12:06 pm
Location: Germany

Re: i have maybe a possible bug in the load

Post by thschrader »

bernd wrote: Tue Jan 31, 2017 12:45 pm Should load constraints be written on analysis type frequency?
In my opinion yes. For example, the first eigenfrequency of a column with
external pressure loading is a little bit lower than without loading.
User avatar
bernd
Veteran
Posts: 12851
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland
Contact:

Re: i have maybe a possible bug in the load

Post by bernd »

Thomas, have you tried on latest FreeCAD. My post is more than a year old!
Post Reply