ElmerGrid 1 2 beam3d Starting program Elmergrid Elmergrid reading in-line arguments Output will be saved to file beam3d. Elmergrid loading data: Loading the geometry from file 'beam3d.grd'. Loading ElmerGrid file version: 210903 Defining the coordinate system (3-DIM). Loading [1] subcell limits in X-direction Loading [1] subcell limits in Y-direction Loading [1] subcell limits in Z-direction Loading material structure End of field Loading boundary conditions Found 4 boundaries Loading minimum number of elements Reached the end of command file Found 1 divisions for grid SetElementDivision The mesh is devided into 1 separate subcells. Created 5 radial and 5 axial divisions, total of 25 elements. Loading ElmerGrid commands from file 'beam3d.grd'. Reached the end of command file Read commands from a file Elmergrid creating and manipulating meshes: 1 cells were created. Numbered 36 knots in 25 4-node elements. Numbering order was and max levelwidth 8. Maximum elementsize is 2.000e-01 and minimum 2.000e-01. 5 element sides between materials -1 and 1 were located to type 1. 5 element sides between materials -2 and 1 were located to type 2. 5 element sides between materials -3 and 1 were located to type 3. 5 element sides between materials -4 and 1 were located to type 4. Created 41 extruded divisions. Elementtype 404 extruded to type 808. Extruded mesh has 1025 elements in 41 levels. Extruded BC 0 of type 1 was created with 205 sides. Extruded BC 1 of type 2 was created with 205 sides. Extruded BC 2 of type 3 was created with 205 sides. Extruded BC 3 of type 4 was created with 205 sides. BC 4 on layer 1 was created with 25 sides. BC 5 on layer 42 was created with 25 sides. Reordering 1512 knots and 1025 elements in 3-dimensions. Ordering with (1*x + 1*y + 1e+04*z) Moving knots to new positions Moving the elements to new positions Moving the parents of the boundary nodes. Indexwidth of the new node order is 51. Created 1025 elements and 1512 nodes by extruding the 2D geometry Elmergrid saving data: Saving mesh in ElmerSolver format to directory beam3d. Saving 1512 coordinates to mesh.nodes. Saving 1025 element topologies to mesh.elements. Saving boundary elements to mesh.boundary. Saving header info to mesh.header. Thank you for using Elmergrid! Send bug reports and feature wishes to peter.raback@csc.fi ../../src/ElmerSolver ELMER SOLVER STARTED AT: 2005/10/07 15:37:58 MAIN: MAIN: =========================================================== MAIN: E L M E R S O L V E R S T A R T I N G, W E L C O M E MAIN: =========================================================== MAIN: MAIN: MAIN: ----------------------- MAIN: Reading Model ... WARNING:: Model Input: Unknown keyword: [model lumping] in section: [solver 1] WARNING:: Model Input: Unknown keyword: [model lumping filename] in section: [solver 1] WARNING:: Model Input: Unknown keyword: [model lumping boundary] in section: [boundary condition 2] WARNING:: Model Input: Unknown keyword: [model lumping] in section: [solver 1] WARNING:: Model Input: Unknown keyword: [model lumping filename] in section: [solver 1] WARNING:: Model Input: Unknown keyword: [model lumping boundary] in section: [boundary condition 2] MAIN: Done MAIN: ----------------------- OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: stress analysis...done. OptimizeBandwidth: Half bandwidth without optimization: 52 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 92 OptimizeBandwidth: Bandwidth optimization rejected, using original ordering. OptimizeBandwidth: --------------------------------------------------------- MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: WARNING:: ListGetLogical: WARNING:: ListGetLogical: Requested property: [Fix Displacement], not found WARNING:: ListGetLogical: StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 1 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 27 % done StressSolve: Assembly: 54 % done StressSolve: Assembly: 82 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 4.06 20 0.1530E+00 40 0.9975E-02 53 0.3516E-08 StressSolve: Result Norm : 1.10449894659951 StressSolve: Relative Change : 2.00000000000000 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 2 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 25 % done StressSolve: Assembly: 51 % done StressSolve: Assembly: 65 % done StressSolve: Assembly: 82 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 3.46 20 0.3976E-01 40 0.6440E-02 48 0.3471E-08 StressSolve: Result Norm : 1.10449894246097 StressSolve: Relative Change : 3.74698367217499E-9 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 3 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 27 % done StressSolve: Assembly: 54 % done StressSolve: Assembly: 82 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 3.42 20 0.1153E+00 40 0.7974E-02 58 0.5973E-08 StressSolve: Result Norm : 0.00333503679524604 StressSolve: Relative Change : 1.98795835167473 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 4 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 26 % done StressSolve: Assembly: 54 % done StressSolve: Assembly: 81 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 3.43 20 0.2361E-01 40 0.3364E-03 54 0.6645E-08 StressSolve: Result Norm : 0.152631647005913 StressSolve: Relative Change : 1.91446796933894 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 5 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 27 % done StressSolve: Assembly: 54 % done StressSolve: Assembly: 82 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 3.43 20 0.9706E-02 40 0.2126E-03 50 0.5305E-08 StressSolve: Result Norm : 0.152631615747434 StressSolve: Relative Change : 2.04796862665448E-7 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 6 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly: 27 % done StressSolve: Assembly: 54 % done StressSolve: Assembly: 81 % done StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Diriclet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done CRS_IncompleteLU: ILU(2) (Real), Starting Factorization: CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 992467 CRS_IncompleteLU: ILU(2) (Real), filling (%) : 347 CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 3.44 20 0.4456E+00 40 0.3791E-03 54 0.9000E-08 StressSolve: Result Norm : 0.0288933278277624 StressSolve: Relative Change : 1.36331994361330 StressSolve: ----------------------------------------- StressSolve: Saving lumped spring matrix to file kmat.dat StressSolve: ----------------------------------------- SolveEquations: (NRM,RELC): ( 0.288933278278E-01 2.00000000000 ) :: stress analysis : *** Elmer Solver: ALL DONE *** Remaining memory: 16 bytes allocated at line 1036 of GeneralUtils.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 362 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 Remaining memory: 156 bytes allocated at line 367 of ElementDescription.f90 ... More segments remain SOLVER TOTAL TIME(CPU,REAL): 98.37 105.16 ELMER SOLVER FINISHED AT: 2005/10/07 15:39:43