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.
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
block 2 volume all
set duplicate block elements on
before you define second block.
This will allow for generation of multiple blocks in cubit.
Ignatios