Welcome to the forums. Please post in English or French.

You are not logged in.

#1 Re: Salome-Meca installation » Problem in Launching Freshly Installed Salome-Meca 2019 / Ubuntu 20.10 » 2021-01-18 06:07:01

Hi all, here is an update of my unfortunate endeavor.

Good news is Salome-Meca can now be launched after a proper driver update via the following terminal command:

$ sudo apt update
$ sudo apt upgrade
$ sudo apt install nvidia-driver-460
$ sudo reboot

Bad news is during meshing the following error occurs:
50847028103_0006b764fc_o.jpg

I can confirmed that package liblapack.so.3 had been installed.
Can anyone please explain what went wrong & how to solve it?

Thanks!

#2 Salome-Meca installation » Problem in Launching Freshly Installed Salome-Meca 2019 / Ubuntu 20.10 » 2021-01-13 02:23:38

hkboondoggle
Replies: 1

Hello, I am having trouble in launching a newly installed SM v2019.0.3-1-universal in my system.

The installation process strictly follows the detailed guide in #20 at the sticky thread _/viewtopic.php?id=24883 (Thanks sb1966!). Yet when launching SM via script, the program is automatically terminated in a flash after finishing the initiation, showing the following error prompt:

============Terminal Output============
asrc@ASRC-CalSer-01:~$ salome_meca/appli_V2019.0.3_universal/./salome 
runSalome running on ASRC-CalSer-01
Searching for a free port for naming service: 2814 - OK
Searching Naming Service  + found in 0.1 seconds 
Searching /Kernel/Session in Naming Service  ++Warning: this type (Study,objref) already exists, it will be ignored.
Warning: this type (SALOME_MED/MEDCouplingFieldDoubleCorbaInterface,objref) already exists, it will be ignored.
****************************************************************
     Warning: module NETGENPLUGIN is improperly configured!
     Module NETGENPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module HYBRIDPLUGIN is improperly configured!
     Module HYBRIDPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module HexoticPLUGIN is improperly configured!
     Module HexoticPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module GHS3DPLUGIN is improperly configured!
     Module GHS3DPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module BLSURFPLUGIN is improperly configured!
     Module BLSURFPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module GMSHPLUGIN is improperly configured!
     Module GMSHPLUGIN will not be available in GUI mode!
****************************************************************
 found in 1.0 seconds 
Start SALOME, elapsed time :   1.1 seconds
asrc@ASRC-CalSer-01:~$ 

For testing things out, I tried restoring the relocated libstdc++.so* files back into ~/salome_meca/V2019.0.3_universal/prerequisites/debianForSalome/lib.
With that, SM could be launched with another set of error prompt below ...

============Terminal Output============
asrc@ASRC-CalSer-01:~$ salome_meca/appli_V2019.0.3_universal/./salome 
runSalome running on ASRC-CalSer-01
Searching for a free port for naming service: 2815 - OK
Searching Naming Service  + found in 0.1 seconds 
Searching /Kernel/Session in Naming Service  +libGL error: MESA-LOADER: failed to open nouveau (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: nouveau
libGL error: MESA-LOADER: failed to open nouveau (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: nouveau
libGL error: MESA-LOADER: failed to open swrast (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: swrast
+Warning: this type (Study,objref) already exists, it will be ignored.
Warning: this type (SALOME_MED/MEDCouplingFieldDoubleCorbaInterface,objref) already exists, it will be ignored.
****************************************************************
     Warning: module NETGENPLUGIN is improperly configured!
     Module NETGENPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module HYBRIDPLUGIN is improperly configured!
     Module HYBRIDPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module HexoticPLUGIN is improperly configured!
     Module HexoticPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module GHS3DPLUGIN is improperly configured!
     Module GHS3DPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module BLSURFPLUGIN is improperly configured!
     Module BLSURFPLUGIN will not be available in GUI mode!
****************************************************************
****************************************************************
     Warning: module GMSHPLUGIN is improperly configured!
     Module GMSHPLUGIN will not be available in GUI mode!
****************************************************************
 found in 1.0 seconds 
Start SALOME, elapsed time :   1.1 seconds

... And giving out fatal error when launching Geometry, or instantly crashed once clicking into Mesh (without error prompt), AsterStudy, or ParaViS:
50830085372_2e7ce82bdb_o.jpg

============Terminal Output in launching AsterStudy============
AsterStudy: Wed-13-08:48:49.275 AsterStudy is activating...
AsterStudy: Wed-13-08:48:49.275 Creating workspace...
AsterStudy: Wed-13-08:48:49.283 Refreshing configuration on localhost...
AsterStudy: Wed-13-08:48:49.896 workspace is ready.
asrc@ASRC-CalSer-01:~$ Terminating SALOME on port 2815...
============Terminal Output in Launching ParaViS============
libGL error: MESA-LOADER: failed to open nouveau (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: nouveau
libGL error: MESA-LOADER: failed to open nouveau (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: nouveau
libGL error: MESA-LOADER: failed to open swrast (search paths /usr/lib/x86_64-linux-gnu/dri:\$${ORIGIN}/dri:/usr/lib/dri)
libGL error: failed to load driver: swrast

Seeing that it might be a GPU driver issue, I tried update the said driver ...

============Terminal Output============
asrc@ASRC-CalSer-01:~$ ubuntu-drivers devices
== /sys/devices/pci0000:00/0000:00:03.1/0000:2b:00.0 ==
modalias : pci:v000010DEd00001D01sv00001043sd000085F4bc03sc00i00
vendor   : NVIDIA Corporation
model    : GP108 [GeForce GT 1030]
driver   : nvidia-driver-450 - distro non-free
driver   : nvidia-driver-440-server - distro non-free
driver   : nvidia-driver-460 - distro non-free recommended
driver   : nvidia-driver-418-server - distro non-free
driver   : nvidia-driver-450-server - distro non-free
driver   : nvidia-driver-390 - distro non-free
driver   : xserver-xorg-video-nouveau - distro free builtin

asrc@ASRC-CalSer-01:~$ sudo ubuntu-drivers autoinstall
[sudo] password for asrc: 
Reading package lists... Done
Building dependency tree       
Reading state information... Done
0 upgraded, 0 newly installed, 0 to remove and 0 not upgraded.
asrc@ASRC-CalSer-01:~$ 

... which didn't help at all as the the aforementioned SM launch issue still persists.

Thus I am now totally clueless in how to solve this. Can anyone please help me out?
Many thanks!

Below are the info of my system should it helps:
50830103667_caa7ae1a77_o.jpg
CPU have passed a stress test in s-tui with a 5-hour session.
RAM sticks have passed Memtest86+ v5.0.1 with 2.5 test cycles in 14 hrs and found 0 error.

#3 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-11-05 02:55:02

sameer21101970 wrote:

It is converging, do below.

Make Sphere (ball)Mesh Fine.
Reduce Time Step to 100

Sammer, how small was the ball mesh size in your setting?
I ended up using this setting:
[ball = 0.1; substrate = 2.5] & number of time steps = 100

But even with this master-to-slave mesh size ratio in 1:2.5, it still gave out convergence issue.
I have attached the log of this latest trial FYI.
In which there was a multiple occasion of <A> <COMPOR2_9> occurred:

!-----------------------------------------------------------------!
   ! <A> <COMPOR2_9>                                                 !
   !                                                                 !
   ! Les déformations deviennent trop grandes : <6.382117E-02>       !
   ! => GROT_GDEP sous COMPORTEMENT n'est plus valide.               !
   !                                                                 !
   ! Pour un calcul en                                               !
   ! grandes déformations                                            !
   ! il faut utiliser GDEF_LOG ou SIMO_MIEHE.                        !
   !                                                                 !
   ! Pour un calcul hyper-élastique, utiliser ELAS_HYPER.            !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!

Is it something related to the convergence issue?

I am unable to test further shrinking the ball mesh size for the time being, as somehow Salome-Meca is giving out fatal error whenever loading up the mesh of all previous case files:
50568174221_aabf5fd653_o.jpg

Any suggestions other than re-installing Salome-Meca?


mib wrote:

maybe this

code-aster.org/doc/default/en/man_u/u2/u2.04.04.pdf (copy and paste in your browser)

from 2.2.3 Case general

to 2.2.7 Quality of the grid

Thanks for the info Mib.
So I suppose this is a hard-coded limitation of CA?
Then this is not a good news for me, because the physics in my case study includes deformation mechanism & local strain hardening effect at the ball / substrate under the same length order.
Are you suggesting that I should look for other simulation tools instead?

#4 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-10-22 08:40:49

Sorry Sameer, but I have to try out your latest suggestion at a later time.
My office desktop, where Code_Aster is installed, is now show BSOD issues. (Thanks Microsoft)

Could you please explain the reason why slave surface mesh has to be finer then that of master?

Also, I find your suggestion of reducing time steps surprising. Because I thought have a smaller time increment could aid the calculation in achieving convergence. Am I mistaken, or just this concept is not applicable in this calculation?

#5 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-10-20 02:16:48

Hello Sameer, thanks for your reply.

However, your suggestion is exactly what I did as in my previous description:

4.) ... and disable "Merge coincident nodes and elements".

So right now, I'm unsure if the latest convergence issue is in fact mesh related or not.

Did you also used compound mesh as well in the reference case you shared on 08/25?

I notice that the mesh size of the substrate is larger than that of the ball in the reference case. Does this factor plays a role in the convergence of the calculation?

#6 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-10-19 03:12:52

I tried to avoid the issue of identical node at the initial contact point where the bottom of the ball meets with the top surface of the substrate through the following steps:

1.) Individually create the mesh of ball and particle respectively

2.) Ensure no existence of duplicating node or mesh element within each individual mesh.

3.) Create groups of node / surface / volume under each individual mesh with their corresponding geometries.

4.) Create a compound mesh from the individual meshes, enable the option "Create groups from input objects", and disable "Merge coincident nodes and elements".

