Coupled Problems - Programming the Finite Element Method - page 422

Ramp
load

1.0

t

t 0

Figure 9.9 Ramp loading

There are 32 equations and the skyline storage is 280

Results at node 21
time load x-disp y-disp porepressure
0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
0.1000E+00 0.2000E+00 0.0000E+00 0.0000E+00 -0.1994E+00
0.2000E+00 0.4000E+00 0.0000E+00 0.0000E+00 -0.3743E+00
0.3000E+00 0.6000E+00 0.0000E+00 0.0000E+00 -0.5125E+00
0.4000E+00 0.8000E+00 0.0000E+00 0.0000E+00 -0.6203E+00
0.5000E+00 0.1000E+01 0.0000E+00 0.0000E+00 -0.7043E+00
0.6000E+00 0.1000E+01 0.0000E+00 0.0000E+00 -0.5703E+00
0.7000E+00 0.1000E+01 0.0000E+00 0.0000E+00 -0.4463E+00
0.8000E+00 0.1000E+01 0.0000E+00 0.0000E+00 -0.3478E+00
0.9000E+00 0.1000E+01 0.0000E+00 0.0000E+00 -0.2709E+00
0.1000E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.2110E+00
0.1100E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1643E+00
0.1200E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1280E+00
0.1300E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.9968E-01
0.1400E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.7764E-01
0.1500E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.6047E-01
0.1600E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.4709E-01
0.1700E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.3668E-01
0.1800E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.2857E-01
0.1900E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.2225E-01
0.2000E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1733E-01
0.2100E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1350E-01
0.2200E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1051E-01
0.2300E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.8187E-02
0.2400E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.6377E-02
0.2500E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.4966E-02
0.2600E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.3868E-02
0.2700E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.3013E-02
0.2800E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.2346E-02
0.2900E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1827E-02
0.3000E+01 0.1000E+01 0.0000E+00 0.0000E+00 -0.1423E-02

Figure 9.10 Results from Program 9.3 example

Turning to the program, the nodal coordinates and steering vector are again provided
by subroutine geom rect with numbering in the x -direction, and the subroutine formke
builds up the ke matrix which is then assembled into the (symmetric) global skyline
matrix kv .