J
DOCUMENTATION AND
USER INSTRUCTIONS OF
PROGRAM `PLATES'
Appendix J contains a brief documentation, user instructions and specimen input and output data les for program PLATES discussed in Chapter 11.
J-1
J-2
J. Documentation and user instructions of program PLATES
J. Documentation and User Instructions of
Program PLATES
J.1 Program Documentation
J.1.1 Overview of the Program
In this section program PLATES which performs the closed form plate analysis
is presented.
J.1.2 Block Structure and Main Routines
This program consists of a main segment and six subroutines as shown in
Figure J.1. The subroutines are now briey described:
FIMA
DATA
INIT
PLATES
COEF
CONS
SUMS
OUTP
FIMA
Fig. 1.1. Flow diagram of program PLATES.
1. MASTER is the main subroutine and it calls all of the other subroutines
in order.
2. FIMA handles the le management. This subroutine open and close the
necessary data les.
2. DATA reads in all the input data values.
3. INIT zeroes all initial values apart from the input values.
4. COEF calculates the coeÆcients P11 ; P12 ; P13 ; etc.
5. CONS calculates the constant values Amn , Bmn and Cmn for a given loading
function.
J.2 Input instructions
J-3
6. SUMS sums the various Fourier series.
7. OUTP prints out the results for each time step in transient dynamic analysis.
J.2 Input Instructions
J.2.1 Glossary of Variable Names
aa,bb
amn,bmn,cmn
bmx,bmy, bmxy
cx,cy,cxy
d11,d22,d12,d21,d66
det
eta,etb
g1
gfs
it
iq
mpon
npon
num
p11,p12,p13, etc.
pi
pz
qx,qy
qmn
s44,s55
title
u,v
w,wux,wuy
xcord,ycord
Length of plate in the x; y directions respectively
unknown Fourier coeÆcients Amn , Bmn and Cmn
Mending moments Mx ; My ; Mxy
Curvatures x;x , y;y , (x;y + y;x )
Flexural rigidities D11 ; D22 ; D12 ; D21 ; D66
determinant of matrix P
x; y coordinates of centroid of load
Shear modication factor
Elastic foundation modulus K
= 1 for non-symmetric loads and
= 2 for symmetric loads
= 1 for uniform loads,
= 2 for concentrated loads,
= 3 for uniform patch loads and
= 4 for cosine patch loads
Maximum number of output points (140)
Number of output points
Number of Fourier terms
CoeÆcients used to evaluate Amn ; Bmn ; Cmn
Load intensity
Shear forces
The term qmn in the Fourier series representation
of the load
Shear rigidities S44 ; S55
Problem title
Dimensions of patch load
w; x ; y
X; Y coordiantes of output points
J-4
J. Documentation and user instructions of program PLATES
J.2.2 Main Structure of Input Data
Input instructions for program PLATES are now provided.
Data set 1: Problem title
| One record (20a4)
cols. variable name
1-80
title
description
Title of the problem
Data set 2: Plate dimensions
cols. variable name
| One record (2f10.5)
description
1-10
aa
side length parallel to x axis
11-20
bb
side length parallel to y axis
Data set 3: Flexural properties
cols. variable name
| One record (4f10.5)
description
1-10
d11
exural rigidity D11
11-20
d22
exural rigidity D22
21-30
d12
exural rigidity D12
31-40
d66
exural rigidity D66
Data set 4: Shear properties
cols. variable name
| One record (2f10.5)
description
1-10
s44
shear rigidity S44
11-20
s55
shear rigidity S55
Data set 5: Control data
cols. variable name
| One record (2i5)
description
1-5
num
number of Fourier terms
6-10
npon
number of output points
J.2 Input instructions
Data set 6: Further plate parameters
cols. variable name
| One record (3f10.5)
description
1-10
pz
load intensity
11-20
gfs
foundation modulus
21-30
g1
shear modication factor
Data set 7: Sampling Points
cols. variable name
| npon records (i5,2f10.5)
description
1-5
jpon
point number
6-15
xcord(.)
x coordinate
16-25
ycord(.)
x coordinate
Data set 8: Loading type
| One record (2i5)
cols. variable name
description
16-20
iq
= 1, uniform loading
= 2, concentrated loading
= 3, uniform loading over rectangular
patch u v
= 4, cosine loading over rectangular
patch u v
21-25
it
= 1, non-symmetric loads
= 2, symmetric loads
Data set 9: Loading position
Omit this data set, if iq < 2
cols. variable name
| One record (2f10.5)
description
1-10
eta
x coordinate of patch parallel to x axis
11-20
etb
x coordinate of patch parallel to x axis
Data set 11: Load patch size
Omit this data set, if iq < 2
cols. variable name
| One record (2f10.5)
description
1-10
u
side length of patch parallel to x axes
11-20
v
side length of patch parallel to x axes
J-5
J-6
J. Documentation and user instructions of program PLATES
J.3 Specimen Input Data File
In this section input and output data les are provided for a numerical example
presented in the text. This information will be of assistance to readers who
wish to run the program contained in the book on their own computer.
J.3.1 Homogeneous Isotropic Rectangular Thin Plates
Problem denition: The analysis of a thin, simply supported, homogeneous
isotropic, square plate is now presented for uniformly distributed loading q .
The plate considered has the following properties: thickness-to-span ratio t=a =
0:1, aspect ratios b=a = 1:0 and D11 = D22 = 1:0, D12 = 0:3, D66 = 0:35,
S44 = S55 = 35000:0, = 0:3 and intensity of load 1.0.
A typical input and output les for P LAT ES is displaced in Table J.1 and
J.2.
Table J.1.
The input data le for isotropic rectangular thin plate.
data set
input data for square uniformly loaded plate (t/l=0.01)
1
1.0
1.0
2
1.0
1.0
0.3
0.35
3
35000.0 35000.0
4
51
3
5
1.0
0.0
1.0
6
1
0.0
0.0
7
2
0.5
0.5
3
0.5
0.0
1
2
8
Table J.2.
The output data le for isotropic rectangular thin plate.
_____
_______ _______ _______ _______
|_____] |
|_____|
|
|______ |______
|
|_____ |
|
|
|______ ______|
Closed Form Solution for Orthotropic Mindlin Plates
Copyright (c) E. Hinton, J. Sienz and M. Ozakca
Swansea, 1988 - 2001
J.3 Specimen input data le
J-7
input data for square uniformly loaded plate (t/l=0.01)
**************************************************
dimension of plate in x direction.(aa)= 0.10000E+01
dimension of plate in y direction.(bb)= 0.10000E+01
flexural rigidity....(d11)=
flexural rigidity....(d22)=
coupling rigidity....(d12)=
torsional rigidity...(d66)=
0.10000000E+01
0.10000000E+01
0.30000000E+00
0.35000000E+00
transverse shear rigidity.....(s44)= 0.35000000E+05
transverse shear rigidity.....(s55)= 0.35000000E+05
shear correction factor.........(g1)= 0.10000000E+01
modulus for elastic foundation.(gfs)= 0.00000000E+00
load intensity..................(pz)= 0.10000000E+01
the no. of terms in fourier series (num)=
the no. of points on plate
(npon)=
51
3
**************************************************
table shows the computed values for
specified x and y coordinates.
for uniform loading.
symmetric loading.
x
y
w
wux
wuy
0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.00000E+00
0.50000E+00 0.50000E+00
0.40645E-02 0.51252E-09 0.51252E-09
J-8
J. Documentation and user instructions of program PLATES
0.50000E+00 0.00000E+00
x
y
0.00000E+00 0.00000E+00 -0.13482E-01
mx
my
mxy
0.00000E+00 0.00000E+00
0.00000E+00 0.00000E+00 -0.32476E-01
0.50000E+00 0.50000E+00
0.47886E-01 0.47886E-01 -0.47205E-16
0.50000E+00 0.00000E+00
0.00000E+00 0.00000E+00 0.11520E-08
x
y
0.00000E+00 0.00000E+00
qx
qy
0.00000E+00 0.00000E+00
0.50000E+00 0.50000E+00 -0.67876E-08 -0.67876E-08
0.50000E+00 0.00000E+00
0.00000E+00 0.33372E+00