Next: 4 Utilities
Up: 3 The domain decomposition method
Prev: 3.1 Introduction
Contents
Main program:
PRSDOM
The two matrices above are stored in a D.S. MUA (one per sub-domain). The matrix of the Neumann problem is placed in an associated array. File SAUV contains a description of the sub-domain interfaces.
The finite element description used and the material constants are done for the whole problem. There is therefore only one D.S. MILI and one D.S. FORC for the collection of sub-domains. Similarly, the list of finite elements used is input in a global fashion. The same finite elements are used for all materials.
The user inputs a generic name: NOM
If the D.S. describes the collection of sub-domains (eg. MILI) its name will be NOM.nameDS (eg. NOM.MILI)
If there is a D.S. per sub-domain (eg. MAIL)
its name will be
NOM.nameDS.0001, NOM.nameDS.0002....
(eg. NOM.MAIL.0001, NOM.MAIL.0002 ...)
The input data structures are the following:
NOPO | called NOM.NOPO.0001 NOM.NOPO.0002 ... | |
MILI | called NOM.MILI | |
FORC | called NOM.FORC |
The output data structures are the following:
B | called NOM.B.0001 NOM.B.0002 ... |
AMAT | called NOM.AMAT.0001 NOM.AMAT.0002 ... |
MUA | called NOM.MUA.0001 NOM.MUA.0002 ... |
SAUV | called NOM.SAUV |
auxiliary file to describe the interfaces |
The data file contains, in sequence:
NOM (A) | generic name of the files |
NDSD (I) | number of sub-domains (sub-structures) |
POBA? (I) | utilization of file POBA |
YES 1 | |
NO 0 | |
NMPOBA (A) | name of file POBA if it is used |
NDIM (I) | dimension |
NBSC (I) | number of curved surfaces |
NOSC (I) | NOSC(NBSC) list of curved surfaces |
if NBSC is non-zero | |
NBLC (I) | number of curved surfaces |
NOLC (I) | NOLC(NBLC) list of curved surfaces |
if NBLC is non-zero | |
NOMBIB (A) | library name (4 characters) |
NTYED (I) | number of straight elements |
LNTYED (A) | list of straight elements (2 words of 4 characters) |
NTYEC (I) | number of curved elements |
LNTYEC (A) | list of curved elements (2 words of 4 characters) |
NBFR (I) | number of triplets describing the boundary conditions |
NOFR (I) | NOFR(NBFR) list of references |
NTYDL(I) | NTYDL(NBFR) list of prescribed d.o.f. |
MNEM (I) | MNEM(NBFR) list of mnemonics (in general VN) |
NTY (I) | variable type (2 single, 5 double precision) |
IOPT (I) | IOPT(4) of THELAS describes the mass, stiffness, RHS, |
stresses | |
0 if the array is not computed | |
1 if not | |
ND (I) | number of d.o.f. per node |
Main program:
SDOMVD
B | called NOM.B.0001 NOM.B.0002 ... | |
AMAT | called NOM.AMAT.0001 NOM.AMAT.0002 ... | |
MUA | called NOM.MUA.0001 NOM.MUA.0002 ... | |
SAUV | called NOM.SAUV auxiliary file describing the interfaces |
Main program:
SDOMVR
B | called NOM.B.0001 NOM.B.0002 ... | |
AMAT | called NOM.AMAT.0001 NOM.AMAT.0002 ... | |
MUA | called NOM.MUA.0001 NOM.MUA.0002 ... | |
SAUV | called NOM.SAUV auxiliary file to describe the interfaces |