Cubit Model error

54 views
Skip to first unread message
Assigned to ignat...@gmail.com by lik...@wp.pl

Calum Macdonald

unread,
Dec 13, 2021, 11:15:19 AM12/13/21
to MoFEM Q&A

Journal File

reset

#volume 1

brick x 10 y 0.5 z 10


#volume 2

brick x 0.5 y 5 z 0.5

move volume 2 x -4 y -2.75 z -4 include_merged

merge volume 1 2

#volume 3

brick x 0.5 y 5 z 0.5

move volume 3 x 4 y -2.75 z -4 include_merged

merge volume 1 3

#volume 4

brick x 0.5 y 5 z 0.5

move volume 4 x 4 y -2.75 z 4 include_merged

merge volume 1 4

#volume 5

brick x 0.5 y 5 z 0.5

move volume 5 x -4 y -2.75 z 4 include_merged

merge volume 1 5


#volume 6

brick x 1 y 0.5 z 1

move volume 6 x -4 y -5.5 z -4 include_merged

merge volume 2 6

#volume 7

brick x 1 y 0.5 z 1

move volume 7 x 4 y -5.5 z -4 include_merged

merge volume 3 7

#volume 8

brick x 1 y 0.5 z 1

move volume 8 x 4 y -5.5 z 4 include_merged

merge volume 4 8

#volume 9

brick x 1 y 0.5 z 1

move volume 9 x -4 y -5.5 z 4 include_merged

merge volume 5 9


#volume 10

brick x 0.5 y 5 z 0.5

move volume 10 x -4 y 2.75 z -4 include_merged

merge volume 1 10

#volume 11

brick x 0.5 y 5 z 0.5

move volume 11 x 4 y 2.75 z -4 include_merged

merge volume 1 11

#volume 12

brick x 0.5 y 5 z 0.5

move volume 12 x 4 y 2.75 z 4 include_merged

merge volume 1 12

#volume 13

brick x 0.5 y 5 z 0.5

move volume 13 x -4 y 2.75 z 4 include_merged

merge volume 1 13


#volume 14

brick x 10 y 0.5 z 10

move volume 14 y 5.25 include_merged

merge volume 14 10

merge volume 14 11

merge volume 14 12

merge volume 14 13

set acisoption integer 'check_level' 70

set overlap maximum gap 0.0005

imprint volume all

merge volume all

set overlap maximum gap 0.0005


undo group begin

volume all scheme tetmesh

mesh volume all

undo group end


block 1 volume all

block 1 name 'MAT_ELASTIC'

block 1 attribute count 2

block 1 attribute index 1 50e6 #young modulus

block 1 attribute index 2 0.1 #poisson ratio


block 2 volume 1

block 2 name "BODY_FORCES"

block 2 attribute count 4

block 2 attribute index 1 2500 #material density

block 2 attribute index 2 0. #constant acceleration in x-direction

block 2 attribute index 3 0. #constant acceleration in y-direction

block 2 attribute index 4 0. #constant acceleration in z-direction


create displacement on volume 7 dof 1 dof 2 dof 3 fix 0

create displacement on volume 8 dof 1 dof 2 dof 3 fix 0

create displacement on volume 9 dof 1 dof 2 dof 3 fix 0

create displacement on volume 6 dof 1 dof 2 dof 3 fix 0


!rm -rf out_*
!{program} \
-my_file Frame.cub \
-ts_dt 0.1 -ts_max_time 4 \
-my_output_prt -10 \
-my_max_post_proc_ref_level 0 \
-my_disp_order 2 \
-my_time_data_file Frame_history.in \
-default_material HOOKE \
-is_linear \
-ts_type alpha \
-ts_alpha_radius 1 2>&1 | tee log


Error

[0] <inform> MoFEM version 0.12.1 (MOAB 5.2.1 Petsc Release Version 3.11.4, Sep, 28, 2019 ) [0] <inform> git commit id 2c75b451207e2bc357f65bbed286fb6664fc69b6 [0] <inform> Local time: 2021-12-13 14:26:51 [0] <inform> UTC time: 2021-12-13 14:26:51 Warning: acis file has sequence numbers! [0] <inform> [MeshsetMng] meshset 12682136550675317183 type BLOCKSET MAT_ELASTICSET msId 1 name MAT_ELASTIC block header: blockCol = 4294967295 blockMat = 0 blockDimension = 3 [0] <inform> [MeshsetMng] meshset 12682136550675317184 type NODESET DISPLACEMENTSET msId 1 [0] <inform> [MeshsetMng] meshset 12682136550675317185 type NODESET DISPLACEMENTSET msId 2 [0] <inform> [MeshsetMng] meshset 12682136550675317186 type NODESET DISPLACEMENTSET msId 3 [0] <inform> [MeshsetMng] meshset 12682136550675317187 type NODESET DISPLACEMENTSET msId 4 [0] <inform> [FieldCore] Add field MESH_NODE_POSITIONS field_id 1 space H1 approximation base AINSWORTH_LEGENDRE_BASE rank 3 meshset 12682136550675317189 [0] <inform> [FieldCore] Add field DISPLACEMENT field_id 2 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317190 [0] <inform> [FECore] Add finite element FORCE_FE [0] <inform> [FECore] Add finite element PRESSURE_FE [0] <inform> [FECore] Add finite element FLUID_PRESSURE_FE [0] <inform> [FieldCore] Add field VELOCITY field_id 4 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317194 [0] <inform> [FieldCore] Add field DOT_DISPLACEMENT field_id 8 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317195 [0] <inform> [FieldCore] Add field DOT_VELOCITY field_id 16 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317196 [0] <inform> [petsc] Block Id 1 Young Modulus 5e+07 Poisson Ration 0.10 Material model HOOKE Nb. of elements 24934 [0] <inform> [FECore] Add finite element ELASTIC [0] <inform> [FECore] Add finite element MASS_ELEMENT [0] <inform> [FECore] Add finite element VELOCITY_ELEMENT [0] <inform> [FECore] Add finite element DAMPER [0] <inform> [FieldCore] Number of dofs 654405 [0] <inform> [FECore] Finite element FORCE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element PRESSURE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element FLUID_PRESSURE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element ELASTIC added. Nb. of elements added 24934 [0] <inform> [FECore] Finite element MASS_ELEMENT added. Nb. of elements added 0 [0] <inform> [FECore] Finite element VELOCITY_ELEMENT added. Nb. of elements added 0 [0] <inform> [FECore] Finite element DAMPER added. Nb. of elements added 0 [0] <inform> [FECore] Number of adjacencies 1246700 [0] <inform> [ProblemCore] Add problem Kuu [0] <inform> [ProblemsManager] FEs ghost dofs on problem Kuu Nb. ghost dof 0 by 0 Nb. local dof 130881 by 130881 [0] <inform> [ProblemCore] Add problem DYNAMICS [0] <inform> [ProblemsManager] FEs ghost dofs on problem DYNAMICS Nb. ghost dof 0 by 0 Nb. local dof 130881 by 130881 [0] <inform> [dynamic] 0 Time 0.00e+00 Elastic energy 0.00e+00 Kinetic Energy 0.00e+00 Total 0.00e+00 [0] <inform> [petsc] 0 TS dt 0.1 time 0. [0] <error> --------------------- MoFEM Error Message --------------------------------------------------- [0] <error> MoFEM version 0.12.1 (MOAB 5.2.1 Petsc Release Version 3.11.4, Sep, 28, 2019 ) [0] <error> MoFEM git commit id 2c75b451207e2bc357f65bbed286fb6664fc69b6 [0] <error> See http://mofem.eng.gla.ac.uk/mofem/html/ guidelines_bug_reporting.html for bug reporting. [0] <error> Write to https://groups.google.com/forum/#!forum/mofem-group to seek help. [0] <error> problem field < VELOCITY > not found, (top tip: check spelling) [0] <error> #1 setOtherGlobalGhostVector() line 572 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/VecManager.cpp [0] <error> #2 problem_basic_method_preProcess() line 1740 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/ProblemsCore.cpp [0] <error> #3 problem_basic_method_preProcess() line 403 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/ProblemsCore.cpp [0] <error> #4 TsSetIFunction() line 79 in /mofem_install/mofem-cephas/mofem/src/petsc/impl/TsCtx.cpp [0] <error> #5 TSComputeIFunction() line 901 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #6 SNESTSFormFunction_Alpha() line 284 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #7 SNESTSFormFunction() line 4735 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #8 SNESComputeFunction() line 2406 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/interface/snes.c [0] <error> #9 SNESSolve_NEWTONLS() line 175 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/impls/ls/ls.c [0] <error> #10 SNESSolve() line 4560 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/interface/snes.c [0] <error> #11 TSAlpha_SNESSolve() line 85 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #12 TSAlpha_Restart() line 124 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #13 TSStep_Alpha() line 185 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #14 TSStep() line 3668 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #15 TSSolve() line 3851 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #16 main() line 907 in /mofem_install/mofem-cephas/mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp #PETSc Option Table entries: -default_material HOOKE -is_linear -ksp_atol 1e-10 -ksp_rtol 1e-10 -ksp_type fgmres -my_disp_order 2 -my_file Frame.cub -my_max_post_proc_ref_level 0 -my_output_prt -10 -my_time_data_file Frame_history.in -pc_factor_mat_solver_type mumps -pc_type lu -snes_atol 1e-7 -snes_linesearch_type basic -snes_max_it 100 -snes_monitor -snes_rtol 1e-7 -snes_type newtonls -ts_alpha_radius 1 -ts_dt 0.1 -ts_max_time 4 -ts_monitor -ts_type alpha #End of PETSc Option Table entries [0] <error> -- MoFEM End of Error Message -- send entire error message to mofem...@googlegroups.com -- -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 102. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them.

