Civil Engineering Reference
In-Depth Information
There are 124 equations
Number of cg iterations to convergence was 64
Node x-disp y-disp z-disp
1 0.0000E+00 0.0000E+00 -0.2246E-01
2 0.1584E-02 0.0000E+00 -0.2255E-01
3 0.3220E-02 0.0000E+00 -0.2333E-01
4 0.0000E+00 0.0000E+00 -0.1849E-01
5 0.1544E-02 0.0000E+00 -0.1884E-01
6 0.0000E+00 0.0000E+00 -0.1443E-01
7 0.7580E-03 0.0000E+00 -0.1435E-01
8 0.1511E-02 0.0000E+00 -0.1411E-01
9 0.0000E+00 0.0000E+00 -0.6164E-02
10 0.2792E-02 0.0000E+00 -0.6430E-02
.
.
.
66 0.0000E+00 0.1572E-02 -0.1028E-03
67 -0.7437E-04 0.1716E-02 -0.5845E-04
68 0.0000E+00 0.0000E+00 0.0000E+00
69 0.0000E+00 0.0000E+00 0.0000E+00
70 0.0000E+00 0.0000E+00 0.0000E+00
The integration point (nip= 1) stresses are:
Element x-coord y-coord z-coord
sig_x sig_y sig_z tau_xy tau_yz tau_zx
1 0.2500E+00 0.5000E+00 -0.5000E+00
-0.2671E-01 -0.1647E+00 -0.9088E+00 0.6148E-02 0.9598E-01 0.4352E-02
2 0.2500E+00 0.5000E+00 -0.1500E+01
0.3986E-01 -0.5316E-01 -0.6298E+00 -0.2139E-02 0.7613E-01 0.4169E-02
3 0.2500E+00 0.1500E+01 -0.5000E+00
-0.2481E-01 -0.1260E+00 -0.1052E+00 0.4842E-02 0.9398E-01 -0.2812E-02
4 0.2500E+00 0.1500E+01 -0.1500E+01
0.2477E-01 -0.8240E-01 -0.2822E+00 -0.3176E-02 0.1214E+00 0.2938E-02
5 0.2500E+00 0.2500E+01 -0.5000E+00
0.3756E-02 0.8610E-02 -0.1470E-01 -0.9130E-03 -0.8730E-02 0.8741E-03
6 0.2500E+00 0.2500E+01 -0.1500E+01
0.7401E-02 -0.4390E-01 -0.2832E-01 -0.5657E-03 0.6083E-01 0.4901E-04
Figure 5.37 Results from Programs 5.5 and 5.6 example
Program 5.6 Three-dimensional strain of an elastic solid using 8-, 14-, or 20-node
brick hexahedra. Mesh numbered in
x
-
z
planes then in the
y
-direction. No global stiff-
ness matrix assembly. Diagonally preconditioned conjugate gradient solver. Vectorised
version.
PROGRAM p56
!-------------------------------------------------------------------------
! Program 5.6 Three-dimensional strain of an elastic solid using
! 8-, 14- or 20-node brick hexahedra. Mesh numbered in x-z
! planes then in the y-direction. No global stiffness matrix
! assembly. Diagonally preconditioned conjugate gradient solver.
! Vectorised version.
!-------------------------------------------------------------------------
USE main; USE geom; IMPLICIT NONE
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::cg_iters,cg_limit,fixed_freedoms,i,iel,k,loaded_nodes,ndim=3, &
ndof,nels,neq,nip,nn,nprops=2,np_types,nod,nodof=3,nr,nst=6,nxe,nye,nze
REAL(iwp)::alpha,beta,big,cg_tol,det,one=1.0_iwp,penalty=1.0e20_iwp,up,
&
zero=0.0_iwp
CHARACTER(LEN=15)::element='hexahedron'; LOGICAL::cg_converged
