00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038 #ifndef MinUnbalDispNorm_h
00039 #define MinUnbalDispNorm_h
00040
00041 #include <StaticIntegrator.h>
00042
00043 #include <LinearSOE.h>
00044 #include <AnalysisModel.h>
00045 #include <FE_Element.h>
00046 #include <Vector.h>
00047
00048 #define SIGN_LAST_STEP 1
00049 #define CHANGE_DETERMINANT 2
00050
00061 class MinUnbalDispNorm : public StaticIntegrator
00062 {
00063 public:
00072 MinUnbalDispNorm(double lambda1, int specNumIterStep,
00073
00074 double dlambda1min, double dlambda1max,
00075 int signFirstStepMethod = SIGN_LAST_STEP);
00076
00080 ~MinUnbalDispNorm();
00081
00082
00086 int newStep(void);
00087
00091 int update(const Vector &deltaU);
00092
00094 int domainChanged(void);
00095
00096
00098 int sendSelf(int commitTag, Channel &theChannel);
00100 int recvSelf(int commitTag, Channel &theChannel,
00101 FEM_ObjectBroker &theBroker);
00102
00104 void Print(ostream &s, int flag =0);
00105
00106 protected:
00107
00108 private:
00110 double dLambda1LastStep;
00111 double specNumIncrStep, numIncrLastStep;
00112
00113 Vector *deltaUhat, *deltaUbar, *deltaU, *deltaUstep;
00114 Vector *phat;
00115
00116 double deltaLambdaStep, currentLambda;
00117 int signLastDeltaLambdaStep;
00118 double dLambda1min, dLambda1max;
00119 double signLastDeterminant;
00121 int signFirstStepMethod;
00122 };
00123
00124 #endif
00125