ignathanas

unread,
Dec 13, 2021, 11:28:14 AM12/13/21
to MoFEM Q&A
Hi Callum,

Thank you for your post. 
One definite error is that you create displacement BCs in the volume whereas they should be prescribed on surfaces.

Ignatios 

Calum Macdonald

unread,
Dec 13, 2021, 11:52:50 AM12/13/21
to MoFEM Q&A
I have changed this however it has given the same error message

Calum Macdonald

unread,
Dec 14, 2021, 12:10:41 PM12/14/21
to MoFEM Q&A

block 1 volume all

block 1 name 'MAT_ELASTIC'

block 1 attribute count 2

block 1 attribute index 1 50e6 #young modulus

block 1 attribute index 2 0.1 #poisson ratio


block 2 volume 1

block 2 name "BODY_FORCES"

block 2 attribute count 4

block 2 attribute index 1 2500 #material density

block 2 attribute index 2 0. #constant acceleration in x-direction

block 2 attribute index 3 0. #constant acceleration in y-direction

block 2 attribute index 4 0. #constant acceleration in z-direction


create displacement on surface 43 dof 1 dof 2 dof 3 fix 0

create displacement on surface 44 dof 1 dof 2 dof 3 fix 0

create displacement on surface 45 dof 1 dof 2 dof 3 fix 0

create displacement on surface 46 dof 1 dof 2 dof 3 fix 0

create displacement on surface 48 dof 1 dof 2 dof 3 fix 0

create displacement on surface 100 dof 1 dof 2 dof 3 fix 0


create displacement on surface 49 dof 1 dof 2 dof 3 fix 0

create displacement on surface 50 dof 1 dof 2 dof 3 fix 0

create displacement on surface 51 dof 1 dof 2 dof 3 fix 0

create displacement on surface 52 dof 1 dof 2 dof 3 fix 0

create displacement on surface 54 dof 1 dof 2 dof 3 fix 0

create displacement on surface 102 dof 1 dof 2 dof 3 fix 0


create displacement on surface 31 dof 1 dof 2 dof 3 fix 0

create displacement on surface 32 dof 1 dof 2 dof 3 fix 0

create displacement on surface 33 dof 1 dof 2 dof 3 fix 0

create displacement on surface 34 dof 1 dof 2 dof 3 fix 0

create displacement on surface 36 dof 1 dof 2 dof 3 fix 0

create displacement on surface 96 dof 1 dof 2 dof 3 fix 0


create displacement on surface 37 dof 1 dof 2 dof 3 fix 0

create displacement on surface 38 dof 1 dof 2 dof 3 fix 0

create displacement on surface 39 dof 1 dof 2 dof 3 fix 0

create displacement on surface 40 dof 1 dof 2 dof 3 fix 0

create displacement on surface 42 dof 1 dof 2 dof 3 fix 0

create displacement on surface 98 dof 1 dof 2 dof 3 fix 0


set acisoption integer 'check_level' 70

set overlap maximum gap 0.0005

imprint volume all

merge volume all

set overlap maximum gap 0.0005


undo group begin

volume all scheme tetmesh

mesh volume all

undo group end


ignathanas

unread,
Dec 14, 2021, 12:15:33 PM12/14/21
to MoFEM Q&A
Hi Callum,

I think BODY_FORCES block has to be set to all volumes rather than only volume 1.

Ignatios

Calum Macdonald

unread,
Dec 14, 2021, 12:27:01 PM12/14/21
to MoFEM Q&A

block 2 volume all


Error

