Statdivr.rdm

PRINTING AND PLOTTING EIGENVECTORS FOR AEROELASTIC STATIC DIVERGENCE, 
statdiva.v707 

M, A. Gockel, November 13, 1998

1. Background

SOL 144 (static aeroelasticity) has a divergence analysis capability as an 
option.  While both eigenvalues and eigenvectors are calculated, only the 
eigenvalues are available as outputs at present.  This alter package allows 
printing and plotting of the eigenvectors.

2. Input

Place 

include 'statdiva.v707'

before the CEND command.  

The alter package imposes some limits on the types of solutions that may be 
requested.  If eigenvectors for more than one Mach number are desired a subcase 
must be prepared for each Mach number that includes a DIVERG = SID command.  
Each DIVERG entry selected must contain a request for only one NROOT (field 3), 
and one Mach number (field 4).  This will result in only one eigenvector per 
subcase.  It will be the lowest one calculated.  More eigensolutions may be 
requested on the EIGC entry, but only the vector for the lowest eigenvalue will 
be output.

If only one divergence solution is required subcases are not necessary.  
However, the limitations on NROOT and Mach number still apply.

3. OUTPUT

Each subcase (eigenvector-subcase pair) is labeled as a static deformation 
vector.  The subcase ID corresponds to the input values.  The eigenvalue is not 
output with the eigenvector.  They are listed earlier in the output.

4. LIMITATIONS

Do not request trim solutions in subcases, only divergence solutions.

Do not request more than one NROOT or Mach number per subcase. 

Element output data is not provided.

5. Example Problem

This problem is a modification of a problem discussed in the Aero Guide under 
the name HA144H.  The full input file is stored under the name statdiv1.dat.  
Selected parts of the input file are listed below

$ file statdiv1.dat   demo statdiva.v707 
. . . [ entries with no changes are deleted]
include 'statdiva.v707' $ print and plot divergence eigenvectors
CEND
title = demo use of statdiva.v707 alter		statdiv1
$TITLE = EXAMPLE HA144B: BAH JET TRANSPORT WING DYNAMIC ANALYSIS HA144H
. . .
$  diverg = 100
cmethod = 100  $ clan
subcase 11  $ add subcase for each mach number.   
  diverg=101 $ one vector per mach number 
subcase 12
  diverg=102
BEGIN BULK
param, post, 0 $ put on xdb for patran plots
. . .
$   use soft bars to draw outlines.  Stiffness is defined on genel
cbar,	13,	13,	1,	3,	0.	0.	1.
=,	24,	=,	2,	4,	==
=,	35	=,	3	5,	==
=,	46,	=,	4,	6,	==
=,	57,	=,	5,	7,	==
=,	68,	=,	6,	8,	==
=,	79,	=,	7,	9,	== $
=,	810,	=,	8,	10,	== $1.e-6,	,	
pbar,	13,	13,	1. $ axial only
mat1,	13,	1.e-6,	,	.3 $ very soft, massless.  for viewing only
. . .
$diverg	100	5	0.0 $ remove old diverg entry
$		NROOT   MACH NUMBER  ONLY ONE OF EACH ALLOWED PER DEIVERGE ENTRY
diverg, 101,	1,	0.0 $ 
diverg,	102,	1,	0.1 
EIGC	100	CLAN	MAX				5
ENDDATA