5f0440fde3
crystallite_integrateStress now raises an error if Fp becomes NaN
2012-10-17 16:19:42 +00:00
31f20e51da
just to be able to tell about the bug I found in the hex slip systems---which I actually forgot to mention during the last commit---I slipped in a syntax error that now is fixed, too...
2012-10-12 20:31:50 +00:00
05c5ea42d7
prepared loop distributions to hold "over shooters"
2012-10-12 18:00:55 +00:00
82a13af474
tidied up the Schmid matrix calculations.
...
added internal check against dilatation in Schmid matrices.
2012-10-12 17:59:50 +00:00
b20c612e25
added error "0" for internal (run time) checks failing
2012-10-12 17:55:23 +00:00
1b2edd7e7d
fixed error in _identity (wrong delta-function)
...
added _exp33 for matrix exponential
added trace33 for matrix trace
2012-10-12 17:54:20 +00:00
324dfda5a2
added comments (doxygen conform) to phenopowerlaw, added warning on specifying h0_sliptwin as it has no effect
2012-10-11 14:49:12 +00:00
e4054c116d
fixed problem with gnu compiler: allocation of variables in constitutive_nonlocal_stateInit was not safe when not using any nonlocal constitution
...
put header of debugging output in constitutive_init at the start of the function in order to make debugging easier
2012-10-09 12:34:57 +00:00
a84e0529c3
make correction of hardening coefficients thread safe
2012-10-05 16:05:51 +00:00
40a2e652e7
introduced material parameter "linetension", which indicates how strong the effect of the line tension on the hardening coefficients is
2012-10-05 14:42:41 +00:00
7cc6108a04
hardening coefficients for fcc are density dependent (line tension effect, Devincre, Kubin, Hoc; 2006; Physical analyses of crystal plasticity by DD simulations)
...
fixed bug with thermal annihilation producing negative density in case of very small dipole density and high climb velocity
2012-10-04 18:08:40 +00:00
80465c0b49
moved option for UTF-8 encoding to spectral interface as this is now called before prec init.
...
if UTF-8 encoding is needed for FEM solver (old compilers might not support it), add the same statement to the corresponding interface routine
2012-10-04 14:22:39 +00:00
e2ea3cdff6
again some changes regarding creation of NaN, ieee-conforming way is not in the complilers yet.
2012-10-02 16:53:03 +00:00
bce8012c89
random distribution of initial dislocation density now really independent of meshsize
2012-10-02 15:26:58 +00:00
d5ce49c471
numerics got some protected statements and is now reading in the new keywords for the solver selection in small letters.
...
spectral solver got cut back facilities + improved output to screen
2012-10-02 15:26:56 +00:00
91b7883c2a
random distribution of initial dislocation density now independent of mesh size
2012-10-02 13:35:34 +00:00
1502a71f0c
added possibility to start with an initial random distribution of dislocation segments with specific overall density
2012-10-02 12:57:24 +00:00
0bcb8f59db
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 12:53:25 +00:00
6301787df4
moved the the variable cutBack from FEsolving.f90 (where it was not used at all) to the interface routines to make responsibility of cut back control clear
...
PLEASE RUN setup_code.py IF YOU DO MARC SIMULATIONS, THE INTERFACE HAS CHANGED
2012-10-02 12:46:58 +00:00
c8b7da432c
set debug_level to be protected as it is read-only for other modules
2012-10-02 09:35:39 +00:00
4c7bcb8df9
modified formated integer output and added error number in case maximum number of cutbacks is reached
2012-10-02 09:30:13 +00:00
9fc46b4950
added /echo/ statements
2012-10-02 08:12:53 +00:00
b99a846503
new scaling feature in numerics (unitlength) did not scale the reference coordinates in the first place; now this is taken care of in mesh.f90, so scaling should work properly
2012-09-25 16:21:58 +00:00
d0fd3c03c5
forgot to check in; belongs to rev 1762
2012-09-24 16:22:25 +00:00
912277d1ed
substituted the call to the flush subroutine with the intrisic flush(6) function
2012-09-24 12:10:28 +00:00
8336d0705c
new numerics parameter "unitlength" that scales the model size in DAMASK; this enables to use very small models in marc without getting any numerical problems for example with contact which does not seem to work when using very small numbers
2012-09-24 06:13:26 +00:00
b587bffdd5
added some 'protected' statements to variables that shouldn't be changed by other modules
2012-09-22 16:49:41 +00:00
8ed69170c5
corrected wrong warning complaining about FFTW string when DAMASK_NUM_THREADS undefined
2012-09-14 16:39:42 +00:00
1cfc467705
added new cut back parameter for spectral solver and instrumented variables for use with doxygen
2012-09-13 09:48:38 +00:00
b9a2890665
fixed bug in loadcase rotation: when not specifying "deg" the rotation was reset to 0,0,0..!
2012-09-06 14:05:28 +00:00
29023cea3c
minor correction of debugging output
2012-09-05 11:19:46 +00:00
330d71864e
fixed bug with stress BC handling
2012-09-05 08:30:15 +00:00
f153866030
introduced factor to control edge contribution to multiplication;
...
dislocation density below a single dislocation per IP considered not significant
2012-09-04 16:56:37 +00:00
22812c9a91
some minor improvement on precision detection: checking only once (in prec and no longer in math and crystallite), added one more 4/8 switch for LAPACK, as there is no single precision FFTW, stopping compilation for spectral method if FLOAT=4
...
new function in IO to print integers without leading zeros, implemented it at some places in the new spectral solver (reporting still needs some serious polishing)
updated preprocessing for documentation to handle precision correctly
2012-08-30 20:26:28 +00:00
299d9e6bc4
fixed bug in seed initialization that didn't allow to redo calculations that started with a random seed. Now the whole random seed array is set to one value, previously only the first one was reported and the remaining ones were random
...
removed old random function used by halton
2012-08-30 14:58:30 +00:00
d80b45c1d9
added relativistic correction of velocity (maximum velocity given by parameter in material.config); probability of thermal activation in velocity law not allowed to be greater than one.
...
output of creation rate for edge jogs now for own system, not collinear.
check for negative densities in dotState now uses absoluteToleranceRho instead of significantRho.
2012-08-30 07:33:13 +00:00
8e3b547cbd
changed numerics.f90 back to the old openmp statement as this doesn't seem to cause the problem and is standard conforming
...
corrected typo in constitutive_j2.f90 that might cause abaqus to crash
now running 20 tests of abaqus in order to have a decent statistic about the crash behavior
improved abaqus_v6.env
2012-08-29 16:16:10 +00:00
84e3f94424
added flush for easier debugging
2012-08-29 15:23:43 +00:00
d83150be1b
now capturing "terminallyIll" for basic solution, should end in a automatic cutback of time steps for the spectral solver
2012-08-29 05:50:42 +00:00
3106b8ce0c
corrected definition of integer type for gfortran
2012-08-29 05:19:35 +00:00
d089c12947
abaqus seems to have similar problems with the openmp lib as marc (use module vs. include file), i hope this fixes the crashes
2012-08-28 19:48:38 +00:00
5e9e8497e6
new solver is now compiling without a PETSc installation, however only the plain basic solver is available then.
...
also removed reporting of PETSc related variables in the case it is not installed/found
2012-08-28 19:19:47 +00:00
0f64289d75
fixed bug causing trouble with ifort 11 (abaqus)
2012-08-28 19:10:54 +00:00
0959ff3299
substituted hand written matrix inversion by LAPACK version with precision selection.
...
also introduced check for inversion into DAMASK_spectral_Utilities.f90 for the stress BC calculation. This part is further improved by using 5% of the reference stiffness to avoid trouble in the fully plastic regime (where the stiffness is underestimated)
Test for Marc 2010 is updated because the new inversion give slightly different results near 0 (order of e-13)
2012-08-28 16:59:45 +00:00
73349d02f5
merged precision info into the file prec.f90 and removed prec_single.f90
...
precision information is now parsed as a makro, defined by the compiler, eg. -DFLOAT=8 for 8 byte long floats or -- in case no values is given -- for the FEM wrapper by the interface files.
Please don't forget to run setup_code.py !!!!!!!!!!!
abaqus_v6.env is adjusted to specify the standard precision (float=8, int=4) and got some extra compilation options, some deprecated warning supressors are removed
Optional Marc modification will follow that allow automatic selection of suitable precision
2012-08-28 16:08:17 +00:00
f3e7b3cec8
added separate output of dipole formation rate for edge and screw
2012-08-27 15:57:31 +00:00
aa380dcb08
added separate output of dipole formation rate for edge and screw
2012-08-27 15:57:22 +00:00
efd1b513cf
fixed bug in divergenceFDM
2012-08-27 15:43:17 +00:00
9214da4465
added missing labels to list of available outputs
2012-08-27 12:41:57 +00:00
96ba5ecae4
moved some more 'mesh related' functions for post processing from math.f90 to mesh.f90
...
f2py functions remaining in math.f90 now uses assumed size arrays in order to have simpler interfaces. This is only working with python 2.7!
changed python pre- and postprocessing scripts.
If you encounter any problems whith core modules, try to remove the old core.so in the lib/damask
2012-08-27 08:04:47 +00:00
Christoph Kords