Civil Engineering Reference
In-Depth Information
This job ran on
32 processors
There are 12465 nodes 225 restrained and 36720 equations
Time after setup is : 0.339999999996507540
Time t cos(omega*t) Displacement Iterations
0.3142E+02 0.9511E+00 0.3469E+03 428
0.6283E+02 0.8090E+00 0.5158E+03 420
0.9425E+02 0.5878E+00 0.4672E+03 419
.
.
.
0.9425E+03 -0.1000E+01 -0.5935E+03 421
0.9739E+03 -0.9511E+00 -0.5696E+03 421
0.1005E+04 -0.8090E+00 -0.4899E+03 422
This analysis took : 562.680000000000291
Figure 12.41
Results from Program 12.9 example
Mesh
No of Processors
Analysis Time(secs)
8 x40 x 8 16 1601
32 562.7
Figure 12.42
Performance statistics: Program 12.9 (IBM SP2)
Program 12.10 Forced vibration analysis of three-dimensional elasto-plastic solid.
Explicit integration in time. Compare Program 11.5.
PROGRAM p1210
!-------------------------------------------------------------------------
! Program 11.7 forced vibration of an elastic-plastic(Von Mises) solid
! using 20-node hexahedral elements; nxe even; viscoplastic strain method
! regular box mesh : lumped mass, explicit integration: parallel version
!-------------------------------------------------------------------------
USE new_library; USE geometry_lib; USE mp_module; USE timing; USE utility
USE global_variables1; USE precision; USE gather_scatter6; IMPLICIT NONE
! ndof,nels,neq,ntot are now global variables - not declared
INTEGER::nxe,nye,nze,nn,nr,nip,nodof=3,nod=20,nst=6,loaded_nodes,nres,
&
i,j,k,jj,iel,ndim=3,nstep ,npri,num_no,no_index_start ,
&
neq_temp,nn_temp, is,it
REAL(iwp)
::aa,bb,cc,rho,dtim,e,v,det,sbary,pload,sigm,f,fnew,fac,
&
volume,sbar,dsbar,lode_theta,real_time
CHARACTER (LEN = 15) :: element = 'hexahedron'
!---------------------------- dynamic arrays------------------------------
REAL(iwp), ALLOCATABLE :: points(:,:),bdylds_pp(:),x1_pp(:),d1x1_pp(:), &
stress(:),pl(:,:),emm(:),d2x1_pp(:),tensor_pp(:,:,:), &
etensor_pp(:,:,:),val(:),dee(:,:),coord(:,:),mm_pp(:),&
jac(:,:),weights(:), der(:,:),deriv(:,:),bee(:,:), &
eld(:),eps(:),sigma(:),bload(:),eload(:),mm_tmp(:,:), &
p_g_co_pp(:,:,:),pmul_pp(:,:),utemp_pp(:,:)
INTEGER, ALLOCATABLE :: rest(:,:), g(:), no(:), num(:),no_local(:),
&
g_num_pp(:,:), g_g_pp(:,:),no_local_temp(:)
!-----------------------input and initialisation--------------------------
timest(1) = elap_time( )
; CALL find_pe_procs(numpe,npes)
