$ DEC/CMS REPLACEMENT HISTORY, Element D200X6.DAT
$ *5    25-MAR-1994 16:07:37 A_BOYADJIAN "68 PLUS/G/ UPDATED FOR USER'S GUIDE JFOR DESIGN SENSITIVITY & OPTIMIZATION"
$ *4    15-SEP-1993 11:59:22 A_BOYADJIAN "68 PLUS/G/ CHANGED FOR MODIFIED USER INTERFACE IN V68 /SSP/EHJ"
$ *3    30-SEP-1992 17:00:30 A_BOYADJIAN "68 PLUS/G/ Changed PARAM,APPC to ANALYSIS = per KLK/SSP"
$ *2    15-OCT-1990 16:14:17 DEC "67 PLUS/G/ MODIFIED TO REDUCE OUTPUT KLK"
$ *1    15-JUN-1990 17:49:32 CMSMGR "66B PLUS/G 66B/ Initial installation of TPL test problems"
$ DEC/CMS REPLACEMENT HISTORY, Element D200X6.DAT
nastran oldq4k
ID MSC, D200X6 $ V68 G_MOORE 24-MAR-1994
TIME  10
SOL 200       $  OPTIMIZATION
INCLUDE 'evds200a.v707'
INCLUDE 'evdsa.v707'
CEND
TITLE    = VIBRATION OF A BEAM.                                D200X6
dsaprt(end=sens)=all
ECHO     = UNSORT
DISP     = ALL
STRESS   = ALL
METHOD   = 1
ANALYSIS = MODES
DESOBJ(MIN) = 1    $ OBJECTIVE FUNCTION DEFINITION
DESSUB =      10   $ CONSTRAINT SET SELECTION
$
BEGIN BULK
$
$------------------------------------------------------------------------
$ ANALYSIS MODEL:
$------------------------------------------------------------------------
$
EIGRL,  1,      ,       ,
GRID,   1,      ,       0.0,    0.0,    -3.0,   ,       123456
GRID,   2,      ,       20.0,   0.0,    -3.0,   ,       2456
GRID,   3,      ,       40.0,   0.0,    -3.0,   ,       2456
GRID,   4,      ,       60.0,   0.0,    -3.0,   ,       2456
GRID,   5,      ,       0.0,    0.0,    3.0,    ,       123456
GRID,   6,      ,       20.0,   0.0,    3.0,    ,       2456
GRID,   7,      ,       40.0,   0.0,    3.0,    ,       2456
GRID,   8,      ,       60.0,   0.0,    3.0,    ,       2456
$
CROD,   1,      201,    5,      6
CROD,   2,      202,    6,      7
CROD,   3,      203,    7,      8
CROD,   7,      201,    1,      2
CROD,   8,      202,    2,      3
CROD,   9,      203,    3,      4
PROD,   201,    1,      1.0,    0.0
PROD,   202,    1,      1.0,    0.0
PROD,   203,    1,      1.0,    0.0
$
CQUAD4, 4,      204,    1,      2,      6,      5
CQUAD4, 5,      205,    2,      3,      7,      6
CQUAD4, 6,      206,    3,      4,      8,      7
PSHELL, 204,    1,      0.2
PSHELL, 205,    1,      0.2
PSHELL, 206,    1,      0.2
$
CONM2,  10,     2,      ,       15.0
CONM2,  11,     3,      ,       15.0
CONM2,  12,     4,      ,       15.0
CONM2,  14,     6,      ,       15.0
CONM2,  15,     7,      ,       15.0
CONM2,  16,     8,      ,       15.0
$
MAT1,   1,      1.03E7, ,       0.3,    0.1
PARAM,  WTMASS, 0.002588
PARAM,  GRDPNT, 1
$
$------------------------------------------------------------------------
$ DESIGN MODEL:
$------------------------------------------------------------------------
$
$...Define the design variables
$DESVAR,ID,     LABEL,  XINIT,  XLB,    XUB,    DELXV
$
DESVAR, 1,      A1,     1.0,    0.01,   100.0
DESVAR, 2,      A2,     1.0,    0.01,   100.0
DESVAR, 3,      A3,     1.0,    0.01,   100.0
DESVAR, 4,      T1,     0.2,    0.001,  10.0
DESVAR, 5,      T2,     0.2,    0.001,  10.0
DESVAR, 6,      T3,     0.2,    0.001,  10.0
$
$...Relate the design variables to analysis model properties
$DVPREL1,ID,    TYPE,   PID,    FID,    PMIN,   PMAX,   C0,     ,       +
$+,     DVIDD1, COEF1,  DVID2,  COEF2,  ...
$
DVPREL1,1,      PROD,   201,    4,      ,       ,       ,       ,       +DP1
+DP1,   1,      1.0
DVPREL1,2,      PROD,   202,    4,      ,       ,       ,       ,       +DP2
+DP2,   2,      1.0
DVPREL1,3,      PROD,   203,    4,      ,       ,       ,       ,       +DP3
+DP3,   3,      1.0
DVPREL1,4,      PSHELL, 204,    4,      ,       ,       ,       ,       +DP4
+DP4,   4,      1.0
DVPREL1,5,      PSHELL, 205,    4,      ,       ,       ,       ,       +DP5
+DP5,   5,      1.0
DVPREL1,6,      PSHELL, 206,    4,      ,       ,       ,       ,       +DP6
+DP6,   6,      1.0
$
$...Identify the analysis responses to be used in the design model
$DRESP1,ID,     LABEL,  RTYPE,  PTYPE,  REGION, ATTA,   ATTB,   ATT1,   +
$+,     ATT2,   ...
$
DRESP1, 1,      W,      WEIGHT
DRESP1, 2,      F1,     FREQ,    ,      ,       1
$
$...Use these responses to define the objective (in case control)
$   and the constraints:
$DCONSTR,DCID,  RID,    LALLOW, UALLOW
DCONSTR,10,     2,      30.                 $ lower bound = 20 Hz
$
$...Optional override of design optimization parameters:
DOPTPRM, IPRINT,2,      DESMAX, 10,     DELP,   0.5,    P1,     1,      +
+,      P2,    15,	DELB	0.01
$
PARAM,ITERATE,YES $ PERFORM ITERATION
PARAM,ITRPRNT,NO $ PRINT AT EACH ITERATION
PARAM,ITMAX,10    $ DO A MAXIMUM OF 10 ITERATIONS
$
PARAM   DSZERO  1.E-2
PARAM, LAMBDAS, -0.8+4  $ HALF-WAY BETWEEN 1st, 2nd ROOTS
$PARAM, LAMBDAS, -1.8+4 $ HALF-WAY BETWEEN 2ND, 3RD ROOTS
$PARAM, LAMBDAS, 0.0 $ NO SHIFT
param,stepsize,0.01
param,eigvecds,1
param,ndv,6
param,nmod,1
param,optexit,4
ENDDATA