📄 testa_pbar5.out
字号:
< 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_pbar5.dat
<<<<---- Data File Echo ---->>>>
TITLE - Four bars in series; applied loads and displ's. Same as testa but using equivalent bars.
NODES - # of nodes: 17; # 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.25
7 7 1.5
8 8 1.75
9 9 2.25
10 10 2.5
11 11 2.75
12 12 3.25
13 13 3.5
14 14 3.75
15 15 4.25
16 16 4.5
17 17 4.75
ELEMENTS - # of elements: 4
# ID elemType propID connectivity
17 1 pBar5 1 1 6 7 8 2
17 2 pBar5 2 2 9 10 11 3
17 3 pBar5 3 3 12 13 14 4
17 4 pBar5 4 4 15 16 17 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 =
130.2910 -181.1640 80.6349 -38.9418 9.1799
-181.1640 440.2116 -375.8730 155.7672 -38.9418
80.6349 -375.8730 590.4762 -375.8730 80.6349
-38.9418 155.7672 -375.8730 440.2116 -181.1640
9.1799 -38.9418 80.6349 -181.1640 130.2910
fe =
0
0
0
0
0
ke =
104.2328 -144.9312 64.5079 -31.1534 7.3439
-144.9312 352.1693 -300.6984 124.6138 -31.1534
64.5079 -300.6984 472.3810 -300.6984 64.5079
-31.1534 124.6138 -300.6984 352.1693 -144.9312
7.3439 -31.1534 64.5079 -144.9312 104.2328
fe =
0
0
0
0
0
ke =
208.4656 -289.8624 129.0159 -62.3069 14.6878
-289.8624 704.3386 -601.3968 249.2275 -62.3069
129.0159 -601.3968 944.7619 -601.3968 129.0159
-62.3069 249.2275 -601.3968 704.3386 -289.8624
14.6878 -62.3069 129.0159 -289.8624 208.4656
fe =
0
0
0
0
0
ke =
52.1164 -72.4656 32.2540 -15.5767 3.6720
-72.4656 176.0847 -150.3492 62.3069 -15.5767
32.2540 -150.3492 236.1905 -150.3492 32.2540
-15.5767 62.3069 -150.3492 176.0847 -72.4656
3.6720 -15.5767 32.2540 -72.4656 52.1164
fe =
0
0
0
0
0
Stiffness Matrix and Load Vector after assembly
kglobal =
Columns 1 through 7
130.2910 9.1799 0 0 0 -181.1640 80.6349
9.1799 234.5238 7.3439 0 0 -38.9418 80.6349
0 7.3439 312.6984 14.6878 0 0 0
0 0 14.6878 260.5820 3.6720 0 0
0 0 0 3.6720 52.1164 0 0
-181.1640 -38.9418 0 0 0 440.2116 -375.8730
80.6349 80.6349 0 0 0 -375.8730 590.4762
-38.9418 -181.1640 0 0 0 155.7672 -375.8730
0 -144.9312 -31.1534 0 0 0 0
0 64.5079 64.5079 0 0 0 0
0 -31.1534 -144.9312 0 0 0 0
0 0 -289.8624 -62.3069 0 0 0
0 0 129.0159 129.0159 0 0 0
0 0 -62.3069 -289.8624 0 0 0
0 0 0 -72.4656 -15.5767 0 0
0 0 0 32.2540 32.2540 0 0
0 0 0 -15.5767 -72.4656 0 0
Columns 8 through 14
-38.9418 0 0 0 0 0 0
-181.1640 -144.9312 64.5079 -31.1534 0 0 0
0 -31.1534 64.5079 -144.9312 -289.8624 129.0159 -62.3069
0 0 0 0 -62.3069 129.0159 -289.8624
0 0 0 0 0 0 0
155.7672 0 0 0 0 0 0
-375.8730 0 0 0 0 0 0
440.2116 0 0 0 0 0 0
0 352.1693 -300.6984 124.6138 0 0 0
0 -300.6984 472.3810 -300.6984 0 0 0
0 124.6138 -300.6984 352.1693 0 0 0
0 0 0 0 704.3386 -601.3968 249.2275
0 0 0 0 -601.3968 944.7619 -601.3968
0 0 0 0 249.2275 -601.3968 704.3386
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
Columns 15 through 17
0 0 0
0 0 0
0 0 0
-72.4656 32.2540 -15.5767
-15.5767 32.2540 -72.4656
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
176.0847 -150.3492 62.3069
-150.3492 236.1905 -150.3492
62.3069 -150.3492 176.0847
load_vector =
0
0
0
0
0
0
0
0
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
0
0
0
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.26
7 0.24
8 0.22
9 0.15179
10 0.10357
11 0.055357
12 0.0017857
13 -0.0035714
14 -0.0089286
15 -0.010714
16 -0.0071429
17 -0.0035714
Strain Energy (assume homogeneous displacement BCs) = 9.471667
ELEMENTAL POSTPROCESSING
ID f_i Sxx(end nodes) energy
1 2 -4.44e-15 -3.55e-15 -2.13e-14 -2 -2 -2 0.08
stress_nd =
-2.0000 -2.0000 -2.0000 -2.0000 -2.0000
ID f_i Sxx(end nodes) energy
2 3.86 3.77e-15 -6.11e-16 -1.33e-15 -3.86 -3.86 -3.86 0.372
stress_nd =
-3.8571 -3.8571 -3.8571 -3.8571 -3.8571
ID f_i Sxx(end nodes) energy
3 0.857 1.11e-16 0 -8.88e-16 -0.857 -0.857 -0.857 0.00918
stress_nd =
-0.8571 -0.8571 -0.8571 -0.8571 -0.8571
ID f_i Sxx(end nodes) energy
4 -0.143 -2.5e-16 1.11e-16 -2.22e-16 0.143 0.143 0.143 0.00102
stress_nd =
0.1429 0.1429 0.1429 0.1429 0.1429
Strain Energy = 0.462143
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -