Martin Diehl
4fb1cb8f87
merged calcmodes, i.e. equilibrium and fulfilling of stress BC is done in one step,
...
made convergence independent of size and resolution,
polishing output in DAMASK_spectral.f90
added function to compute eigenvalues without eigenvectors and function to convert a 3x3 logical to a 9 vector in math.f90
removed obsolete variable in numerics.f90
2011-08-26 14:06:37 +00:00
Martin Diehl
f99bf63397
removed functions added for debugging of divergence calculation to math.f90
...
corrected calculation of stress BC condition. Depending on given BC, the stiffness matrix is reduced and than inversed. Then it is filled with zeros and used for the calculation of the correct change of deformation gradient. All calculation is done using dP/dF
2011-08-10 17:45:37 +00:00
Martin Diehl
34de2e301b
DEBUGGING VERSION for calculation of Divergence (various methods implemented). DO NOT USE. Wait for next update coming in a few minutes
2011-08-10 16:02:13 +00:00
Onur Guevenc
907965272b
added option to build DAMASK_spectral without openMP (make OPENMP=OFF)
...
fixed bug in bc_temperature assignment that was hitting memory.
Temperature is taken from the first loadcase and evolves from there in an adiabatic fashion for the moment. I.e. T-specifications from later loadcases are ignored...
2011-08-03 17:57:28 +00:00
Martin Diehl
b5abd6f203
adding new option OPTIMIZATION with values OFF, DEFENSIVE and AGGRESSIVE to makefile. Polishing of makefile and DAMASK_spectral.f90
2011-08-03 07:54:32 +00:00
Onur Guevenc
64435b8a97
loadcase takes 'temperature' (or 'temp') as input. Previously it was hard-wired.
2011-08-02 13:58:28 +00:00
Martin Diehl
3adb7ab382
corrected makefile, now working again without giving standard values explicitly.
...
did some small modification in order to make it possible to compile with gfortran. Changed NaN=0.0/0.0 to bitwise representation (3 different ways)
2011-08-01 10:11:32 +00:00
Christoph Kords
fb121b1435
* boundary condition masking changed
...
* damper initialized with one
* inversion of Mandelized stiffness tensor does not work, have to use plain tensor
* new functions in math that allow for conversion between Mandel and Plain tensors
2011-07-29 15:57:39 +00:00
Martin Diehl
72d20875de
added some switches and variables to the makefile to make it more flexible
...
DAMASK_spectral.f90 is a "debug version" with a number of different criteria to determine divergence. will be removed later on.
2011-07-25 16:30:21 +00:00
Martin Diehl
09ba92c26e
added and rearranged error messages, polished output and simplified calculation of f depending on fdot or L
...
guessing along former trajectory is now on per default, 'guessreset' and 'dropguessing' switch it off.
2011-07-14 09:37:31 +00:00
Martin Diehl
56340fd487
changed input of loadcase. Now fdot (time derivative) can be used instead of velocity gradient. Velocity gradient needs to have each line fully or not at all defined, as for other loadcases the stress BCs are not known in advance. Also added the possibility to keep guessing along trajectory of former loadcase.
...
changed back to use the compliance of initial linear material behavior.
added counter of non-converged steps
renamed compiler flags in makefile
2011-07-13 16:33:12 +00:00
Martin Diehl
a08eaecf4e
polishing (names of variables, spaces, comments etc.)
...
renamed mpie to DAMASK in IO.f90 error message
2011-07-07 15:27:35 +00:00
Martin Diehl
a561ef1cf5
corrected calculation of divergence in Fourier space, removed normalization of normdyad (was useless), now using correct compliance for calculation of stress BC.
2011-07-07 10:03:55 +00:00
Franz Roters
18a5841bc5
keyword freq for frequency of output writing was not evaluated when reading loadcase file
2011-07-07 09:18:05 +00:00
Onur Guevenc
e889b69bec
Ability to define loglinear time increments by using 'logincs' or 'logsteps' instead of 'incs' in loadcase file is added.
2011-07-06 13:10:18 +00:00
Philip Eisenlohr
314adae1bb
spectral solver can now write result at reduced output frequencies: use 'f', 'freq', or 'frequency' keyword in loadcase definition.
...
example: ... time 1 incs 100 freq 10 will write output at times 0.0, 0.1, 0.2, ... , i.e., every ten increments only.
2011-06-15 17:48:14 +00:00
Philip Eisenlohr
4a694fa7fd
changed wildcard letter from '#' to '*' now consistent with new IO comment parsing.
...
fixed memory bug with bc_maskvector.
some brushing up here and there...
2011-06-06 15:20:28 +00:00
Christoph Kords
ae4c8fa2d8
renaming of "mpieNumThreadsInt" to "DAMAKS_NumThreadsInt" was forgotten here
2011-06-06 08:51:07 +00:00
Christoph Kords
78f74e8991
forgot to use "mesh_ipCenterOfGravity"
2011-05-26 09:23:13 +00:00
Christoph Kords
cac45cff96
ip coordinates are now updated every cycle; this is needed for the nonlocal internal stress fields
...
* Marc: node displacements are added to initial node coordinates (mesh_node0) to get current node positions (mesh_node), then ip coordinates are deduced
* Abaqus: ip coordinates are directly updated, no update of node coordinates!
* Spectral: for the moment no update of either ip or node coordinates! passing only dummy values with initial ip coordinates
2011-05-24 15:57:59 +00:00
Philip Eisenlohr
fa98133f8e
renamed interface from "mpie_interface" to "DAMASK_interface"
2011-05-11 17:01:03 +00:00
Philip Eisenlohr
e5a2d829b0
new substructure with 'include' and 'config' directories
...
renaming of mpie_... to DAMASK for main usersubroutines
extension of element outputs from 5 to 8 digits (FFT!!!)
2011-05-11 16:38:45 +00:00