testa_pbar2.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_pbar2.dat


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

TITLE - Four bars in series; applied loads and displ's. Same as testa but using equivalent bars.
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    pBar2        1        1        2
  5     2    pBar2        2        2        3
  5     3    pBar2        3        3        4
  5     4    pBar2        4        4        5
PROPERTIES - # of properties: 4
  #   ID     propType         properties
  1    1    MatpBar            25             1             0             0
  2    2    MatpBar            20             1             0             0
  3    3    MatpBar            40             1             0             0
  4    4    MatpBar            10             1             0             0
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

ke =

    25   -25
   -25    25


fe =

     0
     0


ke =

    20   -20
   -20    20


fe =

     0
     0


ke =

    40   -40
   -40    40


fe =

     0
     0


ke =

    10   -10
   -10    10


fe =

     0
     0



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         f_i                 Sxx(end nodes)          energy
   1        2       -2          -2      -2       0.08

stress_nd =

   -2.0000   -2.0000

  ID         f_i                 Sxx(end nodes)          energy
   2     3.86    -3.86       -3.86   -3.86      0.372

stress_nd =

   -3.8571   -3.8571

  ID         f_i                 Sxx(end nodes)          energy
   3    0.857   -0.857      -0.857  -0.857    0.00918

stress_nd =

   -0.8571   -0.8571

  ID         f_i                 Sxx(end nodes)          energy
   4   -0.143    0.143       0.143   0.143    0.00102

stress_nd =

    0.1429    0.1429


Strain Energy = 0.462143