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 [2] subcell limits in X-direction Loading [2] subcell limits in Y-direction Loading [2] 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. 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. 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 (v 5.3.2) STARTED AT: 2007/05/16 12:33:37 MAIN: MAIN: ========================================== MAIN: E L M E R S O L V E R S T A R T I N G MAIN: Library version: 5.3.2 MAIN: ========================================== MAIN: MAIN: MAIN: ----------------------- MAIN: Reading Model ... MAIN: Done MAIN: ----------------------- OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: stress analysis...done. OptimizeBandwidth: Half bandwidth without optimization: 44 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: StressSolve: Using six pure forces and moments to compute the spring matrix StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 1 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 1.10449894671367 StressSolve: Relative Change : 2.00000000000000 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 2 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 1.10449894670985 StressSolve: Relative Change : 3.458430947261155E-012 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 3 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 3.335036805817577E-003 StressSolve: Relative Change : 1.98795835168286 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 4 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 0.152631638307863 StressSolve: Relative Change : 1.91446796430361 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 5 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 0.152631638308373 StressSolve: Relative Change : 3.337434393813648E-012 StressSolve: StressSolve: StressSolve: ------------------------------------- StressSolve: DISPLACEMENT SOLVER ITERATION 6 StressSolve: ------------------------------------- StressSolve: StressSolve: Starting assembly... StressSolve: Assembly done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StressSolve: Set boundaries done StressSolve: Result Norm : 2.889332771941279E-002 StressSolve: Relative Change : 1.36332002475097 StressSolve: ----------------------------------------- StressSolve: Saving lumped elastic spring to file kmat.dat StressSolve: ----------------------------------------- SolveEquations: (NRM,RELC): ( 0.288933277194E-01 2.00000000000 ) :: stress analysis : *** Elmer Solver: ALL DONE *** SOLVER TOTAL TIME(CPU,REAL): 2.15 2.27 ELMER SOLVER FINISHED AT: 2007/05/16 12:33:40