While the error prompt of <A> <CALCULEL5_49>  under LIAISON_MAIL no longer occur, the calculation is still plagued by the convergence issue where the time increment kept on decreasing and finally trigger the following:

!-----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ADAPTATION_11>                                                       !
   !                                                                                   !
   !  The value of the time step selected   1.862645169939e-13 is lower than PAS_MINI. !
   !-----------------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
   !                                                                 !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
   !                                                                 !
   !    If an ABNORMAL_ABORT occured, that is not reproductible      !
   !    on official Linux version, or if you need some help to       !
   !    understand the message above, please report issues at:       !
   !                                                                 !
   !                      support@simulease.com                      !
   !                                                                 !
   !-----------------------------------------------------------------!

Any suggestions / ideas in handling the latest issue? Thanks!
As usual, the comm, log & med files of this latest trial are attached FYI.

@sammer: Would you mind showing me the details of the meshing procedure in your reference case?

#7 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-10-15 04:32:09

Grateful if anyone can assist.

After detail reviewing of the calculation log, one of the possible cause of the singularity matrix problem is identified with the following error prompt during the setup of adhesion relationship in AFFE_CHAR_MECA:

Création du fichier au format MED 3.3.1.
   
   !------------------------------------------------------------------------------!
   ! <A> <CALCULEL5_49>                                                                      !
   !                                                                                                             !
   !  LIAISON_MAIL:                                                                                 !
   !  The linear relation intended to eliminate the slave node N10 is a tautology !
   !  because the mesh Master in with                                                                     !
   ! respect to this node has this same node in its connectivity.                           !
   !  It is thus not written.                                                                                        !
   !                                                                                                             !
   !                                                                                                             !
   ! This is a warning. If you do not understand the meaning of this    !
   !  warning, you can obtain unexpected results!                                  !
   !------------------------------------------------------------------------------!

Node N10 appears to be the initial contact point where the bottom of the ball meets with the top surface of the substrate at T = 0. Is this a problem mesh creation? And if so, how should this be solved?



Moreover, there seems to be other hiccups in the mesh I am using.
This is reflected by the multiple (5 times) warnings during AFFE_MODELE:

!----------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                                       !
   !                                                                                                                            !
   !  - > the mesh M155 carries an edge finite element, but it does not border  !
   !  any element having" a rigidity".                                                                     !
   !                                                                                                                            !
   !  - > Risks &                                                                                                        !
   ! advices:                                                                                                             !
   !  That can involve problems of" null pivot" at the time of the resolution.      !
   !  If the resolution of the linear systems does                                                   !
   ! not pose problems, you                                                                                    !
   !  can be unaware of this message.                                                                    !
   !  If not, check the definition of model (AFFE_MODELE) by avoiding the use  !
   !                                                                                                                            !
   ! of operand TOUT='OUI'.                                                                                   !
   !                                                                                                                            !
   !                                                                                                                            !
   ! This is a warning. If you do not understand the meaning of this                   !
   !  warning, you can obtain unexpected results!                                                !
   !----------------------------------------------------------------------------------------!

What does this error suppose to mean?
Mesh elements in question (M155 ~ M159) locate at -x edge of the top surface of the substrate.
The whole mesh has no duplicating edge or mesh element.

#8 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-09-28 07:05:35

Hello, could some one please help me out?

As stated in my previous post, I am facing issue of singular matrix whenever repeating the calculation using a finer mesh.
The mesh is entirely NetGen 1D-2D-3D, and has no coincident nodes / elements. I have included the dump of mesh info FYI.

The comm, log & med files are attached in my previous post above.
I am running code aster version 14.4.0 on Windows-10-10.0.17134-SP0.

Many thanks.

#9 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-09-04 03:14:33

@sammer: Thanks for sharing the your reference case! It did converge nicely. I am trying to replicate the settings to my existing study for having a finer mesh resolution on the substrate surface, though it turns out having issue of singular matrix (the comm, log & med files are attached FYI). Would you mind giving it a look? Thanks!

@mueller: Thanks for your input! I did try out "GROT_GDEP" which successfully launched the calculation. Do you mean that whether "SIMO_MIEHE" is applicable or not depends on the settings of "DEFI_MATERIAU"?

#10 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-08-25 04:17:41

Hi sameer, sorry for the late replay. The lock down of the campus screws up my work flow.

I have re-done the geometry & meshing based on your recommendation (as shown in attached figure below). However, it is now reporting an issue in "DEFI_CONTACT":

!--------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <CONTACT2_13>                                                                                          !
   !                                                                                                                    !
   ! Contact methods with a grid.                                                                                       !
   ! The contact zone number 1 contains 1 nodes common to surfaces Masters and slaves.                                  !
   ! Check the                                                                                                          !
   ! definition of your contact surfaces or inform one of key words SANS_NOEUD/SANS_GROUP_NO/SANS_MAILLE/SANS_GROUP_MA. !
   !--------------------------------------------------------------------------------------------------------------------!

I tried adding this additional settings in "DEFI_CONTACT":
SANS_GROUP_NO=('Particle_Face', )

But then it led to another fatal error when initializing the "DYNA_NON_LINE":

.. __stg1_cmd9:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0012 Concept of the type:  evol_noli
  # ------------------------------------------------------------------------------------------
  solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         COMPORTEMENT=_F(TOUT='OUI',
                                         RELATION='VMIS_JOHN_COOK',
                                         DEFORMATION='SIMO_MIEHE',
                                         RESI_CPLAN_RELA=1.E-06,
                                         RESI_INTE_RELA=1.E-06,
                                         ITER_INTE_MAXI=20,
                                         ITER_CPLAN_MAXI=1,
                                         SYME_MATR_TANG='OUI',
                                         PARM_THETA=1.0,),
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F(LIST_INST=timeline,
                                      PRECISION=1.E-06,),
                         SCHEMA_TEMPS=_F(SCHEMA='HHT',
                                         ALPHA=-0.6,
                                         MODI_EQUI='NON',
                                         FORMULATION='DEPLACEMENT',
                                         COEF_MASS_SHIFT=0.0,),
                         SOLVEUR=_F(METHODE='MUMPS',
                                    RENUM='AUTO',
                                    NPREC=8,
                                    ELIM_LAGR='LAGR2',
                                    STOP_SINGULIER='OUI',
                                    TYPE_RESOL='AUTO',
                                    ACCELERATION='AUTO',
                                    LOW_RANK_SEUIL=0.0,
                                    PRETRAITEMENTS='AUTO',
                                    POSTTRAITEMENTS='AUTO',
                                    PCENT_PIVOT=20,
                                    RESI_RELA=-1.0,
                                    GESTION_MEMOIRE='AUTO',
                                    FILTRAGE_MATRICE=-1.0,
                                    MIXER_PRECISION='NON',
                                    MATR_DISTRIBUEE='NON',),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F(REAC_INCR=1,
                                   MATRICE='TANGENTE',
                                   REAC_ITER=1,
                                   REAC_ITER_ELAS=0,
                                   MATR_RIGI_SYME='NON',),
                         CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                        ITER_GLOB_ELAS=25,
                                        ARRET='OUI',),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F(PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INFO=1,)

../../../bibc/supervis/aster_module.c 2343 : Echec lors de la creation du comportement (lccree/create) !
EXECUTION_CODE_ASTER_EXIT_4830=0
<INFO> Code_Aster run ended, diagnostic : <F>_ABNORMAL_ABORT

Any suggestions / ideas in handling the latest issue? Thanks!
As usual, the comm, log & med files of this latest trial are attached FYI.

50265984777_4da9714359_o.jpg

#11 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-08-03 08:03:44

sameer, would you mind sharing the details of your test run?

'Cause I have managed to get 'Liaison_Mail' implemented in my latest trial, but it is still giving out  <MECANONLINE9_3> under 'VMIS_JOHN_COOK'.

Also, switching DEFORMATION of behavior of DYNA_NON_LINE from 'GROT_GDEP' to 'SIMO_MIEHE' still instantly leads to <F>_ABNORMAL_ABORT with EXECUTION_CODE_ASTER_EXIT_8627=0.

#12 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-07-31 09:25:00

Hello sameer, thank you for your reply.
May I ask a few questions about the 'Liaison_Mail' approach:

1. I was bumped with the following problem when trying out 'Liaison_Mail'. Can you please explain what the matter is?

.. __stg1_cmd6:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0009 Concept of the type:  char_meca
  # ------------------------------------------------------------------------------------------
  BC_Fix = AFFE_CHAR_MECA(MODELE=model1,
                          DDL_IMPO=_F(GROUP_NO=('Substrate_Bottom', ),
                                      DX=0.0,
                                      DY=0.0,
                                      DZ=0.0,),
                          LIAISON_MAIL=_F(GROUP_MA_MAIT=('Substrate_Top', ),
                                          GROUP_MA_ESCL=('Particle_Face', ),
                                          TYPE_RACCORD='MASSIF',
                                          ELIM_MULT='NON',),
                          VERI_NORM='OUI',
                          DOUBLE_LAGRANGE='OUI',
                          INFO=1,)

   
   !-----------------------------------------------------------------------------!
   ! <EXCEPTION> <CALCULEL4_55>                                                  !
   !                                                                                                          !
   !  Pas de mailles à projeter ou en correspondance.                            !
  !                                                                                                          !
   !  Dans le cas de l'utilisation de AFFE_CHAR_MECA / LIAISON_MAIL, les mailles !
   ! maîtres                                                                                        !
   !  doivent avoir la même dimension que l'espace de modélisation :             !
   !  - mailles volumiques pour un modèle 3D                                     !
   !  - mailles                                                                                           !
   ! surfaciques pour un modèle 2D                                                     !
   !-----------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
  !                                                                                                          !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
  !                                                                                                          !
   !    If an ABNORMAL_ABORT occured, that is not reproductible      !
   !    on official Linux version, or if you need some help to              !
   !    understand the message above, please report issues at:       !
  !                                                                                                          !
   !                      support@simulease.com                                         !
  !                                                                                                          !
   !-----------------------------------------------------------------!

Destruction of the concept  BC_Fix.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR


2. Would you please elaborate on your suggestion "Apply Velocity or acceleration to ball" ? 'Cause I have already set the initial velocity of the ball in the following command:

init_vel = CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
                        OPERATION='AFFE',
                        MODELE=model1,
                        AFFE=_F(GROUP_NO=('Particle_Solid', ),
                                NOM_CMP=('DX', 'DY', 'DZ'),
                                VALE=(0.0, 0.0, -1.765),),
                        INFO=1,)

solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         COMPORTEMENT=_F( ... ),
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F( ... ),
                         SCHEMA_TEMPS=_F( ... ),
                         SOLVEUR=_F( ... ),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F( ... ),
                         CONVERGENCE=_F( ... ),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F( ... ),
                         INFO=1,)


3. How does setting up 'Liaison_Mail' help tackling the convergence issue of 'DYNA_NON_LINE' with 'VMIS_JOHN_COOK' enabled?

Thanks again!

#14 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-07-17 06:52:22

Thanks Konyaro for your inspection.

