J2PlaneStrain.hGo to the documentation of this file.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.6 $ 00017 // $Date: 2006/08/04 18:18:37 $ 00018 // $Source: /usr/local/cvs/OpenSees/SRC/material/nD/J2PlaneStrain.h,v $ 00019 00020 #ifndef J2PlaneStrain_h 00021 #define J2PlaneStrain_h 00022 00023 // Written: Ed "C++" Love 00024 // 00025 // J2PlaneStrain isotropic hardening material class 00026 // 00027 // Elastic Model 00028 // sigma = K*trace(epsilion_elastic) + (2*G)*dev(epsilon_elastic) 00029 // 00030 // Yield Function 00031 // phi(sigma,q) = || dev(sigma) || - sqrt(2/3)*q(xi) 00032 // 00033 // Saturation Isotropic Hardening with linear term 00034 // q(xi) = simga_0 + (sigma_infty - sigma_0)*exp(-delta*xi) + H*xi 00035 // 00036 // Flow Rules 00037 // \dot{epsilon_p} = gamma * d_phi/d_sigma 00038 // \dot{xi} = -gamma * d_phi/d_q 00039 // 00040 // Linear Viscosity 00041 // gamma = phi / eta ( if phi > 0 ) 00042 // 00043 // Backward Euler Integration Routine 00044 // Yield condition enforced at time n+1 00045 // 00046 // Send strains in following format : 00047 // 00048 // strain_vec = { eps_00 00049 // eps_11 00050 // 2 eps_01 } <--- note the 2 00051 // 00052 // set eta := 0 for rate independent case 00053 // 00054 00055 #include <stdio.h> 00056 #include <stdlib.h> 00057 #include <math.h> 00058 00059 #include <Vector.h> 00060 #include <Matrix.h> 00061 00062 #include <J2Plasticity.h> 00063 00064 class J2PlaneStrain : public J2Plasticity { 00065 00066 //-------------------Declarations------------------------------- 00067 00068 public : 00069 00070 //null constructor 00071 J2PlaneStrain( ) ; 00072 00073 //full constructor 00074 J2PlaneStrain( int tag, 00075 double K, 00076 double G, 00077 double yield0, 00078 double yield_infty, 00079 double d, 00080 double H, 00081 double viscosity = 0 ) ; 00082 00083 00084 //elastic constructor 00085 J2PlaneStrain( int tag, double K, double G ) ; 00086 00087 //destructor 00088 ~J2PlaneStrain( ) ; 00089 00090 //make a clone of this material 00091 NDMaterial* getCopy( ) ; 00092 00093 //send back type of material 00094 const char* getType( ) const ; 00095 00096 const char *getClassType(void) const {return "J2PlaneStrain";}; 00097 00098 //send back order of strain in vector form 00099 int getOrder( ) const ; 00100 00101 00102 //get the strain and integrate plasticity equations 00103 int setTrialStrain( const Vector &strain_from_element) ; 00104 00105 //unused trial strain functions 00106 int setTrialStrain( const Vector &v, const Vector &r ) ; 00107 int setTrialStrainIncr( const Vector &v ) ; 00108 int setTrialStrainIncr( const Vector &v, const Vector &r ) ; 00109 00110 //send back the strain 00111 const Vector& getStrain( ) ; 00112 00113 //send back the stress 00114 const Vector& getStress( ) ; 00115 00116 //send back the tangent 00117 const Matrix& getTangent( ) ; 00118 const Matrix& getInitialTangent( ) ; 00119 00120 //this is mike's problem 00121 int setTrialStrain(const Tensor &v) ; 00122 int setTrialStrain(const Tensor &v, const Tensor &r) ; 00123 int setTrialStrainIncr(const Tensor &v) ; 00124 int setTrialStrainIncr(const Tensor &v, const Tensor &r) ; 00125 const Tensor& getTangentTensor( ) ; 00126 00127 //jeremic@ucdavis.edu 22jan2001 const Tensor& getStressTensor( ) ; 00128 //jeremic@ucdavis.edu 22jan2001 const Tensor& getStrainTensor( ) ; 00129 00130 //swap history variables 00131 int commitState( ) ; 00132 int revertToLastCommit( ) ; 00133 int revertToStart( ) ; 00134 00135 //sending and receiving 00136 int sendSelf(int commitTag, Channel &theChannel) ; 00137 int recvSelf(int commitTag, Channel &theChannel, 00138 FEM_ObjectBroker &theBroker ) ; 00139 00140 00141 private : 00142 00143 //static vectors and matrices 00144 static Vector strain_vec ; //strain in vector notation 00145 static Vector stress_vec ; //stress in vector notation 00146 static Matrix tangent_matrix ; //material tangent in matrix notation 00147 } ; //end of J2PlaneStrain declarations 00148 00149 #endif
Generated on Mon Oct 23 15:05:14 2006 for OpenSees by |
opensees-support @ berkeley.edu
Supported by the National Science Foundation
©2006, UC Regents