教程5.13g加密所属官方ws

1、WORKSHOP 12EXTERNAL RESPONSE (DRESP3)Objective:Learn how to set up an external response optimization task using the same box design problem used in Workshop 2. You will write Fortran routines to define the volume and surface area instead of using DEQATN entries.Suggested Exercise StepsGenerate an In

2、put File for MSC.NASTRAN users.Compare the Input File Modify DRESP3 Server Routine, r3sgrtCompare the modified server routine, r3sgrtModify server routine to define volume and surface areaCompare modified server routine, r3svaldBuild the response serverSubmit Input file for MSC.Nastran for analysisV

3、iew ResultsModel DescriptionObjective Function: Minimize the material surface area.S = 2 * ( W * H + L * H + 2.0 * W * L )Design Variables: W (Width), H (Height), L (Length)Constraints:Volume: V = H * W * L 2.0Initial DesignOptimum Design H = 1.5 H = 2.0W = 1.5 W = 1.0 L = 1.5 L = 1.0 S = 18.0 S = 1

4、2.0 V = 3.375 V = 2.0 Following the instructions given in the input file, wkshp12_fillin.dat to complete the exercise. The following input file is provided as a starting point.1. Generating an Input File for MSC.NASTRAN Users$ Define an external response group using a CONNECT entry. Notice that$ the

5、 group name must be referenced by the subsequent DRESP3 entriesCONNECT TIME 5diag 8,56 $SOL 200CENDTITLE = BOX MATERIAL MINIMIZATION WITH A VOLUME CONSTRAINT D200X8DESOBJ = 2DESSUB = 1SUBTITLE = INITIAL DESIGN: W=1.5, L=1.5, AND H=1.5$DISPLACEMENT = ALLSTRESS = NONESPC = 1ANALYSIS = STATICSLOAD = 1B

6、EGIN BULK$ DUMMY STRUCTURAL MODEL (NOT USED)FORCE 1 2 1.0E7 1.0 0.0 0.0GRID 1 0.0 0.0 0.0GRID 2 1.0 0.0 0.0MAT1 1 1.0E7 0.3 0.1CROD 1 1 1 2PROD 1 1 1.0SPC1 1 23456 2SPC1 1 123456 1$ DUMMY DESIGN PROPERTY RELATIONS$DVPREL11PROD1431.0$ BOX DESIGN PROBLEM DEFINITION$DESVAR DESVAR 2 L 1.5 0.1 10.0DESVAR

7、 3 H 1.5 0.1 10.0$ Define volume and surface area using DRESP3 entriesDRESP3DRESP3DCONSTR 1 1 2.0 1.0E35DSCREEN DRESP3 -100.0 1DOPTPRM DESMAX 10delp1.0gmax1.e-16ctmin1.e-16ENDDATA2. Compare the Input FileThe completed input file is shown below:CONNECT DRESP3 WKSHOP EXTRESPTIME 5diag 8,56 $SOL 200CEN

8、DTITLE = BOX MATERIAL MINIMIZATION WITH A VOLUME CONSTRAINT D200X8DESOBJ = 2DESSUB = 1SUBTITLE = INITIAL DESIGN: W=1.5, L=1.5, AND H=1.5$DISPLACEMENT = ALLSTRESS = NONESPC = 1ANALYSIS = STATICSLOAD = 1BEGIN BULK$ DUMMY STRUCTURAL MODEL (NOT USED)$FORCE 1 2 1.0E7 1.0 0.0 0.0GRID 1 0.0 0.0 0.0GRID 2 1

9、.0 0.0 0.0MAT1 1 1.0E7 0.3 0.1CROD 1 1 1 2PROD 1 1 1.0SPC1 1 23456 2SPC1 1 123456 1$ DUMMY DESIGN PROPERTY RELATIONS$DVPREL11PROD1431.0$ BOX DESIGN PROBLEM DEFINITION$DESVAR 1 W 1.5 0.1 10.0DESVAR 2 L 1.5 0.1 10.0DESVAR 3 H 1.5 0.1 10.0DRESP3 1 VOLUME WKSHOP VOLUME DESVAR 1 2 3DRESP3 2 SURFACE WKSHO

10、P SURFACE DESVAR 1 2 3DCONSTR 1 1 2.0 1.0E35DSCREEN DRESP3 -100.0 1DOPTPRM DESMAX 10delp1.0gmax1.e-16ctmin1.e-16ENDDATA3. Modify DRESP3 Server Routine, R3SGRTModify DRESP3 server routine, r3grt to define response types. SUBROUTINE R3SGRT(GRPID,TYPNAM,ERROR)C 09/18/00;14:24:09 cdp 2001t1 na001130C Up

11、date dr3srv per SHZC C PURPOSE: VERIFY THE EXTERNAL RESPONSE TYPEC ArgumentsC GRPID input integer - Group IDC TYPNAM input character*8 - Name of external response typeC ERROR input/output integer -error code for the call.CC MethodC Match the user input: typnam with the list of availableC external re

12、sponse types. If no match is found, set error code.C Called byC R3CGRTC NOTE:C The writer of this routine is responsible to specifyC NTYPES and R3TYPE.C C VARIABLES PASSED IN INTEGER GRPID, ERROR CHARACTER*8 TYPNAMC LOCAL VARIABLES INTEGER BADTYPCC declaration for character array R3TYPE and initiali

13、ze character C array R3TYPE with corresponding R3TYPE from DRESP3s DATA BADTYP/7554/ ERROR = 0 DO 100 ITYPE = 1, NTYPES IF (TYPNAM .EQ. R3TYPE(ITYPE) THEN GOTO 200 END IF100 CONTINUE ERROR = BADTYP200 CONTINUE RETURN END4. Compare the Server Routine, R3SGRTThe listing of server routine, r3grt is sho

14、wn here: SUBROUTINE R3SGRT(GRPID,TYPNAM,ERROR)C 09/18/00;14:24:09 cdp 2001t1 na001130C Update dr3srv per SHZC CC PURPOSE: VERIFY THE EXTERNAL RESPONSE TYPECC GRPID input integer - Group IDC TYPNAM input character*8 - Name of external response typeC ERROR input/output integer -error code for the call

15、.CC MethodC Match the user input: typnam with the list of availableC external response types. If no match is found, set error code.CC Called byC R3CGRTCC NOTE:C The writer of this routine is responsible to specifyC NTYPES and R3TYPE.C CC VARIABLES PASSED INC INTEGER GRPID, ERROR CHARACTER*8 TYPNAMCC

16、 LOCAL VARIABLESC INTEGER NTYPES, BADTYP PARAMETER(NTYPES=2) CHARACTER*8 R3TYPE(NTYPES)C DATA BADTYP/7554/ DATA R3TYPE/VOLUME,SURFACE/ ERROR = 0 DO 100 ITYPE = 1, NTYPES IF (TYPNAM .EQ. R3TYPE(ITYPE) THEN GOTO 200 END IF100 CONTINUE ERROR = BADTYP200 CONTINUEC RETURN END Modify DRESP3 server routine

17、, r3vald to define volume and surface area.5. Modify Server Routine to Define Volume and Surface Area SUBROUTINE R3SVALD(GRPID,TYPNAM, . NITEMS,ARGLIS, . NSIZE, ARGVAL, . NWRDA8,ARGCHR, . DR3VAL,ERROR)C CC PURPOSE: COMPUTE THE EXTERNAL RESPONSECC GRPID input integer - Group idC TYPNAM input characte

18、r*8 - Name of external response typeC NITEMS input integer - dimension of array ARGLISC NSIZE input integer - dimension of array ARGVALC NWRDA8 input integer - dimension of character array ARGCHRC ARGLIS input integer - array of no. of items for each argumentC ARGVAL input double - array of all argu

19、ment valuesC ARGCHR input character*8 - array of charactersC DR3VAL output double - value of the external responseC ERROR input/output integer -error code for the call.C 0 = print error messagesC 1 = do not print error messages.CC Portion of routine is omittedCC LOCAL VARIABLESC INTEGER BADTYPC decl

20、are user defined variables C DATA BADTYP /7554/CC THE USER-SUPPLIED EQUATION TO DEFINE THE EXTERNAL RESPONSEC DEQATN 1 F(W,L,H)=W*L*HC DEQATN 2 F(W,L,H)=2.0*(W*H+L*H+2.0*W*L)C ERROR = 0CC define equations of volume and surface areaC RETURN END6. Compare the modified Server Routine, R3SVALDThe listin

21、g of server routine, r3vald is shown here. Notice there are many ways to modify this routine: SUBROUTINE R3SVALD(GRPID,TYPNAM, . NITEMS,ARGLIS, . NSIZE, ARGVAL, . NWRDA8,ARGCHR, . DR3VAL,ERROR)C 07/02/01;11:30:11 cdp V2001 nxf01800C Update per SHZC CC PURPOSE: COMPUTE THE EXTERNAL RESPONSECC GRPID i

22、nput integer - Group idC TYPNAM input character*8 - Name of external response typeC NITEMS input integer - dimension of array ARGLISC NSIZE input integer - dimension of array ARGVALC NWRDA8 input integer - dimension of character array ARGCHRC ARGLIS input integer - array of no. of items for each arg

23、umentC ARGVAL input double - array of all argument valuesC ARGCHR input character*8 - array of charactersC DR3VAL output double - value of the external responseC ERROR input/output integer -error code for the call.C 0 = print error messagesC 1 = do not print error messages.CC MethodC A)Set up variou

24、s parameters from the argument listC B)Evaluate the external response based on the given TYPNAMC C)Return BADTYP error if TYPNAM is not matched here.CC Called byC SendR3svaldC C6. Compare the modified Server Routine (Cont.)CC VARIABLES PASSED INC CHARACTER*8 TYPNAM, ARGCHR(NWRDA8) INTEGER GRPID, NIT

25、EMS, NSIZE, ARGLIS(NITEMS), ERROR DOUBLE PRECISION ARGVAL(NSIZE), DR3VALCC LOCAL VARIABLESC INTEGER BADTYP DOUBLE PRECISION W,L,HC DATA BADTYP /7554/CC THE USER-SUPPLIED EQUATION TO DEFINE THE EXTERNAL RESPONSEC DEQATN 1 F(W,L,H)=W*L*HC DEQATN 2 F(W,L,H)=2.0*(W*H+L*H+2.0*W*L)C ERROR = 0CC set up par

26、ameters for various argument itemsC W = argval(1) L = argval(2) H = argval(3) IF (TYPNAM .EQ. VOLUME) THEN VOL = W*L*H dr3val = VOL ELSE IF (TYPNAM .EQ. SURFACE) THEN SURF = 2.D0*(W*H+L*H+2.D0*W*L) DR3VAL = SURF ELSE ERROR = BADTYP END IF RETURN END Build the response server using the following comm

27、and: msc2001 ./dr3srv build Generate a connection evaluator file, called wkshp12_eval that associates the external response group with the actual server program. The full content of file, wkshp12_eval is shown here. Submit the Nastran job by adding a new keyword, gmconn where gmconn refers to the ev

28、aluator file name. 7. Build Response ServerEXTRESP, - , /home/wkshop/server/dr3serv To submit the MSC.Nastran .dat file, find an available UNIX shell window and at the command prompt enter Nastran fn scr=yes gmconn=wkshp12_eval8. Submit Input File to MSC.Nastran For AnalysisWhen the run is completed

29、, edit the wkshp12.f06 file and search for the word FATAL. If no matches exist, search for the word WARNING. Determine whether existing WARNING messages indicate modeling errors.While still editing wkshp12.f06, search for the word “H I S T O R Y”. Compare the results obtained in the .f06 file with the following:9. View Results * S U M M A R Y O F D E S I G N C Y C L E H I S T O R Y * (HARD CONVERGENCE ACHIEVED)