hexagon logo

How to extract stiffness at a node in f06 file or .pch file for a particular loadstep ?

How to extract stiffness at a node in f06 file or .pch file for a particular loadstep ?
Parents
  • The differential stiffness is the stiffness that results from including the higher-order terms (non-linear terms) of the strain-displacement relations. This is developed in non-linear solutions (from small and large displacement effects), as well as, linear buckling and pre-stressed normal modes (from small displacement effects).
     
    For linear statics, you can get the stiffness for any given node using the alter found in the following article:
     
     
    Here is the alter:
     
    MALTER 'MALTER.*KGG.*BGG.*MGG.*K4GG'      $
    MATMOD  EQEXINS,USET0,SILS,CASECC,,/VECX,/17//1 $
    PARTN  KGG,VECX,VECX/,,,KEEP/0/0        $
    DIAGONAL KEEP/KEEPD//              $  
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'PRINT K for Case Control, user-defined GRIDs'$
    MESSAGE //'  '                 $
    MESSAGE //'Diagonal ONLY'            $
    MATPRN  KEEPD//                 $
    MESSAGE //'  '                 $
    MESSAGE //'Full ONLY'              $
    MATPRN  KEEP//                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    ENDALTER                     $    
     
    This requires you at add the following to your Case Control portion of the Nastran input file:
     
    $ define the grids that you want the stiffness for in this SET, for example, Grid 30 will be printed
    set 99=30
    PARTN=99
     
    I'll attach a sample input file. Note that this same ALTER and SET/PARTN will work for most solution sequences, such as SOL101 (linear statics), SOL103 (normal modes) and SOL105 (linear buckling) and many others... but it is always working with the KGG stiffness matrix, which is prior to any reductions.
     
     
     
Reply
  • The differential stiffness is the stiffness that results from including the higher-order terms (non-linear terms) of the strain-displacement relations. This is developed in non-linear solutions (from small and large displacement effects), as well as, linear buckling and pre-stressed normal modes (from small displacement effects).
     
    For linear statics, you can get the stiffness for any given node using the alter found in the following article:
     
     
    Here is the alter:
     
    MALTER 'MALTER.*KGG.*BGG.*MGG.*K4GG'      $
    MATMOD  EQEXINS,USET0,SILS,CASECC,,/VECX,/17//1 $
    PARTN  KGG,VECX,VECX/,,,KEEP/0/0        $
    DIAGONAL KEEP/KEEPD//              $  
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'PRINT K for Case Control, user-defined GRIDs'$
    MESSAGE //'  '                 $
    MESSAGE //'Diagonal ONLY'            $
    MATPRN  KEEPD//                 $
    MESSAGE //'  '                 $
    MESSAGE //'Full ONLY'              $
    MATPRN  KEEP//                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    MESSAGE //'  '                 $
    ENDALTER                     $    
     
    This requires you at add the following to your Case Control portion of the Nastran input file:
     
    $ define the grids that you want the stiffness for in this SET, for example, Grid 30 will be printed
    set 99=30
    PARTN=99
     
    I'll attach a sample input file. Note that this same ALTER and SET/PARTN will work for most solution sequences, such as SOL101 (linear statics), SOL103 (normal modes) and SOL105 (linear buckling) and many others... but it is always working with the KGG stiffness matrix, which is prior to any reductions.
     
     
     
Children
No Data