ArcLength1.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 ** Developed by: ** 00015 ** Frank McKenna (fmckenna@ce.berkeley.edu) ** 00016 ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** 00017 ** Filip C. Filippou (filippou@ce.berkeley.edu) ** 00018 ** ** 00019 ** ****************************************************************** */ 00020 00021 // $Revision: 1.2 $ 00022 // $Date: 2003/02/14 23:00:46 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/integrator/ArcLength1.h,v $ 00024 00025 00026 // File: ~/analysis/integrator/ArcLength1.h 00027 // 00028 // Written: fmk 00029 // Created: 07/98 00030 // Revision: A 00031 // 00032 // Description: This file contains the class definition for ArcLength1. 00033 // ArcLength1 is an algorithmic class for perfroming a static analysis 00034 // using the arc length scheme, that is within a load step the follwing 00035 // constraint is enforced: 00036 // i=1 delta U^T delta U + alpha^2 delta lambda^2 = delta s^2 00037 // i>1 dU^T delta U + alpha^2 dLambda delta lambda = 0 00038 // where dU is change in nodal displacements for step, dLambda is 00039 // change in applied load and ArcLength1 is a control parameter. 00040 // 00041 // What: "@(#) ArcLength1.h, revA" 00042 00043 #ifndef ArcLength1_h 00044 #define ArcLength1_h 00045 00046 #include <StaticIntegrator.h> 00047 00048 class LinearSOE; 00049 class AnalysisModel; 00050 class FE_Element; 00051 class Vector; 00052 00053 class ArcLength1 : public StaticIntegrator 00054 { 00055 public: 00056 ArcLength1(double arcLength, double alpha = 1.0); 00057 00058 ~ArcLength1(); 00059 00060 int newStep(void); 00061 int update(const Vector &deltaU); 00062 int domainChanged(void); 00063 00064 int sendSelf(int commitTag, Channel &theChannel); 00065 int recvSelf(int commitTag, Channel &theChannel, 00066 FEM_ObjectBroker &theBroker); 00067 00068 void Print(OPS_Stream &s, int flag =0); 00069 00070 protected: 00071 00072 private: 00073 double arcLength2; 00074 double alpha2; 00075 Vector *deltaUhat, *deltaUbar, *deltaU, *deltaUstep; 00076 Vector *phat; // the reference load vector 00077 double deltaLambdaStep, currentLambda; 00078 int signLastDeltaLambdaStep; 00079 }; 00080 00081 #endif 00082 |