Method setDomain

SRC/element/shell/ShellMITC4.cpp:359–419

set domain

//set domain
void  ShellMITC4::setDomain( Domain *theDomain ) 
{  
  int i, j ;
  static Vector eig(3) ;
  static Matrix ddMembrane(3,3) ;

  //node pointers
  for ( i = 0; i < 4; i++ ) {
     nodePointers[i] = theDomain->getNode( connectedExternalNodes(i) ) ;
     if (nodePointers[i] == 0) {
       opserr << "ShellMITC4::setDomain - no node " << connectedExternalNodes(i);
       opserr << " exists in the model\n";
     }
     const Vector &nodeDisp=nodePointers[i]->getTrialDisp();
     if (nodeDisp.Size() != 6) {
       opserr << "ShellMITC4::setDomain - node " << connectedExternalNodes(i);
       opserr << " NEEDS 6 dof - GARBAGE RESULTS or SEGMENTATION FAULT WILL FOLLOW\n";
     }       
     if (!m_initialzed) {
         init_disp[i][0] = nodeDisp(0);
         init_disp[i][1] = nodeDisp(1);
         init_disp[i][2] = nodeDisp(2);
         init_disp[i][3] = nodeDisp(3);
         init_disp[i][4] = nodeDisp(4);
         init_disp[i][5] = nodeDisp(5);
     }
  }

  //compute drilling stiffness penalty parameter
  const Matrix &dd = materialPointers[0]->getInitialTangent( ) ;

  //assemble ddMembrane ;
  for ( i = 0; i < 3; i++ ) {
      for ( j = 0; j < 3; j++ )
         ddMembrane(i,j) = dd(i,j) ;
  } //end for i 

  //eigenvalues of ddMembrane
  eig = LovelyEig( ddMembrane ) ;
  
  //set ktt 
  //Ktt = dd(2,2) ;  //shear modulus 
  Ktt = min( eig(2), min( eig(0), eig(1) ) ) ;
  //Ktt = dd(2,2);

  //basis vectors and local coordinates
  computeBasis( ) ;

  for (i = 0; i < 4; i++)
  {
    if (theDamping[i] && theDamping[i]->setDomain(theDomain, 8)) {
      opserr << "ShellMITC4::setDomain -- Error initializing damping\n";
      exit(-1);
    }
  }

  this->DomainComponent::setDomain(theDomain);