search

calculix / CalculiX-Examples

CalculiX examples by Prof. Martin Kraska from Brandenburg University of Applied Sciences. Excellent starting point to master parametric modelling with CGX and CCX.
MIT License
299 stars 92 forks source link

CCX: kinematic coupling shows wrong stress distribution #18

Open mkraska opened 7 years ago

mkraska commented 7 years ago

Example Linear/StressConc1

Force plot and stress distribution indicate that forces are applied to nodes outside the specified surface

mkraska commented 4 years ago

Problem still present in 2.15. see example page

tournesol59 commented 3 years ago

Hi I have tried to get a mesh for the Test/Joints hinge example. I did it with cgx 2.15 but I was not able to work it with ccx the job was pending after writing on the console the number of nodes. I did many other of your problems with ccx but this one not. Would you have an idea to get mire messages?

mkraska commented 3 years ago

Can you please attach the problem files and the console output of ccx? Also some info on your OS would be helpful. It is not clear to me what actually failed: CGX or CCX?

tournesol59 commented 3 years ago

thank you for the quick reply. I was on travel. Here the replies. I did not understand why there is not an include input of dpin.inc and kpin.inc in solve.inp since it seems (as I understood) the idea of coupling to two ref center node the circular surfaces of the hinge) F Peugny

Here is the output of $ccx solve2 (my OS is Windows I used a precompiled version 2.14 of CalculiX found in bConverged...)

ERROR reading SOLID SECTION: SOLID SECTION can not be used for beam or shell elements Faulty element: 1 ERROR in gen3delem: first thickness in node 1 of element 1 is zero

Then $ccx solve (I invoke it first but I could not understand why an include pin.inc was called since it is writen with pre2.fbl --> I commented finally the line *include, input=pin.inc when using $ccx solve)

You are using an executable made on pon., 30 kwi 2018 09:46:04 The numbers below are estimated upper bounds number of:

 nodes:         70657
 elements:          1071
 one-dimensional elements:         3
 two-dimensional elements:         1068
 integration points per elements:         50
 degree of freedom per nodes:         3
 layers per element:         2

 distributed facial loads:         20
 distributed volumetric loads:         0
 concentrated loads:         0
 single point constraint:         78003
 multiple point constraints:         128521
 terms in all multiple point constraints:        822529
 tie constraints:         0
 dependent nodes tied by cyclic constraints:         0
 dépendent nodes in pre-tension constraints:        0

 sets:         11
 terms in all sets:         6508
 materials:         1
 constants per material and temperature:          2
 temperature points per material:         1
 plastic data points per material:         0

 orientation:         2136
 amplitudes:         3
 data points in all amplitudes:         3
 print requests:         0
 transformations:          0
 property cards:         0

ERROR reading SOLID SECTION: SOLID SECTION can not be used for beam or shell elements Faulty element: 1 ERROR in gen3delem: first thickness in node 1 of element 1 is zero

On 9/19/20, mkraska notifications@github.com wrote:

Can you please attach the problem files and the console output of ccx? Also some info on your OS would be helpful. It is not clear to me what actually failed: CGX or CCX?

-- You are receiving this because you commented. Reply to this email directly or view it on GitHub: https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-695359281

-- Frédéric Peugny 06 50 32 92 74 Résidence les Berges du Midi, Bat. A, N°19 25 rue du Bosc 31150 FENOUILLET France

mkraska commented 3 years ago

When running the simulation, the control script renames either dpin.inc or kpin.inc to pin.inc. There is an include command for pin.inc in solve.inp. did you execute the CGX scripts according to the instruction? Does test.py run?

----- Am 25. Sep 2020 um 22:11 schrieb tournesol59 notifications@github.com:

thank you for the quick reply. I was on travel. Here the replies. I did not understand why there is not an include input of dpin.inc and kpin.inc in solve.inp since it seems (as I understood) the idea of coupling to two ref center node the circular surfaces of the hinge) F Peugny