I have shorten the calculation time from 10 to 3.24 (<DVP_1> would happen at 3.26), and the calculation does converge well.
Although the result is physically incorrect again with the ball re-bounces backward after the collision, whereas the reference I'm following [1] has the ball undergoes plastic deformation and stays attached to the substrate.

The material properties and initial velocity of the ball in my code-aster simulation follow the presented data in [1] and are the normalized with respect to the density of copper, the speed of sound, the diameter of the ball, and the room temperature.

After inspection, I found that the RELATION of behavior of DYNA_NON_LINE by default is 'ELASTIC'.
Thinking that this might be the reason of the discrepancy of the calculated result, I changed it to 'VMIS_JOHN_COOK'.
Other changes included halving the time increment to 5e-3, and enabling Deformation in CLAC_CHAMP.

Now it appears to be having trouble in convergence, as time increment automatically down-sizing itself to the order of -16 once the ball impacts the substrate &  'VMIS_JOHN_COOK' kicks into action, which ultimately leads to <MECANONLINE9_3>.
The comm, log & med files of this latest trial are attached FYI.

I also tried switching DEFORMATION of behavior of DYNA_NON_LINE from 'GROT_GDEP' to 'SIMO_MIEHE', as high strain rate deformation is expected upon the collision.
The calculation refused to launch and instantly gave out <F>_ABNORMAL_ABORT with EXECUTION_CODE_ASTER_EXIT_8627=0.

Any suggestions / ideas in handling the latest issue?

[1] T. Schmidt, F. Gärtner, H. Assadi, et al., "Development of a generalized parameter window for cold spray deposition", Acta Materialia, 2006, 54, pp 729-742.

#15 Re: Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-07-03 09:40:59

Thinking that the problem might be caused by the stretching of the mesh, I have replace the whole mesh set up with an uniform setting shown below:

50071164116_1540215690_o.jpg

However, the same exception still occurs.
The comm, log & med files of this latest trial are attached FYI.

#16 Code_Aster usage » Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE) » 2020-06-29 06:13:27

hkboondoggle
Replies: 29

Hi all,

I am trying to simulate an impact scenario with a ball travelling on negative z-axis direction and being about to collide on a static plate.

With the aid of konyaro, Johannes, mecour, and others (I really appreciate it guys!), I have set up the case which details can be found in the attached zip file.
It includes the .py script for geometry & mesh generation (.med too large), the .comm file and .log of the problematic Run case.

The following error happened as the Run Case was being launched:

.. __stg1_cmd9:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0012 Concept of the type:  evol_noli
  # ------------------------------------------------------------------------------------------
  solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F(LIST_INST=timeline,
                                      PRECISION=1.E-06,),
                         SCHEMA_TEMPS=_F(SCHEMA='HHT',
                                         ALPHA=-0.3,
                                         MODI_EQUI='NON',
                                         FORMULATION='DEPLACEMENT',
                                         COEF_MASS_SHIFT=0.0,),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F(REAC_INCR=1,
                                   MATRICE='TANGENTE',
                                   REAC_ITER=1,
                                   REAC_ITER_ELAS=0,
                                   MATR_RIGI_SYME='NON',),
                         CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                        ITER_GLOB_ELAS=25,
                                        ARRET='OUI',),
                         SOLVEUR=_F(RENUM='AUTO',
                                    NPREC=8,
                                    ELIM_LAGR='LAGR2',
                                    STOP_SINGULIER='OUI',
                                    TYPE_RESOL='AUTO',
                                    ACCELERATION='AUTO',
                                    LOW_RANK_SEUIL=0.0,
                                    PRETRAITEMENTS='AUTO',
                                    POSTTRAITEMENTS='AUTO',
                                    PCENT_PIVOT=20,
                                    RESI_RELA=-1.0,
                                    GESTION_MEMOIRE='AUTO',
                                    FILTRAGE_MATRICE=-1.0,
                                    MIXER_PRECISION='NON',
                                    MATR_DISTRIBUEE='NON',
                                    METHODE='MUMPS',),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F(PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INFO=1,)


Comme vous n'avez pas défini explicitement le comportement, tout le modèle est supposé élastique en petites perturbations.
   Liste des comportements
   Affecté sur 186954 éléments
     Relation                             : ELAS
     Déformation                          : PETIT
     Nombre total de variables internes   : 1
            V1 : VIDE
  Le système linéaire à résoudre a 734008 degrés de liberté:
   - 734008 sont des degrés de liberté physiques
     (ils sont portés
par 244274 noeuds du maillage)
   - 0 sont les couples de paramètres de Lagrange associés
     aux 0 relations linéaires dualisées.
  La matrice est de taille 734008 équations.
  Elle contient 31072249 termes non nuls si elle est symétrique et 61410490 termes non
nuls si elle n'est pas symétrique.
  Soit un taux de remplissage de   0.011 %.
    Lecture de l'état initial
      Il n'y a pas d'état initial défini. On prend un état initial nul.
      Le champ <DEPL> est initialisé a zéro
      Le champ <SIEF_ELGA> est initialisé a zéro
      Le champ <VARI_ELGA> est initialisé a zéro
   
   !-----------------------------------------------------------------!
   ! <A> <ALGELINE7_20>                                              !
   !                                                                 !
   !  On ne peut pas remplir la composante DX du noeud numéro N1     !
   !  pour le champ "&&NMETL2.CHAMP.CONV"                            !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------------!
   ! <A> <MECANONLINE_2>                                                                    !
   !                                                                                        !
   !  Lors de la recopie du champ &&NMETL2.CHAMP.CONVER donné dans ETAT_INIT de la commande !
   !  STAT_NON_LINE vers le champ                                                           !
   ! &&NMCH1P.VITMOI, certaines composantes de &&NMCH1P.VITMOI                              !
   !  ont du être mises à zéro.                                                             !
   !                                                                                        !
   !  Ce problème survient lorsque le champ donné                                           !
   ! dans ETAT_INIT ne comporte                                                             !
   !  pas assez de composantes, on complète donc par des zéros.                             !
   !                                                                                        !
   !                                                                                        !
   ! This is a warning. If you do not understand the meaning of this                        !
   !  warning, you can obtain unexpected results!                                           !
   !----------------------------------------------------------------------------------------!
   
      Le champ <VITE> est lu dans ETAT_INIT, par un champ donné explicitement
      Le champ <ACCE> est initialisé a zéro

One considers an acceleration initial.

The initial state does not have acceleration given.
It is computed.


Filing of the initial state

Filing of the fields
    Field stored  DEPL at time  0.000000000000e+00 for the sequence number  0
    Field stored  SIEF_ELGA at time  0.000000000000e+00 for the sequence number  0
    Field stored  VARI_ELGA at time  0.000000000000e+00 for the sequence number  0
    Field stored  COMPORTEMENT at time  0.000000000000e+00 for the sequence number  0
    Field stored  VITE at time  0.000000000000e+00 for the sequence number  0
    Field stored  ACCE at time  0.000000000000e+00 for the sequence number  0
    Field stored  CONT_NOEU at time  0.000000000000e+00 for the sequence number  0
   Il y a 0 points initialement en contact et 0 points exclus.
-----------------------------------------------------------------------------------------------------------------------------------------

Time of computation:   4.000000000000e-03
-----------------------------------------------------------------------------------------------------------------------------------------
|     CONTACT    |     NEWTON     |     RESIDU     |     RESIDU     |     OPTION     |     CONTACT    |     CONTACT    |     CONTACT    |
|    BCL. GEOM.  |    ITERATION   |     RELATIF    |     ABSOLU     |   ASSEMBLAGE   |   NEWTON GENE  |    PRESSURE    |     CRITERE    |
|    ITERATION   |                | RESI_GLOB_RELA | RESI_GLOB_MAXI |                |   VARI. CONT.  |    ERROR       |    VALEUR      |
-----------------------------------------------------------------------------------------------------------------------------------------
Traceback returned by pdbhelp:
[0] print_trace_
[1] utmess_core_
[2] utmess_
[3] assert_
[4] mmcaln_
[5] mmvppe_
[6] te0365_
[7] te0000_
[8] calcul_
[9] nmelcv_
[10] nmfocc_
[11] nmprta_
[12] nmpred_
[13] nmnewt_
[14] op0070_
[15] execop_
[16] expass_
[17] initaster_fonctions
[18] PyMethodDef_RawFastCallKeywords
[19] PyMethodDef_RawFastCallKeywords
[20] PyEval_EvalFrameDefault
[21] PyMethodDef_RawFastCallKeywords
[22] PyEval_EvalFrameDefault
[23] PyMethodDef_RawFastCallKeywords
[24] PyEval_EvalFrameDefault
   
   !-----------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <DVP_1>                                                                           !
   !                                                                                               !
   ! Program error.                                                                                !
   !                                                                                               !
   ! Condition not met:                                                                            !
   !     ASTER_FALSE                                                                               !
   ! File                                                                                          !
   ! /home/siavelis/PROJETS/aster/windows/build/codeaster-src/bibfor/cont_elem/mmcaln.F90, line 63 !
   !                                                                                               !
   ! --------------------------------------------                                                  !
   ! Contexte du message :                                                                         !
   !    Option         : CHAR_MECA_CONT                                                            !
   !    Type d'élément : COT6T6                                                                    !
   !                                                                                               !
   ! Maillage       : mesh                                                                         !
   !    Maille         : XXX                                                                       !
   !    Type de maille : TRIA66                                                                    !
   !    Cette maille appartient aux groupes de mailles suivants                                    !
   ! :                                                                                             !
   !                                                                                               !
   !    Position du centre de gravité de la maille :                                               !
   !       x=nan y=0.000000 z=0.000000                                                             !
   !                                                                                               !
   !                                                                                               !
   !                                                                                               !
   ! Il y a probablement une erreur dans la programmation.                                         !
   ! Veuillez contacter votre assistance technique.                                                !
   !-----------------------------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
   !                                                                 !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
   !                                                                 !
   !    If an ABNORMAL_ABORT occured, that is not reproductible      !
   !    on official Linux version, or if you need some help to       !
   !    understand the message above, please report issues at:       !
   !                                                                 !
   !                      support@simulease.com                      !
   !                                                                 !
   !-----------------------------------------------------------------!

Could you please advice in dealing this problem?
Thanks!

I have attached the mesh info below FYI:

50056298003_4ee4aa89f8_o.jpg

50057117222_c947c2115b_o.jpg

#17 Introduce yourself / Présentez vous » Greetings from Hong Kong » 2020-06-22 04:59:12

hkboondoggle
Replies: 0

Hi all,

I'm hkboondoggle & I'm a PhD candidate.

I'm attempting to simulate a spray particle colliding & attaching itself to a substrate with code_aster. High strain rate deformation on both entities are expected upon the collision.

Nice to meet you!

Regards,
hkboondoggle

#18 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-22 04:52:20

Wow, I wasn't expecting to receive such attention to my help request throughout the weekend.
Thanks for all of your enthusiastic feedback every one! I really appreciate it!

First of all, I have to admit that I'm a complete doolie for the Salome-Meca / Code_aster platform.
So during the geometry creation & meshing stage, I basically followed this video tutorial:
youtube .com / watch?v=eOvjusbQ_RE

Now, with the disclaimer done, may I inquire the following:

1.

konyaro wrote:

there are a few mistakes in your model and the main issue is that you're meshing solids AND faces. You must simply define groups in the geometry and not duplicate faces.

Konyaro, could you please clarify? Because aren't face mesh is required for "DEFI_CONTACT" ?

Or do you mean having one should assign the contact rule to the 3D elements of the ball & substrate instead?
I was approaching this way but encounter problems (see: viewtopic.php?id=24974). You did mention that "I need surfaces or edges in 3D"  for "DEFI_CONTACT" in that post. So it bears the question: how do one obtain surfaces or edges in 3D?


2.

AsterO'dactyle wrote:

For these kind of element: probably degenerated "flat" elements.

konyaro wrote:

you mesh has quadrangles only on the surface. I corrected a few points in your comm files.

I just wanted to point out that one may think that a Netgen mesh with the option "Quad-dominated" produces a nice hexahedral mesh. This is not the case, the elements are pyramids connected to tetrahedrons.

Given that, should I use Netgen ONLY for meshing the ball? Or should I avoid using Netgen at all?


@Anirudh: I'm attempting to simulate a spray particle colliding & attaching itself to a substrate. High strain rate deformation on both entities are expected upon the collision. I'll introduce myself there right away.


Many thanks guys!

#19 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-19 09:23:23

Hi mecour, appreciate for your feedback!

I've merged the nodes (tolerance = 5e-4) as well as the mesh and gave it a shot.

Yet, the same problem of <EXCEPTION> <FACTOR_78> still persists. Any thoughts?

The Message Log of the latest run is attached FYI.

Thanks!

#20 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-18 10:22:00

Hello, could some one please help me out?

The links of the mesh file are attached FYI in case it helps the 'dissection'.
I am using mega.nz to host them since they exceed the size limit of attachment.

Mesh file: https: // mega.nz / file / HcoXHIxD#tp3AeMwVUqwFjUy2hEwiFzFT2v4aAWEU92Wb7nuR334  (38.7MB)

I am running code aster version 14.4.0 on Windows-10-10.0.17134-SP0.

Many thanks.

#21 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-17 06:58:36

Johannes, would the exception be mesh related as well?
I'm wondering this because the log of the study has multiple warning about the mesh which I don't quite understand:

.. __stg1_cmd0:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0002 Concept of the type:  maillage_sdaster
  # ------------------------------------------------------------------------------------------
  mesh = LIRE_MAILLAGE(UNITE=20,
                       INFO_MED=1,
                       FORMAT='MED',
                       VERI_MAIL=_F(VERIF='OUI',
                                    APLAT=1.E-3,),
                       INFO=1,)

Vérification du maillage.
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M22 et M1 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M23 et M2 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M24 et M3 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M25 et M4 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M26 et M5 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !--------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                    !
   !                                                                                      !
   ! The warning message  MODELISA8_13 was emitted 5 time, it will not be displayed more. !
   !                                                                                      !
   !                                                                                      !
   ! This is a warning. If you do not understand the meaning of this                      !
   !  warning, you can obtain unexpected results!                                         !
   !--------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------!
   ! <A> <MODELISA4_8>                                                             !
   !                                                                               !
   !  - > Phase of checking of the mesh: presence of double meshes (or triple,')   !
   !  The multiple meshes are meshes of different names                            !
   ! which have same connectivity                                                  !
   !  (they lean on the same nodes).                                               !
   !                                                                               !
   !  - > Risks & advices:                                                         !
   !  The risk is to modelize 2 times (or more)                                    !
   ! space. One can for example have                                               !
   !  a model 2 times too heavy or 2 times too rigid.                              !
   !  Note: the meshes concerned are printed in the                                !
   ! file" message".                                                               !
   !  On this mesh, it is imprudent to affect quantities with key word TOUT='OUI'. !
   !                                                                               !
   !                                                                               !
   ! This is a warning. If you do not understand the meaning of this               !
   !  warning, you can obtain unexpected results!                                  !
   !-------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M19392 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000429, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M20423 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000138, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M20942 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000007, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M21441 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000005, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M22848 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000695, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !--------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                    !
   !                                                                                      !
   ! The warning message  MODELISA8_15 was emitted 5 time, it will not be displayed more. !
   !                                                                                      !
   !                                                                                      !
   ! This is a warning. If you do not understand the meaning of this                      !
   !  warning, you can obtain unexpected results!                                         !
   !--------------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------!
   ! <A> <MODELISA4_9>                                                    !
   !                                                                      !
   !  - > Phase of checking of the mesh: presence of degenerated meshes.  !
   !  The relationship between the smallest edge on largest is            !
   ! lower than 0,001                                                     !
   !                                                                      !
   !  - > Risks & advices:                                                !
   !  Check your mesh. The presence of such meshes can lead to            !
   !  problems of convergence and harm                                    !
   ! the quality of the results.                                          !
   !                                                                      !
   !                                                                      !
   ! This is a warning. If you do not understand the meaning of this      !
   !  warning, you can obtain unexpected results!                         !
   !----------------------------------------------------------------------!
   

------------ MAILLAGE mesh     - IMPRESSIONS NIVEAU  1 ------------

MA-16-JUIN-2020 09:31:03                                                       

NOMBRE DE NOEUDS                      519665

NOMBRE DE MAILLES                     365355
                              SEG3                  898
                              TRIA6                  28
                              QUAD8                8281
                              TETRA10            347867
                              PYRAM13              8281

NOMBRE DE GROUPES DE NOEUDS                6
                              GR_1_Particle_Surface             2341
                              Particle_Surface                  2341
                              GR_1_Substrate_Surface            4961
                              Substrate_Surface                 4961
                              GR_1_Substrate_Bottom             4961
                              Substrate_Bottom                  4961

NOMBRE DE GROUPES DE MAILLES               4
                              Substrate                       336925
                              Particle                         19223
                              Substrate_Shell                   7520
                              Particle_Shell                     789

--------------------------------------------------------------------------------

# Memory (Mo):    216.99 /    162.55 /    163.38 /    107.86 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0002 user+syst:          2.77S (syst:          0.16S, Elaps:          0.00S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd1:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0004 Concept of the type:  modele_sdaster
  # ------------------------------------------------------------------------------------------
  model1 = AFFE_MODELE(MAILLAGE=mesh,
                       AFFE=_F(TOUT='OUI',
                               PHENOMENE='MECANIQUE',
                               MODELISATION=('3D', ),),
                       INFO=1,
                       DISTRIBUTION=_F(PARTITIONNEUR='METIS',
                                       METHODE='SOUS_DOMAINE',),
                       VERI_JACOBIEN='OUI',)


Sur les 365355 mailles du maillage mesh, on a demandé l'affectation de 365355, on a pu en affecter 365355
Modélisation     Type maille  Élément fini     Nombre
                  SEG3         MECA_ARETE3      898
                  TRIA6        MECA_FACE6       28
                  QUAD8        MECA_FACE8       8281
3D               PYRAM13      MECA_PYRAM13     8281
3D               TETRA10      MECA_TETRA10     347867
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9245 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9292 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9588 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9591 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9661 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                   !
   !                                                                                     !
   ! The warning message  MAILLAGE1_1 was emitted 5 time, it will not be displayed more. !
   !                                                                                     !
   !                                                                                     !
   ! This is a warning. If you do not understand the meaning of this                     !
   !  warning, you can obtain unexpected results!                                        !
   !-------------------------------------------------------------------------------------!

.. __stg1_cmd7:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0010 Concept of the type:  char_contact
  # ------------------------------------------------------------------------------------------
  contact = DEFI_CONTACT(MODELE=model1,
                         STOP_INTERP='NON',
                         LISSAGE='NON',
                         VERI_NORM='OUI',
                         ITER_GEOM_MAXI=10,
                         RESI_GEOM=1.E-2,
                         ALGO_RESO_GEOM='POINT_FIXE',
                         REAC_GEOM='AUTOMATIQUE',
                         ITER_CONT_MULT=4,
                         STOP_SINGULIER='OUI',
                         NB_RESOL=10,
                         ITER_GCP_MAXI=0,
                         RECH_LINEAIRE='ADMISSIBLE',
                         PRE_COND='SANS',
                         ZONE=_F(APPARIEMENT='MAIT_ESCL',
                                 GROUP_MA_MAIT='Substrate_Shell',
                                 GROUP_MA_ESCL='Particle_Shell',
                                 NORMALE='MAIT_ESCL',
                                 VECT_MAIT='AUTO',
                                 VECT_ESCL='AUTO',
                                 GLISSIERE='NON',
                                 TYPE_PROJECTION='ORTHOGONALE',
                                 DIST_POUTRE='NON',
                                 DIST_COQUE='NON',
                                 DIST_APPA=-1.0,
                                 TOLE_PROJ_EXT=0.5,
                                 RESOLUTION='OUI',
                                 ALGO_CONT='CONTRAINTE',),
                         FORMULATION='DISCRETE',
                         FROTTEMENT='SANS',)

T
T
   
   !------------------------------------------------------------------------------------------------------------------------------------!
   ! <A> <CONTACT_8>                                                                                                                    !
   !                                                                                                                                    !
   ! Formulation discrète du contact.                                                                                                   !
   !  -> Il y a des éléments de type QUAD8 sur les surfaces de contact. Ces éléments ne permettent pas                                  !
   ! de respecter exactement la condition de contact.                                                                                   !
   !     Afin d'empêcher une pénétration intempestive des surfaces, on a procédé à des                                                  !
   ! liaisons cinématiques (LIAISON_DDL) entre les noeuds milieux et les noeuds sommets, sur les deux surfaces (maître et esclave).     !
   !                                                                                                                                    !
   !                                                                                                                                    !
   ! Risques et conseils :                                                                                                              !
   !    - Ces liaisons supplémentaires peuvent provoquer des incompatibilités avec les conditions limites en                            !
   ! particulier dans le cas de symétries, ce qui se traduira par une matrice singulière.                                               !
   !      Dans ce cas, il est possible de                                                                                               !
   ! n'appliquer les conditions aux limites concernées que sur les noeuds sommets (on pourra utiliser la commande DEFI_GROUP pour créer !
   ! le groupe de noeuds adéquat).                                                                                                      !
   !    - Il est toutefois préférable d'utiliser sur les surfaces de contact des éléments de type QUAD9.                                !
   ! Pour cela, transformez les éléments volumiques HEXA20 en HEXA27 ou PENTA15 en PENTA18 (par exemple avec la commande                !
   ! CREA_MAILLAGE).                                                                                                                    !
   !                                                                                                                                    !
   !                                                                                                                                    !
   ! This is a warning. If you do not understand the meaning of this                                                                    !
   !  warning, you can obtain unexpected results!                                                                                       !
   !------------------------------------------------------------------------------------------------------------------------------------!

