Moduul PAREST: programmeurs handleiding

(1)

Citation for published version (APA):

Courage, W. M. G., & Hendriks, M. A. N. (1988). Moduul PAREST: programmeurs handleiding. (DCT rapporten; Vol. 1988.063). Technische Universiteit Eindhoven.

Document status and date: Gepubliceerd: 01/01/1988

Document Version:

Uitgevers PDF, ook bekend als Version of Record

Please check the document version of this publication:

• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.

• The final author version and the galley proof are versions of the publication after peer review.

• The final published version features the final layout of the paper including the volume, issue and page numbers.

Link to publication

General rights

Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain

• You may freely distribute the URL identifying the publication in the public portal.

If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:

www.tue.nl/taverne

Take down policy

If you believe that this document breaches copyright please contact us at:

openaccess@tue.nl

providing details and we will investigate your claim.

(2)

WFW 88.063

door Wim Courage en

Max

Hendriks

november 1988

T.U. Eindhoven

Afdeling Werktuigbouwkunde Vakgroep WFW

(3)

(4)

INHOUD O

.

INLEIDING

...

1

.

INPUT SEGMENTEN

...

1.1 Segment EXPERI

...

1.2 Segment ESTIMA

...

2

.

MODUUL PAREST

...

2.1 FILOS items

...

2.2 Segmenten

...

2.2.1 Segment ITCHCK

...

2.2.2 Segment ITNINI

...

2.2.3 Segment KALMAN

...

2.2.4.Segment NLCHCK

...

2.2.5 Segment NLINIT

...

2.2.6 Segment PARDEV

...

2.2.7 Segment PEINIT

...

2.2.8 Segment STCHCK

...

2.2.9 Segment STPINI

...

LITERATUUR

...

APPENDIX A: De jobstream

...

APPENDIX B: Filos items met bijbehorende segmenten

...

5 6 6 8 9 9 11 12 13 14 16 17 18 20 22 23 24 25 26

(5)

(6)

O. INLEIDING

PAREST is een moduul bedoeld voor het schatten van materiaal parameters ( PARa- meter ESTimation ). Het is geimplementeerd in het eindige elementen pakket DIANA. Deze programmeurshandleiding beschrijft het moduul PAREST en twee aan moduul INPUT toegevoegde segmenten en moet gezien worden als een uitbreiding van de DIANA programmeurshandleiding (Witte, 1988). Tevens is er een gebruikershandleiding voor PAREST verschenen (Hendriks, 1988) als uitbreiding van de DIANA gebruikershandleiding (Witte, 1988).

PAREST is gebaseerd op het idee het (extended) KALMAN-filter t e gebruiken als

parameter-schattingstechniek (Oomens, 1988; Hendriks, 1988). Dit is een iteratieve

techniek hetgeen tot uiting komt in het herhaaldelijk uitvoeren van complete eindige elementen sommen, met naar verwachting steeds betere waarden voor de materiaal par amet ers.

Deze eindige elementen berekeningen worden uitgevoerd met behulp van de andere modulen van het DIANA-pakket. Hiermee wordt bereikt dat in principe alle

modelmogelijkheden van DIANA t e gebruiken zijn bij het schatten van de parameters. Voor een toelichting op de theorie en de uitwerking van het KALMAN-filter wordt verwezen naar de gebruikershandleiding (Hendriks, 1988).

In hoofdstuk 1 worden de twee aan moduul INPUT toegevoegde segmenten beschreven, met voor elk segment:

- een omschrijving van de functie van het segment,

- een korte beschrijving van de routines, de besturing en de gegevensstroom van en naar de filosfile.

In hoofdstuk 2 wordt moduul PAREST beschreven. Hiervan wordt vermeld:

- een korte omschrijving van de taak van het moduul,

- een opsomming van de gegevens van de FILOS-file die in het moduul worden

- een opsomming van de segmenten waaruit het moduul bestaat, met voor elk gebruikt (gelezen en/of geschreven),

segment een omschrijving van de functie en een korte beschrijving van de routines, de besturing en de gegevensstroom van en naar de FILOS-file.

In appendix A wordt een schematische weergave gegeven van de jobstream zoals deze door de command interpreter van moduul PAREST wordt gegenereerd. In appendix B tenslotte wordt per FILOS-item uit moduul PAREST aangegeven in welke segmenten deze een rol spelen.

(7)

1.

INPUT SEGMENTEN

