testc.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 =

testc.dat


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

TITLE - Four bar truss (Chandrupatla pg 106)
NODES - # of nodes: 4; # of dimensions: 2
  #      ID                x                 y              z
  1      1                 0             0
  2      2                40             0
  3      3                40            30
  4      4                 0            30
ELEMENTS - # of elements: 4
  #     ID   elemType   propID   connectivity
  4     1    truss2d        1        1        2
  4     2    truss2d        1        2        3
  4     3    truss2d        1        1        3
  4     4    truss2d        1        3        4
PROPERTIES - # of properties: 1
  #   ID     propType         properties
  1    1    MatTruss      2.95e+07             1
CONSTRAINTS - # of constraints: 5
  #       nodeID    dof     value
  1         1         1              0
  2         1         2              0
  3         2         2              0
  4         4         1              0
  5         4         2              0
LOADS - # of loads: 2
  #       nodeID    dof     value
  1         2         1          20000
  2         3         2         -25000


Stiffness Matrix and Load Vector after assembly
kglobal =

   1.0e+06 *

  Columns 1 through 7 

    1.1151    0.2832   -0.7375         0   -0.3776   -0.2832         0
    0.2832    0.2124         0         0   -0.2832   -0.2124         0
   -0.7375         0    0.7375         0         0         0         0
         0         0         0    0.9833         0   -0.9833         0
   -0.3776   -0.2832         0         0    1.1151    0.2832   -0.7375
   -0.2832   -0.2124         0   -0.9833    0.2832    1.1957         0
         0         0         0         0   -0.7375         0    0.7375
         0         0         0         0         0         0         0

  Column 8 

         0
         0
         0
         0
         0
         0
         0
         0


load_vector =

     0
     0
     0
     0
     0
     0
     0
     0


Load Vector after assembly of point loads
load_vector =

           0
           0
       20000
           0
           0
      -25000
           0
           0


Constrained Stiffness Matrix and Load Vector
STIFFNESS MATRIX
r/c        1        2        3        4        5        6        7        8
  1  1.1e+06        0        0        0        0        0        0        0
  2        0  2.1e+05        0        0        0        0        0        0
  3        0        0  7.4e+05        0        0        0        0        0
  4        0        0        0  9.8e+05        0        0        0        0
  5        0        0        0        0  1.1e+06  2.8e+05        0        0
  6        0        0        0        0  2.8e+05  1.2e+06        0        0
  7        0        0        0        0        0        0  7.4e+05        0
  8        0        0        0        0        0        0        0        1

LOAD VECTOR
Row     Value
  1              0
  2              0
  3          20000
  4              0
  5              0
  6         -25000
  7              0
  8              0

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

SOLUTION VECTOR
 dof#      value
   1              0
   2              0
   3       0.027119
   4              0
   5      0.0056497
   6      -0.022246
   7              0
   8              0

Strain Energy (assume homogeneous displacement BCs) = 549.258475
ELEMENTAL POSTPROCESSING
  ID     f1_x          f1_y          f2_x          f2_y      magnitude      stress   strain energy
   1       -20000            0        20000            0        20000        20000       271.19
  ID     f1_x          f1_y          f2_x          f2_y      magnitude      stress   strain energy
   2            0        21875            0       -21875        21875       -21875       243.31
  ID     f1_x          f1_y          f2_x          f2_y      magnitude      stress   strain energy
   3       4166.7         3125      -4166.7        -3125       5208.3      -5208.3       22.989
  ID     f1_x          f1_y          f2_x          f2_y      magnitude      stress   strain energy
   4       4166.7            0      -4166.7            0       4166.7       4166.7        11.77

Strain Energy = 549.258475