#22 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-16 04:18:06

Hi Johannes, Thanks for your feedback.

I've corrected the the initial velocity field as instructed.

init_vel = CREA_CHAMP(identifier='8:1',
                      AFFE=_F(GROUP_MA=('Particle_Shell', 'Particle'),
                              GROUP_NO=('Particle_Surface', ),
                              NOM_CMP=('DX', 'DY', 'DZ'),
                              VALE=(0.0, 0.0, -1.765)),
                      MODELE=model1,
                      OPERATION='AFFE',
                      TYPE_CHAM='NOEU_DEPL_R')

However, the same problem of <EXCEPTION> <FACTOR_78> still persists. Could you please advice?

Thanks!

#23 Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE) » 2020-06-15 07:07:29

hkboondoggle
Replies: 15

Hi all,

I am new to code_aster, and I am trying to simulate an impact scenario with a ball travelling on negative z-axis direction and being about to collide on a static plate.

This is code of the case setting:

DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(identifier='0:1',
                     UNITE=20)

model1 = AFFE_MODELE(identifier='1:1',
                     AFFE=_F(MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE',
                             TOUT='OUI'),
                     MAILLAGE=mesh)

Cu_L_ISO = DEFI_MATERIAU(identifier='2:1',
                         ECRO_COOK=_F(A=0.0631779,
                                      B=0.283815,
                                      C=0.025,
                                      EPSP0=1e-06,
                                      M_PUIS=1.09,
                                      N_PUIS=0.31,
                                      TMELT=4.708333,
                                      TROOM=1.0),
                         ELAS=_F(ALPHA=5.903e-08,
                                 COEF_AMOR=0.01,
                                 E=106.92,
                                 NU=0.31,
                                 RHO=1.0),
                         THER=_F(LAMBDA=0.006422,
                                 RHO_CP=0.959))

fieldma0 = AFFE_MATERIAU(identifier='3:1',
                         AFFE=_F(MATER=(Cu_L_ISO, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model1)

listr = DEFI_LIST_REEL(identifier='4:1',
                       DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=10.0,
                                     NOMBRE=2500),
                       TITRE='Time_step')

timeline = DEFI_LIST_INST(identifier='5:1',
                          DEFI_LIST=_F(LIST_INST=listr),
                          METHODE='MANUEL')

BC_Fix = AFFE_CHAR_MECA(identifier='6:1',
                        DDL_IMPO=_F(DX=0.0,
                                    DY=0.0,
                                    DZ=0.0,
                                    GROUP_NO=('Substrate_Bottom', )),
                        MODELE=model1)

contact = DEFI_CONTACT(identifier='7:1',
                       MODELE=model1,
                       ZONE=_F(APPARIEMENT='MAIT_ESCL',
                               GROUP_MA_ESCL=('Particle_Shell', ),
                               GROUP_MA_MAIT=('Substrate_Shell', ),
                               NORMALE='MAIT_ESCL',
                               VECT_ESCL='AUTO',
                               VECT_MAIT='AUTO'))

init_vel = CREA_CHAMP(identifier='8:1',
                      AFFE=_F(GROUP_MA=('Particle', ),
                              GROUP_NO=('Particle_Surface', ),
                              NOM_CMP=('DZ', ),
                              VALE=(-1.765, )),
                      MODELE=model1,
                      OPERATION='AFFE',
                      TYPE_CHAM='NOEU_DEPL_R')

solver = DYNA_NON_LINE(identifier='9:1',
                       CHAM_MATER=fieldma0,
                       CONTACT=contact,
                       ETAT_INIT=_F(VITE=init_vel),
                       INCREMENT=_F(LIST_INST=timeline,
                                    NUME_INST_FIN=10,
                                    NUME_INST_INIT=0),
                       MODELE=model1,
                       SCHEMA_TEMPS=_F(FORMULATION='DEPLACEMENT',
                                       MODI_EQUI='OUI',
                                       SCHEMA='HHT'))
FIN()

However, during execution of the study, the following error happened:

.. __stg1_cmd9:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0012 Concept of the type:  evol_noli
  # ------------------------------------------------------------------------------------------
  solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F(LIST_INST=timeline,
                                      NUME_INST_INIT=0,
                                      NUME_INST_FIN=10,
                                      PRECISION=1.E-06,),
                         SCHEMA_TEMPS=_F(SCHEMA='HHT',
                                         MODI_EQUI='OUI',
                                         ALPHA=-0.3,
                                         FORMULATION='DEPLACEMENT',
                                         COEF_MASS_SHIFT=0.0,),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F(REAC_INCR=1,
                                   MATRICE='TANGENTE',
                                   REAC_ITER=1,
                                   REAC_ITER_ELAS=0,
                                   MATR_RIGI_SYME='NON',),
                         CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                        ITER_GLOB_ELAS=25,
                                        ARRET='OUI',),
                         SOLVEUR=_F(RENUM='AUTO',
                                    NPREC=8,
                                    ELIM_LAGR='LAGR2',
                                    STOP_SINGULIER='OUI',
                                    TYPE_RESOL='AUTO',
                                    ACCELERATION='AUTO',
                                    LOW_RANK_SEUIL=0.0,
                                    PRETRAITEMENTS='AUTO',
                                    POSTTRAITEMENTS='AUTO',
                                    PCENT_PIVOT=20,
                                    RESI_RELA=-1.0,
                                    GESTION_MEMOIRE='AUTO',
                                    FILTRAGE_MATRICE=-1.0,
                                    MIXER_PRECISION='NON',
                                    MATR_DISTRIBUEE='NON',
                                    METHODE='MUMPS',),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F(PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INFO=1,)


Comme vous n'avez pas défini explicitement le comportement, tout le modèle est supposé élastique en petites perturbations.
   Liste des comportements
   Affecté sur 365355 éléments
     Relation                             : ELAS
     Déformation                          : PETIT
     Nombre total de variables internes   : 1
            V1 : VIDE
  Le système linéaire à résoudre a 1658619 degrés de liberté:
   - 1558995 sont des degrés de liberté physiques
     (ils sont
portés par 519665 noeuds du maillage)
   - 99624 sont les couples de paramètres de Lagrange associés
     aux 49812 relations
