教程5.13g加密所属

1、WORKSHOP 3FRAME STRUCTURE DESIGN WITH DIFFERENT SCALINGS10,000 lbs51A3, I3A1, I1A2, I210,000 lbs51A3, I3A1, I1A2, I2ObjectivesDesign the following frame to match the prescribed displacements at grid points 1 and 5.This workshop also illustrates the importance of scaling in optimization.Suggested Exe

2、rcise StepsGenerate an input file for Part A for MSC.Nastran Users.Compare the MSC.Nastran Input File for Part A.Submit Input File to MSC.Nastran for analysis.View ResultsXY Plots of Objective Function and Design Variable Function for Part A.Generate an input file for Part B for MSC.Nastran Users.Su

3、bmit Input File to MSC.Nastran for analysis.View ResultsXY Plots of Objective Function and Design Variable Function for Part B.Generate an input file for Part C for MSC.Nastran Users.Submit Input File to MSC.Nastran for analysis.View ResultsXY Plots of Objective Function and Design Variable Function

4、 for Part C.Model DescriptionPart AObjective Function: Match the displacements at grid points 1 and 5 to the values shown in the table using the following formulation:Design Variables: A1, A2, A3, I1, I2, I3Part BRedo the problem by normalizing the displacement using the foolowing formulation for th

5、e objective function:Model Description (Cont.)Part CRedo the problem by normalizing the design variables and the displacements using the following formulation for the objective function:whereand 1. Generate an Input File for MSC.NASTRAN Users$wkshp3a.dat$ID MSC, D200X10TIME 5SOL 200CENDTITLE = DISPL

6、ACEMENT MATCHING EXERCISESUBTITLE = PORTAL FRAME SUBJECT TO HORIZONTAL LOAD NO SCALINGECHO = SORTSPC = 10SET 10 = 1,5DISP = 10$ CASE CONTROL FOR ANALYSIS TYPE AND OBJECTIVE FUNCTION$ANALYSIS = .DESOBJ = .$SUBCASE 1 LABEL = LOAD CONDITION 1 LOAD = 1BEGIN BULKPARAM,POST,-1PARAM AUTOSPC YESPARAM GRDPNT

7、 0$ Analysis Model Description$MAT1 1 69000.0 0.30000 2.798-6SPC1 10 123456 9 13FORCE 1 5 10000.0 1.00000 0.00000 0.00000$GRID 1 500.0 1000.0 0.00000 345GRID 2 250.0 1000.0 0.00000 345GRID 3 0.000 1000.0 0.00000 345GRID 4 -250.0 1000.0 0.00000 345GRID 5 -500.0 1000.0 0.00000 345GRID 6 -500.0 750.00

8、0.00000 345GRID 7 -500.0 500.00 0.00000 345GRID 8 -500.0 250.00 0.00000 345vGRID 9 -500.0 0.0000 0.00000GRID 10 500.0 750.00 0.00000 345GRID 11 500.0 500.00 0.00000 345GRID 12 500.0 250.00 0.00000 345GRID 13 500.00 0.00000 0.00000Generate an input file using the data from pages 3-5 and 3-6 Part A. U

9、se the following input file as a starting point.CBAR 1 1 1 2 0.000 1.00000 0.000CBAR 2 1 2 3 0.000 1.00000 0.000CBAR 3 1 3 4 0.000 1.00000 0.000CBAR 4 1 4 5 0.000 1.00000 0.000CBAR 5 2 5 6 1.00000 0.000 0.000CBAR 6 2 6 7 1.00000 0.000 0.000CBAR 7 2 7 8 1.00000 0.000 0.000CBAR 8 2 8 9 1.00000 0.000 0

10、.000CBAR 9 3 1 10 1.00000 0.000 0.000CBAR 10 3 10 11 1.00000 0.000 0.000CBAR 11 3 11 12 1.00000 0.000 0.000CBAR 12 3 12 13 1.00000 0.000 0.000PBAR 1 1 100. 1.+7 1.+7 1.+7PBAR 2 1 50. 1.+6 1.+6 1.+7PBAR 3 1 150. 0.3+7 3.+6 1.+7$ $ Design Model Description$ Design Variables (A, I) DESVAR.$ Design Vari

