CCX: kinematic coupling shows wrong stress distribution #18
Comments
|
Problem still present in 2.15. see example page |
|
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? |
|
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? |
|
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:
#18 (comment)
--
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
|
|
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 <[email protected]>:
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:
> #18 (comment)
--
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, [
#18 (comment) |
view it on GitHub ] , or [
https://github.com/notifications/unsubscribe-auth/ADDC4DV2JBJKTJFKSPJXUXTSHT2P5ANCNFSM4DKEDOSQ
| unsubscribe ] .
|
|
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 <[email protected]> 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:
> >
#18 (comment)
> --
> 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, [
>
#18 (comment)
|
> 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
<#18 (comment)>,
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
|
|
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:
>
>
>
> > >
> #18 (comment)
>
>
>
>
>
>
>
> > --
>
>
>
> > 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, [
>
>
>
> >
> #18 (comment)
> |
>
>
>
> > 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
> <#18 (comment)>,
> 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
|
|
I don't see any attached files. Where do I find them? |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Example Linear/StressConc1
Force plot and stress distribution indicate that forces are applied to nodes outside the specified surface
The text was updated successfully, but these errors were encountered: