Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/tdegeus/gmatelastoplasticfinitestrainsimo
Simo elasto-plastic model
https://github.com/tdegeus/gmatelastoplasticfinitestrainsimo
Last synced: 11 days ago
JSON representation
Simo elasto-plastic model
- Host: GitHub
- URL: https://github.com/tdegeus/gmatelastoplasticfinitestrainsimo
- Owner: tdegeus
- License: mit
- Created: 2018-11-11T14:34:38.000Z (almost 6 years ago)
- Default Branch: main
- Last Pushed: 2022-10-15T14:36:33.000Z (about 2 years ago)
- Last Synced: 2024-04-16T06:53:13.848Z (7 months ago)
- Language: C++
- Size: 929 KB
- Stars: 1
- Watchers: 3
- Forks: 1
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# GMatElastoPlasticFiniteStrainSimo
[](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo/actions)
[](https://tdegeus.github.io/GMatElastoPlasticFiniteStrainSimo)
[](https://anaconda.org/conda-forge/gmatelastoplasticfinitestrainsimo)
[](https://anaconda.org/conda-forge/python-gmatelastoplasticfinitestrainsimo)
Simo elasto-plastic model.
An overview of the theory can be found in `docs/theory/readme.tex`
conveniently compiled to this [PDF](docs/theory/readme.pdf).
# Disclaimer
This library is free to use under the
[MIT license](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo/blob/master/LICENSE).
Any additions are very much appreciated, in terms of suggested functionality, code,
documentation, testimonials, word-of-mouth advertisement, etc.
Bug reports or feature requests can be filed on
[GitHub](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo).
As always, the code comes with no guarantee.
None of the developers can be held responsible for possible mistakes.
Download:
[.zip file](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo/zipball/master) |
[.tar.gz file](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo/tarball/master).
(c - [MIT](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo/blob/master/LICENSE))
T.W.J. de Geus (Tom) | [email protected] | www.geus.me |
[github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo](https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo)
# Python implementation
## Partial example
```python
import GMatElastoPlasticFiniteStrainSimo.Cartesian3d as GMat
import GMatTensor.Cartesian3d as tensor
shape = [...]
K = np.empty(shape)
G = np.empty(shape)
...
GMat.ElasticXd model(K, G);
...
F = tensor.ArrayXd(shape).I2
...
model.F = F
print(model.Sig)
```
## Installation
### Using conda
```bash
conda install -c conda-forge python-gmatelastoplasticfinitestrainsimo
```
Note that *xsimd* and hardware optimisations are **not enabled**.
To enable them you have to compile on your system, as is discussed next.
### From source
> You need *xtensor*, *xtensor-python* and optionally *xsimd* as prerequisites.
> The easiest is to use *conda* to get the prerequisites:
>
> ```bash
> conda install -c conda-forge xtensor xsimd xtensor-python
> ```
>
> If you then compile and install with the same environment you should be good to go.
> Otherwise, a bit of manual labour might be needed to treat the dependencies.
```bash
git checkout https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo.git
cd GMatElastoPlasticFiniteStrainSimo
# Only if you want to use hardware optimisation:
export SKBUILD_CONFIGURE_OPTIONS="-DUSE_SIMD=1"
python -m pip install . -v
```
# C++ implementation
## Partial example
```cpp
#include
#include
namespace GMat = GMatElastoPlasticFiniteStrainSimo::Cartesian3d;
namespace TEnsor = GMatTensor::Cartesian3d;
int main()
{
static const size_t rank = ...;
xt::xtensor K = ...;
xt::xtensor G = ...;
GMat::ElasticXd model(K, G);
...
xt::xtensor F = Tensor::ArrayX2(K.shape()).I2();
...
// all necessary computation are done at this point
model.set_F(F);
...
// get reference to stress
auto Sig = model.Sig();
return 0;
}
```
## Installation
### Using conda
```bash
conda install -c conda-forge gmatelastoplasticfinitestrainsimo
```
### From source
```bash
git checkout https://github.com/tdegeus/GMatElastoPlasticFiniteStrainSimo.git
cd GMatElastoPlasticFiniteStrainSimo
cmake -Bbuild
cd build
cmake --install .
```
## Compiling
## Using CMake
### Example
Your `CMakeLists.txt` can be as follows
```cmake
cmake_minimum_required(VERSION 3.1)
project(example)
find_package(GMatElastoPlasticFiniteStrainSimo REQUIRED)
add_executable(example example.cpp)
target_link_libraries(example PRIVATE GMatElastoPlasticFiniteStrainSimo)
```
### Targets
The following targets are available:
* `GMatElastoPlasticFiniteStrainSimo`
Includes the library and its dependencies.
* `GMatElastoPlasticFiniteStrainSimo::assert`
Enables IO-assertions by defining `GMATELASTOPLASTICFINITESTRAINSIMO_ENABLE_ASSERT`.
* `GMatElastoPlasticFiniteStrainSimo::debug`
Enables assertions of all dependencies.
* `GMatElastoPlasticFiniteStrainSimo::compiler_warings`
Enables compiler warnings (generic).
### Optimisation
It is advised to think about compiler optimisation and enabling *xsimd*.
Using *CMake* this can be done using the `xtensor::optimize` and `xtensor::use_xsimd` targets.
The above example then becomes:
```cmake
cmake_minimum_required(VERSION 3.1)
project(example)
find_package(GMatElastoPlasticFiniteStrainSimo REQUIRED)
find_package(xtensor REQUIRED)
find_package(xsimd REQUIRED)
add_executable(example example.cpp)
target_link_libraries(example PRIVATE
GMatElastoPlasticFiniteStrainSimo
xtensor::optimize
xtensor::use_xsimd)
```
See the [documentation of xtensor](https://xtensor.readthedocs.io/en/latest/).
## By hand
Presuming that the compiler is `c++`, compile using:
```
c++ -I/path/to/GMatElastoPlasticFiniteStrainSimo/include ...
```
Note that you have to take care of the *xtensor* dependency, the C++ version, optimisation,
enabling *xsimd*, ...
## Using pkg-config
Presuming that the compiler is `c++`, compile using:
```
c++ `pkg-config --cflags GMatElastoPlasticFiniteStrainSimo` ...
```
Note that you have to take care of the *xtensor* dependency, the C++ version, optimization,
enabling *xsimd*, ...
# References / Credits
+ The model is described in
*M.G.D. Geers (2004).
Finite strain logarithmic hyperelasto-plasticity with softening:
a strongly non-local implicit gradient framework.
Computer Methods in Applied Mechanics and Engineering, 193(30–32), 3377–3401,
[doi: 10.1016/j.cma.2003.07.014](https://doi.org/10.1016/j.cma.2003.07.014)*.
+ [xtensor](https://github.com/QuantStack/xtensor) is used under the hood.
+ The eigenvalue and eigenvector decomposition is done by an algorithm due to
[Joachim Kopp](https://www.mpi-hd.mpg.de/personalhomes/globes/3x3/).
See also:
*J. Kopp (2008).
Efficient numerical diagonalization of hermitian 3x3 matrices.
Int. J. Mod. Phys. C. 19:523-548,
arXiv: [physics/0610206](https://arxiv.org/abs/physics/0610206)*.
# Upgrading instructions
## Upgrading to >v0.3.*
The individual material point and the array of material points was fully integrated.
In addition, the number of copies was reduced.
### C++
There is only a single class `Elastic`. It's functions where renamed:
* `.setDefGrad(...)` -> `.set_F(...)`
* `.Stress()` -> `.Sig()` (now returns a reference).
* `.stress(...)`: deprecated.
* `.Tangent()` -> `.C()` (now returns a reference).
* `.tangent(...)`: deprecated.
### Python
There is only a single class `Elastic`. It's functions are converted to properties:
* `.setDefGrad(...)` -> `.F = ...`
* `.Stress()` -> `.Sig` (now returns a reference).
* `.stress(...)`: deprecated.
* `.Tangent()` -> `.C` (now returns a reference).
* `.tangent(...)`: deprecated.
## Upgrading to >v0.2.*
`xtensor_fixed` was completely deprecated in v0.2.0, as were the type aliases
`Tensor2` and `Tensor4`.
Please update your code as follows:
* `Tensor2` -> `xt::xtensor`.
* `Tensor4` -> `xt::xtensor`.
**Tip:** Used `auto` as return type as much as possible.
This simplifies implementation, and renders is less subjective to library
return type changes.
Compared to v0.1.0, v0.2.0 has some generalisations and efficiency updates.
This requires the following changes:
* `Matrix` has been generalised to `Array`. Practically this requires changing:
- `Matrix` to `Array<2>` in C++.
- `Matrix` to `Array2d` in Python.
Note that `Array1d`, `Array3d`, are also available.
* `Array.check` ->
```cpp
if (xt::any(xt::equal(array.type(), Type::Unset))) {
throw std::runtime_error("Please set all points");
}
```
Note however that it is no longer required to set all points,
unset points are filled-up with zeros.
* Strain is now stored as a member.
Functions like `stress` now return the state based on the last specified deformation gradient,
specified using `setDefGrad(F)`. This leads to the following changes:
- `stress`: no argument.
- `tangent`: no argument, single return value (no longer returns stress).
# Change-log
## v0.3.0
Complete API overhaul.
## v0.2.1
* Using scikit-build, setuptools_scm, xtensor-python (#21)
* CMake clean-up (#21)
## v0.2.0
Compared to v0.1.0, v0.2.0 has some generalisations and efficiency updates.
This requires the following changes:
* `Matrix` has been generalised to `Array`. Practically this requires changing:
- `Matrix` to `Array<2>` in C++.
- `Matrix` to `Array2d` in Python.
Note that `Array1d`, `Array3d`, are also available.
* `Array` now sets zeros for all `Type::Unset` points.
The function `check` is deprecated accordingly.
* Strain is now stored as a member.
Functions like `stress` now return the state based on the last specified deformation gradient,
specified using `setDefGrad(F)`. This leads to the following changes:
- `stress`: no argument.
- `tangent`: no argument, single return value (no longer returns stress).
* Tensor operations are now provided centrally in the GMat eco-system,
by GMatTensor.