11、able to Properties Relation DVPREL1.$ Design Responses DRESP1.$ Equation to describe Objective function $ The goal is to match the prescribed displacements.$ Define a table containing the prescribed displacements.DRESP2.$DEQATN.$DTABLE.$ Screening and Optimization Control Data DSCREEN,DISP,-0.5,1$23

12、45678910DOPTPRMDESMAX20DELP.8DPMIN.1ENDDATA2. Compare the MSC.NASTRAN Input File$ soln3a.datID MSC, D200X10TIME 5SOL 200CENDTITLE = DISPLACEMENT MATCHING EXERCISESUBTITLE = PORTAL FRAME SUBJECT TO HORIZONTAL LOAD NO SCALINGECHO = SORTSPC = 10SET 10 = 1,5DISP = 10ANALYSIS = STATICSDESOBJ = 700SUBCASE

13、 1 LABEL = LOAD CONDITION 1 LOAD = 1BEGIN BULKPARAM,POST,-1PARAM AUTOSPC YESPARAM GRDPNT 0$ Analysis Model DescriptionMAT1 1 69000.0 0.30000 2.798-6SPC1 10 123456 9 13FORCE 1 5 10000.0 1.00000 0.00000 0.00000$GRID 1 500.0 1000.0 0.00000 345GRID 2 250.0 1000.0 0.00000 345GRID 3 0.000 1000.0 0.00000 3

14、45GRID 4 -250.0 1000.0 0.00000 345GRID 5 -500.0 1000.0 0.00000 345GRID 6 -500.0 750.00 0.00000 345GRID 7 -500.0 500.00 0.00000 345GRID 8 -500.0 250.00 0.00000 345GRID 9 -500.0 0.0000 0.00000GRID 10 500.0 750.00 0.00000 345GRID 11 500.0 500.00 0.00000 345GRID 12 500.0 250.00 0.00000 345GRID 13 500.00

15、 0.00000 0.00000CBAR 1 1 1 2 0.000 1.00000 0.000CBAR 2 1 2 3 0.000 1.00000 0.000CBAR 3 1 3 4 0.000 1.00000 0.000CBAR 4 1 4 5 0.000 1.00000 0.000CBAR 5 2 5 6 1.00000 0.000 0.000CBAR 6 2 6 7 1.00000 0.000 0.000CBAR 7 2 7 8 1.00000 0.000 0.000CBAR 8 2 8 9 1.00000 0.000 0.000CBAR 9 3 1 10 1.00000 0.000

16、0.000CBAR 10 3 10 11 1.00000 0.000 0.000CBAR 11 3 11 12 1.00000 0.000 0.000CBAR 12 3 12 13 1.00000 0.000 0.000PBAR 1 1 100. 1.+7 1.+7 1.+7PBAR 2 1 50. 1.+6 1.+6 1.+7PBAR 3 1 150. 0.3+7 3.+6 1.+7$ Design Model Description$ Design Variables DESVAR,10,AREAHB,100.0,0.05,5000.0DESVAR,20,AREAC1,50.0,0.05,

17、5000.0DESVAR,30,AREAC2,150.0,0.05,5000.0DESVAR,40,IHB,10.E6,0.1,1.0E9DESVAR,50,IC1,1.E6,0.1,1.0E9DESVAR,60,IC2,3.E6,0.1,1.0E9$ Design Variable to Properties Relation DVPREL1,100,PBAR,1,4,5.0,0.0,+DVPA+DVPA,10,1.0DVPREL1,200,PBAR,2,4,5.0,0.0,+DVPB+DVPB,20,1.0DVPREL1,300,PBAR,3,4,5.0,0.0,+DVPC+DVPC,30

18、,1.0DVPREL1,400,PBAR,1,5,1.0,0.0,+DVPD+DVPD,40,1.0DVPREL1,500,PBAR,2,5,1.0,0.0,+DVPE+DVPE,50,1.0DVPREL1,600,PBAR,3,5,1.0,0.0,+DVPF+DVPF,60,1.0$The completed MSC.NASTRAN input file for PART A is shown below:(Contd) 2. Compare the MSC.NASTRAN Input File (Cont.)$ Design Responses DRESP1,100,HD1,DISP,10

