lfssolver.cpp

来自「Finite element program for mechanical pr」· C++ 代码 · 共 49 行

#include "lfssolver.h"
#include "global.h"
#include "globmat.h"
#include "loadcase.h"
#include "gmatrix.h"
#include "mechprint.h"
#include "flsubdom.h"
#include "elemswitch.h"

/**
   function solves linear problem with floating subdomains
   it may be used for pull-out tests
   
   JK, 5.11.2005
*/
void solve_linear_floating_subdomains ()
{
  long i,lcid;
  double *lhs,*rhs;
  //flsubdom *fsd;
  
  lcid=0;
  
  Fsd = new flsubdom;
  
  lhs=Lsrs->give_lhs (lcid);
  rhs=Lsrs->give_rhs (lcid);
  
  mefel_right_hand_side (lcid,rhs);
  
  stiffness_matrix (lcid);
  
  Fsd->initialization (Ndofm,Mp->ense,rhs);
  
  Fsd->solve_lin_alg_system (lhs,rhs);
  
  print_init(-1, "wt");    
  for (i=0;i<Lsrs->nlc;i++){
    compute_ipstrains (i);
    compute_ipstresses (i);
    compute_req_val (i);
    print_step(i, 0, 0.0, NULL);    
  }
  print_close();
  
  
}