![]() |
FemProcess
The GeMA Fem Process Plugin
|
Basic class for solving Non linear FEM problems. More...
#include <gmpFemNLSolver.h>
Public Member Functions | |
GmpFemNLSolver (GmElementMesh *, GmSimulationData *, const QList< GmpFemPhysics * > &, GmNumSolver *, const GmpFemSolverOptions &, const GmLogCategory &) | |
void | enableFastUpdate (bool mode) |
virtual bool | init () |
Prepares the solver for assembling matrices by creating the assembler object and allocating the needed memory. This operations only need to be done once if the solver is used in a loop. More... | |
virtual bool | runStep (double dt, double *newt, bool *conv, double *err, int *niter, bool eval=false) |
virtual bool | geostatic (QString svId="u", QString gaId="Em") |
The geostatic implementation (femNLSolver.geostatic) This function manages the geostatic process. | |
![]() | |
GmpFemSolver (GmElementMesh *mesh, GmSimulationData *simulation, const QList< GmpFemPhysics * > &physics, GmNumSolver *solver, const GmpFemSolverOptions &options, const GmLogCategory &logger) | |
Constructor. Expects to receive as parameters the mesh we are acting upon, the list of physics objects that will cooperate to create the global stiffness matrix and the solver used to solve the resulting linear system. More... | |
virtual | ~GmpFemSolver (void) |
Destructor. | |
bool | run () |
Execute the process. Returns true on success. | |
bool | calcLinearResidual (double *rnorm, double *maxNodeDiff, double *avgNodeDiff) |
When solving a non linear system by repeated iterations, this function aims to calculate the residual of the actual solution. More... | |
void | setMatrixCombinerObject (const GmpFemAssemblerMatrixCombiner *combiner) |
Sets the combiner object that will be used by the assembler to merge elemental data from multiple matrices (eg. C and K) into an equivalent matrix. More... | |
void | setVectorCombinerObject (const GmpFemAssemblerVectorCombiner *combiner) |
Sets the combiner objects that will be used by the assembler to merge elemental data from multiple vectors (eg. Fe and Fi) into an equivalent vector. More... | |
bool | update (int mode, QString &err) |
Informs the solver of an external change in model dofs, elements or BCs that must be forwarded to the assembler. See the documentation of GmpFemAssembler::update() for more details. | |
virtual bool | stateAboutToBeSaved (GmStateDump *state) |
virtual bool | stateSaved (GmStateDump *state) |
virtual bool | stateAboutToBeLoaded (GmStateDump *state) |
virtual bool | stateLoaded (GmStateDump *state) |
Informs physics objects that a load operation was completed giving them an opportunity to adjust its state in case the loaded data is not exactly what the internal state needs to be. It will also call derivedResults() so that previous results can be recovered. In that call, the nonLinearSolver parameter is obtained from the stateLoadedCalcDerivedResultsNonLinear() virtual function. | |
![]() | |
virtual const QMetaObject * | metaObject () const const |
virtual void * | qt_metacast (const char *) |
virtual int | qt_metacall (QMetaObject::Call, int, void **) |
QObject (QObject *parent) | |
virtual bool | event (QEvent *e) |
virtual bool | eventFilter (QObject *watched, QEvent *event) |
QString | objectName () const const |
void | setObjectName (const QString &name) |
bool | isWidgetType () const const |
bool | isWindowType () const const |
bool | signalsBlocked () const const |
bool | blockSignals (bool block) |
QThread * | thread () const const |
void | moveToThread (QThread *targetThread) |
int | startTimer (int interval, Qt::TimerType timerType) |
int | startTimer (std::chrono::milliseconds time, Qt::TimerType timerType) |
void | killTimer (int id) |
T | findChild (const QString &name, Qt::FindChildOptions options) const const |
QList< T > | findChildren (const QString &name, Qt::FindChildOptions options) const const |
QList< T > | findChildren (const QRegExp ®Exp, Qt::FindChildOptions options) const const |
QList< T > | findChildren (const QRegularExpression &re, Qt::FindChildOptions options) const const |
const QObjectList & | children () const const |
void | setParent (QObject *parent) |
void | installEventFilter (QObject *filterObj) |
void | removeEventFilter (QObject *obj) |
QMetaObject::Connection | connect (const QObject *sender, const char *signal, const char *method, Qt::ConnectionType type) const const |
bool | disconnect (const char *signal, const QObject *receiver, const char *method) const const |
bool | disconnect (const QObject *receiver, const char *method) const const |
void | dumpObjectTree () |
void | dumpObjectInfo () |
void | dumpObjectTree () const const |
void | dumpObjectInfo () const const |
bool | setProperty (const char *name, const QVariant &value) |
QVariant | property (const char *name) const const |
QList< QByteArray > | dynamicPropertyNames () const const |
void | destroyed (QObject *obj) |
void | objectNameChanged (const QString &objectName) |
QObject * | parent () const const |
bool | inherits (const char *className) const const |
void | deleteLater () |
Protected Member Functions | |
virtual bool | addStateItemsToGroup (GmStateDump *state, int groupId) |
Adds to 'state' the data items that should be saved for this FEM process. Should probably be overriden by derived classes to add their own requirements. More... | |
virtual bool | fillStateControlMapData (QVariantMap *map) |
Method called for filling the state dump control regustered in addStateItemsToGroup() | |
virtual bool | stateControlMapDataLoaded (QVariantMap *map) |
Method called when the state dump control registered in addStateItemsToGroup() has been loaded. | |
virtual bool | stateLoadedCalcDerivedResultsNonLinear () |
Returns the value of the "nonLinearSolver" parameter that will be passed to the calcDerivedResults() call from the standard implementation of stateLoaded() | |
virtual bool | fillStateControlMapDataStaticSolver (QVariantMap *map) |
Virtual method to fill the map with the needed continuation method data. | |
virtual bool | stateControlMapDataLoadedStaticSolver (QVariantMap *map) |
Virtual method to Load the scalar attributes needed by continuation method. | |
virtual bool | fillStateControlMapDataTransientSolver (QVariantMap *map) |
Virtual method to load the scalar attributes needed by transient solver. | |
virtual bool | stateControlMapDataLoadedTransientSolver (QVariantMap *map) |
Virtual method to fill the map with the needed transient solver data. | |
virtual bool | initElementSets (GmNumSolver *) |
Helpper function used to initialize the matrix / vector sets with the definition of the needed types and also the relationships of those sets with the assembler. More... | |
virtual void | setElementSets (bool, bool, bool, bool, bool) |
virtual bool | applyFixedBoundaryConditions (int *) |
Apply Dirichlet (fixed) boundary conditions to the global equation system. More... | |
virtual bool | IntialFixedBoundaryConditions (int *) |
virtual bool | applyNonzeroFixedBCsbyModification (int *) |
virtual bool | isDynamic () const |
virtual bool | isNonlinear () const |
virtual bool | convergence () const |
analysis | |
virtual bool | convergenceN (GmVector &_erL, GmVector &_erD) const |
virtual bool | convergenceDC (GmVector &_er) const |
virtual void | update () |
virtual void | restore () |
virtual bool | LoadError (GmVector &_er) const |
virtual bool | DispError (GmVector &_er) const |
virtual bool | VelocityError (GmVector &_er) const |
virtual bool | FluxError (GmVector &_er) const |
virtual bool | FluxMagnitude (GmVector &_qa) const |
virtual bool | MaxErr (GmVector, double &ermax, double &tolmax) |
Maximun erro in the current itarion. | |
virtual bool | MaxVal (GmVector &array, GmVector) const |
Maximun value of the array. More... | |
virtual bool | PhygroupName (int &kgroup, QString var) const |
Physical group. More... | |
virtual bool | divergingIncrement (double ermax, double &ermax_1, double &ermax_2) |
diverging Increment | |
virtual bool | adaptativeTimeAbqs (bool conv, int iter, double dt, double &newdt, double errmax, int inc, bool iter4) |
Adaptative time increment based on abaqus std. | |
virtual bool | adaptativeTime (bool conv, int iter, double dt, double &newdt, double err, double tol) |
Adaptative time increment. | |
virtual void | saveInternalForceAttributes (GmVector &dv) |
Save internal force attibutes such as Reaction and Volume flux. | |
virtual void | setLoad () |
virtual void | setTime () |
virtual void | setState () |
virtual void | setTransientSystem (double) |
evaluates transient system | |
virtual void | setTransientSystem2 (double) |
virtual GmpFemPhysics::FemResultType | updateResidualForceAndStiffnessMatrix (double l) |
Updates the stiffness matrix. More... | |
virtual GmpFemPhysics::FemResultType | fillStiffnessMatrixBroyden () |
Evaluates the stiffness matrix using Broyden's method. More... | |
virtual GmpFemPhysics::FemResultType | fillStiffnessMatrixBFGS () |
Evaluates the stiffness matrix using BFGS method. More... | |
virtual double | evalLocalError (GmVector, GmVector) |
calc local error | |
virtual double | evalTimeIncrementSize (int, double, double, GmVector &) |
calc time increment size | |
virtual double | loadPredictor () |
Load predictor strategies. | |
virtual double | loadCorrector () |
Load corrector strategy. | |
virtual bool | staticLinearSolver (double step, double *newt, bool *conv, double *err, int *niter, bool eval) |
virtual bool | staticSolver (double step, double *newt, bool *conv, double *err, int *niter, bool eval) |
The NLSolver.stepStaticNL() static solver implementation This function expects to receive as parameters: 1) The current step 2) A pointer for the new load factor -> returns new load factor. | |
virtual bool | transientSolver (double dt, double *newt, bool *conv, double *err, int *niter, bool eval) |
The transient nonlinear implementation (femNLSolver.transient) This function manages the transient step. | |
virtual bool | transientLinearSolver (double) |
virtual bool | iterativeSolver (double dt, double *newt, bool *conv, double *err, int *niter, bool eval) |
The transient nonlinear implementation (femNLSolver.transient) This function manages the transient step. More... | |
virtual bool | stepStandardDc (double dt, double *newt, bool *conv, double *err, int *niter, bool eval) |
The transient nonlinear implementation (femNLSolver.transient) This function manages the transient step. More... | |
virtual bool | stepTransientNonlinearAut (double dt, double *newt) |
virtual bool | stepNonlinearAut (double dt, double *newt) |
![]() | |
bool | initSolver (GmpFemAssembler::FixedDofMode assemblerMode, bool assemblerReverseMapping, bool enableFastUpdate) |
Basic implementation of the init function receiving as parameters the configuration options that are sent to the assembler. More... | |
virtual bool | initResultAttributes (QString prefix) |
Helpper function used to register the set of result attributes managed by the fem solver. To avoid name clashes, all of them should prepend the given prefix (usually "fem") to the attribute name. More... | |
virtual bool | cleanup () |
Dealocates memory and sets allocated resources to NULL. As a convenience, returns false. | |
GmpFemPhysics::FemResultType | prepareMatrices (bool skipFixedBcs=false) |
Auxilliary function used to fill the K, and f matrices / vectors along with any other matrices/vectors stored in _matSet and _vecSet. More... | |
GmpFemPhysics::FemResultType | saveElementData (int iter) |
Auxilliary function used to save to the configured file data from the set of element matrices and vectors. In this process, physics are called again to fill element matrices and vectors (according to the active matrices and vectors in _vecSet and _matSet) but the results are NOT added to the global matrices. | |
bool | solveLinearSystem (bool xFilled) |
Solves the linear system K.x = f taking K and f from the single equivalent matrix/vector or from matrix K / vector f, depending on the assembler configuration. Results are saved in _x. Prints the used matrices / result as configured in _printOptions. More... | |
double | timeConvert (double val) |
Auxiliary function created for being used by derived classes that converts a time value given in the current simulation time unit and returns that value converted to the physics expected unit. | |
GmpFemPhysics::FemResultType | traverseElements () |
Fills the stiffness matrix _K and the force vector _f by traversing elements asking physics for the local element matrix / force vector. More... | |
GmpFemPhysics::FemResultType | traverseBoundaryElements () |
Simillar to traverseElements(), this function steps through each element that belongs to an edge or face boundary condition, asking physics for additional contributions to the stiffness matrix _K and the force vector _f. | |
GmpFemPhysics::FemResultType | traverseExternalLoads () |
Simillar to traverseElements(), this function steps through each element associated to an external load force, asking physics for additional contributions to the force vector _f. | |
GmpFemPhysics::FemResultType | traverseContactBoundaries () |
Simillar to traverseElements(), this function steps through each possible contact from the given contact boundary conditions and asks the physics for additional contributions to the solver matrices and vectors. This is different from the other element based loops since the returned matrices are not restricted to a single element, but can combine nodes from several ones. | |
GmpFemPhysics::FemResultType | traverseElementsForSaving (FILE *f, int iter) |
Similar to traverseElements() but instead of adding the elements to the assembler, saves the element data to the given file. More... | |
GmpFemPhysics::FemResultType | fillElementData (const GmElement *e, int physIndex) |
Helper function used by traverseElements() to get data for a single element calling phys->fillElementData() | |
GmpFemPhysics::FemResultType | fillElementBoundaryData (const GmElement *e, int physIndex, const GmBoundaryCondition *bc, int bcIndex, int bcListIndex, int border, const GmCellBoundary *b) |
Helper function used by traverseBoundaryElements() to get data for a single element calling phys->fillElementDataForBc() | |
GmpFemPhysics::FemResultType | fillContactData (GmMatrixDof &localDofMap, const GmContactBoundaryCondition *cbc, int physIndex, int index1, int index2) |
Helper function used by traverseContactBoundaries() to get data for a single contact pair calling phys->fillContactData() | |
bool | addFixedForces () |
Adds fixed nodal forces, as seen by each physics, to the global vector _f. | |
bool | collectFixedBcs (QVector< bool > &fixedRows, QVector< double > &fixedValues, QList< int > &dofIndex) |
Fills the set of vectors received as parameters with the complete set of fixed boundary conditions, as seen by the full set of physics in use, checking for possible conflicts in conditions. More... | |
virtual void | collectGlobalContactPairs (const GmContactBoundaryCondition *cbc, QList< QPair< int, int > > &contactPairs) |
Global search rule used to define contact surfaces. Should be implemented to allow global contact conditions, filling contactPairs with the indices in cbc of the in-contact surfaces. | |
void | printElementData (const GmElement *e, int ndof, const int *dofMapping, const GmpFemPhysics *p, const GmBoundaryCondition *bc=NULL, const GmContactBoundaryCondition *cbc=NULL, int loadId1=-1, int id2=-1) |
Print element dof mapping, local matrices and vectors depending on the current solver print options. More... | |
void | printGlobalData (bool linear) |
Print global matrices and vectors, either after assembling or before solving the linear system. | |
![]() | |
QObject * | sender () const const |
int | senderSignalIndex () const const |
int | receivers (const char *signal) const const |
bool | isSignalConnected (const QMetaMethod &signal) const const |
virtual void | timerEvent (QTimerEvent *event) |
virtual void | childEvent (QChildEvent *event) |
virtual void | customEvent (QEvent *event) |
virtual void | connectNotify (const QMetaMethod &signal) |
virtual void | disconnectNotify (const QMetaMethod &signal) |
Protected Attributes | |
double | _dl |
double | _dl0 |
double | _ddl |
double | _Lold |
double | _tin |
double | _attempt = 0.0 |
double | _dt |
double | _oldIter |
double | _hlast |
int | _ninc = 0 |
counter of consecutive converged increments | |
int | _incT = 0 |
counter of consecutive increments, i, and if no cut-back within those increments | |
bool | _iter4old = false |
corresponding to True If no more than IG iterations are required in two consecutive increments | |
GmVector | _v |
GmVector | _a |
GmVector | _dx |
GmVector | _dv |
GmVector | _xold |
GmVector | _vold |
GmVector | _dxt |
GmVector | _ddx |
GmVector | _ddv |
GmVector | _ddxt |
GmVector | _ddxr |
GmVector | _FeGs |
GmVector | _erL |
GmVector | _erD |
GmVector | _r0 |
GmVector | _qta |
An overall time-averaged value of the typical flux for field α so far during this step including the current increment. | |
GmVector | _qtamax |
The time-averaged value of the largest flux corresponding to the field αduring this step, excluding the current increment. | |
QList< int > | _fixedDof |
bool | _fastUpdate |
![]() | |
int | _solverId |
A unique index for this solver object used by the result attribute auto prefix feature. | |
GmElementMesh * | _mesh |
Mesh with the nodes and cells data. | |
GmSimulationData * | _simulation |
Simulation Data. | |
QList< GmpFemPhysics * > | _physicsList |
List of physics that will cooperate to generate the system response. | |
GmNumSolver * | _solver |
Numeric solver that will be used to solve the set of linear equations. | |
const GmLogCategory & | _logger |
Logger for eventual messages. | |
QMutex | _printMutex |
Mutex for serializing calls to printElementData when traversing elements in parallel. | |
bool | _enabledWarn [NumDisabledWarnings] |
Vector storing which warnings are enabled. | |
GmpFemSolverOptions | _solverOptions |
The set of solver options. | |
UnitConverter * | _timeConv |
Unit converter from the currentTimeUnit to the physics unit. Used only by derived solvers. | |
GmpFemAssembler * | _assembler |
Assembler used for adding element contributions to vectors and matrices. | |
const GmpFemAssemblerMatrixCombiner * | _matCombiner |
The combiner object used with the assembler to merge data from multiple elemental matrices. | |
const GmpFemAssemblerVectorCombiner * | _vecCombiner |
The combiner object used with the assembler to merge data from multiple elemental vectors. | |
GmpFemVectorSet | _vecSet |
Set of global/element vectors filled by physics. | |
GmpFemMatrixSet | _matSet |
Set of global/element matrices filled by physics. | |
GmVector | _x |
State vector. | |
GmVector | _r |
Residual vector. | |
unsigned | _preCount |
Number of times that prepareMatrices() was called. | |
unsigned | _runCount |
Number of times that run() was called. | |
QVariantMap | _controlMap |
Variant map used for storing dump control data for the solver. | |
int | _iterResAttr |
The registered id for the "fem.iter" result attribute. | |
int | _iterResErrAttr |
The registered id for the "fem.iterErr" result attribute. | |
Additional Inherited Members | |
![]() | |
enum | PrintOptions { PRINT_ELEMENT_MATRICES = 0x001, PRINT_ELEMENT_VECTORS = 0x002, PRINT_ELEMENT_DOF_MAPPING = 0x004, PRINT_EQ_MATRIX = 0x008, PRINT_EQ_VECTOR = 0x010, PRINT_GLOBAL_MATRICES = 0x020, PRINT_GLOBAL_VECTORS = 0x040, PRINT_LINEAR_MATRIX = 0x080, PRINT_LINEAR_VECTOR = 0x100, PRINT_LINEAR_RESULT = 0x200 } |
Flag options used to compose the printOptions constructor parameter. More... | |
![]() | |
QString | tr (const char *sourceText, const char *disambiguation, int n) |
QString | trUtf8 (const char *sourceText, const char *disambiguation, int n) |
QMetaObject::Connection | connect (const QObject *sender, const char *signal, const QObject *receiver, const char *method, Qt::ConnectionType type) |
QMetaObject::Connection | connect (const QObject *sender, const QMetaMethod &signal, const QObject *receiver, const QMetaMethod &method, Qt::ConnectionType type) |
QMetaObject::Connection | connect (const QObject *sender, PointerToMemberFunction signal, const QObject *receiver, PointerToMemberFunction method, Qt::ConnectionType type) |
QMetaObject::Connection | connect (const QObject *sender, PointerToMemberFunction signal, Functor functor) |
QMetaObject::Connection | connect (const QObject *sender, PointerToMemberFunction signal, const QObject *context, Functor functor, Qt::ConnectionType type) |
bool | disconnect (const QObject *sender, const char *signal, const QObject *receiver, const char *method) |
bool | disconnect (const QObject *sender, const QMetaMethod &signal, const QObject *receiver, const QMetaMethod &method) |
bool | disconnect (const QMetaObject::Connection &connection) |
bool | disconnect (const QObject *sender, PointerToMemberFunction signal, const QObject *receiver, PointerToMemberFunction method) |
![]() | |
enum | DisabledWarnings { PrescribedForceForInvalidDof, PrescribedForceForFixedDof, FixedBcForInvalidDof, ConflictingFixedBcValue, NumDisabledWarnings } |
Enums describing the set of warnings that can be disabled by simulation options. More... | |
![]() | |
virtual void | meshChanged () |
Slot called when the underlying mesh has been changed. More... | |
![]() | |
objectName | |
Basic class for solving Non linear FEM problems.
|
protectedvirtual |
Adds to 'state' the data items that should be saved for this FEM process. Should probably be overriden by derived classes to add their own requirements.
Auxiliar variant map dump item that notifies us for loading and saving data to the map
< Our "father" solver object
Reimplemented from GmpFemSolver.
|
protectedvirtual |
Apply Dirichlet (fixed) boundary conditions to the global equation system.
Apply boundary conditions, clearing rows and columns in the stiffness matrix and adjusting the global force vector as needed.
Follows the procedures described by Fellipa in sections 3.6.1 and 3.6.2 (pages 3-12 and 3-13). Our algorithm does the same thing as the one presented in figure 21.5 (page 21-7), but in a different way (hopefully more clear).
Returns in numFixed the number of degrees of freedom "removed" from the system.
Reimplemented from GmpFemSolver.
|
protectedvirtual |
Evaluates the stiffness matrix using BFGS method.
Evaluates stiffness matrix and internal forces using BFGS method.
|
protectedvirtual |
Evaluates the stiffness matrix using Broyden's method.
Evaluates stiffness matrix and internal forces using Broyden's method.
|
virtual |
Prepares the solver for assembling matrices by creating the assembler object and allocating the needed memory. This operations only need to be done once if the solver is used in a loop.
When overriding this function, if you need to pass to the assmebler different configurations from the default (REMOVE_DOF & false), consider calling initSolver() in your reimplementation instead of init().
Reimplemented from GmpFemSolver.
|
protectedvirtual |
Helpper function used to initialize the matrix / vector sets with the definition of the needed types and also the relationships of those sets with the assembler.
In a nutshell, this function is responsible for callint the initTypes() methods from matrix and vector sets. It is also its responsibility to call other configuration functions such as GmpFemMatrixSet::setTransposedDofVector(), GmpFemVectorSet::addTransposedDofVector() and GmpFemVectorSet::addSavedDofVector() as needed to drive the method that will be used by the assembler.
This is the function to be reimplemented when a different solver driver needs to ask elements for additional matrices / vectors.
Returns false on errors.
Reimplemented from GmpFemSolver.
|
protectedvirtual |
The transient nonlinear implementation (femNLSolver.transient) This function manages the transient step.
Maximun value of the array.
Maximun value.
|
protectedvirtual |
Physical group.
Maximun value.
|
protectedvirtual |
The transient nonlinear implementation (femNLSolver.transient) This function manages the transient step.
|
protectedvirtual |
Updates the stiffness matrix.
Updates stiffness matrix for static nonlinear solver.