Voor een uitgebreide beschrijving van moduul INPUT wordt verwezen naar de DIANA-programmeurs handleiding (Witte, 1988). Hier worden slechts de twee toegevoegde segmenten EXPERI en ESTIMA besproken.

1.1 Segment EXPERI Pibd: /src/diana_parest/in/exii

Segment EXPERI interpreteert de tabel 'EXPERI' uit de invoerfile. Deze bevat de data over het experiment, te weten: de coordinaten van de meetpunten en voor elke tijd-

belastingstap de verplaatsingen en de covarianties van de verplaatsingen in de meetpunten. Routines: Pad: ex1 i

/

Ivli exll/lvll exi i /h l 2 exii/lv12 exil/lvl2 exi i /h l 3 ex1 i /lv13 ex1 i / I d 3 Routine mainex exper excoor exdisp excova comat dcmat rpmat Doel hoofd routine

verwerking van tabel 'EXPERI' subtabel coordinates

subtabel displacements subtabel covariances volle cov. matrix diagonaal cov. matrix

herhaling van andere tijd stap Besturing:

Het JCL record voor segment EXPERI wordt gegenereerd door aanroep van

JCLXEQ (zie progr. handleiding

,

Witte 1988, hfdst. 8) met mod = 2HIN

segexe = 6HEX11

ircsgp Segment stuur parameters?

<=

O : Nee

> O : Ja, Recordnummer van SF.SGMINF (output van sub- routine SEGPAR)

Recordnummer van SF.SGMINF, vult COMMON /SET/ Recordnummer van SF.SGMINF, vult COMMON /READ/

(8)

Genevenss t room Input SF.EX11 lfn=49 Output VL.EXNDIM int VL.EXNMSR int VL.EXNPNT int

AR.EXMPTR int (NSTEP,EXNPNT,EXNDIM) AR.EXSPTR int

AR.PNTFLG int SF.COVARD lfn=19 SF.EXCOOR lfn=23 SF.EXDISP lfn=25

Experimental data input.

Max. number of exp. values per experimental point.

Number of experiment al dat a values. Number of experimental points. Storage pointer for exp. displ.

Starting address pointer of exp. displ. Flag for (non)active exp. point. Covariances of exp. displacements. Coordinates of experimental points. Exp. displ. values in exp. points.

(9)

1.2 Segment ESTIMA Pad: /src/diana~arest/in/esii

Segment ESTIMA interpreteert de tabel 'ESTIMA' uit de invoerfile. Deze bevat de data over de te schatten parameters, te weten: de parameters, de covarianties van de parameters en voor elke tijd of belastingstap de covarianties van het model.

Rout ines: Pad: esli/lvli esli/lvll es1 i /lv12 es1 i /lv12 es1 i /lv12 eslí/lvl3 esll/lvl3 esíï/lvl3 Routine maines estim espara escova tabrcd comat dcmat rpmat Doel hoofd routine

verwerking van tabel 'ESTIMA' subt abel parameters

process par. subtable record. subtabel covariances par. or model volle cov. matrix

diagonaal cov. matrix

herhaling van andere tijd stap Besturing :

Het JCL record voor segment ESTIMA wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

Witte 1988, hfdst. 8) met

mod segexe ircsgp jclpar jclparYI 2 = 2HIN = 6HES11

Segment stuur Parameters?

<=

O : Nee

> O : Ja, Recordnummer van SF.SGMINF (output van sub- routine SEGPAR)

Gegevensstroom Input

Recordnummer van SF.SGMINF, vult COMMON /SET/ Recordnummer van SF.SGMINF, vult COMMON /READ/

SF.ES 11 output lfn=49 VL.NKPAR int SF.COVARM lfn=20 SF.COVARP lfn=21 SF. PRMTRS lfn=3 1

Estimation data input.

Number of parameters io estimate. Covariances of model.

Covariances of parameters. Info of parameters to estimate.

(10)

2. MODUUL PAREST

Het DIANA pakket is modulair opgebouwd waarbij elk moduul een bepaalde taak voor een eindige elementen berekening uitvoert. Elk moduul doet dit door aanroep van zijn eigen segmenten, ofwel programma's, waaruit het is opgebouwd.

PAREST is een moduul voor het schatten van parameters. Het is gebaseerd op het

Kalman-filter en maakt gebruik van herhaaldelijk uit te voeren eindige elementen sommen. Deze sommen worden uitgevoerd met behulp van andere modulen uit het DIANA-pakket. Omdat PAREST dus niet alleen eigen segmenten aanroept, maar ook segmenten uit andere modulen, en in feite deze modulen volledig uitvoert, kan PAREST gezien worden als een soort super moduul.

2.1 FILOS items

Moduul PAREST maakt gebruik van segmenten uit de modulen ELASSE, ELMAT, LOADS, SOLVE en NONLIN. Voor de in deze segmenten gebruikte FILOS-items wordt verwezen naar de DIANA programmeurshandleiding (Witte, 1988).

De volgende gegevens van de FILOS-file worden expliciet in moduul PAREST gebruikt. ty.naam specificatie Omschrijving

VL.EXNDIM int VL.EXNMSR int VL.EXNPNT int VL.IKITER int VL.IKMOD int VL.IKNORM rea VL.IKPAR int VL.IKSTEP int VL.KURSTP int VL.MBASTY int VL.MEL int VL.MNODL int VL.NKITER int VL.NKMOD int VL.NKPAR int VL.NKSTEP int VL.NNODE int VL.RPARSV rea

Max. number of exp. values per experiment al point.

Number of experimental data values. per time/load step.

Number of experimental points. Number of K alman it er at ion.

Counter for repeated use of Kalman matrix.

User value of norm based on the difference between calculated and exp. displacements.

Counter for (non)linear analysis to be executed, if IKPAR>O it is also the number of the parameter for which a gradient is to be calculated. Number of Kalman step.

Number of load/time step executed. Number of different basis types. Highest element number.

M. number of nodes per element. M. nmbr of Kalman iter. per step. Number of repeated use of Kalman matrix (mod. kalman).

Number of parameters to estimate. M. number of Kalman steps. Highest node number.

(11)

AR.ELMTYP int MEL)

AR.EXMPTR int

I

NSTEP ,EXNPNT,EXNDIM) AR.EXSPTR int

AR.IKPNTR int AR.KAMTRX AR.KASTEP

AR.NODBAS int MBASTY,NNODE) AR.PNTFLG int

i

EXNPNT)

SF.BASDOF lfn=16 SF.BASOR1 lfn=17 SF.BOUNDS lfn=18 SF.COVARD lfn=19 SF.COVARM lfn=20 SF.COVARP lfn=21 SF.DISPLA lfn=22 SF.EXC0 O

R

lfn= 23 SF. EXDEVI lfn= 24 SF.EXDISP lfn=25 SF.EXELC0 lfn=26 SF.EXNLD1 lfn=27 SF.PEOUTC lfn=28 SF.PEOUTP lfn=29 SF.PEOUTR lfn=30 SF.PRMTRS lfn=31 SF.SICOPY lfn=.. SF.TOTDIS lfn=22 DS.MAT lfn= 11,12 DS.NOD lfn= 13,14 DS.TYP lfn= 15

Element type numbers.

Storage pointer for exp. displ.

Starting address pointer of exp. displ. Pointer for steps and iterations.

Kalman gain matrix.

Corresponding load/time step for Kalman step.

Basisnumbers per node

Flag for (non)active exp. point. System d.0.f. numbers per basis Orientations of syst. d.0.f. per basis. Bounds for parameter values.

Covariances of exp. displacements. Covariances of model.

Covariances of parameters. Displacements of linear analysis. Coordinates of experimental points. Disp1.-par. gradients in exp. points. Exp. displ. values in exp. points. Loc. coordinates of exp. points and corresponding element numbers. Calc. displ. values of exp. points. Report of covariances for parameters. Report of parameter values.

Residual norms of executed calculations.

Info on parameters t o be estimated. Contains a copy of SF.SGMINF to make several restarts of nonlin. Displacements of nonlinear analysis. Materials data structure.

Node nmbrs and coord. per element. Element type data structure.

(12)

2.2 Segmenten

Moduul PAREST bevat de volgende segmenten: Sepment Toep ass ing.

ITCHCK To indicate a legal termination of a KALMAN iteration loop. ITNINI

KALMAN

Initialization for each iteration of a Kalman step. Calculation of new parametervalues and covariances.

NLCHCK To indicate that the next load/time step has to be executed within a Kalman iteration.

NLINIT Initialization of a (non)linear calculation, pertubation of a material parameter in accordance with IKPAR.

To execute a modified Kalman calculation for IKPAR<O, or to store the displacements for IKPAR=O, or

to calculate disp1.-par. gradients for IKPAR>O. Init ia1 dat a for parameter estimation,

evaluation of experimental points. PARDEV

PEINIT

