Scribd
Upload
833 views4 pages

Format of The Stiffness Matrix and Load Vector Output in Nastran PCH File

The document describes the output format of stiffness matrices and load vectors from Nastran and a Fortran code for a structural analysis problem. It shows the 8x8 stiffness matrix with elements corresponding to displacements of each of 4 nodes in x and y directions. It also shows the load vector with loads in x and y directions applied at each node. Comments explain that Nastran only outputs the symmetric part of the stiffness matrix in the .pch file and that the bar element has 2 degrees of freedom per node for planar structures.

Uploaded by

kspetsnaz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
833 views4 pages

Format of The Stiffness Matrix and Load Vector Output in Nastran PCH File

The document describes the output format of stiffness matrices and load vectors from Nastran and a Fortran code for a structural analysis problem. It shows the 8x8 stiffness matrix with elements corresponding to displacements of each of 4 nodes in x and y directions. It also shows the load vector with loads in x and y directions applied at each node. Comments explain that Nastran only outputs the symmetric part of the stiffness matrix in the .pch file and that the bar element has 2 degrees of freedom per node for planar structures.

Uploaded by

kspetsnaz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 4

y, v(m)

Element numbers

3 (0.15,0.26)

F3x
(2)

(5)

(3)

4 (0.5,0)
1 (0,0)

(1)

(4)

x, u (m)

F4y

2 (0.15,0)

Only symmetric part of stiffness matrices are output in .pch file. For the bar element we have 2 dof per node for planar structure. Therefore, each
node has x and y displacement dof.

(1,1)
(1,1)

(1,1)

(1,1)
(1,1)

(1,1)
(1,1)

(1,1)

(1,2)
(1,2)
(1,2)
(1,2)
(1,2)
(1,2)
(1,2)

( 2,1)
( 2,1)
( 2,1)
( 2,1)
( 2,1)
( 2,1)

( 2,2)
( 2,2)
( 2,2)
( 2,2)
( 2,2)

(3,1)
(3,1)
(3,1)
(3,1)

(3,2)
(3,2)
(3,2)

( 4,1)
( 4,1)

(1,1)
(1,2)

( 2,1)

( 2,2)
(3,1)

(
3
,
2
)

( 4,1)

( 4,4) ( 4,2)

( Node 1, xdof
( Node 1, ydof

( Node 2, xdof

( Node 2, ydof

( Node 3, xdof
( Node 3, ydof

( Node 4, xdof
( Node 4, ydof

1)
2)
1)

2)

1)
2)

1)
2)

Stiffness matrix output of Nastran in .pch file

DMIG
DMIG*
*
DMIG*
*
*
DMIG*
*
*
DMIG*
*
DMIG*
*
*
*
DMIG*
*
*
*
*
*
DMIG*
*
*
*
*

KAAX
KAAX

DMIG*
*
*
*
*

KAAX

KAAX

2
3
3
4

2
0
1
1
1 2.099611628D+07
1
2
1 4.037713011D+06
2 6.998702552D+06
2
1
1-1.866666647D+07
1 2.666666638D+07
2
2
2 1.076923065D+07
3
1
1-2.329449814D+06
2-4.037713011D+06
1 6.467739991D+06
3
2
1-4.037713011D+06
2-6.998702552D+06
2-1.076923065D+07
1 9.635545937D+05
2 2.005159375D+07
4
1
1-7.999999915D+06
1-4.138290177D+06
2 3.074158417D+06
1 1.213829009D+07

8 (order of the stiffness matrix. 8X8)


(first row up to the diagonal element=(1,1)
(element of stiffness matrix multiplying the dof (1,1))
(second row up to the diagonal element=(1,2)
(element of stiffness matrix multiplying the dof (1,1))
(element of stiffness matrix multiplying the dof (1,2))
(third row up to the diagonal element=(2,1)
(element of stiffness matrix multiplying the dof (1,1))
(element of stiffness matrix multiplying the dof (2,1))
(fourth row up to the diagonal element=(2,2)
(element of stiffness matrix multiplying the dof (2,2))
(fifth row up to the diagonal element=(3,1)
(element of stiffness matrix multiplying the dof (1,1))
(element of stiffness matrix multiplying the dof (1,2))
(element of stiffness matrix multiplying the dof (3,1))
(sixth row up to the diagonal element=(3,2)
(element of stiffness matrix multiplying the dof (1,1))
(element of stiffness matrix multiplying the dof (1,2))
(element of stiffness matrix multiplying the dof (2,2))
(element of stiffness matrix multiplying the dof (3,1))
(element of stiffness matrix multiplying the dof (3,2))
(seventh row up to the diagonal element=(4,1)
(element of stiffness matrix multiplying the dof (2,1))
(element of stiffness matrix multiplying the dof (3,1))
(element of stiffness matrix multiplying the dof (3,2))
(element of stiffness matrix multiplying the dof (4,1))

3
3
4
4

4
2
1 3.074158417D+06
2-2.283660538D+06
1-3.074158417D+06
2 2.283660538D+06

(eight row up to the diagonal element=(4,2)


(element of stiffness matrix multiplying the
(element of stiffness matrix multiplying the
(element of stiffness matrix multiplying the
(element of stiffness matrix multiplying the

1
1
1

KAAX
KAAX
KAAX

KAAX

KAAX

1
2
2
1
1
3
1
1
2
3
3

dof
dof
dof
dof

(3,1))
(3,2))
(4,1))
(4,2))

Stiffness matrix output of the fortran code


0.20996E+08 0.40377E+07 -0.18667E+08 0.00000E+00 -0.23294E+07 -0.40377E+07 0.00000E+00 0.00000E+00
0.40377E+07 0.69987E+07 0.00000E+00 0.00000E+00 -0.40377E+07 -0.69987E+07 0.00000E+00 0.00000E+00
-0.18667E+08 0.00000E+00 0.26667E+08 0.00000E+00 0.00000E+00 0.00000E+00 -0.80000E+07 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00 0.10769E+08 0.00000E+00 -0.10769E+08 0.00000E+00 0.00000E+00
-0.23294E+07 -0.40377E+07 0.00000E+00 0.00000E+00 0.64677E+07 0.96355E+06 -0.41383E+07 0.30742E+07
-0.40377E+07 -0.69987E+07 0.00000E+00 -0.10769E+08 0.96355E+06 0.20052E+08 0.30742E+07 -0.22837E+07
0.00000E+00 0.00000E+00 -0.80000E+07 0.00000E+00 -0.41383E+07 0.30742E+07 0.12138E+08 -0.30742E+07
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.30742E+07 -0.22837E+07 -0.30742E+07 0.22837E+07

Load vector output of Nastran


DMIG*
*
*
*
*
*
*
*
*

PAX
1
1
2
2
3
3
4
4

1
0
1-8.692391441D+03
2-5.020649943D+03
1 9.094947018D-13
2-5.796402259D+03
1-1.676312023D+03
2 1.427217856D+04
1 1.044870346D+04
2-3.806709905D+03

(load
(load
(load
(load
(load
(load
(load
(load

acting
acting
acting
acting
acting
acting
acting
acting

on
on
on
on
on
on
on
on

node
node
node
node
node
node
node
node

Load vector output of the fortran code


-0.86924E+04
-0.50207E+04
0.00000E+00
-0.57964E+04
-0.16763E+04
0.14272E+05

1
1
2
2
3
3
4
4

in
in
in
in
in
in
in
in

x
y
x
y
x
y
x
y

direction:1)
direction:2)
direction:1)
direction:2)
direction:1)
direction:2)
direction:1)
direction:2)

0.10449E+05
-0.38067E+04

You might also like