Here is the output of $ccx solve2 (my OS is Windows I used a precompiled version 2.14 of CalculiX found in bConverged...)

ERROR reading SOLID SECTION: SOLID SECTION can not be used for beam or shell elements Faulty element: 1 ERROR in gen3delem: first thickness in node 1 of element 1 is zero

Then $ccx solve (I invoke it first but I could not understand why an include pin.inc was called since it is writen with pre2.fbl --> I commented finally the line *include, input=pin.inc when using $ccx solve)

You are using an executable made on pon., 30 kwi 2018 09:46:04 The numbers below are estimated upper bounds number of:

nodes: 70657 elements: 1071 one-dimensional elements: 3 two-dimensional elements: 1068 integration points per elements: 50 degree of freedom per nodes: 3 layers per element: 2

distributed facial loads: 20 distributed volumetric loads: 0 concentrated loads: 0 single point constraint: 78003 multiple point constraints: 128521 terms in all multiple point constraints: 822529 tie constraints: 0 dependent nodes tied by cyclic constraints: 0 dépendent nodes in pre-tension constraints: 0

sets: 11 terms in all sets: 6508 materials: 1 constants per material and temperature: 2 temperature points per material: 1 plastic data points per material: 0

orientation: 2136 amplitudes: 3 data points in all amplitudes: 3 print requests: 0 transformations: 0 property cards: 0

ERROR reading SOLID SECTION: SOLID SECTION can not be used for beam or shell elements Faulty element: 1 ERROR in gen3delem: first thickness in node 1 of element 1 is zero

On 9/19/20, mkraska notifications@github.com wrote:

Can you please attach the problem files and the console output of ccx? Also some info on your OS would be helpful. It is not clear to me what actually failed: CGX or CCX?

-- You are receiving this because you commented. Reply to this email directly or view it on GitHub: https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-695359281

-- Frédéric Peugny 06 50 32 92 74 Résidence les Berges du Midi, Bat. A, N°19 25 rue du Bosc 31150 FENOUILLET France

— You are receiving this because you authored the thread. Reply to this email directly, [ https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-699131878 | view it on GitHub ] , or [ https://github.com/notifications/unsubscribe-auth/ADDC4DV2JBJKTJFKSPJXUXTSHT2P5ANCNFSM4DKEDOSQ | unsubscribe ] .

tournesol59 commented 3 years ago

No unfortunately I did not notice test.py under Windows I have problems with sys commands but I have another version on fedora that I will try. Thank you a lot FPeugny Le ven. 25 sept. 2020 à 23:15, mkraska notifications@github.com a écrit :

When running the simulation, the control script renames either dpin.inc or kpin.inc to pin.inc. There is an include command for pin.inc in solve.inp.

did you execute the CGX scripts according to the instruction? Does test.py run?

----- Am 25. Sep 2020 um 22:11 schrieb tournesol59 < notifications@github.com>:

thank you for the quick reply. I was on travel. Here the replies. I

did not understand why there is not an include input of dpin.inc and

kpin.inc in solve.inp since it seems (as I understood) the idea of

coupling to two ref center node the circular surfaces of the hinge)

F Peugny

Here is the output of

$ccx solve2

(my OS is Windows I used a precompiled version 2.14 of CalculiX found

in bConverged...)

ERROR reading SOLID SECTION: *SOLID SECTION can

not be used for beam or shell elements

Faulty element: 1

*ERROR in gen3delem: first thickness

in node 1 of element 1

is zero

Then

$ccx solve

(I invoke it first but I could not understand why an include pin.inc

was called since it is writen with pre2.fbl --> I commented finally

the line *include, input=pin.inc when using

$ccx solve)

You are using an executable made on pon., 30 kwi 2018 09:46:04

The numbers below are estimated upper bounds

number of:

nodes: 70657

elements: 1071

one-dimensional elements: 3

two-dimensional elements: 1068

integration points per elements: 50

degree of freedom per nodes: 3

layers per element: 2

distributed facial loads: 20

distributed volumetric loads: 0

concentrated loads: 0

single point constraint: 78003

multiple point constraints: 128521

terms in all multiple point constraints: 822529

tie constraints: 0

dependent nodes tied by cyclic constraints: 0

dépendent nodes in pre-tension constraints: 0

sets: 11

terms in all sets: 6508

materials: 1

constants per material and temperature: 2

temperature points per material: 1

plastic data points per material: 0

orientation: 2136

amplitudes: 3

data points in all amplitudes: 3

print requests: 0

transformations: 0

property cards: 0

ERROR reading SOLID SECTION: *SOLID SECTION can

not be used for beam or shell elements

Faulty element: 1

*ERROR in gen3delem: first thickness

in node 1 of element 1

is zero

On 9/19/20, mkraska notifications@github.com wrote:

Can you please attach the problem files and the console output of ccx? Also

some info on your OS would be helpful. It is not clear to me what actually

failed: CGX or CCX?

--

You are receiving this because you commented.

Reply to this email directly or view it on GitHub:

https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-695359281

--

Frédéric Peugny

06 50 32 92 74

Résidence les Berges du Midi, Bat. A, N°19

25 rue du Bosc https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g> 31150 FENOUILLET https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g> France https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

—

You are receiving this because you authored the thread.

Reply to this email directly, [

https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-699131878 |

view it on GitHub ] , or [

https://github.com/notifications/unsubscribe-auth/ADDC4DV2JBJKTJFKSPJXUXTSHT2P5ANCNFSM4DKEDOSQ

| unsubscribe ] .

—

You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-699158895, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADDYQV4BIX6GQ4IGLEXKR6TSHUB7PANCNFSM4DKEDOSQ .

-- Frédéric Peugny 06 50 32 92 74 Résidence les Berges du Midi, Bat. A, N°19 25 rue du Bosc 31150 FENOUILLET France

tournesol59 commented 2 years ago

Hello, Thank you for your reply. I did not manage to solve the problem a year ago. Nethertheless I have concentrated myself on modeling a contact and I tried another model on your site mkraska: Contact/Shell0. I think I encountered a similar problem as in "Joints" because my PYTHON27 is broken and I had to manually run the run.fbd file several time. I changed only the pc-ns.inp file to swep the frequency calculation, taking a force and a boundary condition into account in the static analysis. All in all the messages I get from the run.fbd reads that no convergence is reached after some iterations. (Example below) I know that these are few information from my side perhaps the best way to inform you is to take my attached modification and simply run tar xvf Shell0.tar cd Shell0 cgx -c run.fbd

(perhaps that will evoke you an idea of troubleshooting.. I have also attached an image of what I got for this check) Thanks Frederic (Stuttgart fellow) ^

and below are the screen message from ccx:

iteration 10

Determining the structure of the matrix: number of equations 930 number of nonzero lower triangular matrix elements 45696

Using up to 1 cpu(s) for the stress calculation.

Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

Factoring the system of equations using the unsymmetric spooles solver Using 1 cpu for spooles.

Using up to 1 cpu(s) for the stress calculation.

average force= 263722819289.723724 time avg. forc= 263722819289.723724 largest residual force= 3015411567852.941895 in node 107 and dof 2 largest increment of disp= 1.566085e+05 largest correction to disp= 1.006406e+03 in node 52 and dof 2

estimated number of iterations till convergence = 17 too slow convergence; the increment size is decreased to 7.812500e-03 the increment is reattempted

*ERROR: too many cutbacks best solution and residuals are in the frd file

On Sat, Sep 26, 2020 at 1:05 AM Frederic Peugny @.***> wrote:

No unfortunately I did not notice test.py under Windows I have problems with sys commands but I have another version on fedora that I will try. Thank you a lot FPeugny Le ven. 25 sept. 2020 à 23:15, mkraska @.***> a écrit :

When running the simulation, the control script renames either dpin.inc or kpin.inc to pin.inc. There is an include command for pin.inc in solve.inp.

did you execute the CGX scripts according to the instruction? Does test.py run?

----- Am 25. Sep 2020 um 22:11 schrieb tournesol59 < @.***>:

thank you for the quick reply. I was on travel. Here the replies. I

did not understand why there is not an include input of dpin.inc and

kpin.inc in solve.inp since it seems (as I understood) the idea of

coupling to two ref center node the circular surfaces of the hinge)

F Peugny

Here is the output of

$ccx solve2

(my OS is Windows I used a precompiled version 2.14 of CalculiX found

in bConverged...)

ERROR reading SOLID SECTION: *SOLID SECTION can

not be used for beam or shell elements

Faulty element: 1

*ERROR in gen3delem: first thickness

in node 1 of element 1

is zero

Then

$ccx solve

(I invoke it first but I could not understand why an include pin.inc

was called since it is writen with pre2.fbl --> I commented finally

the line *include, input=pin.inc when using

$ccx solve)

You are using an executable made on pon., 30 kwi 2018 09:46:04

The numbers below are estimated upper bounds

number of:

nodes: 70657

elements: 1071

one-dimensional elements: 3

two-dimensional elements: 1068

integration points per elements: 50

degree of freedom per nodes: 3

layers per element: 2

distributed facial loads: 20

distributed volumetric loads: 0

concentrated loads: 0

single point constraint: 78003

multiple point constraints: 128521

terms in all multiple point constraints: 822529

tie constraints: 0

dependent nodes tied by cyclic constraints: 0

dépendent nodes in pre-tension constraints: 0

sets: 11

terms in all sets: 6508

materials: 1

constants per material and temperature: 2

temperature points per material: 1

plastic data points per material: 0

orientation: 2136

amplitudes: 3

data points in all amplitudes: 3

print requests: 0

transformations: 0

property cards: 0

ERROR reading SOLID SECTION: *SOLID SECTION can

not be used for beam or shell elements

Faulty element: 1

*ERROR in gen3delem: first thickness

in node 1 of element 1

is zero

On 9/19/20, mkraska @.***> wrote:

Can you please attach the problem files and the console output of ccx? Also

some info on your OS would be helpful. It is not clear to me what actually

failed: CGX or CCX?

--

You are receiving this because you commented.

Reply to this email directly or view it on GitHub:

https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-695359281

--

Frédéric Peugny

06 50 32 92 74

Résidence les Berges du Midi, Bat. A, N°19

25 rue du Bosc https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g> 31150 FENOUILLET https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g> France https://www.google.com/maps/search/25+rue+du+Bosc+%0D%0D%0A+31150+FENOUILLET+%0D%0D%0A+France?entry=gmail&source=g

—

You are receiving this because you authored the thread.

Reply to this email directly, [

https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-699131878 |

view it on GitHub ] , or [

https://github.com/notifications/unsubscribe-auth/ADDC4DV2JBJKTJFKSPJXUXTSHT2P5ANCNFSM4DKEDOSQ

| unsubscribe ] .

—

You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/mkraska/CalculiX-Examples/issues/18#issuecomment-699158895, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADDYQV4BIX6GQ4IGLEXKR6TSHUB7PANCNFSM4DKEDOSQ .

-- Frédéric Peugny 06 50 32 92 74 Résidence les Berges du Midi, Bat. A, N°19 25 rue du Bosc 31150 FENOUILLET France

-- Frédéric Peugny 06 50 32 92 74 Résidence les Berges du Midi, Bat. A, N°19 25 rue du Bosc 31150 FENOUILLET France

mkraska commented 2 years ago

I don't see any attached files. Where do I find them?