D
DOCUMENTATION AND
USER INSTRUCTIONS OF
PROGRAM `SANOPT-P'
Appendix D contains the documentation for the program SANOPT-P
developed in Chapters 4-9. It gives some information about the general
features of the programs. It also presents detailed input instructions and
the specimen input data les.
D-1
D-2
D. Documentation and user instructions of program SANOPT-P
D. Documentation and User Instructions of
Program SANOPT-P
D.1 Program Documentation
D.1.1 Overview of the Program
Program SANOPT-P deals with the post processing of results obtained with
the linear elastic stress analysis and shape optimization of shells of revolution
and also prismatic shells of rectangular and curved planform which are supported on diaphragms at two opposite edges with the other two edges arbitrarily restrained. It plots the results of SANOPT-S and SANOPT-F programs.
D.1.2 Block Structure and Main Routines
A modular approach is adopted, in that separate subroutines are employed to
perform the various plotting operation. Each module in turn is composed of
one or more subroutines relevant only to its own needs and, in some cases, of
subroutines which are common to several modules. The main segment of the
program structure for SANOPT-P is given in Figure D.1. Only the important
subroutines are shown. The modules are described in relation to their general
functions as follows:
1. FIMA handles the le management. This subroutine open and close the
necessary data les.
2. DATAALL reads the plotting parameters.
3. DATAPOST reads input data including, control data, element connectivity, nodal point coordinates, boundary conditions, element thicknesses,
deections and stresses. This subroutine reads the static analysis results.
4. DATAmode reads input data including, control data, element connectivity,
coordinates of points, boundary conditions, element thicknesses, modal
amplitudes of each harmonic. This subroutine reads the free vibration or
buckling analysis results.
5. POSTPLO set up the plotting parameters.
6. PLOT MESHPO1 plots the generated mesh with the xed nodes and the
loads with point loads marked (static analysis results).
7. PLOT MESHPAPO plots the generated mesh with the xed nodes and
the required parameters (static analysis results).
8. PLOT MESHPO2 plots the generated mesh with the xed nodes and the
loads with point loads marked (static analysis results).
9. PLOT MESHMOP1 plots the generated mesh with the xed nodes and
the required parameters (free vibration or buckling analysis results).
D.1 Program documentation
D-3
FIMA
DATAALL
DATAPOST
DATAMODE
FIMA
POSTPLO
SANOPT-P
PLOT MESHPO1
PLOT MESHPAPO
PLOT MESHPO2
PLOT MESHPAPO
PLOT MESHMOP1
PLOT MESHMOP2
Fig. D.1.
Main block diagram of program SANOPT-P.
10. PLOT MESHMOP2 plots the generated mesh with the xed nodes and
the required parameters (free vibration or buckling analysis results).
D.1.3 File Structures
Program SANOPT-S uses the following les for input and output:
unit 5 (from keyboard) reads name of the les and some plotting control
parameters.
unit 6 (print screen) some basic data of preprocessing.
unit 10 reads the result le of analysis programs SANOPT-S and SANOPTF. This les include control data, element connectivity, coordinates of points,
boundary conditions, and element thicknesses It also reads stresses and deections for static and modal amplitudes for free and buckling problems.
D.1.4 Main Dimensions and Limitations of the Program
The program at present can handle 1000 elements/strips and 1000 points.
D-4
D. Documentation and user instructions of program SANOPT-P
D.2 Input Instructions
D.2.1 Main Structure of Input Data File for Static Problems
In the following section user instructions for preparing the input data are
presented for static problems. Note that this data le is automatically prepared
by SANOPT-S. The program reads following information from unit 10.
Data set 1: Plotting control parameter
cols.
variable name
| One record (*)
description
ijob
type of the problem
= 1, axisymmetric static
= 2, nite strip static
= 3, axisymmetric dynamic
= 4, nite strip dynamic
iinit
plotting control parameter
= 1, use default values
6= 1, user dened plotting variables
N.B: default values can be changed interactively.
Data set 2: Title of problem
| One record (a80)
cols.
variable name
description
1-80
title
Title of the problem to be plotted
Data set 3: Control parameters
cols.
| One record (16i5)
variable name
description
npoib
Number of nodal point in geometry
neleb
Number of element
nnode
Number of nodes in each element or strip
nnbon
Number of boundary conditions
nnlon
Number of point load applied to structure
npara
Number of parameters to be plotted
= 9 for static problems (see Data Set 1.9)
D.2 Input instructions
D-5
Data set 4: Nodal point coordinates
One record for each nodal point. Total of npoib records (*)
cols.
variable name
jpoin
coorb(jpoin,1)
coorg(jpoin,2)
description
point number
x1 -coordinate of nodal point
x2 -coordinate of nodal point
Data set 5: Nodal displacement data
One record for each nodal point. Total of npoib records (*)
cols.
variable name
jpoin
displ(jpoin,1)
displ(jpoin,2)
description
point number
u-displacement of nodal point
w-displacement of nodal point
Data set 6: Element connectivity data
One record for each element. Total of neleb records (*)
cols.
variable name
description
jelem
element number
lnodb(jelem,1)
the rst node number of. the element
..
..
..
..
..
.
lnodb(jelem,inode) the nth node number of the element
Data set 7: Fixity condition data
One record for each restrained point. Total of nnbog records (*)
cols.
variable name
nboug(iboun)
icode(iboun,1)
icode(iboun,2)
icode(iboun,3))
description
Nodal point number at which boundary condition is specied
Constraint on u-displacement,
=0, no constraint
=1, constraint
Constraint on v -displacement
Constraint on w-displacement
Data set 8: Loading data
One record for each point load. Total of nnlon records (*)
cols.
variable name
nloan(iload)
rload(iload,1)
rloag(iload,2)
rloag(iload,3)
description
Key point number at which load is applied
Value of point load in x-direction
Value of point load in y -direction
Value of point load in x-direction
D-6
D. Documentation and user instructions of program SANOPT-P
D.2.2 Main Structure of Input Data File for Free Vibration and
Buckling Problems
In the following section user instructions for preparing the input data are presented for free vibration and buckling problems. Note that this data le is
automatically prepared by SANOPT-F. The program reads following information from unit 10.
Data set 1: Plotting control parameter
cols.
variable name
| One record (*)
description
ijob
type of the problem
= 1, axisymmetric static
= 2, nite strip static
= 3, axisymmetric dynamic
= 4, nite strip dynamic
iinit
plotting control parameter
= 1, use default values
6= 1, user dened plotting variables
N.B: default values can be changed interactively.
Data set 2: Title of problem
| One record (a80)
cols.
variable name
description
1-80
title
Title of the problem to be plotted
Data set 3: Control parameters
cols.
| One record (16i5)
variable name
description
npoib
Number of nodal point in geometry
neleb
Number of element
nnode
Number of nodes in each element or strip
nnbon
Number of boundary conditions
npara
Number of mode (or harmonic) to be plotted
= 9 for static problems (see Data Set 1.9)
D.2 Input instructions
D-7
Data set 4: Nodal point coordinates
One record for each nodal point. Total of npoib records (*)
cols.
variable name
jpoin
description
point number
coorb(jpoin,1)
x1 -coordinate of nodal point
coorg(jpoin,2)
x2 -coordinate of nodal point
Data set 5: Element connectivity data
One record for each element. Total of neleb records (*)
cols.
variable name
jelem
description
element number
lnodb(jelem,1)
the rst node number of. the element
..
..
..
..
..
.
lnodb(jelem,inode) the nth node number of the element
Data set 6: Element thicknesses data
One record for each element. Total of neleb records (*)
cols.
variable name
jelem
description
element number
ethic(jelem,1)
thickness of the rst node
..
.. of the element
..
..
..
..
ethic(jelem,inode) thickness of the the nth node of the element
Data set 7: Fixity condition data
One record for each restrained point. Total of nnbog records (*)
cols.
variable name
nboug(iboun)
description
Nodal point number at which boundary condition is specied
icode(iboun,1)
Constraint on u-displacement,
=0, no constraint
=1, constraint
icode(iboun,2)
Constraint on v -displacement
icode(iboun,3))
Constraint on w-displacement
D-8
D. Documentation and user instructions of program SANOPT-P
Data set 8: Parameters for plotting
One record for each point. Total of npoibnpara records (*)
cols.
variable name
description
jpoin
Point number
jpara
Parameter number to be plotted
paras(jpoin,jpara,idime) u- value of nodal point
paras(jpoin,jpara,idime) w-value of nodal point
D.3 Specimen Data Files
In this section output graphic les are provided for a numerical example
presented in the book. Input data les are produced by analyses programs
SANOPT-S and SANOPT-F. The SANOPT-P input result le is not provided for that reason. This information will be of assistance to readers who
wish to run the program contained in the book on their own computer. For
clarity of the manual, presentation is limited to one static and one optimization
examples.
D.3.1 Cylindrical Shell Roof Subjected to Self Weight Loading
The cylindrical shell roof analysed in Chapter 5 is considered for post processing plotting example. It is supported by diaphragms at each curved end and is
free along the other two straight edges. The output of the program for static
problems are given in Figures D.2-D.10.
The output postscript le for cylindrical shell roof subjected to self weight
loading: displacement
Fig. D.2.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Membrane force Ny
Fig. D.3.
D.3 Specimen data les
D-9
The output postscript le for cylindrical shell roof subjected to self weight
loading: Membrane force N`y
Fig. D.4.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Bending moment M`
Fig. D.5.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Bending moment My
Fig. D.6.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Bending moment M`y
Fig. D.7.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Shear force Q`
Fig. D.8.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Shear force Qy
Fig. D.9.
The output postscript le for cylindrical shell roof subjected to self weight
loading: Equivalent stress eq
Fig. D.10.
D-10
D. Documentation and user instructions of program SANOPT-P
D.3.2 Thin Circular Plate
The thin circular plate analysed in Chapter 7 is considered for post processing
plotting example. The modal amplitudes associated with various frequencies
are given in Figures D.11-D.15.
The output postscript le for modal amplitudes associated with various
frequencies of thin circular plate: 1st harmonic
Fig. D.11.
The output postscript le for modal amplitudes associated with various
frequencies of thin circular plate: 2nd harmonic
Fig. D.12.
D.3 Specimen data les
D-11
The output postscript le for modal amplitudes associated with various
frequencies of thin circular plate: 3rd harmonic
Fig. D.13.
The output postscript le for modal amplitudes associated with various
frequencies of thin circular plate: 4th harmonic
Fig. D.14.
D-12
D. Documentation and user instructions of program SANOPT-P
The output postscript le for modal amplitudes associated with various
frequencies of thin circular plate: 5th harmonic
Fig. D.15.