My title

00001 /* ****************************************************************** **
00002 **    OpenSees - Open System for Earthquake Engineering Simulation    **
00003 **          Pacific Earthquake Engineering Research Center            **
00004 **                                                                    **
00005 **                                                                    **
00006 ** (C) Copyright 1999, The Regents of the University of California    **
00007 ** All Rights Reserved.                                               **
00008 **                                                                    **
00009 ** Commercial use of this program without express permission of the   **
00010 ** University of California, Berkeley, is strictly prohibited.  See   **
00011 ** file 'COPYRIGHT'  in main directory for information on usage and   **
00012 ** redistribution,  and for a DISCLAIMER OF ALL WARRANTIES.           **
00013 **                                                                    **
00014 ** ****************************************************************** */
00015                                                                         
00016 // $Revision: 1.2 $
00017 // $Date: 2001/01/23 08:46:28 $
00018 // $Source: /usr/local/cvs/OpenSees/SRC/material/nD/J2AxiSymm.h,v $
00019 
00020 // Written: Ed "C++" Love
00021 //
00022 // J2AxiSymmetric isotropic hardening material class
00023 // 
00024 //  Elastic Model
00025 //  sigma = K*trace(epsilion_elastic) + (2*G)*dev(epsilon_elastic)
00026 //
00027 //  Yield Function
00028 //  phi(sigma,q) = || dev(sigma) ||  - sqrt(2/3)*q(xi) 
00029 //
00030 //  Saturation Isotropic Hardening with linear term
00031 //  q(xi) = simga_0 + (sigma_infty - sigma_0)*exp(-delta*xi) + H*xi 
00032 //
00033 //  Flow Rules
00034 //  \dot{epsilon_p} =  gamma * d_phi/d_sigma
00035 //  \dot{xi}        = -gamma * d_phi/d_q 
00036 //
00037 //  Linear Viscosity 
00038 //  gamma = phi / eta  ( if phi > 0 ) 
00039 //
00040 //  Backward Euler Integration Routine 
00041 //  Yield condition enforced at time n+1 
00042 //
00043 //  Send strains in following format :
00044 // 
00045 //     strain_vec = {   eps_00
00046 //                      eps_11
00047 //                 eps_22         
00048 //                    2 eps_01   }   <--- note the 2
00049 // 
00050 //  set eta := 0 for rate independent case
00051 //
00052 
00053 #include <iostream.h>
00054 #include <stdio.h> 
00055 #include <stdlib.h> 
00056 #include <math.h> 
00057 
00058 #include <Vector.h>
00059 #include <Matrix.h>
00060 
00061 #include <J2Plasticity.h>
00062 
00063 
00064 
00065 
00067 class J2AxiSymm : public J2Plasticity {
00068 
00070 
00071   public : 
00072 
00073   //null constructor
00074   J2AxiSymm( ) ;
00075 
00077   J2AxiSymm(   int    tag, 
00078                    double K,
00079                    double G,
00080                    double yield0,
00081                    double yield_infty,
00082                    double d,
00083                    double H,
00084                    double viscosity = 0 ) ;
00085 
00086 
00088   J2AxiSymm( int tag, double K, double G ) ;
00089 
00091   ~J2AxiSymm( ) ;
00092 
00094   NDMaterial* getCopy( ) ;
00095 
00097   const char* getType( ) const ;
00098 
00100   int getOrder( ) const ;
00101 
00103   int setTrialStrain( const Vector &strain_from_element) ;
00104 
00106   int setTrialStrain( const Vector &v, const Vector &r ) ;
00108   int setTrialStrainIncr( const Vector &v ) ;
00110   int setTrialStrainIncr( const Vector &v, const Vector &r ) ;
00111 
00113   const Vector& getStrain( ) ;
00114 
00116   const Vector& getStress( ) ;
00117 
00119   const Matrix& getTangent( ) ;
00120 
00122   int setTrialStrain(const Tensor &v) ;
00124   int setTrialStrain(const Tensor &v, const Tensor &r) ;    
00126   int setTrialStrainIncr(const Tensor &v) ;
00128   int setTrialStrainIncr(const Tensor &v, const Tensor &r) ;
00130   const Tensor& getTangentTensor( ) ;
00131 //jeremic@ucdavis.edu 22jan2001  const Tensor& getStressTensor( ) ;
00132 //jeremic@ucdavis.edu 22jan2001  const Tensor& getStrainTensor( ) ;  
00133 
00135   int sendSelf(int commitTag, Channel &theChannel) ;  
00137   int recvSelf(int commitTag, Channel &theChannel, 
00138                FEM_ObjectBroker &theBroker ) ;
00139   
00141   private :
00142 
00143   //static vectors and matrices
00144   static Vector strain_vec ;     
00145   static Vector stress_vec ;     
00146   static Matrix tangent_matrix ; //material tangent in matrix notation
00147 
00148           
00149 } ; //end of J2AxiSymm declarations
00150 
00151 
00152 ÿ