INIT MASTER(S)
NASTRAN SYSTEM(442)=-1, SYSTEM(319)=1
ID Femap,Femap
SOL SEMODES
TIME 1000
CEND
$  Benchmark Supercritical Wing (BSCW) Faux model$
$  Generated by Jennifer Heeg, 2014 in support of the
$  AePW-2 analysis effort
$
$
$  Note:  a separate .bdf has to be generated to use with the DDF program
$  that maps the mode shapes to the aerodynamic grid.  That .bdf has formatted
$  cards rather than the free format used here.

  TITLE = BSCW Faux aeroelastic model
  Autospc(noprint)
  ECHO = sort(except grid, cquad4, ctria3, CHEXA, CTETRA, CBEAM, CPENTA,+
         CORD2R, PSOLID, RBE2 )
  DISPLACEMENT = ALL
  SPCFORCE(PLOT) = ALL
$  ESE(PLOT) = ALL
$
  METHOD = 1
$
BEGIN BULK
$
PARAM,  POST,    -1
PARAM,  OGEOM,   NO
PARAM,  AUTOSPC,YES
PARAM,  GRDPNT,   2
PARAM,  NOFISR,   1
EIGRL         1                      20       0                    MASS
$  May 2, 2014:
$  The published properties of the PAPA mounted wing are not relative to the 30%chord
$  location, but rather relative to the midchord
$  to correct this error, all grid points at the 30% chord are being moved to the midchord location
$  CORRECTED VERSION:  GRID 2 is at the wing root at the center of rotation(x/c = 0.5, 16 inch chord --> 8 inches from leading edge)
$  GRID 2 is at the wing root at the center of rotation(x/c = 0.3, 16 inch chord = 4.8 inches from leading edge)
$  GRId 1000 is the ground point that the springs will attach to

GRID,1000,,8.00,0.0,0.0,,123456

$  adding grid points to define a rigid wing
$  matrix of 3 points chordwise by 5 points spanwise.
$  chord-wise points are at leading edge, rotational axis and trailing edge
$  spanwise points are at root, 30% 60 % span, 95% span,and wing tip



GRID,1,,0.00,0.00,0.00,,
GRID,2,,8.00,0.00,0.00,,1246
GRID,3,,16.00,0.00,0.00,,
GRID,4,,0.00,9.60,0.00,,
GRID,5,,8.00,9.60,0.00,,
GRID,6,,16.00,9.60,0.00,,
GRID,7,,0.00,19.20,0.00,,
GRID,8,,8.00,19.20,0.00,,
GRID,9,,16.00,19.20,0.00,,
GRID,10,,0.00,30.40,0.00,,
GRID,11,,8.00,30.40,0.00,,
GRID,12,,16.00,30.40,0.00,,
GRID,13,,0.00,32.00,0.00,,
GRID,14,,8.00,32.00,0.00,,
GRID,15,,16.00,32.00,0.00,,


CQUAD4,3001,3000,1,2,5,4
CQUAD4,3002,3000,2,3,6,5
CQUAD4,3003,3000,4,5,8,7
CQUAD4,3004,3000,5,6,9,8
CQUAD4,3005,3000,7,8,11,10
CQUAD4,3006,3000,8,9,12,11
CQUAD4,3007,3000,10,11,14,13
CQUAD4,3008,3000,11,12,15,14



PSHELL,3000,3000,1.0,3000,,3000


$  celas1 1000 is the vertical displacement linear spring
$  celas1 1001 is the torsional rotational linear spring
CELAS1,1000,1000,2,3,1000,3
CELAS1,1001,1001,2,5,1000,5
PELAS,1000,219.75
PELAS,1001,35568.

$  leading edge points
CONM2,3000,1,,0.0357
CONM2,3004,4,,0.0357
CONM2,3007,7,,0.0357
CONM2,3010,10,,0.0357
CONM2,3013,13,,0.0357
$  rotation axis points
CONM2,3002,2,,0.0290,,,,,+3000
+3000,,,10.476
CONM2,3008,8,,0.0290
CONM2,3005,5,,0.0290
CONM2,3011,11,,0.0290
CONM2,3014,14,,0.0290

$  trailing edge points
CONM2,3003,3,,0.0357
CONM2,3006,6,,0.0357
CONM2,3009,9,,0.0357
CONM2,3012,12,,0.0357
CONM2,3015,15,,0.0357



$  very stiff light material that I may use in future renditions
$MAT1,3000,2.85+19,,0.3,1.098-15
MAT1,3000,2.85+10,,0.3,1.098-12
$
$
ENDDATA bb7aeb05
$