testa_pbar3.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_pbar3.dat


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

TITLE - Four bars in series; applied loads and displ's. Same as testa but using equivalent bars.
NODES - # of nodes: 9; # of dimensions: 1
  #      ID                x                 y              z
  1      1                 1
  2      2                 2
  3      3                 3
  4      4                 4
  5      5                 5
  6      6               1.5
  7      7               2.5
  8      8               3.5
  9      9               4.5
ELEMENTS - # of elements: 4
  #     ID   elemType   propID   connectivity
  9     1    pBar3        1        1        6        2
  9     2    pBar3        2        2        7        3
  9     3    pBar3        3        3        8        4
  9     4    pBar3        4        4        9        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 =

   58.3333  -66.6667    8.3333
  -66.6667  133.3333  -66.6667
    8.3333  -66.6667   58.3333


fe =

     0
     0
     0


ke =

   46.6667  -53.3333    6.6667
  -53.3333  106.6667  -53.3333
    6.6667  -53.3333   46.6667


fe =

     0
     0
     0


ke =

   93.3333 -106.6667   13.3333
 -106.6667  213.3333 -106.6667
   13.3333 -106.6667   93.3333


fe =

     0
     0
     0


ke =

   23.3333  -26.6667    3.3333
  -26.6667   53.3333  -26.6667
    3.3333  -26.6667   23.3333


fe =

     0
     0
     0



Stiffness Matrix and Load Vector after assembly
kglobal =

  Columns 1 through 7 

   58.3333    8.3333         0         0         0  -66.6667         0
    8.3333  105.0000    6.6667         0         0  -66.6667  -53.3333
         0    6.6667  140.0000   13.3333         0         0  -53.3333
         0         0   13.3333  116.6667    3.3333         0         0
         0         0         0    3.3333   23.3333         0         0
  -66.6667  -66.6667         0         0         0  133.3333         0
         0  -53.3333  -53.3333         0         0         0  106.6667
         0         0 -106.6667 -106.6667         0         0         0
         0         0         0  -26.6667  -26.6667         0         0

  Columns 8 through 9 

         0         0
         0         0
 -106.6667         0
 -106.6667  -26.6667
         0  -26.6667
         0         0
         0         0
  213.3333         0
         0   53.3333


load_vector =

     0
     0
     0
     0
     0
     0
     0
     0
     0


Load Vector after assembly of point loads
load_vector =

     2
     0
    -3
    -1
     0
     0
     0
     0
     0


Constrained Stiffness Matrix and Load Vector

 Stiffness and force terms only printed for small problems

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

SOLUTION VECTOR
 dof#      value
   1           0.28
   2            0.2
   3      0.0071429
   4      -0.014286
   5              0
   6           0.24
   7        0.10357
   8     -0.0035714
   9     -0.0071429

Strain Energy (assume homogeneous displacement BCs) = 4.290714
ELEMENTAL POSTPROCESSING
  ID         f_i                 Sxx(end nodes)          energy
   1        2        0       -2          -2      -2       0.08

stress_nd =

   -2.0000   -2.0000   -2.0000

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

stress_nd =

   -3.8571   -3.8571   -3.8571

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

stress_nd =

   -0.8571   -0.8571   -0.8571

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

stress_nd =

    0.1429    0.1429    0.1429


Strain Energy = 0.462143