STCHCK Indication of the next Kalman step to be executed, Print out of estimation report to output file.

Initialization for each Kalman step. STPINI

Voor een schematische weergave van de jobstream, zoals gegenereerd door de command interpreter, wordt verwezen naar appendix A.

(13)

2.2.1 Segment ITCHCK - Pad: /src/diana_paxest/pe/icïï

Segment ITCHCK ( ITeration CHeCK ) checked op beeindiging van een Kalman iteratie loop. Routines Pad: ici i /Ml ici i /lvli ici i /lvli Routine maic ovicin ovicma Doel hoofd routine initialisatie

Kalman iteratie check Besturinn:

Het JCL record voor segment ITCHCK wordt gegenereerd door aanroep van

JCLXEQ (zie progr. handleiding

,

Witte 1988, hfdst. 8) met mod segexe ircsgp = 2HPE Genevensstroom = 6HIC11

Segment stuur parameters? = -1 : Nee VL.IKITER int VL.IKNORM rea VL.IKSTEP int VL.NKITER int SF.PEOUTR lfn=30 input input input input input out put

(14)

2.2.2 Segment ITNINI

-

Pad: /src/&ana_parest/pe/iiíí

Segment ITNINI ( ITeratioN INItialisation ) zorgt voor de initialisatie bij elke iteratie van een Kalman stap.

Routines

Pad: Routine

ii i i

/

lvl i

iii i /lvli oviiin iii i /Mi oviima

maii

Doel

hoofd routine initialisatie

reset van filos items Besturing:

Het JCL record voor segment ITNINI wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod = 2HPE segexe = 6HII11

ircsgp Segment stuur parameters? =

-1

: Nee

Gegevenss t room

VLIKITER int input

VL.IKMOD int output

VL.IKPAR int output

VL.IKSTEP int input VL.NKITER int input VL.NKMOD int input

(15)

2.2.3 Segment KALMAN - Pad: /src/diana_paxest/pe/kalî

Segment KALMAN berekent voor elke Kalman iteratie voor elke Kalman stap nieuwe parameters en covarianties.

Routines Pad: ka1 l/lvll ka1 1 / h l l kall/lvll ka1 1 /lv12 kall/lvl2 ka1 1 /lv12 ka1 l/lv12 ka1 1 / h l 2 kall/lvl2 Routine maka ovkain ovkama kprecp kalman pldlt kcorcp kcorpm putpar Doel hoofd routine initialisatie Kalman evaluatie

predictor value van par. covariantie berekening van Kalman gain matrix matrix inversie

corrector value van par. covariantie nieuwe mat. parameters

nieuwe mat. par. op DS.MAT schrijven Besturing :

Het JCL record voor segment KALMAN wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod segexe ircsgp

= 2HPE = 6HKA11

Segment stuur parameters? = -1 : Nee Gegevensst room VL.EXNDIM int VL.EXNMSR int VL.EXNPNT int VL.IKITER int VL.IKSTEP int VL.NKITER int VL.NKMOD int VL.NKPAR int VL.NKSTEP int AR.EXMPTR int AR.EXSPTR int AR.IKNPTR int AR.KAMTRX rea input input input input input input input input input input input

input out put output

(16)

AR.KASTEP int AR.PNTFLG int SF.BOUNDS lfn=18 SF.COVARD lfn=19 SF.COVARM lfn=20 SF.COVARP lfn=21 SF.EXDEV1 lfn=24 SF.EXDISP lfn=25 SF.EXNLD1 lfn=27

SF.PEOUTC lfn=28 out put SF.PEOUTP lfn=29 output SF.PEOUTR lfn=30 output SF.PRMTRS lfn= 3 1 input

DS.MAT lfn=11,12 input output input input input input input input output input input input

(17)

2.2.4 Segment NLCHCK - Pad: /src/ûïana_paxest/pe/ncíí

Segment NLCHCK ( NonLin CHeCK ) checkt of voor een niet lineaire som de belastingstap behorende bij de Kalman stap uitgevoerd is.

Routines Pad: ncii/lvll ncíi/lvll ncll/lvll Routine manc ovncin ovncma Doel hoofd routine initialisatie

Kalman step check Besturing:

Het JCL record voor segment NLCHCK wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod segexe ircsgp

= 2HPE = 6HNC11

Segment stuur parameters? =

-1

: Nee Gegevenss t room VL.IKSTEP int VL.KURSTP int AR.KASTEP int AR.PNTFLG int input input input input

(18)

2.2.5 Segment

