help to output the stiffness matrix in global coordinates

[kirin]

hello

I followed the older thread(

http://www.eng-tips.com/viewthread.cfm?qid=62686#p...)

but I can't output the matrix.

My steps

1) From the solution sequence (DIAG 14) I found the first DMAP module that could give me results on MGG and KGG

DMAP-DMAP INSTRUCTION
NO.
1 SUBDMAP SEMODES $ SOL 103 - SUPERELEMENT NORMAL MODES
2 $ IFP, PHASE0
2 TYPE DB AXIC,PVT,BULK,GEOM1Q,TIMSIZ,K4JJ,YG,RMG,CASESX,USET0,GEOM1S,
VIEWTB,GPSNTS,BGPDTVU,IBULK,EHT,EHTA,CASERS,PVTBC,EDOM,DRG,
MEDGE,MFACE,GDNTAB,PELSETS,ELEMVOL,PVAL,ERROR,TEMPALL,ERRORG,
$ PHASE 1
KJJ,MJJ,BJJ,USET,GOT,KAA,DYNAMICS,DYNAMIC,KGG,PCDBS,AGG,

… ETC …

Then at line 96 I found an instruction containing the input (GPLS, USET0, SILS) and output data blocks (MGG, KGG) I need:

96 CALL PHASE1DR,
SLIST ,EMAP ,DITS ,MATPOOLS,EDTS ,
DYNAMICS,GEOM4S ,EPTS ,CASECC ,MPTS ,
CASES ,GPLS ,EQEXINS ,CASESX ,CSTMS ,
BGPDTS ,SILS ,ECTS ,MAPS ,PVTS ,
,VGFS ,GDNTAB ,CONTACT ,DEQATN ,
DEQIND ,GPSNTS ,ECTAS ,MEDGE ,EHT ,
RMG ,USET0 ,YG ,SLT ,ETT ,
,EPTA ,EHTA ,SCSTM ,PVTBC ,
, ,EQEXIN , , ,
, ,AG0 ,GEQMAP , ,
,SPCPART , , , ,
DRG ,VGQ ,XSOP2DIR,GEOM1QS , ,
,GEOM2S ,GEOM4QS ,SELIST ,POSTCDB ,
PCDB ,IBULK , , , ,
MMCDB ,EXTDB , , ,MP1 ,
MPD1 ,MP3 , ,EDOM ,YMU ,
ACTNODS ,YSD /
EST ,GPECT ,VELEM ,KELM ,KDICT ,
MELM ,MDICT ,KJJ ,MJJ ,BJJ ,
K4JJ , , ,KGG ,GM ,
USET ,KSF ,KFS ,KSS , ,
KAA ,MPHVZ ,LOO , ,GOT ,
KLL ,LLL ,ULL ,DAR ,MGG ,
BGG ,K4GG ,MAA ,CMLAMA1 ,MR ,
MLR ,MLL ,MRR ,CMPHA ,CMPHO ,
CMLAMA ,CEE ,GOQ ,MEE ,MAR ,
MEA ,GGE ,BAA ,K4AA , ,
, , , , ,
,UOX ,ESTL ,ESTNL ,AGG ,
BTOPO ,KMG ,MMG ,MSF ,BMG ,
BSF ,K4MG ,K4SF , ,ORSEQ ,
LRSEQ ,VGFD ,DITID , , ,
, , , , ,
,BK4ZZ ,DKFF , ,DMFF ,
, , ,MSS ,BSS ,
K4SS ,BELM ,BDICT , ,MP2SLAB ,
, , , , ,

I included the following DMAP alter in the executive control deck:

$
COMPILE SEMODES $
ALTER 96 $
MATGPR GPLS,USET0,SILS,MGG//'G' $
MATGPR GPLS,USET0,SILS,KGG//'G' $
$

After all steps ,I analysis it

But I could not output the matrices results in the .f06 file