Civil Engineering Reference
In-Depth Information
There are 30 equations
Result at node 18
time load x-disp y-disp cg iters
0.0000E+00 0.1000E+01 0.0000E+00 0.0000E+00
0.1000E+01 0.9553E+00 0.7363E+00 0.2167E+01 17
0.2000E+01 0.8253E+00 0.2231E+01 0.5226E+01 17
0.3000E+01 0.6216E+00 0.3252E+01 0.7398E+01 17
0.4000E+01 0.3624E+00 0.3987E+01 0.1046E+02 17
0.5000E+01 0.7074E-01 0.4832E+01 0.1398E+02 16
0.6000E+01 -0.2272E+00 0.5272E+01 0.1678E+02 16
0.7000E+01 -0.5048E+00 0.5086E+01 0.1839E+02 16
0.8000E+01 -0.7374E+00 0.4370E+01 0.1894E+02 17
0.9000E+01 -0.9041E+00 0.3255E+01 0.1834E+02 17
0.1000E+02 -0.9900E+00 0.1789E+01 0.1650E+02 17
0.1100E+02 -0.9875E+00 0.8842E-01 0.1349E+02 18
0.1200E+02 -0.8968E+00 -0.1674E+01 0.9574E+01 18
0.1300E+02 -0.7259E+00 -0.3373E+01 0.5082E+01 18
0.1400E+02 -0.4903E+00 -0.4916E+01 0.3731E+00 17
0.1500E+02 -0.2108E+00 -0.6151E+01 -0.4127E+01 17
0.1600E+02 0.8750E-01 -0.6859E+01 -0.7907E+01 17
0.1700E+02 0.3780E+00 -0.6839E+01 -0.1050E+02 17
0.1800E+02 0.6347E+00 -0.5986E+01 -0.1157E+02 18
0.1900E+02 0.8347E+00 -0.4304E+01 -0.1096E+02 18
0.2000E+02 0.9602E+00 -0.1902E+01 -0.8706E+01 18
Figure 11.17
Results from Program 11.6 example
Program 11.7 Forced vibration analysis of an elastic-plastic (von Mises) solid in
plane strain using rectangular 8-node quadrilateral elements. Lumped mass. Mesh
numbered in
x
-or
y
-direction. Explicit time integration.
PROGRAM p117
!-------------------------------------------------------------------------
! Program 11.7 Forced vibration analysis of an elastic-plastic (Von Mises)
! solid in plane strain using rectangular 8-node quadrilateral
! elements. Lumped mass. Mesh numbered in x- or y-direction.
! Explicit time integration.
!-------------------------------------------------------------------------
USE main; USE geom; IMPLICIT NONE
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::i,j,iel,k,loaded_nodes,ndim=2,ndof=16,nels,neq,nip=4,nn,nod=8, &
nodof=2,npri,nprops=4,np_types,nr,nres,nst=4,nstep,nxe,nye
REAL(iwp)::area,det,dsbar,dtim,f,fac,fmax,fnew,lode_theta,one=1.0_iwp,
&
pt5=0.5_iwp,sigm,time,zero=0.0_iwp
CHARACTER(LEN=15)::element='quadrilateral'
!-----------------------dynamic arrays------------------------------------
INTEGER,ALLOCATABLE::etype(:),g(:),g_g(:,:),g_num(:,:),nf(:,:),node(:), &
num(:)
REAL(iwp),ALLOCATABLE::bdylds(:),bee(:,:),bload(:),coord(:,:),dee(:,:), &
der(:,:),deriv(:,:),d1x1(:),d2x1(:),eld(:),eload(:),eps(:),
&
etensor(:,:,:),diag(:),g_coord(:,:),jac(:,:),mm(:,:),pl(:,:),
&
points(:,:),prop(:,:),sigma(:),stress(:),tensor(:,:,:),val(:,:),
&
weights(:),x1(:),x_coords(:),y_coords(:)
!-----------------------input and initialisation--------------------------
OPEN(10,FILE='fe95.dat'); OPEN(11,FILE='fe95.res')
READ(10,*)nxe,nye,np_types; CALL mesh_size(element,nod,nels,nn,nxe,nye)
ALLOCATE(nf(nodof,nn),points(nip,ndim),weights(nip),g_coord(ndim,nn),
&
