OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
tf1.C
Go to the documentation of this file.
1 /*
2  *
3  * ##### ##### ###### ###### ### ###
4  * ## ## ## ## ## ## ## ### ##
5  * ## ## ## ## #### #### ## # ##
6  * ## ## ## ## ## ## ## ##
7  * ## ## ## ## ## ## ## ##
8  * ##### ##### ## ###### ## ##
9  *
10  *
11  * OOFEM : Object Oriented Finite Element Code
12  *
13  * Copyright (C) 1993 - 2013 Borek Patzak
14  *
15  *
16  *
17  * Czech Technical University, Faculty of Civil Engineering,
18  * Department of Structural Mechanics, 166 29 Prague, Czech Republic
19  *
20  * This library is free software; you can redistribute it and/or
21  * modify it under the terms of the GNU Lesser General Public
22  * License as published by the Free Software Foundation; either
23  * version 2.1 of the License, or (at your option) any later version.
24  *
25  * This program is distributed in the hope that it will be useful,
26  * but WITHOUT ANY WARRANTY; without even the implied warranty of
27  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
28  * Lesser General Public License for more details.
29  *
30  * You should have received a copy of the GNU Lesser General Public
31  * License along with this library; if not, write to the Free Software
32  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
33  */
34 
35 #include "tf1.h"
36 #include "timestep.h"
37 #include "mathfem.h"
38 #include "classfactory.h"
39 
40 namespace oofem {
42 
43 void
44 TF1 :: computeValueAt(FloatArray &answer, TimeStep *tStep, const FloatArray &coords, ValueModeType mode)
45 // Returns the value of the receiver at time and given position respecting the mode.
46 {
47  FloatArray cd(3);
48  double t, k, c, h, result;
49  answer.resize(1);
50 
51  k = 2.0e-4;
52  c = 1.0;
53  h = 200.0;
54 
55 
56  t = tStep->giveTargetTime();
57  for ( int i = 1; i <= coords.giveSize(); i++ ) {
58  cd.at(i) = coords.at(i);
59  }
60 
61  result = -1.e-5 - k *macbra( atan ( c *t + cd.at ( 2 ) - h ) );
62 
63  if ( ( mode == VM_Incremental ) && ( !tStep->isTheFirstStep() ) ) {
64  t = tStep->giveTargetTime() - tStep->giveTimeIncrement();
65  result -= -1.e-5 - k *macbra( atan ( c *t + cd.at ( 2 ) - h ) );
66  }
67 
68  answer.at(1) = result;
69 }
70 
73 {
74  return IRRT_OK;
75 }
76 } // end namespace oofem
virtual void computeValueAt(FloatArray &answer, TimeStep *tStep, const FloatArray &coords, ValueModeType mode)
Computes components values of load at given point - global coordinates (coordinates given)...
Definition: tf1.C:44
REGISTER_BoundaryCondition(BoundaryCondition)
double & at(int i)
Coefficient access function.
Definition: floatarray.h:131
ValueModeType
Type representing the mode of UnknownType or CharType, or similar types.
Definition: valuemodetype.h:78
double macbra(double x)
Returns the positive part of given float.
Definition: mathfem.h:115
double giveTargetTime()
Returns target time.
Definition: timestep.h:146
double giveTimeIncrement()
Returns solution step associated time increment.
Definition: timestep.h:150
bool isTheFirstStep()
Check if receiver is first step.
Definition: timestep.C:134
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: tf1.C:72
Class representing vector of real numbers.
Definition: floatarray.h:82
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
Class representing the general Input Record.
Definition: inputrecord.h:101
int giveSize() const
Returns the size of receiver.
Definition: floatarray.h:218
the oofem namespace is to define a context or scope in which all oofem names are defined.
Class representing solution step.
Definition: timestep.h:80
void resize(int s)
Resizes receiver towards requested size.
Definition: floatarray.C:631

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Jan 2 2018 20:07:31 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011