Het programmapakket MARC en geometrisch niet-lineaire berekeningen

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

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.

Het prograninriapakket

WFV 56.030

MARC

en geometrisch niet-lineaire berekeningen.

P.A.A. van Hoogstraten

Deel A van het verslag van stage 4k805

juli 1986

-0.2- o. I 0.2 O. 3 1. I 2. I 2. I 2. P 6.1 4.c 4.1 4.6 4*? 4. I 5.1 5.1 5.

g

P

c 52

e

n'

P

N

[ S I

I -

1

e-

1

I - ]

I-

1

cp

62.9.1.)

(2. I . S.>

IF

=

IR-

(2.1. IO)

-

2.5-

4'1

2

=x,

42.2. S.>

(3.2.1)

f

-3.4- (3.3.1.9 (3.3.2. C3.33.)

Q

3.3.

e.>

43.3.5.)

3 Ql

(x,-x,p)-

-

-

24

-

A,

+

1-23

toe

: -Y

x,

=

A,

=

A,

Q

3.3. d. ) 4 3.3.7. ) ( 3 . 3 . 8 . ) (3.3. I f . )

(3.3. i%)

(3.3. re)

(3.3.15)

-

4.3-

4.3.

t

Y

x

3

o

_{@ &}

-4.9-

krek belaskinq -

-

bereùenincf

-5.2- LhCTIC!TEITXZODdLUC 2.E5 'Es1 NDyAf6T)ûûRcNEDE i. EGINLWTE I . iUARSCONTRACTIE~AICIENT 8.3 .EES FILE ( = I ) 0. 'HEDRETICCH TOTAL (=I 1

'HEMIETISCY UPDATED (=i ) 0.

LASTIClrEITcMoDilLuC 2.E5

EG1WiAET)WRSNEüE I.

EGIVLENGTE l I

W ~ ~ O Y T M C T I E ~ F F ICIENT Q z 3

EEC FILE <=I> 0. ~

HES?ET!SCli TOTAL ( = l ) 0. HEORETICCH UPDATED (=i ) Q.

~~

Fig. 5.3.

~ A S T I C I ~ E I T S ~ D & ~ 2 .E

EG!VDWARSDO~SNEDE I .

EGINLENGTE I .

vlA!EZM\ITRAfT I ECMFF IC1 ENT .3 EEC FILE ( = I > a.

HEORETECH TOTAL ( = I 1 @.

HECRETISCH UPDATED (=i > 0 =

1-1

I

I

--

m 1.18 1.B

Fig.5.u. :

-5.I-

AXIALE KRACHT (TOTAL , V . V E-I

* 40 . i 0 - .m -.ia -.a -.a - ~ 4 8 I70 .sa Í.m Í . i a 1.23 1.3 LACTICITEITSMODuLüS 2.E5 EE1NDWARSI)~CNEE 1 . EE? NLENGTE 1.

W AECONTRACT I EMIEFFICI ENT a. 3

EEC FILE (=I ) 0.

HEDRETISCH TUTAL (=i 1 a.

HEOF?ET!XX UPDATED (=1) 0.

L-cJ

5.5. ..AÇTICiTEITSMOD'uC 2 .E5 ZINDUAWOClffSNEDE i . EGIVLENGÏE 1 . W&WONTRACT iEC!EFF!CIENT 8.3 EEC FILE ( = I ) 0. HEORETISCH TOTAL ( = I ) 0. HEOREÏISCH UPDATED (-1 ) 0.

AXIALE %AWT (UFDATED,V.S> i-I ... _AST!C!TEIWlQDULUC 2 . E5 3IiUDYARCDDOESNEDE I . IGINLEffiTE 1 . IARSCONTRACTIEMIEFF ICIEVT 8.3 ?ES FILE ( = I > 0* iEDRETICCH TOTAL ( = 1 ) 0- 4EORETI”H WDATED (=í ) 0.

-5.6.- 1 3 % - 'm -1 .m LASTICITEITÇNODWUC 2.E EF!NDWAEDLMRCNEDE EGINLENGTE WA~DNTRACTIECMFF EIENT 0.3

EEC FILE (=i)

HEDRETISCH UPDATED (=1) @.

I

LASTIC?TE1?WYiLi& 2 .E5 EGINDUARSDDOKVEDE i . EG1 NLENGTE I . WARSCONTRACTIECMFFICIENT 0.3 EES FILE ( = i > a. HEMIETISCH TDTK : = I ) a. HEORETISCH uWATED (=1 ) 0.

-6.1-

e =

E A A ,

A

=

A,

A,

-

2v

A,

- -

N

-

A;"" EA0

HEMZETIÇCH TOTAL (=1 > 8 .

HEORETISM UPDATED (=I ) B.

HEORETIÇCH UPDATED (=i )

_ASTICITEITCMODULdS 2 .E5 3INDidARSDDJRSNEDE 1. :GINLENGTE 1 . IARSCONTRACTIECOEFFICIENT 0.3 fEC FILE (=1) 0. E O R E T I W TOTAL (=! ) 8. IEORETICCH UPDATED (=I ) 8 .

O

Bijlage E : subroutine :TEKEN, file GRFILE

-

R.i

-

TITLE MARC.TREKSTAAF.VRU

TITLE voorgeschreven belasting

T 1 TLE

TITLE elementtype 9

TITLE

COMMENT parameter cards

updated lagrangian formulering

. . .

SIZING,10000, : w f A r u ; m L c

ELEMENTS

,

9 : elemefit keuze

ALL POINTS : n.v.L ; zie uoorbeelJ 3

LARGE DISP : bevelrehing met y o & v e v p l a 4 ~ ~ ; f i g @ n

UPDATE END

COMMENT model definition cards

C ~ N N E C T ~ ~ ~ T Y

'

b e t d e s C Q Y d : Y/JCfGted /'gvqngi4n fOOTtVLllC?Ylh9

zolider : to+,/ ,, I, 1,9,1,2i COORDINATES 1,0.0,0.0,0.Q, 2,1.0,0.0,0.0, GEOMETRY 1 1 1

.o,

1 CONTROL 5013, 0.Q1,0.50,0.10,0.05,0,0.2 POST 0, r r Q i

1

i POINT LOAD 1, 1000

,

2 , PROPERTY

1 ,

2E5, O. 3 i i , ~ O U N ~ A R Y C O ~ ~ ~ T I O N S 2, O.O,Q.O,O.O, 1 TO 3

TITLE

COMMENT parameter cards

SIZINC,l0000, ELEMENTS,S, ALL POINTS LARGE DISP UPDATE, END

COMNIENT model definition cards

CONNECTIVITY

. . .

1,0.0,0.0,0.0, 2,1.0,0.0,0.0r 1 , I

.o,

1 CONTROL 50r3i 0.01,0.50,0.10,0.05,0,0.2 POST Or i r0,lt PROPERTY 1, 2E5 i O. 3

,

1 ,

BOUNDARY CONDITIONS 21 0 . 0 , 0 . 0 , 0 . 0 , GEOMETRY 1 TO 3 1, 0.005,0.0i0.0, 1 to 3 2, END OPTION

~ O ~ M ~ N T historv definition cards AUTO LOAD

50 1

PROPORTIONAL I N C ~ E M E N T I,l.O,

-

A.3-

TITLE MARC.TREKSTAAF.VEU

TITLE voorgeschreven belasting

TITLE updated lagrangian formulering

TITLE elementtype 21

TITLE

COMMENT parameter cards

SIZING,30000, ELEMENTS,21, ALL POINTS LARGE DISP UPDATE END

C O ~ M E N T model definition cards

CONNECTIVITY

. . .

1,21

,

1,2 3

,

4

,

5

,

6 7

,

8

,

9, 10, 11,12

,

13 14, 15,16,17,18,19,20, : 2 ' ~ keoopp6fig .& 20 COORDINATES : k n o o p p w L coördihqtefi 1,0.0,0.0,0.0, 2,1.0,0.0,0.0, 3,1.0,1.0,0.0, 4,0.0,1.0,0.0, 5,Q.0,0.0,1.0, 6,1.0,0.0,2.0, 711.011.0r'l.0, 8,0.0,2.0,1.0, 9,0.5,0.0,0.0, 10,1.0,0.5,0.0, 11,0.5,1.0,0.0, 12,Q.0,0.5,0.OI 13,0.5,0.0,1.0, 14,1.0,0.5,1.0, 15,0.5,1.0,1.0, 16,0.0,0.5,1.0, 17,0.0,0.0,0.5, 18,1.0,0.0,0.5, 19,1.0,1.0,0.5, 20,0.0,1.0,0.5, CONTROL 5 0 1 3 , 0.01,0.50,Q.10,0.05,~,0.2 POST 0, i I Q , 1 1

1, - 1 . OE3 PROPERTY 1 1 2E5 ,O. 3 I 8,16,20, knooppunt 8,16,20 /I /, END OPTION

T I T L E ~ ~ ~ ~ . ~ R E ~ S T ~ A F . W U

T I T L E voorgesclireven verplaatsing

T I T L E updated lagrangian formulering

T I T L E elementtype 21

T I T L E

COMMENT parameter cards

SIZING,30000,

ALL POINTS

LARGE DISF

UPDATE END

COMMENT model definition cards

COMNECTIVITY

. . .

EIJMENTS I 2 I l,O.O,O.O,O.Qr 2,4.0,0.0,0.0, 3,1.0r1.0,0-O, 4 , 0 . 0 , 4 .O,O.O, 5,0.O,O.O,l.Or 6,1.0,0.0,1.0, 7,1.0,1.0~1.0, 8,0.0,1.0,1.0, 9,0.5,0.0,0.0, 10,1.0,0.5,0.0, 11,0.5,1.0,0.0, 12,0.0,0.5,0.0, 13,0.5,0.0r1.0r 14,1.0,0.5,1.0, 15,0.5,1.0,1.0, I 6 , O . O , G . 5 , l . O , 17,0.0,0.0,0.5, 18,1.0,0.0,0.5, 19,1.0,1.0,0.5, 20,0.0,1.0,0.5, 5 0 , 3 , 0 . 0 1 ~ 0 . 5 0 ~ 0 . 1 0 ~ 0 ~ 0 5 ~ 0 ~ 0 ~ 2 POST CONTROL Or i iOt1 r

183' 1284, 0.0,o.

o,

112, 0 . 0 ,

1 ,

8,l6,ZO, 0.005, 4 1 2,3,6,7,18,14,18,19, EMD OPTION C O M ~ ~ ~ T h i s t o r y definition cards AUTO LOAD 50 f P R O F O R ~ I O ~ ~ ~ I ~ C ~ E M ~ ~ T l I 1 . O I C O ~ T I ~ U ~ 5 , 4 7 1

- B i -

Het aanroepen van het prosrammapakket Marc

Bij het gebruik van het programmapakket Narc op de Prime F750 afdelings- computer, moet de invoerfile op 'POOL' staan. Door het intoetsen van 'MARC' kan de rekenrun worden gestart. Na dit commando gebeurt het volgende:

MARC

Info about MARC can be found in MARC)DOC.

( 1 ) Inputfile: xxxxx

(2) Outputfile (xxxxx.MRCOUT): xxxxx.MRCOUT

( 3 ) Queue (NIGNT): y w v y

(4) Marc version (Kl):

<cr>

[JOB rev 1 9 . 4 1

( 5 ) Your job, #jobnr, was submitted to queue yyyyy.

Bome=<RLAD>POO&>

. . .

(Het onderstreepte moet door de gebruiker worden ingetoetst)

Er wordt een invoerfile gevraagd.

Er wordt een naam voor de uitvoerfile ( xxxxx.mrcout )

voorgesteld. Deze kan worden geaccepteerd d.m.v. een return of er kan een andere naam worden ingevoerd.

Er wordt gevraagd in welke queue de berekening moet worden geplaatst. Be mogelijkheden zijn :

FAST : Wanneer niet te veel andere rekenjobs draaien, voert Marc kleine berekeningen direct uit. MARC! : Idem voor middellange berekeningen

NIGNT: Voor grote rekenpartijen. Deze worden ' s nachts uitgevoerd.

Er wordt gevraagd met welke versie van MARC moet worden gerekend. D.m.v. een <cr> wordt de voorgestelde versie (Kl) geaccepteerd. Momenteel is dit de enig beschikbare versie.

Als alle voorgaande vragen correct zijn beantwoord volgt tenslotte de melding ( 5 ) . De berekening wordt als j o b

worden bij het verwerken van de zg. POST16-files tot files van een bruikbaar formaat. Het starten van de filetranslator gebeurt d.m.v. het intoetsen van het commando 'FT':

(Het onderstreepte moet door de gebruiker worden ingetoetst.)

-

FT

ENTER TRANSLATOR#

M - MARC-POST TO SDRC-UNIVERSAL

F - FAST, SIMPLE STRIPPED VERSION OF M

S - SDR~-UN~VERSAL TO MARC-INPUT

P - NARC-POST TO NARC-POST

(BINARY TO FORMATTED AND V.V. )

I

-

IDEAS-UNIVERSAL TO MARC-INPUT Q - QUIT

ENTER TRANS LA TOR^ &i

: start van de filetranslator

: kies voor filetranslator M MARC-POST to UNIVERSAL translator, version 2.2 feb. ' 8 4

if you don't know what t o do, type M (for menu)

or : for help (if there is any)

INPUTFILE (POSTIS) : <CT>

HEADER

C O ~ N E C ~ I V ~ T Y ~ O O R D ~ N A T E ~ LOADSTEP

( 1 )

: <CI>

LOADSTEP = I TINE = 0.000 INC =

NODAL DATA

BINARY FILE (Y/N) ( U ) : < C S >

CREATED SET 1 CREATED SET 2 ENTER ~ O M ~ A N D : &i RM

-

READJIARC DC - ~ ~ ~ A - S ~ L E C T I O N W U

-

~ R ~ T ~ - U N I V E R S A L Q

-

QUIT EMTER COPINAND: READ-MARC ENTER ~ O P I N ~ ~ D : ~ NF

-

MEW-FILE US

-

USE-LOADSTEP W

-

WATCH ST

-

STATUS : POST16 is de invoerfile

,

: dit is een binaire file

: begin bij belastingsincrement 1

O FREQUENCY = 0.000

O SUBIMC =

: nu wordt het menu getoond

-

c.2- RE AD-PI ARC ENTER COMMAND: READ-NARC,WATCH ENTER COMMAND:

a

CR

-

CREATE GE - GENERATE SP

-

SPECIAL LI - LIST WR

-

WRITE ST

-

STATUS R~~D-~ARC,WATCH ENTER C O M ~ ~ N D : CR

VARIABLE TYPE (END-OF-LIST):M

E

-

END-OF-LIST T - TIME

I - INC~EMENT

F - FREQUENCY

EN - ELEPIE~TAL DATA AT NODE

EI - DATA AT INT-POINT N

-

MODAL DATA

SI?

-

SPARE

s

-

s m m c

VARIABLE TYPE (ENDOF-LIST):@ N O D E L A ~ ~ L ( 1 ) :

1

NODAL DATA- TEM : ( C K >

o :

1 : 2 : 3 : 4 : 5 : 6 : NONE FIRST DISFLACE~ENT C O M F O ~ ~ N T SECOEJU D I ~ ~ L A C ~ ~ E N T C O M ~ O N ~ N T THIRD DIS~LACEPIENT C O M P O ~ E ~ T FIRST REACTION CON~ONENT

SECOND REACTION ~ O ~ F O N E N T THIRD REACTION COMPON~NT enter a number ( 1 ) : j,

: kies hieruit WATCH

: nu wordt eexst de gewenste data

: geselecteerd; in dit geval : knooppuntsverplaatsingen en : knooppuntskrachten

.

( op deze wijze wordt, alle gewenste data

: alleen knooppuntsdaka nodig!

achtereenvolgens geselecteerd 1

met ST wordt gecontroleerd of

-NODAL DATA, NODE 2 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE 3 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE 4 FIRST REACTION COMPONENT -NODAL DATA, NODE 5 FIRST REACTION COMPONENT -NODAL DATA, NODE 6 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE 7 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE 8 FIRST REACTION C O M P O N ~ ~ T -NODAL DATA, NODE 10 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE I 2 FIRST REACTION COMPONENT -NODAL DATA, NODE 1 4 FIRST DISPLACEMENT COMPONENT -NODAL DATA, NODE 16 FIRST REACTION COMPONE~T -NODAL DATA, NODE 17 FIRST REACTION COMPONENT

-NODAL DATA, NODE 18 FIRST DISPLAC~MENT C O ~ P ~ N E N T

-NODAL DATA, NODE 19 FIRST ~ I S ~ L A C ~ M E N T COMPONENT : deze serie is gebruikt

-NOBAL DATA, BODE 20 FIRST REACTION COMPONENT : voor eiementtype 24

R E ~ D - ~ A ~ C , W A T C ~

ENTER C O ~ ~ A N D : SE : tenslotte wordt de file met

: de gekozen data gegenereerd,

LOADSTEP ( 1 ) : A

ENTER COMMAND:

O ~ T P U ~ F I L E (HISTORY) : (czewenste naam)

NEW FILE OM (Y/N) ( Y ) : < c ~ >

~ E A D - M A R C , ~ A T C ~

EN TE^ C O ~ A N D : &

: en weggeschreven.

: de filetranslator wordt nu

-

D.1-

C PROGRAMMA NIDDEL: STAGE1

c

_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _

C

C TREKSTAAF 21 EN OPTELLEN VAN ACHT REACTIEKRACHTEN AAN IN~LEMMING C VAM DIT ELEMENT.

MJDDELEN VAN ACHT KNOOPPUNTSVERPLAATSINGEN VAN UITEINDE ELEMENT-

C C INVOER : C C UITVOER : C C

c

C

e

REAL* 4 REAL" 4

FORMAAT MAX. 16*500 (POSTl6-FILE VAN MARCUITVOER M.B.V. FILETRANSLATOR M TOT DIT FORMAAT BEWERKEN).

FILE MMM

,

FORMAAT MAX. 2*500

EERSTE KOLOM BEVAT 50171 VAN DE ~EACTIEKRACHTEN (INKL~MMING) TWEEDE KOLOM BEVAT GEMIDDELDE KNOOPPUNTSVERPLAATSINGEN

IN XI-XI VAN HET VRIJE UITEINDE.

€i( 16,500) F(2,500) INTEGER*4 FNAME(20) ,NH, I CAVE H,NH C C $INSERT SYSCO~>A$KEYS

10 CALL TNOWA ('GEEF READ (1,201 FNANE 20 FORMAT (20A4) IF ( .MOT. OPEN$A NH= 1 ~ ~ L E N A A M :',INTC(16)) (A$WEADrFNA~E,INTS(80),INTS49))) GOTO 10

30 READ (13,40,ENB=50,ER~=50) (H(IfNH)fI=lrj6)

F f

1

,

NH) = H ( 1 NH) +H ( 4

,

NH) +€I( 5

,

MH) +H ( 8 i NK) +Hf 10, NH) +H( 12 NH) + F(2,NH) = ( H ( 2 , N H ) t M ( 3 , N H ) t H ( 6 , N H ) + H ~ ? ~ N H ~ + H ( 9 , N H ) + H ( ~ ~ ~ N ~ ~ ~ NH=NM-t 1 GOTO 30 ~H(l3~N€i)tH(16~~H) &€i(14,NH)+H(15,NH))la.0 40 FORRAT (6El3.5) 50 CALL ~ ~ (A$WRITr ' ~ ~ ~ M 1 , % N T S ( 3 ) , I N T C ( I O ) ~ ~ ~ $ A NH =NM-I DO 100 I = 1,NH WRITE (14g40) F(lrI),F(2,I) 100 CONTINUE

CALL CLOS$A (INTSIlO)) CALL CLOS$A (INTS(9)) CALL EXIT

C C C C C C C C C C C C

COMMON/ IFNAMES/ IN(20), IHELP(20)

DATA IN / 'WFW)', 'HOOG', ' S T ) S ' , 'TAGE', 'I>RE', 'WERK', END

& 'ING)', 'GRFI', 'LE ' , 11"' ' / , IHELP f 20"' ' i

S U B R O U ~ I N ~ FUNCT : ~ E W ~ R K T EEN FILE MET BEREKENDE WAARDEN

EN PLOT DEZE AL DAN NIET MET HET T ~ ~ ~ R E ~ I S C H E RES~LTAAT

DE NAAM FUNCT IS VERP~ICHT ! .

_ _ _ _ _ - _ - - -

S U B ~ O U ~ I ~ ~ ~ FUNCT ( X , Y, N, OX, GEG, IY) LOGICAL OK

INTEGER"4 N, IY

REAL*4 X(I), Y(1), GEGtl)

B ~ T E ~ E N I S VAN DE LOKALE V ~ R I ~ B ~ L E N : C C GEG C C C C C C C

c

IY C C C C C C C C C

c

C C C

-IN- BEVAT G ~ G ~ V E N S VOOR ~EREKENING VAN DE KROMME GEGfIf = DE ELASTICITEI~SMODULU~

GEGf2) = DE B~GINDWARSD~ORSN~DE

GEG(3) = DE B E G ~ N ~ E N G T E

GEG(4) = BE DWARSCONTRACTIE~ODUL~S GEG(5) = LEES FILE, VRAAG NAAM

GEGC6) = TOTAL LAGR. GEG(7) = UPDATED LAGR. GETEKEND MOET WORDEN. IY=I: VE~LE~GIMGSFACTOR -

-IN- KEUZE PARA ME TE^ OM ABM TE GEVEN WELKE GRbnFIEK

AXIALE KRACHT BIJ UPDATED ~ ~ ~ C H R I J V ~ N G EN V ~ O R ~ E S C H R E V E ~ BELASTING IY=2: V~RLE~GINGSFACTO~ -

AXIALE KRACHT BIJ UPDATED BESCHRIJVING

EN VOORGESC~REVEN V E R ~ ~ h A T S I N ~ IY=3: VERLENGINGSFACTOR

-

AXIALE KRACHT BIJ TOTAL BESCHRIJVING EN VOORGESC~REV~N BELASTING

IY=4: VERLENGINGSFACTOR

-

AXIALE KRACHT BIJ TOTAL BESCHRIJVING

-

E . 2 -

C N -OUT- AANTAL PUNTEN IN GRAFIEK

C! OX -OUT- .TRUE. OM PLAATJE TE TEKENEN C X -OUT- BEVAT X-COORDINATEN

C P -OUT- BEVAT Y-COORDINATEN

c

C

____________________---

C DE FILE MET DE ~EREKENINGSRESULTATEN MAG MAXIMAAL HET FORMAAT

C 2

*

500 HEBBEN.

C DE EERSTE KOLOM MOET BE RESULTERENDE REAC~~EKRACHT AAN DE C INKLEMMING VAN DE TREKSTAA~ BEVATTEN. ( = - I * BELASTING )

C DE TWEEDE KOLOM MOET DE GEMIDDELDE KNOOPPUNTSVERPLAATSING IN XI-RI

C VAN HET VRIJE UITEINDE BEVATTEN.

c . . .

C

c

LOGICAL LEES,TOTAL,U~DATE REAL*4 E, NU, A, LO, H(2,500) SAVE H,NH C $ ~ ~ SU~&OM>A$K~YS S ~ ~ T INTEGER*4 FNAME (2018 NHf

1

10 20 30 40 50 60 70 E = GEG(1) A = GEGt21 L O = GEG(3) NU = GEGf4)

LEES = (GEG(5) .EQ. 1.01

GEGCS) = 0.0

TOTAL = (GEG(6) .EQ. 1.0) UPDATE = (GEG(7) .EQ. 1.0)

IF (.NOT. LEES) SOT0 60

CALL TNOUA ('GEEF ~ILENAAM : I ,INTS(16))

READ (1,20) F N ~ ~ ~

FORHAT 620A4)

IF ( .NOT. OPEN$A ( A$READ,FNAME,I~TS(80),INT~(9))) GOT8 10

NH= 'I

READ (13,40,EMD=50rERR=50) K(l,NH), H(2,NH)

NH = NH 4- 1

GOT0 30

FORMAT (2E13.5)

CALL CLOS$A ( INTS(9))

NH = NH - 1 DO 70 I =

1 ,

NH

X(1) = H ( 2 , I ) / LO 4- 1,o

YfI) = H(1,I) / (-1.O"E) / A

IF (TOTAL) Y(I)=O.5*X(I)*((X(I)**2)-1) IF (UPDATE) U(I)=ALOG(X(I))*(X(I)**(-2*NU)) CONTINUE R=NH OK = .TRUE. RETURN END

GEGEVENS ELASTICITEITS~ODULUS BEGINDWARSDOORSNEDE BEGINLENGTE DWARSCONTRACTIECOEFFPCIENT LEES FILE ( = I ) THEORETISCH TOTAL (=1) THEORETISCH UPDATED (=1) 2E5 1.0 1.0 0 . 3 1 O O

AXIALE KRACHT (UPDATED, V. B) [

-3

V E R L E N G ~ N G S F A C T ~ R

[-I

AXIALE KRACHT (UPDATED,V.V)

[-I

V E ~ ~ ~ N G I N G ~ F A C T O R

[-I

AXIALE KRACHT (TOTAL ,V.B

[-I

V ~ R ~ ~ ~ G

[-I

~ ~ G ~ ~ ~ C T ~ ~

AXIALE KRACHT (TOTAL , V . V [ - ]

~ ~ ~ L ~ ~ G ~

6-3

~ G ~ F ~ C T ~ R

0 0 0 0 0 0

-F. I C C C C C C C C

c

C C

c

c

C

PROGRAMMA OgKEER: STAGE1

OWKEREN VAN INVOERFILES VOOR HET PROGRAMMA TEKEN INVOER VOOR OMKEER: MAX 2

*

500

___-__---__---

EERSTE KOLOM: RESULTERENDE REACTIEKRACHT INKLEMMING TWEEDE KOLOM: GEMIDDELDE K N O O ~ P U N T S V € R ~ L A A T ~ I N G ~ N

IN XI-RI VAM VRIJE UITEINDE.

UITVOER : VOLGORDE (VAN INCREMENTEM) OMGEDRAAID, NAAM: I'TTT"

REAL*4 H(2,500)

IWTEGER*4 FMAME(20) ,NH, I, J

SAVE H,NH

$INSERT SYSCûM>A$KEYS

I O CALL TNOUA ('GEEF F~LENAAM :',INTS(16)) 20 FORMAT (20A4)

READ (1,20) NA^^

IF ( .NOT. OPEN$A ~ A $ R E A D , F N A M E ~ I ~ T S ~ 8 0 ) ~ 1 ~ T S ~ g ) ) ) GOT0 10

NH= 1

NH=NWt 4

GOT0 30

30 READ (13,40,END=50,ERR=50) H(l,NH), H(2,NH)

40 FORMAT (2E13.5)

50 CALL OPEN$A (A$WRIT,'TTT',INTS(3),fNTS(10))

NW = NH-1

DO 100 I = 1,NH

J=NHti-f

WRITE (14,401 M(I,J),H(2,J) 100 CONTINUE

CALL CLOS$A {INTSCIO)) -CAL% CLOS$& (INTS(9)E

CALL EXIT END