Subversion Repositories OpenSees

#ifndef RAFourSteelPCPlaneStress_h

#define RAFourSteelPCPlaneStress_h

// File: RAFourSteelPCPlaneStress.h

//

// Written: ALaskar

// Created: 2007.09

//

// Description: This file contains the class definition for 

// RAFourSteelPCPlaneStress material. 

// RC material will contain steel in four directions

// For Detailed explanation of the model, please refer to the book

// entitled "Unified Theory of Concrete Structures,"

// by Thomas T.C. Hsu and Y.L. Mo, John Wiley & Sons, April 2010.

#include <NDMaterial.h>

#include <UniaxialMaterial.h>

#include <Matrix.h>

#include <Vector.h>

#include <ID.h>

#include <Tensor.h>

class RAFourSteelPCPlaneStress : public NDMaterial

{

  public:

  RAFourSteelPCPlaneStress ( int      tag,

                             double   RHO,

                             UniaxialMaterial *t1,

                             UniaxialMaterial *t2,

                             UniaxialMaterial *s1,

                             UniaxialMaterial *s2,

                             UniaxialMaterial *c1,

                             UniaxialMaterial *c2,

                             double   ANGLE1,

                             double   ANGLE2,

                             double   ANGLE3,

                             double   ANGLE4,

                             double   ROU1,

                             double   ROU2,

                             double   ROU3,

                             double   ROU4,

                             double     PSTRAIN1,

                             double     PSTRAIN2,

                             double   FPC,

                             double   FY1,

                             double     FY2,

                             double   E,

                             double   EPSC0);                                    

  RAFourSteelPCPlaneStress();

  ~RAFourSteelPCPlaneStress();                           

  double getRho(void);

  int setTrialStrain(const Vector &v); // the one is used

  int setTrialStrain(const Vector &v, const Vector &r);

  int setTrialStrainIncr(const Vector &v);

  int setTrialStrainIncr(const Vector &v, const Vector &r);

  const Matrix &getTangent(void);

  const Matrix &getInitialTangent(void) {return this->getTangent();};

  const Vector &getStress(void);

  const Vector &getStrain(void);

  const Vector &getCommittedStress(void);

  const Vector &getCommittedStrain(void);    

  int commitState(void);

  int revertToLastCommit(void);

  int revertToStart(void);

  NDMaterial *getCopy(void);

  NDMaterial *getCopy(const char *type);

  void Print(OPS_Stream &s, int flag = 0);

  int sendSelf(int commitTag, Channel &theChannel);

  int recvSelf(int commitTag, Channel &theChannel, FEM_ObjectBroker &theBroker);

  const char *getType(void) const { return "PlaneStress"; };

  int getOrder(void) const { return 3;};

 protected:

 private:

  double   rho;

  UniaxialMaterial **theMaterial; // pointer of the materials 

  Response **theResponses;

  double   angle1;    // angel of the first steel layer to x coordinate 

  double   angle2;    // angel of the second steel layer to x coordinate

  double   angle3;    // angel of the third steel layer to x coordinate 

  double   angle4;    // angel of the forth steel layer to x coordinate

  double   rou1;      // steel ratio of the first steel layer

  double   rou2;      // steel ratio of the second steel layer

  double   rou3;      // steel ratio of the third steel layer

  double   rou4;      // steel ratio of the forth steel layer

  double        pstrain1;       // prestressing strain in first tendon layer

  double        pstrain2;       // prestressing strain in second tendon layer

  double   fpc;       // compressive strength of the concrete

  double   fy1;       // yield stress of the bare tendons

  double         fy2;           // yield stress of the bare steel bars

  double   E0;        // young's modulus of the steel

  double   epsc0;     // compressive strain of the concrete

  double   Tstress[3];  // Trial stresses

  double   lastStress[3];  // Last committed stresses, added for x, k

  int      steelStatus;  // check if steel yield, 0 not yield, 1 yield

  int      dirStatus;    // check if principle direction has exceed 90 degree, 1 yes, 0 no

  double   citaStrain;      // principle strain direction

  double   citaStress;     // principle stress direction

  double   miu12;        // Hsu/Zhu ratio

  double   miu21;        // Hsu/Zhu ratio

  double   G12;

  // for damgage factor D=1-0.4*epslonC'/epslon0; epslon0=0.002

  // Trial values

  int TOneReverseStatus;         // Trial reverse status for concrete One, 1 reverse, 0 no

  double TOneNowMaxComStrain;

  double TOneLastMaxComStrain;

  int TTwoReverseStatus;         // Trial reverse status for concrete Two, 1 reverse, 0 no

  double TTwoNowMaxComStrain;

  double TTwoLastMaxComStrain;

  // Converged values

  int COneReverseStatus;         // Converged reverse status for concrete One, 1 reverse, 0 no

  double COneNowMaxComStrain;

  double COneLastMaxComStrain;

  int CTwoReverseStatus;         // Converged reverse status for concrete Two, 1 reverse, 0 no

  double CTwoNowMaxComStrain;

  double CTwoLastMaxComStrain;

  double DDOne; // damage factor for concrete One

  double DDTwo; // damage factor for concrete Two

  double tt1;

  double tt2;

  double xxx;

  double kkk;

  double beta1;

  double beta2;

  Vector strain_vec;

  Vector stress_vec;   

  Matrix tangent_matrix;

  int determineTrialStress(void);

};

#endif