testa.out

matlab实例

                              < M A T L A B >
                  Copyright 1984-2002 The MathWorks, Inc.
                      Version 6.5.0.180913a Release 13
                                Jun 18 2002

  Using Toolbox Path Cache.  Type "help toolbox_path_cache" for more info.
 
  To get started, type one of these: helpwin, helpdesk, or demo.
  For product information, visit www.mathworks.com.
 
*********************************************************************
* FEMjr v 2.0 - Finite Element Method junior - An FEM framework     *
*********************************************************************
Enter name of input data file: 
inputfile =

testa.dat


<<<<---- Data File Echo ---->>>>

TITLE - Four springs in series; applied loads and displ's  (Allaire pg 85)
NODES - # of nodes: 5; # of dimensions: 1
  #      ID                x                 y              z
  1      1                 1
  2      2                 2
  3      3                 3
  4      4                 4
  5      5                 5
ELEMENTS - # of elements: 4
  #     ID   elemType   propID   connectivity
  5     1    spring1d        1        1        2
  5     2    spring1d        2        2        3
  5     3    spring1d        3        3        4
  5     4    spring1d        4        4        5
PROPERTIES - # of properties: 4
  #   ID     propType         properties
  1    1    MatSpring            25
  2    2    MatSpring            20
  3    3    MatSpring            40
  4    4    MatSpring            10
CONSTRAINTS - # of constraints: 2
  #       nodeID    dof     value
  1         2         1            0.2
  2         5         1              0
LOADS - # of loads: 3
  #       nodeID    dof     value
  1         1         1              2
  2         3         1             -3
  3         4         1             -1


Stiffness Matrix and Load Vector after assembly
kglobal =

    25   -25     0     0     0
   -25    45   -20     0     0
     0   -20    60   -40     0
     0     0   -40    50   -10
     0     0     0   -10    10


load_vector =

     0
     0
     0
     0
     0


Load Vector after assembly of point loads
load_vector =

     2
     0
    -3
    -1
     0


Constrained Stiffness Matrix and Load Vector
STIFFNESS MATRIX
r/c        1        2        3        4        5        6        7        8
  1       25        0        0        0        0
  2        0       45        0        0        0
  3        0        0       60      -40        0
  4        0        0      -40       50        0
  5        0        0        0        0       10

LOAD VECTOR
Row     Value
  1              7
  2              9
  3              1
  4             -1
  5              0

<<<<---- Solution Results ---->>>>

SOLUTION VECTOR
 dof#      value
   1           0.28
   2            0.2
   3      0.0071429
   4      -0.014286
   5              0

Strain Energy (assume homogeneous displacement BCs) = 1.890714
ELEMENTAL POSTPROCESSING
   ID      f1           f2     energy
    1          2         -2        0.08
   ID      f1           f2     energy
    2       3.86      -3.86       0.372
   ID      f1           f2     energy
    3      0.857     -0.857     0.00918
   ID      f1           f2     energy
    4     -0.143      0.143     0.00102

Strain Energy = 0.462143