linéaires dualisées.
  La matrice est de taille 1658619 équations.
  Elle contient 67789758 termes non nuls si elle est symétrique et 133920897 termes
non nuls si elle n'est pas symétrique.
  Soit un taux de remplissage de   0.005 %.
    Lecture de l'état initial
      Il n'y a pas d'état initial défini. On prend un état initial nul.
      Le champ <DEPL> est initialisé a zéro
      Le champ <SIEF_ELGA> est initialisé a zéro
      Le champ <VARI_ELGA> est initialisé a zéro
   
   !-----------------------------------------------------------------!
   ! <A> <ALGELINE7_20>                                                 !
   !                                                                                         !
   !  On ne peut pas remplir la composante DX du noeud numéro N1     !
   !  pour le champ "&&NMETL2.CHAMP.CONV"                            !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------------!
   ! <A> <MECANONLINE_2>                                                                    !
   !                                                                                        !
   !  Lors de la recopie du champ &&NMETL2.CHAMP.CONVER donné dans ETAT_INIT de la commande !
   !  STAT_NON_LINE vers le champ                                                           !
   ! &&NMCH1P.VITMOI, certaines composantes de &&NMCH1P.VITMOI                              !
   !  ont du être mises à zéro.                                                             !
   !                                                                                        !
   !  Ce problème survient lorsque le champ donné                                           !
   ! dans ETAT_INIT ne comporte                                                             !
   !  pas assez de composantes, on complète donc par des zéros.                             !
   !                                                                                        !
   !                                                                                        !
   ! This is a warning. If you do not understand the meaning of this                        !
   !  warning, you can obtain unexpected results!                                           !
   !----------------------------------------------------------------------------------------!
   
      Le champ <VITE> est lu dans ETAT_INIT, par un champ donné explicitement
      Le champ <ACCE> est initialisé a zéro

One considers an acceleration initial.

The initial state does not have acceleration given.
It is computed.

   
   !---------------------------------------------------------------!
   ! <EXCEPTION> <FACTOR_78>                                       !
   !                                                               !
   ! Solver MUMPS:                                                 !
   !  Overflow of the terms of the matrix. One detected at least:  !
   !  - 593685 too small terms,                                    !
   !  - 5 too large terms.                                         !
   !---------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
   !                                                                 !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
   !                                                                 !
   !    If an ABNORMAL_ABORT occured, that is not reproductible      !
   !    on official Linux version, or if you need some help to       !
   !    understand the message above, please report issues at:       !
   !                                                                 !
   !                      support@simulease.com                      !
   !                                                                 !
   !-----------------------------------------------------------------!


Criterion of destruction of the file ( 1.00%) associated with the base VOLATILE exceeded 7.89%
Nombre of records used:  4965

Occupied disc volume:  3879 Mo.
Maximum number of records:  62914


  Ouverture en écriture du fichier ./vola.1

Validation of the concept  solver.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR


Could you please explain the problem, and show me a solution to it?
Thanks!

#24 Re: Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <CONTACT2_11> (DYNA_NON_LINE) » 2020-06-08 02:40:56

Hello, could some one please help me out?

The links of the case files are attached FYI in case it helps the 'dissection'.
I am using mega.nz to host them since they exceed the size limit of attachment.

Mesh file: https: // mega.nz / file / HcoXHIxD#tp3AeMwVUqwFjUy2hEwiFzFT2v4aAWEU92Wb7nuR334  (38.7MB)
RunCase result: https: // mega.nz / file / 3U5HGSrL#ptV3rdkg6oiUyeMDuBosT1OIVH3Bz0EFBxu5EoLU69M  (47.7MB)

I am running code aster version 14.4.0 on Windows-10-10.0.17134-SP0.

Many thanks.

#25 Code_Aster usage » [Solved]Advice needed on <EXCEPTION> <CONTACT2_11> (DYNA_NON_LINE) » 2020-06-05 08:09:50

hkboondoggle
Replies: 2

Hi all,

I am new to code_aster, and I am trying to simulate an impact scenario with a ball travelling on negative z-axis direction and being about to collide on a static plate.

I have set up the contact rule as follows:
49972752751_0e2bdfa8fa_o.jpg
49972752661_c399b114f6_o.jpg
49973012952_c2a7ec17b3_o.jpg
49973013002_67fe9fbaf2_o.jpg
49972234013_bc2cdd3f38_o.jpg


However, during execution of the study, the following error happened:

.. __stg1_cmd7:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0010 Concept of the type:  char_contact
  # ------------------------------------------------------------------------------------------
  contact = DEFI_CONTACT(MODELE=model1,
                         STOP_INTERP='NON',
                         LISSAGE='NON',
                         VERI_NORM='OUI',
                         ITER_GEOM_MAXI=10,
                         RESI_GEOM=1.E-2,
                         ALGO_RESO_GEOM='POINT_FIXE',
                         REAC_GEOM='AUTOMATIQUE',
                         ITER_CONT_MULT=4,
                         STOP_SINGULIER='OUI',
                         NB_RESOL=10,
                         ITER_GCP_MAXI=0,
                         RECH_LINEAIRE='ADMISSIBLE',
                         PRE_COND='SANS',
                         ZONE=_F(APPARIEMENT='MAIT_ESCL',
                                 GROUP_MA_MAIT='Substrate',
                                 GROUP_MA_ESCL='Particle',
                                 NORMALE='MAIT_ESCL',
                                 VECT_MAIT='AUTO',
                                 VECT_ESCL='AUTO',
                                 GLISSIERE='NON',
                                 TYPE_PROJECTION='ORTHOGONALE',
                                 DIST_POUTRE='NON',
                                 DIST_COQUE='NON',
                                 DIST_APPA=-1.0,
                                 TOLE_PROJ_EXT=0.5,
                                 RESOLUTION='OUI',
                                 ALGO_CONT='CONTRAINTE',),
                         FORMULATION='DISCRETE',
                         FROTTEMENT='SANS',)

T
T
   
   !------------------------------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <CONTACT2_11>                                                                                                                                          !
   !                                                                                                                                                                                                 !
   ! At least one of the meshes of contact which you defined is of dimension 3, but the dimension of your problem is:  3.  !
   ! This mesh is                                                                                                                                                                            !
   ! thus not an edge mesh. There must be an error in your setting in data.                                                                            !
   !                                                                                                                                                                                                 !
   ! Advice:                                                                                                                                                                                    !
   ! Check your AFFE_MODELE and the type of your                                                                                                                   !
   ! meshes in the definition of contact surfaces.                                                                                                                       !
   !------------------------------------------------------------------------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------------------!
   !                                                                                                            !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
   !                                                                                                            !
   !    If an ABNORMAL_ABORT occured, that is not reproductible       !
   !    on official Linux version, or if you need some help to                  !
   !    understand the message above, please report issues at:            !
   !                                                                                                            !
   !                      support@simulease.com                                              !
   !                                                                                                            !
   !-----------------------------------------------------------------------------!

Destruction of the concept  contact.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR


Could you please explain the problem, and show me a solution to it?
Thanks!