[0] <inform> MoFEM version 0.12.1 (MOAB 5.2.1 Petsc Release Version 3.11.4, Sep, 28, 2019 ) [0] <inform> git commit id 2c75b451207e2bc357f65bbed286fb6664fc69b6 [0] <inform> Local time: 2021-12-14 17:17:40 [0] <inform> UTC time: 2021-12-14 17:17:40 Warning: acis file has sequence numbers! [0] <inform> [MeshsetMng] meshset 12682136550675317183 type BLOCKSET MAT_ELASTICSET msId 1 name MAT_ELASTIC block header: blockCol = 4294967295 blockMat = 0 blockDimension = 3 [0] <inform> [MeshsetMng] meshset 12682136550675317184 type NODESET DISPLACEMENTSET msId 1 [0] <inform> [MeshsetMng] meshset 12682136550675317185 type NODESET DISPLACEMENTSET msId 2 [0] <inform> [MeshsetMng] meshset 12682136550675317186 type NODESET DISPLACEMENTSET msId 3 [0] <inform> [MeshsetMng] meshset 12682136550675317187 type NODESET DISPLACEMENTSET msId 4 [0] <inform> [MeshsetMng] meshset 12682136550675317188 type NODESET DISPLACEMENTSET msId 5 [0] <inform> [MeshsetMng] meshset 12682136550675317189 type NODESET DISPLACEMENTSET msId 6 [0] <inform> [MeshsetMng] meshset 12682136550675317190 type NODESET DISPLACEMENTSET msId 7 [0] <inform> [MeshsetMng] meshset 12682136550675317191 type NODESET DISPLACEMENTSET msId 8 [0] <inform> [MeshsetMng] meshset 12682136550675317192 type NODESET DISPLACEMENTSET msId 9 [0] <inform> [MeshsetMng] meshset 12682136550675317193 type NODESET DISPLACEMENTSET msId 10 [0] <inform> [MeshsetMng] meshset 12682136550675317194 type NODESET DISPLACEMENTSET msId 11 [0] <inform> [MeshsetMng] meshset 12682136550675317195 type NODESET DISPLACEMENTSET msId 12 [0] <inform> [MeshsetMng] meshset 12682136550675317196 type NODESET DISPLACEMENTSET msId 13 [0] <inform> [MeshsetMng] meshset 12682136550675317197 type NODESET DISPLACEMENTSET msId 14 [0] <inform> [MeshsetMng] meshset 12682136550675317198 type NODESET DISPLACEMENTSET msId 15 [0] <inform> [MeshsetMng] meshset 12682136550675317199 type NODESET DISPLACEMENTSET msId 16 [0] <inform> [MeshsetMng] meshset 12682136550675317200 type NODESET DISPLACEMENTSET msId 17 [0] <inform> [MeshsetMng] meshset 12682136550675317201 type NODESET DISPLACEMENTSET msId 18 [0] <inform> [MeshsetMng] meshset 12682136550675317202 type NODESET DISPLACEMENTSET msId 19 [0] <inform> [MeshsetMng] meshset 12682136550675317203 type NODESET DISPLACEMENTSET msId 20 [0] <inform> [FieldCore] Add field MESH_NODE_POSITIONS field_id 1 space H1 approximation base AINSWORTH_LEGENDRE_BASE rank 3 meshset 12682136550675317205 [0] <inform> [FieldCore] Add field DISPLACEMENT field_id 2 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317206 [0] <inform> [FECore] Add finite element FORCE_FE [0] <inform> [FECore] Add finite element PRESSURE_FE [0] <inform> [FECore] Add finite element FLUID_PRESSURE_FE [0] <inform> [FieldCore] Add field VELOCITY field_id 4 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317210 [0] <inform> [FieldCore] Add field DOT_DISPLACEMENT field_id 8 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317211 [0] <inform> [FieldCore] Add field DOT_VELOCITY field_id 16 space H1 approximation base AINSWORTH_BERNSTEIN_BEZIER_BASE rank 3 meshset 12682136550675317212 [0] <inform> [petsc] Block Id 1 Young Modulus 5e+07 Poisson Ration 0.10 Material model HOOKE Nb. of elements 24934 [0] <inform> [FECore] Add finite element ELASTIC [0] <inform> [FECore] Add finite element MASS_ELEMENT [0] <inform> [FECore] Add finite element VELOCITY_ELEMENT [0] <inform> [FECore] Add finite element DAMPER [0] <inform> [FieldCore] Number of dofs 654405 [0] <inform> [FECore] Finite element FORCE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element PRESSURE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element FLUID_PRESSURE_FE added. Nb. of elements added 0 [0] <inform> [FECore] Finite element ELASTIC added. Nb. of elements added 24934 [0] <inform> [FECore] Finite element MASS_ELEMENT added. Nb. of elements added 0 [0] <inform> [FECore] Finite element VELOCITY_ELEMENT added. Nb. of elements added 0 [0] <inform> [FECore] Finite element DAMPER added. Nb. of elements added 0 [0] <inform> [FECore] Number of adjacencies 1246700 [0] <inform> [ProblemCore] Add problem Kuu [0] <inform> [ProblemsManager] FEs ghost dofs on problem Kuu Nb. ghost dof 0 by 0 Nb. local dof 130881 by 130881 [0] <inform> [ProblemCore] Add problem DYNAMICS [0] <inform> [ProblemsManager] FEs ghost dofs on problem DYNAMICS Nb. ghost dof 0 by 0 Nb. local dof 130881 by 130881 [0] <inform> [dynamic] 0 Time 0.00e+00 Elastic energy 0.00e+00 Kinetic Energy 0.00e+00 Total 0.00e+00 [0] <inform> [petsc] 0 TS dt 0.1 time 0. [0] <error> --------------------- MoFEM Error Message --------------------------------------------------- [0] <error> MoFEM version 0.12.1 (MOAB 5.2.1 Petsc Release Version 3.11.4, Sep, 28, 2019 ) [0] <error> MoFEM git commit id 2c75b451207e2bc357f65bbed286fb6664fc69b6 [0] <error> See http://mofem.eng.gla.ac.uk/mofem/html/ guidelines_bug_reporting.html for bug reporting. [0] <error> Write to https://groups.google.com/forum/#!forum/mofem-group to seek help. [0] <error> problem field < VELOCITY > not found, (top tip: check spelling) [0] <error> #1 setOtherGlobalGhostVector() line 572 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/VecManager.cpp [0] <error> #2 problem_basic_method_preProcess() line 1740 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/ProblemsCore.cpp [0] <error> #3 problem_basic_method_preProcess() line 403 in /mofem_install/mofem-cephas/mofem/src/interfaces/impl/ProblemsCore.cpp [0] <error> #4 TsSetIFunction() line 79 in /mofem_install/mofem-cephas/mofem/src/petsc/impl/TsCtx.cpp [0] <error> #5 TSComputeIFunction() line 901 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #6 SNESTSFormFunction_Alpha() line 284 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #7 SNESTSFormFunction() line 4735 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #8 SNESComputeFunction() line 2406 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/interface/snes.c [0] <error> #9 SNESSolve_NEWTONLS() line 175 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/impls/ls/ls.c [0] <error> #10 SNESSolve() line 4560 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/snes/interface/snes.c [0] <error> #11 TSAlpha_SNESSolve() line 85 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #12 TSAlpha_Restart() line 124 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #13 TSStep_Alpha() line 185 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/impls/implicit/alpha/alpha1.c [0] <error> #14 TSStep() line 3668 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #15 TSSolve() line 3851 in /tmp/root/spack-stage/spack-stage-petsc-3.11.4-vg3wvej7t3tymxaiywmjotoq52th6dqb/spack-src/src/ts/interface/ts.c [0] <error> #16 main() line 907 in /mofem_install/mofem-cephas/mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp #PETSc Option Table entries: -default_material HOOKE -is_linear -ksp_atol 1e-10 -ksp_rtol 1e-10 -ksp_type fgmres -my_disp_order 2 -my_file Frame.cub -my_max_post_proc_ref_level 0 -my_output_prt -10 -my_time_data_file Frame_history.in -pc_factor_mat_solver_type mumps -pc_type lu -snes_atol 1e-7 -snes_linesearch_type basic -snes_max_it 100 -snes_monitor -snes_rtol 1e-7 -snes_type newtonls -ts_alpha_radius 1 -ts_dt 0.1 -ts_max_time 4 -ts_monitor -ts_type alpha #End of PETSc Option Table entries [0] <error> -- MoFEM End of Error Message -- send entire error message to mofem...@googlegroups.com -- -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 102. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them.

ignathanas

unread,
Dec 14, 2021, 12:45:02 PM12/14/21
to MoFEM Q&A
Hi Callum,

Please add also the line:

set duplicate block elements on

before you define second block.

This will allow for generation of multiple blocks in cubit.


Ignatios


Reply all
Reply to author
Forward
0 new messages