$ beigp1   from rlh April 10, 2000
NASTRAN NLINES=9999
TIME	999
SOL	106
DIAG	8
Include 'beigpa.v707' $
CEND
TITLE = SMALL UNSYMMETRIC BUCKLING PROBLEM        beigp1
Subtitle = large enough to use clan
$ PCRIT = (PI/2)**2 * EI / L**2 = 246.74...
  SPC	= 200
$ try two subcases, per rlh suggestion
subcase 1
  load = 150
  NLPAR	= 400
  DISPL	= ALL
  OLOAD	= ALL
  NLSTR	= NONE
Subcase 2
  LOAD	= 300
  NLPAR	= 400
  PARAM	BUCKLE 2
  METHO	= 500
  DISPL	= ALL
  OLOAD	= ALL
  NLSTR	= NONE
BEGIN BULK
GRDSET							345
GRID	1		0.
=	*1	=	*1.
$=9 $ 4x size of model
=39
CBEAM	1	100	1	2		1.
=	*1	=	*1	*1	=	=
$=8
=38
PBEAM	100	100	.01	.001	.001		.002
MAT1	100	1.+7		.3
SPC1	200	126	1
$FORCE	300	11		200.	-1.
FORCE1	150	41	100.	11	10
FORCE1	300	41	200.	11	10
NLPARM,	400	2,	,	iter,	1,	,	,	YES
$NLPARM	400	2		AUTO				YES
EIGRL	500			1
PARAM	LGDISP	1
ENDDATA

COMPILE MODERS
ALTER 'READ *KXX,MXX'(,-1) $
PARAML	MXX//'TRAILER'/1/S,N,WD1 $
PARAML	MXX//'TRAILER'/2/S,N,WD2 $
PARAML	MXX//'TRAILER'/3/S,N,WD3 $
PARAML	MXX//'TRAILER'/4/S,N,WD4 $
PARAML	MXX//'TRAILER'/5/S,N,WD5 $
MESSAGE	//' MXX MATRIX SYMMETRY CHECK'/WD1/WD2/WD3/WD4/WD5 $
TRNSP	MXX/TMXX $
ADD5	MXX,TMXX,,,/UMXX/.5/-.5 $
MATMOD	MXX,UMXX,,,,/MSCR,/20 $
MATGPR,	,,,MXX//'H'/'H' $
MATGEN,	/IXX/1/V,N,WD1 $
READ    MXX,IXX,,,,,,,,,,,,/LAMM,PHIM,MM,OEIGM,,/
        'MODES'/S,N,NEIGV/1///-1/'HOU'////99 $
OFP     OEIGM,LAMM// $
$
PARAML	KXX//'TRAILER'/1/S,N,WD1 $
PARAML	KXX//'TRAILER'/2/S,N,WD2 $
PARAML	KXX//'TRAILER'/3/S,N,WD3 $
PARAML	KXX//'TRAILER'/4/S,N,WD4 $
PARAML	KXX//'TRAILER'/5/S,N,WD5 $
MESSAGE	//' KXX MATRIX SYMMETRY CHECK'/WD1/WD2/WD3/WD4/WD5 $
TRNSP	KXX/TKXX $
ADD5	KXX,TKXX,,,/UKXX/.5/-.5 $
MATMOD	KXX,UKXX,,,,/KSCR,/20 $
MATGPR,	,,,KXX//'H'/'H' $
READ    KXX,IXX,,,,,,,,,,,,/LAMK,PHIK,MK,OEIGK,,/
        'MODES'/S,N,NEIGV/1///-1/'HOU'////99 $
OFP     OEIGK,LAMK// $