19、,1,1DRESP1,200,VD1,DISP,11,2,1DRESP1,300,RD1,DISP,12,6,1DRESP1,400,HD5,DISP,20,1,5DRESP1,500,VD5,DISP,21,2,5DRESP1,600,RD5,DISP,22,6,5$ Equation to describe Objective function$ These equations scale the objective to get it away from zero.$DRESP2,700,ERROR,10,+DRESP2A+DRESP2A,DTABLE,T1,T2,T3,T4,T5,T6

20、,+DRESP2B+DRESP2B,DRESP1,100,200,300,400,500,600$DEQATN 10 OBJ(T1,T2,T3,T4,T5,T6,U1,U2,U3,U4,U5,U6) + = ( (U1-T1)*2 + (U2-T2)*2 + (U3-T3)*2 + + (U4-T4)*2 + (U5-T5)*2 + (U6-T6)*2 )/1.e-3$ 2 3 4 5 6 7 8 9 10DTABLE* T1 9.1449 E-01 T2 -9.56418E-02 +* T3 -9.7538 E-04 T4 1.068138 +* T5 9.56418E-02 T6 -1.0

21、012 E-03 +*$ Screening and Optimization Control Data DSCREEN,DISP,-0.5,1$ 2 3 4 5 6 7 8 9 10DOPTPRM DESMAX 20 DELP .8 DPMIN .1ENDDATAThe completed MSC.NASTRAN input file for PART A is shown below (Cont.): To submit the MSC.Nastran .dat file, find an available UNIX shell window and at the command pro

22、mpt enter nastran wkshp3a scr=yes. Monitor the run using the UNIX ps command:3. Submit Input File to MSC.Nastran For Analysis * 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)(SOFT CONVERGENCE ACHIEVED) NUMBER OF FINITE ELEMENT ANALYSES COMPLETED 10 NUMBER OF OPTIM

23、IZATIONS W.R.T. APPROXIMATE MODELS 9OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY CYCLE OBJECTIVE FROMOBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE NUMBER APPROXIMATE EXACT OF OF OPTIMIZATION ANALYSIS APPROXIMATION CONSTRAINT INITIAL 2.595013E+04 N/A 1 1.938923E+04 1.916654E+04 1.161876E-02 N/A 2 1.59339

24、1E+04 1.586977E+04 4.041748E-03 N/A 3 1.423345E+04 1.421479E+04 1.312590E-03 N/A 4 6.164483E+03 5.614997E+03 9.786058E-02 N/A 5 7.664992E+02 7.488561E+02 2.356002E-02 N/A 6 1.450883E+02 1.304158E+02 1.125052E-01 N/A 7 1.244497E+02 1.244324E+02 1.391819E-04 N/A 8 1.222592E+02 1.222557E+02 2.895602E-0

25、5 N/A 9 1.222557E+02 1.222557E+02 0.000000E+00 N/A When the run is completed, edit the wkshp3a.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 wkshp3a.f06, search for th

26、e word “H I S T O R Y”. Compare the results obtained in the .f06 file with the following:4. View Results4. View Results (Cont.) DESIGN VARIABLE HISTORY INTERNAL | EXTERNAL | | DV. ID. | DV. ID. | LABEL | INITIAL : 1 : 2 : 3 : 4 : 5 : 1 | 10 | AREAHB | 1.0000E+02 : 1.7973E+02 : 3.2351E+02 : 5.8310E+0

27、2 : 1.0496E+03 : 1.8946E+03 : 2 | 20 | AREAC1 | 5.0000E+01 : 9.0000E+01 : 1.6200E+02 : 2.9074E+02 : 5.2333E+02 : 9.4199E+02 : 3 | 30 | AREAC2 | 1.5000E+02 : 2.7000E+02 : 4.8600E+02 : 8.7480E+02 : 1.5746E+03 : 2.8344E+03 : 4 | 40 | IHB | 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 1.0038E+07

