$ FILE CBAR1G.DAT - CREATE CRAIG_BAMPTON MODELS AND OTM
$
SOL 103
TIME 10
diag 8,14,56
INCLUDE 'alter1ga.v707' $ ALTER TO CREATE BOUNDARY MATRICES
CEND
TITLE = CREATE EXTERNAL SUPERELEMENT MODEL OF S.E. 10          CBAR1G
SUBTITLE = TEST RESIDUAL VECTORS
SEALL = 10
$SEDR = 10
SPC = 1
$
$ loadset and second subcase included for generation of residual vectors
$
LOADSET = 99
SUPER =ALL
$
$ place output requests for OTM quantities before first subcase
$
  DISP = ALL
  FORCE = ALL
stress = all
SUBCASE 1
  LABEL = CMS OF SUPERELEMENT 10
  METHOD=101
$ FOLLOWING ADDED TO CALCULATE MODAL DAMPING FOR THE SUPERELEMENT
  SDAMP = 10
$
$ second subcase no longer necessary for residual vectors, see 
$ bulk data section for DTI,resmeth and eigr 2000 - now used instead
$
$SUBCASE 2 $ only need of residual vectors are to be calculated
$  LABEL = EIGENVALUE SOLUTION TO ORTHOGONALIZE RESIDUAL VECTORS
$  METHOD = 2
BEGIN BULK
dti,resmeth,1,cr,1,dw,5,5,1,+resmeth
+resmeth,2000
eigr,2000,agiv,,,,2000
param,drmh,1
PARAM,GRDPNT,0
$
$ SET FIXEDB TO -1 TO GET COMPONENT MODES DATA RECOVERY
$  NECESSARY FOR CREATION OF OTM
$
PARAM,FIXEDB,-1
$
$ SET MAKEMAT TO SEID FOR MATRIX GENERATION
$
PARAM,MAKEMAT,10
$
EIGR,101,MGIV,0.,100.,,4
$
$ EIGR FOR RESIDUAL VECTOR OTRHOGONALIZATION
$ ASK FOR ALL MODES OF REDUCED PROBLEM
$  THIS WILL GIVE THE SAME MODES AS FOUND ORIGINALLY
$  PLUS ONE 'HIGH-FREQUENCY' MODE FOR EACH STATIC
$  LOADING VECTOR USED
$
EIGR,2,MGIV,,,,100
$
$ LSEQ and loadings for use in creating residual vectors
$  (the loadings should be the static distribution of forces
$   moments, etc which will be applied dynamically - this is
$   not required, but helps to improve the accuracy)
$
PARAM,RESvec,yes 
LSEQ,99,5,10
LSEQ,99,6,20
FORCE,10,8,,1.,1.,0.,0.
FORCE,20,7,,1.,0.,1.,1.
$
$ ADD MODAL DAMPING FOR S.E. 10
$
PARAM,HFREQ,99999.
TABDMP1,10,
,0.,.01,100.,.01,endt
$
$ ADD MODAL COORDINATES FOR S.E. 10
$  (it is preferred to use GRID points or a multiple of
$   6 SPOINTS to allow for multi-level assemblies)
$
SPOINT,1001,THRU,1012
SEQSET1,10,0,1001,THRU,1012
$
$ DEFINE SEUSET U2 FOR INTERIOR POINTS WITH POTENTIAL LOADINGS
$  AND WHICH WILL BE USED FOR VISUALIZING THE SUPERELEMENT WHEN
$  IT IS ATTACHED TO THE ASSEMBLY
$
DEFUSET,U2,U2
SEUSET,10,U2,8,123456
SEUSET,10,U2,9,123456
SEUSET,10,U2,10,123456
$
GRID,1,,0.,0.,0.
=,*(1),=,*(10.),==
=(8)
GRID,11,,0.,1.,0.,,123456
GRID,12,,0.,1.,0.,,123456
CBAR,1,1,1,2,11
=,*(1),=,*(1),*(1),==
=(2)
CBAR,5,1,5,6,12
=,*(1),=,*(1),*(1),==
=(3)
$
$ DEFINE SUPERELEMENT
$
SESET,10,6,THRU,10
SESET,10,12
PBAR,1,1,1.,10.,10.,10.
,1.,1.,-1.,1.,-1.,-1.,1.,-1.
PARAM,AUTOSPC,YES
MAT1,1,30.+6,,.3,.283
SPC1,1,123456,1
$ add rigid element to test data recovery with MPC equations
rbar,99,6,7,123456
ENDDATA