NLINIT

-

Pad:

/src/dianajarest/pe/niíi

Segment NLINIT ( NonLin INITialisation ) voert de benodigde initialisaties voor een (niet) lineaire berekening uit. Bovendien vindt pertubatie van een parameter plaats,

overeenkomstig de waarde van IKPAR. Routines Pad: nil i /lvll niil/lvll ni 1 1 /lvll niii/lvli ni i i /lv12 niil/lvl2 Routine mani ovniin ovnima sgbgcp kprmod putpar Doel hoofd routine ini tialisat ie

originele data terug op SF.SGMINF. par amet er pertubat ie

par amet er pertubat ie

nieuwe parameters op DS.MAT Besturing:

Het JCL record voor segment NLINIT wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod segexe ircsgp

= 2HPE

= 6 ~ ~ 1 1 1

Segment stuur parameters? = -1 : Nee Gegevensstroom VL.IKMOD int VL.IKPAR int VL.MEL int VL.NKMOD int VL.NKPAR int SF.PRMTRS lfn=31 SFSICOPY lfn=.. DS.MAT lfn= 11,12 input input input input input input input input output

(19)

2.2.6 Segment PARDEV - Pak /src/diana_paxest/p/pdll

Segment PARDEV ( PARtial DEriVatives ) voert afhankelijk van IKPAR de volgende taken uit:

IKPAR

<

O IKPAR = O IKPAR

>

O

Modified Kalman berekening

opslaan van berekende verplaatsingen berekening van disp1.-param. afgeleiden. Routines Pad: pdl 1 /lvll pdll/lvll pd 1 1

/

lvl 1 pdll/lvl2 pdl l/lv13 pdll/pdllib pdll/lvl2 p d l l /lv13 p d l l /lv12 p d l l /lv12 pdl 1 /lv13 Routine mapd ovpdin ovpdma pddisp pdnwpa pdstdi pddevi pddisw putpar putpar

...

Doel hoofd routine ini t ialisat ie verpl./afgeleiden/param. berekening verplaat singen element switch element afh. routines modified Kalman

nieuwe mat. par. op DS-MAT schrijven opslaan van verplaatsingen

partiele afgeleiden

mat. par. op DS.MAT schrijven Besturing:

Het JCL record voor segment PARDEV wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod segexe ircsgp jclpar( i) Gegevensst room = 2HPE = 6HPD11

Segment stuur parameters? = -1 : Nee VL.EXNDIM int VL.EXNMSR int VL.EXNPNT int VL.IKITER int VL.IKMOD int VL.IKNORM rea VL.IKPAR int = O : nonlinear analysis = 1 : linear analysis input input input input input input input output

(20)

VL.IKSTEP int VL.MBASTY int VL.MEL int VL.MNODL int VL.NKITER int VL.NKMOD int VL.NKPAR int VL.NKSTEP int VL.NNODE int VL.RPARSV rea AR.ELMTYP int AR.EXMPTR int AR.EXSPTR int AR.IKNPTR int AR.KAMTRX rea AR.KASTEP int AR.NQDBAS int AR.PNTFLG int SF.BASDOF lfn=16 SF.BASOR1 lfn=17 SF.BOUNDS lfn=18 SF.DISPLA lfn=22 SF.EXDEV1 lfn=24 SF. EXD IS P lfn= 25 SF. EXELCO lfn= 26 SF. EXNLDI lfn= 27 SF.PEOUTP lfn=29 SF.PEOUTR lfn=30 SF.PRMTRS lfn=31 SF.TOTDIS lfn=22 input input input input input input input input input input input input input input input input input input input input input input input input input input input out put out put out put

DS.MAT lfn=11,12 input output DS.NOD lfn=13,14 input output

(21)

2.2.7 Segment PEINIT - Pad: /src/ûiana_parest/pe/piil

Segment PEINIT ( ParEst INITialisation ) zorgt voor de initialisatie bij een PAREST berekening. Bovendien worden de exp. punten geevalueerd, d.w.z. er wordt gekeken in welk element ze liggen en lokale coordinaten worden berekend.

Routines Pad: pil i /lvll pil i/lvlí piií/lvli piil/lvll pi 1 i /lv12 pil i/lv13 pil i/lv13 pil l/pillib Routine mapi ovpiin ovpima sgbgcb msrpnt exldst xelcor

...

Doel hoofd routine initialisatie

save inhoud van SF.SGMINF evaluatie van exp. punten evaluatie van exp. punten dicht sbijzi jnd element lokale coordinaten element afh. routines Besturing:

Het JCL record voor segment PEINIT wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod = 2HPE segexe = 6HPI11

ircsgp Segment stuur parameters? = -1 : Nee

Recordnummer op SF.SGMINF met Kalman step data Recordnummer op SF.SGMINF met IKNORM value = NKITER = NKMOD Genevensstroom VL.EXNMSR int VL.EXNPNT int VL.IKNORM rea VL.IKSTEP int VL.MEL int VL.MNODL int VL.NKITER int VL.NKMOD int VL.NKPAR int VL.NKSTEP int VL.NNODE int input input input input input input output out put out put output out put

(22)

AR.ELMTYP int AR.IKNPTR int AR.KAMTRX rea AR.KASTEP int AR.PNTFLG int SF.EXCOOR lfn=23 SF.EXDEV1 lfn=24 SF.EXELC0 lfn=26 SF.EXNLD1 lfn=27 SF.PEOUTC lfn=28 SF.PEOUTP lfn=29 SF.PEOUTR lfn=30 SF. SIC OPY DS.NOD lfn=13,14 input DS.TYP lfn= 15 input lfn=

. .

input output out put output output output out put out put out put out put output input input

(23)

2.2.8 Segment STCHCK -

P d

/src/ûïana_pa.rest/pe/sci i

Segment STCHCK ( STep CHeCK ) kijkt of een volgende Kalman stap uitgevoerd moet worden. Zo niet dan wordt het programma beeindigt met een print out van een estimation report naar de uitvoerfile.

Routines Pad: Routine scl l/lvll masc scl l/lvll ovscin scl l/lvll ovscma scl l/lv12 peout Doel hoofd routine initialisatie

volgende Kalman stap

/

estimation report estimation report

Besturing:

Het JCL record voor segment STCHCK wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod = 2HPE segexe = 6HSC11

ircsgp Segment stuur parameters? = -1 : Nee Genevensstroom VL.IKSTEP VL.NKITER VL.NKMOD VL.NKPAR VL.NKSTEP AR.IKPNTR AR.KASTEP SF.PEOUTC SF.PEOUTP SF.PEOUTR SF.PRMTRS int int int int int int int lfn=28 lfn=29 lfn=30 lfn=31 input output input input input input input input input input input input

(24)

2.2.9 Segment

STPINI

-

Pad:

/src/ûiana_paxest/pe/siii

Segment STPINI ( STeP INItialisation ) initialiseert voor elke Kalman stap Routines Pad: si i i / h l i si i i /lvli sii i /lvli Routine masi ovsiin ovsima Doel hoofd routine ini t ialisat ie reset filos items Besturing :

Het JCL record voor segment STPINI wordt gegenereerd door aanroep van JCLXEQ (zie progr. handleiding

,

mod segexe ircsgp = 2HPE Gegevensst room = 6HSI11

Segment stuur parameters? =

-1

: Nee VL.IKITER int VL.IKSTEP int VL.NKSTEP int output input input

(25)

LITERATUUR

Hendriks, M., ter publikatie, T.U. Eindhoven, 1988.

Hendriks, M., Moduul PAREST; gebruikershandleiding, T.U. Eindhoven, WFW rapport nr. 88461,1988. Oomens, C., ter publikatie,

T.U. Eindhoven, 1988.

Witte, F. de, Handleidingen DIANA, deel iii, gebruikershandleiding, deel iv, programmeurs handleiding,

T.N.0

- I.B.B.C., Rijswijk, 1988.

(26)

{ segment STPINI }

for IKITER = 1 until NKITER do begin

{ segment ITNINI }

for IKPAR = IKBGN until IKEND do

APPENDIX A

Hieronder volgt een schematische weergave van de jobstream voor moduul PAREST, zoals gegenereerd door de command interpret or.

NKSTEP NKITER IKBGN IKEND

._

.-

._

.-

._

.-

begin

max. number of Kalman steps }

max. number of Kalman iterations per step }

equal zero for no modified Kalman, equal -NKMOD for modified Kalman,

NKMOD = max. number of modified Kalman executions per iteration }

NKPAR = number of parameters to estimate } { equal NKPAR,

segment NLINIT }

ELASSE segments ELMAT segments LOAD segments }

if {linear analysis} then

I

{SOLVE segments }

else end

{ segment PARDEV }

{ NONLIN segments

+

segment NLCHCK }

" i

I

end

r

segment ITCHCK segment KALMAN } }

end