28、: 1.0363E+07 : 5 | 50 | IC1 | 1.0000E+06 : 1.0000E+06 : 1.0000E+06 : 1.0000E+06 : 1.8000E+06 : 3.2401E+06 : 6 | 60 | IC2 | 3.0000E+06 : 3.0000E+06 : 3.0000E+06 : 3.0000E+06 : 3.6348E+06 : 6.5427E+06 : INTERNAL | EXTERNAL | | DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 : 1 | 10 | AREAHB | 1.8

29、946E+03 : 9.2515E+02 : 8.5986E+02 : 8.5986E+02 : 2 | 20 | AREAC1 | 1.6956E+03 : 3.0520E+03 : 5.0000E+03 : 5.0000E+03 : 3 | 30 | AREAC2 | 5.0000E+03 : 5.0000E+03 : 5.0000E+03 : 5.0000E+03 : 4 | 40 | IHB | 1.0694E+07 : 1.0694E+07 : 1.0694E+07 : 1.0694E+07 : 5 | 50 | IC1 | 5.0565E+06 : 5.0566E+06 : 5.0

30、567E+06 : 5.0567E+06 : 6 | 60 | IC2 | 7.9001E+06 : 7.9002E+06 : 7.9003E+06 : 7.9003E+06 :5. XY Plots of Design ResultsMoment of InertiaFigure 3.1 Objective FunctionFigure 3.2 Design Variables (Moment of Inertia)5. XY Plots of Design Results (Cont.)Figure 3.3 Design Variables (Area)6. Generate an Inp

31、ut File for MSC.NASTRAN Users$ soln3b.dat$ID MSC, D200X10TIME 5SOL 200CENDTITLE = DISPLACEMENT MATCHING EXERCISESUBTITLE = PORTAL FRAME SUBJECT TO HORIZONTAL LOAD PARTIAL SCALINGECHO = SORTSPC = 10SET 10 = 1,5DISP = 10ANALYSIS = STATICSDESOBJ = 700SUBCASE 1 LABEL = LOAD CONDITION 1 LOAD = 1BEGIN BUL

32、KPARAM,POST,-1PARAM AUTOSPC YESPARAM GRDPNT 0$ Analysis Model Description$MAT1 1 69000.0 0.30000 2.798-6SPC1 10 123456 9 13FORCE 1 5 10000.0 1.00000 0.00000 0.00000$GRID 1 500.0 1000.0 0.00000 345GRID 2 250.0 1000.0 0.00000 345GRID 3 0.000 1000.0 0.00000 345GRID 4 -250.0 1000.0 0.00000 345GRID 5 -50

33、0.0 1000.0 0.00000 345GRID 6 -500.0 750.00 0.00000 345GRID 7 -500.0 500.00 0.00000 345GRID 8 -500.0 250.00 0.00000 345GRID 9 -500.0 0.0000 0.00000GRID 10 500.0 750.00 0.00000 345GRID 11 500.0 500.00 0.00000 345GRID 12 500.0 250.00 0.00000 345GRID 13 500.00 0.00000 0.00000The completed MSC.NASTRAN in

34、put file for PART B is shown below:(Cont.)CBAR 1 1 1 2 0.000 1.00000 0.000CBAR 2 1 2 3 0.000 1.00000 0.000CBAR 3 1 3 4 0.000 1.00000 0.000CBAR 4 1 4 5 0.000 1.00000 0.000CBAR 5 2 5 6 1.00000 0.000 0.000CBAR 6 2 6 7 1.00000 0.000 0.000CBAR 7 2 7 8 1.00000 0.000 0.000CBAR 8 2 8 9 1.00000 0.000 0.000CB

35、AR 9 3 1 10 1.00000 0.000 0.000CBAR 10 3 10 11 1.00000 0.000 0.000CBAR 11 3 11 12 1.00000 0.000 0.000CBAR 12 3 12 13 1.00000 0.000 0.000PBAR 1 1 100. 1.+7 1.+7 1.+7PBAR 2 1 50. 1.+6 1.+6 1.+7PBAR 3 1 150. 0.3+7 3.+6 1.+7$ Design Model Description$ Design Variables DESVAR,10,AREAHB,100.0,0.05,5000.0D

