UpdateStress.cpp File Reference

#include <limits>
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/DEM/Bond/Beam.hpp"
#include "GranOO3/DEM/SupportShape.hpp"
#include "GranOO3/Core/Sensor.hpp"
#include "GranOO3/Core/Problem.hpp"
#include "GranOO3/Physic/Time.hpp"
#include "GranOO3/Math/UniRandom.hpp"
#include "GranOO3/PlugIn/UpdateStress.hpp"
#include "GranOO3/DEM/Bond/PlasticBeam.hpp"

Namespaces
	GranOO3
	GranOO3::PlugIn
	This namespace group the entire PlugIn library.

Functions
	GranOO3::PlugIn::REGISTER_GRANOO_PLUGIN_DESC (UpdateStress, desc)

Variables
static const char *	desc

Variable Documentation

desc

const char* desc

static

Initial value:

= R"V0G0N(
 
This PlugIn computes the stress for each specified `DEM::Element` with a 
`DEM::Stress` property. Note that if the `DEM::Stress` properties do not exist 
the plugin will automatically create these properties.
 
The Stress can be computed in two ways : 
- `VIRIAL` which uses the virial stress method [1,2]
- `LSQ` method based on least square method (experimental)
 
 
### Info ###
+ Process on         : `DEM::Element` set with `DEM::Stress` property
+ Multithreaded      : No
+ Recommended section: Processing 
 
 
### Required param ###
Method (string, must be "VIRIAL" or "LSQ")
:    Choose here the method that will be used to compute the stress
 
### Optional param ###
NeighbourLevel(int, default=1)
:    Choose the level number of neighbour to compute the virial stress. 
     The default value (one) is a good compromise between accuracy 
     and speed of computation. You can choose the zero value to speed up 
     the comutation. Therefore results are less accurate.
 
Periodic(bool, default=False)
:    Manage periodic boundaries or not. If this the case, the stress is computed 
     over the periodic boundaries.
 
UpdatePrincipalStress (bool, default=True)
:    Tell here if you want to update the principal virial stress. It 
     requires to compute the eigen values of the tensile stress tensor matrix. 
     This step is very time consuming but it is necessary for failure computation.
 
 
IncludeVelocity (bool, default=True)
:    Include the velocity in the computatio of the equivalent virial stress tensor 
     This option is only available for the "VIRIAL" method
 
VolumeFraction (float, default=auto)
:    By default, the volumic fraction is automatically computed. You can force it 
     by setting a value here. Note that the volumic fraction is very important for 
     the accuracy of the virial stress computation. This option is only available for 
     the "VIRIAL" method
 
Periodic (bool, default=false)
:    manage periodic boundaries or not... default is false, it means that you don't need 
     to manage periodic boundaries (free boundaries) 
 
### Example ###
```
<UPDATE-STRESS Method="Virial" VolumeFraction"0.68" NeighbourLevel="1" />
```
 
### Simulation example ###
`00011_{lattice}{dem}{mechanical}{brittle}{beam}_TENSILE-TEST/` and other brittle examples
 
 
> [1]: https://en.wikipedia.org/wiki/Virial_stress
> [2]: Andre, D., Jebahi, M., Iordanoff, I., Charles, J. L., & Neauport, J. (2013). 
> Using the discrete element method to simulate brittle fracture in the indentation of a silica glass 
> with a blunt indenter. 
> Computer Methods in Applied Mechanics and Engineering, 265, 136-147.
)V0G0N"