36、ESVAR,20,AREAC1,50.0,0.05,5000.0DESVAR,30,AREAC2,150.0,0.05,5000.0DESVAR,40,IHB,10.E6,0.1,1.0E9DESVAR,50,IC1,1.E6,0.1,1.0E9DESVAR,60,IC2,3.E6,0.1,1.0E9$ Design Variable to Properties Relation $DVPREL1,100,PBAR,1,4,5.0,0.0,+DVPA+DVPA,10,1.0DVPREL1,200,PBAR,2,4,5.0,0.0,+DVPB+DVPB,20,1.0DVPREL1,300,PBA

37、R,3,4,5.0,0.0,+DVPC+DVPC,30,1.0DVPREL1,400,PBAR,1,5,1.0,0.0,+DVPD+DVPD,40,1.0DVPREL1,500,PBAR,2,5,1.0,0.0,+DVPE+DVPE,50,1.0DVPREL1,600,PBAR,3,5,1.0,0.0,+DVPF+DVPF,60,1.0The completed MSC.NASTRAN input file for PART B is shown below (Cont.):6. Generate an Input File for MSC.NASTRAN Users (Cont.)(Cont

38、.)$ Design Responses DRESP1,100,HD1,DISP,10,1,1DRESP1,200,VD1,DISP,11,2,1DRESP1,300,RD1,DISP,12,6,1DRESP1,400,HD5,DISP,20,1,5DRESP1,500,VD5,DISP,21,2,5DRESP1,600,RD5,DISP,22,6,5$ Equation to describe Objective function$DRESP2,700,ERROR,10,+DRESP2A+DRESP2A,DTABLE,T1,T2,T3,T4,T5,T6,+DRESP2B+DRESP2B,DR

39、ESP1,100,200,300,400,500,600$ These equations scale the objective to get it away from zero.$234567812345678DEQATN 10 OBJ(T1,T2,T3,T4,T5,T6,U1,U2,U3,U4,U5,U6) + =( (1.-U1/T1)*2 + (1.-U2/T2)*2 + (1.-U3/T3)*2 + + (1.-U4/T4)*2 + (1.-U5/T5)*2 + (1.-U6/T6)*2 )/1.e-3$2345678910DTABLE*T19.1449E-01T2-9.56418

40、E-02+*T3-9.7538E-04T41.068138+*T59.56418E-02T6-1.0012E-03+*$ Screening and Optimization Control Data DSCREEN,DISP,-0.5,1$ 2 3 4 5 6 7 8 9 10DOPTPRM DESMAX 20 DELP .8 DPMIN .1ENDDATA6. Generate an Input File for MSC.NASTRAN Users (Cont.)The completed MSC.NASTRAN input file for PART B is shown below (

41、Cont.): To submit the MSC.Nastran .dat file, find an available UNIX shell window and at the command prompt enter nastran wkshp3b scr=yes. Monitor the run using the UNIX ps command:7. Submit Input File to MSC.Nastran For Analysis * 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 CONVERGE

42、NCE ACHIEVED)(SOFT CONVERGENCE ACHIEVED) NUMBER OF FINITE ELEMENT ANALYSES COMPLETED 7 NUMBER OF OPTIMIZATIONS W.R.T. APPROXIMATE MODELS 6OBJECTIVE AND MAXIMUM CONSTRAINT HISTORY CYCLE OBJECTIVE FROMOBJECTIVE FROM FRACTIONAL ERROR MAXIMUM VALUE NUMBER APPROXIMATE EXACT OF OF OPTIMIZATION ANALYSIS AP

43、PROXIMATION CONSTRAINT INITIAL 1.865449E+05 N/A 1 6.939589E+04 6.605289E+04 5.061096E-02 N/A 2 2.876406E+04 2.817723E+04 2.082625E-02 N/A 3 1.769462E+04 1.761205E+04 4.688396E-03 N/A 4 1.516328E+04 1.515656E+04 4.434832E-04 N/A 5 1.484496E+04 1.484487E+04 6.052176E-06 N/A 6 1.484487E+04 1.484487E+04

44、 0.000000E+00 N/A When the run is completed, edit the wkshp3b.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 wkshp3b.f06, search for the word “H I S T O R Y”. Compare t

45、he results obtained in the .f06 file with the following:8. View Results8. View Results (Cont.) 1 DISPLACEMENT MATCHING EXERCISE JANUARY 30, 2003 MSC.NASTRAN 4/ 9/01 PAGE 40 PORTAL FRAME SUBJECT TO HORIZONTAL LOAD PARTIAL SCALING 0 DESIGN VARIABLE HISTORY INTERNAL | EXTERNAL | | DV. ID. | DV. ID. | L

46、ABEL | INITIAL : 1 : 2 : 3 : 4 : 5 : 1 | 10 | AREAHB | 1.0000E+02 : 1.8000E+02 : 3.2379E+02 : 5.8282E+02 : 1.0491E+03 : 1.1104E+03 : 2 | 20 | AREAC1 | 5.0000E+01 : 8.9982E+01 : 1.6189E+02 : 2.9140E+02 : 5.2451E+02 : 8.5592E+02 : 3 | 30 | AREAC2 | 1.5000E+02 : 2.7000E+02 : 4.8600E+02 : 8.7480E+02 : 8

47、.9565E+02 : 8.8249E+02 : 4 | 40 | IHB | 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 1.0000E+07 : 5 | 50 | IC1 | 1.0000E+06 : 1.0000E+06 : 1.0000E+06 : 1.0000E+06 : 1.0000E+06 : 1.0001E+06 : 6 | 60 | IC2 | 3.0000E+06 : 3.0000E+06 : 3.0000E+06 : 3.0000E+06 : 3.0000E+06 : 3.0000E+0

48、6 : INTERNAL | EXTERNAL | | DV. ID. | DV. ID. | LABEL | 6 : 7 : 8 : 9 : 10 : 11 : 1 | 10 | AREAHB | 1.1104E+03 : 2 | 20 | AREAC1 | 8.5592E+02 : 3 | 30 | AREAC2 | 8.8249E+02 : 4 | 40 | IHB | 1.0000E+07 : 5 | 50 | IC1 | 1.0001E+06 : 6 | 60 | IC2 | 3.0000E+06 :9. XY Plot of Design ResultsObjective Func

49、tionAreaFigure 3.4 Objective FunctionFigure 3.5 Design Variables (Area)9. XY Plot of Design Results (Cont.)Moment of InertiaFigure 3.6 Design Variables (Moment of Inertia)10. Generate an Input File for Part C$ soln3c.datID MSC, D200X10TIME 5SOL 200CENDTITLE = DISPLACEMENT MATCHING EXERCISESUBTITLE =

50、 PORTAL FRAME SUBJECT TO HORIZONTAL LOAD FULL SCALINGECHO = SORTSPC = 10SET 10 = 1,5DISP = 10ANALYSIS = STATICSDESOBJ = 700SUBCASE 1 LABEL = LOAD CONDITION 1 LOAD = 1BEGIN BULKPARAM,POST,-1PARAM AUTOSPC YESPARAM GRDPNT 0$ Analysis Model DescriptionMAT1 1 69000.0 0.30000 2.798-6SPC1 10 123456 9 13FOR

51、CE 1 5 10000.0 1.00000 0.00000 0.00000GRID 1 500.0 1000.0 0.00000 345GRID 2 250.0 1000.0 0.00000 345GRID 3 0.000 1000.0 0.00000 345GRID 4 -250.0 1000.0 0.00000 345GRID 5 -500.0 1000.0 0.00000 345GRID 6 -500.0 750.00 0.00000 345GRID 7 -500.0 500.00 0.00000 345GRID 8 -500.0 250.00 0.00000 345GRID 9 -5

52、00.0 0.0000 0.00000GRID 10 500.0 750.00 0.00000 345GRID 11 500.0 500.00 0.00000 345GRID 12 500.0 250.00 0.00000 345GRID 13 500.00 0.00000 0.00000Generate an input file called wkshp3c.dat using wkshp3b.dat as a starting point. The completed MSC.Nastran input file for PART C is shown below.(Cont.)10.

53、Generate an Input File for Part C (Cont.)CBAR 1 1 1 2 0.000 1.00000 0.000CBAR 2 1 2 3 0.000 1.00000 0.000CBAR 3 1 3 4 0.000 1.00000 0.000CBAR 4 1 4 5 0.000 1.00000 0.000CBAR 5 2 5 6 1.00000 0.000 0.000CBAR 6 2 6 7 1.00000 0.000 0.000CBAR 7 2 7 8 1.00000 0.000 0.000CBAR 8 2 8 9 1.00000 0.000 0.000CBA

54、R 9 3 1 10 1.00000 0.000 0.000CBAR 10 3 10 11 1.00000 0.000 0.000CBAR 11 3 11 12 1.00000 0.000 0.000CBAR 12 3 12 13 1.00000 0.000 0.000PBAR 1 1 100. 1.+7 1.+7 1.+7PBAR 2 1 50. 1.+6 1.+6 1.+7PBAR 3 1 150. 0.3+7 3.+6 1.+7$ Design Model Description$ Design Variables $ Note these are unit design variabl

55、es$ You will have to scale final dvs by coeff.$DESVAR,10,AREAHB,1.0,0.05,5000.0DESVAR,20,AREAC1,1.0,0.05,5000.0DESVAR,30,AREAC2,1.0,0.05,5000.0DESVAR,40,IHB,1.0,0.1,1.0E9DESVAR,50,IC1,1.0,0.1,1.0E9DESVAR,60,IC2,1.0,0.1,1.0E9$ Design Variable to Properties Relation DVPREL1,100,PBAR,1,4,5.0,0.0,+DVPA+

56、DVPA,10,100.0DVPREL1,200,PBAR,2,4,5.0,0.0,+DVPB+DVPB,20,50.0DVPREL1,300,PBAR,3,4,5.0,0.0,+DVPC+DVPC,30,150.0DVPREL1,400,PBAR,1,5,1.0,0.0,+DVPD+DVPD,40,10.0E6DVPREL1,500,PBAR,2,5,1.0,0.0,+DVPE+DVPE,50,1.0E6DVPREL1,600,PBAR,3,5,1.0,0.0,+DVPF+DVPF,60,3.0E6$(Cont.)10. Generate an Input File for Part C (

57、Cont.)$ Design Responses DRESP1,100,HD1,DISP,10,1,1DRESP1,200,VD1,DISP,11,2,1DRESP1,300,RD1,DISP,12,6,1DRESP1,400,HD5,DISP,20,1,5DRESP1,500,VD5,DISP,21,2,5DRESP1,600,RD5,DISP,22,6,5$ Equation to describe Objective function$DRESP2,700,ERROR,10,+DRESP2A+DRESP2A,DTABLE,T1,T2,T3,T4,T5,T6,+DRESP2B+DRESP2

58、B,DRESP1,100,200,300,400,500,600$ These equations scale the objective to get it away from zero.$DEQATN 10 OBJ(T1,T2,T3,T4,T5,T6,U1,U2,U3,U4,U5,U6) + =( (1.-U1/T1)*2 + (1.-U2/T2)*2 + (1.-U3/T3)*2 + + (1.-U4/T4)*2 + (1.-U5/T5)*2 + (1.-U6/T6)*2 )/1.e-3$ 2 3 4 5 6 7 8 9 10DTABLE* T1 9.1449 E-01 T2 -9.56

59、418E-02 +* T3 -9.7538 E-04 T4 1.068138 +* T5 9.56418E-02 T6 -1.0012 E-03 +*$ Screening and Optimization Control Data DSCREEN,DISP,-0.5,1$ 2 3 4 5 6 7 8 9 10$DOPTPRM DESMAX 20$ TRY THE ABOVE DOPTRM AND CHECK OUT DISPLACEMENTSDOPTPRM DESMAX 20 DELP .8 DPMIN .1ENDDATA To submit the MSC.Nastran .dat fil

60、e, find an available UNIX shell window and at the command prompt enter nastran wkshp3c scr=yes. Monitor the run using the UNIX ps command:11. Submit Input File to MSC.Nastran For Analysis * 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)(SOFT CONVERGENCE ACHIEVED)