28259b2c46
switched dNeq and dEq to relative tolerance, removed single precision (makes things complicated
...
and was never used anyway)
2016-03-20 23:20:58 +01:00
97b52f60e7
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-03-20 22:28:49 +01:00
9704eb6926
small correction to variable declaration layout
2016-03-15 17:23:01 -04:00
654e58faf2
integer overflow due do missing brackets, wrong forwarding for large files
2016-03-14 22:30:55 +01:00
cf21d5ad49
did not compile without these statements
2016-03-13 19:01:01 +01:00
52ba6e19a0
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-03-12 17:09:18 +01:00
e84e90448e
Merge branch 'development' into CCodeUse
2016-03-11 21:38:39 +01:00
b9b490d02d
almost working
2016-03-11 20:59:14 +01:00
63fe0c6b91
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-03-09 17:55:54 +01:00
ed1948689f
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-03-09 11:40:31 -05:00
20fe46a030
corrected logical mask for reporting stress in load cases
2016-03-09 11:40:12 -05:00
849d487ab6
autogenerated file not needed in repository
2016-03-09 16:08:42 +01:00
976803ba0b
corrected error discovered by Haiming Zhang
2016-03-09 15:36:11 +01:00
564ff68fbb
first ideas
2016-03-08 23:32:18 +01:00
0c1fe236b4
fixed issue with multiple allocation of param%outputID (was done every line)
2016-03-07 17:13:36 -05:00
83ccd7c982
introduced (in)equality comparison for double precision
2016-03-05 22:25:13 +01:00
54d0a0a704
some important lines got lost during last commit
2016-03-01 10:58:52 +01:00
a7665bdab9
removed empty line remaining from old ID string
2016-02-29 14:26:06 +01:00
3f93eb9878
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-02-26 22:11:42 +01:00
c7ab5a9396
simplified and got MPI Heidelberg solution for eigenvalues/vectors back
2016-02-26 21:52:54 +01:00
ca3e1f0da0
further simpilifcations
2016-02-26 20:06:36 +01:00
17e75a1e0b
some more simplifications
2016-02-26 17:51:34 +01:00
76b67e88eb
re-indroduced special spectral decomposition for 33 tensors
2016-02-26 16:35:55 +01:00
5d0900ee2e
plasticity test (phenoplus) working again with changed polar decomposition
2016-02-26 15:36:24 +01:00
8419e293d0
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-02-22 15:47:14 +01:00
c87625b30c
wrong state handling
2016-02-22 15:45:18 +01:00
2bfc070797
clarified function description
2016-02-08 17:33:17 +01:00
fc1a21eb78
ignore marc aliases
2016-02-04 19:55:44 +01:00
6bd17b3e4e
integrated new version handling
2016-02-03 09:52:11 +01:00
ebe7707ac7
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
...
Conflicts:
code/math.f90
2016-02-03 09:15:15 +01:00
34b21cb278
LAPACK version as backup when analytic eigenvalues fail
2016-02-03 09:10:37 +01:00
08150564cf
added info like in spectral solver
2016-02-02 23:34:42 +01:00
a56f720e36
LAPACK version as backup when analytic eigenvalues fail
2016-02-02 18:59:04 +01:00
dc1e8f9def
clearifying names
2016-02-02 13:23:45 +01:00
15e1c7edd8
invariants for symmetric matrix
2016-02-02 09:44:51 +01:00
b28e70e36a
old version removed in rev 6896521bf2, seems to be more stable
2016-02-02 09:09:07 +01:00
bf04ee60f0
polishing
2016-02-02 08:45:47 +01:00
eaf9b41b7a
rewrote to follow formulas in paper, still having problems
2016-02-02 08:44:57 +01:00
2a86eef778
Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development
2016-02-01 21:42:50 +01:00
c7f1b677cc
needed cleaning and copyright
2016-02-01 21:37:27 +01:00
f13bdddc6c
wrong state alias access
2016-02-01 11:41:39 +01:00
1f5c54e5c1
analytic expression for eigenvalues of 3x3 matrix
2016-02-01 07:48:42 +01:00
5dc0f28179
analytical variant not stable enough (spectral decomposition of 0-Matrix)
2016-01-31 17:49:56 +01:00
9bc05acba2
re-introduced special eigenvalues routine for 3x3 matrices
2016-01-31 12:25:26 +01:00
66731575af
more prominent place for config files
2016-01-28 23:19:26 +01:00
e298d2ebaa
new version handling was broken
2016-01-28 08:56:12 +00:00
f2c7eea161
git compatible way of version handling for Fortran and Python
2016-01-27 16:48:27 +00:00
6e4239bb7c
removed Id
2016-01-27 14:06:21 +00:00
cdf09d6822
Bug fix for v4780 so it works for non transformation case
2016-01-26 13:20:45 +00:00
8fbbfaa106
Rotate elastic constants (fcc to hex)
2016-01-25 14:08:38 +00:00
288ab15ef9
more precise file names
2016-01-24 21:05:36 +00:00
a727f2b043
output related changes in last commit do not work since constitutive.f90 needs size and name in array notation
2016-01-22 08:13:05 +00:00
45c7bfa93d
introduced pointer aliases for state, abstol and dotate and type structure for input parameters
2016-01-22 01:08:36 +00:00
17b27271a5
init was missing (only important for restart)
2016-01-21 09:06:25 +00:00
eda03ed3e0
need to write out initial results for new CPFEM2
2016-01-20 16:19:05 +00:00
143266dbd4
missing variables from FESolving
2016-01-18 16:52:18 +00:00
072fa58027
fixed missing variable name change
2016-01-18 16:13:27 +00:00
a46fcc53b7
corrected a typo.
2016-01-18 09:07:48 +00:00
7ebd5e29e3
line too long
2016-01-17 19:42:24 +00:00
12653bf1a4
more rename g -> c/ipc, wrong ipc iteration for crystallite_integrateStateRK4 and single run fixed
2016-01-17 17:56:24 +00:00
b10774aae9
unused variable
2016-01-17 15:06:35 +00:00
cfd6579024
split CPFEM module into one for spectral solver (also suitable for PETSc FEM) and one for commercial FEM codes as all the ping pong handling, cut back detection etc. is not needed for spectral.
...
Needs certainly a better name, or might be superfluous at all if crystallite, constitutive, and homogenization take care of their data for restart write/read and forwarding
2016-01-17 15:03:54 +00:00
0799570a03
renamed g -> c/ipc
2016-01-17 14:50:33 +00:00
1caec2e608
calc done is used only internally
2016-01-17 14:15:38 +00:00
3d2187cc07
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 13:29:42 +00:00
ff840ae8ad
remove superfluous code and commented
2016-01-17 11:14:06 +00:00
11df75dfb2
added documentation and changed some names
2016-01-16 17:27:19 +00:00
83a6d00c60
remove superfluous variables
2016-01-16 07:06:34 +00:00
208196d213
introduced pointers as names. thx to Chuanlai for doing the work
2016-01-15 15:26:24 +00:00
117c8edfbd
corrected comment
2016-01-15 14:03:24 +00:00
684992bf50
exchanged mappingConstitutive array for more descriptive "phaseAt and phasememberAt"
2016-01-15 00:19:44 +00:00
4b10e4792e
adjusted function for polar decomposition to actual need (no return of U)
2016-01-12 21:39:31 +00:00
f090a1b216
generalized dimension of calls to lapack, for performance reason special 3x3 variants will follow
2016-01-12 11:00:23 +00:00
2eafefe652
removed norm functions from math in favor of intrinsic and simplified (mostly by using existing functions, merge intrinsic and array constructors)
2016-01-10 13:34:26 +00:00
519cd29c6f
using norm2 function as L2 norm (vector, math_norm3) and Frobenius norm (tensor, math_norm33)
2016-01-09 16:01:30 +00:00
53d07215a9
reverting back to converged Lp as the initial guess for next iteration
2016-01-09 15:03:18 +00:00
885b9095c4
option not supported by slightly older gfortran versions. thx to Luv for pointing out
2016-01-09 12:33:33 +00:00
5378dbdcc9
cleaned
2016-01-09 12:12:31 +00:00
b5ab94a0f1
added new plasticity module
2016-01-09 12:12:05 +00:00
7ba031adb3
fixed wrong Li and dLi_dTstar_3333 initialization when Tstar is zero.
2016-01-09 03:41:56 +00:00
a30b596e7f
combined J2 and I2 isotropic mechanical response. Default is identical to old plastic_j2. /dilatation/ flag adds (spherical) Li calculation based on hydrostatic pressure.
2016-01-08 19:45:20 +00:00
6cf92913a3
use current Lp for next iteration instead of converged Lp inherited from end of former increment as starting guess
2016-01-08 19:01:53 +00:00
1d1672ebfc
added math_spherical33 to calculate hydrostatic part of a 2nd order tensor
2016-01-08 18:57:37 +00:00
cbb6b7a689
adding fpe report for gfortran, some compiler opions are now used for linking (Debug + Optimization)
2016-01-07 12:58:19 +00:00
1857e47f75
div by zero in output
2016-01-07 11:48:30 +00:00
86885cb2ba
wrong initialization order, thanks to N Grilli for reporting the bug
2016-01-06 16:46:37 +00:00
268c7f359c
wrong index in output
2015-12-23 21:15:49 +00:00
74086fd0fb
correct initialization
2015-12-23 16:13:08 +00:00
7127403bd9
usused variables
2015-12-23 11:25:25 +00:00
dfe09f11f5
cleaner handling of NaN
2015-12-22 10:03:15 +00:00
ee47c9302a
small readability improvement
2015-12-22 09:22:57 +00:00
e127bd8db4
should reflect name of parameter
2015-12-21 07:27:21 +00:00
2a2f558416
update physics in phenoplus module
2015-12-17 16:11:13 +00:00
5f37ade722
followed some suggestions from gfortran
2015-12-16 22:42:08 +00:00
5901af1212
removing code fragments after introducing pointer state
2015-12-16 06:15:21 +00:00
9baec1a0ff
get rid of overly long file names 1/2
2015-12-15 20:45:54 +00:00
c710709b1a
example Makefile typically does not work, deleted
2015-12-15 14:34:03 +00:00
37aa559bd4
fixed missing structure type in material.config, small changes in the other code
2015-12-15 14:30:17 +00:00
88999c847a
reverted handling of NaN in converged reason, seems to be a compiler bug
2015-12-14 20:04:59 +00:00
dc990a88d6
unused variables
2015-12-14 19:48:11 +00:00
3a7f4bf43a
added new set of accelerated spectral methods based on derivative approximations. use spectral_derivative to set the kind of derivative you like. valid options are:
...
‘continuous’: computed continuous derivatives (standard spectral method)
‘central_difference’: central difference approximation of derivatives
‘fwbw_difference’: forward difference for gradient and backward difference divergence
2015-12-14 18:12:09 +00:00
da913ec3c7
bug in dot state, thx to Haiming for pointing out
2015-12-11 18:36:58 +00:00
76fe37e496
commented and respected compiler warning
2015-12-08 18:10:06 +00:00
4174e6f9e2
removed dislokmc, use disloucla for tungsten simulations
2015-12-05 22:07:03 +00:00
6ecbd3ae56
makefile had problems with static flag, mentat spectral Box was still using table.croak
2015-12-02 20:02:19 +00:00
0e76ef85c1
checking return code in util.exececute, cleaning makefile, deleted dummy docs
2015-12-02 13:21:17 +00:00
22a9d65825
spectral solver accepting NaN in SNES converged reason as this might happen for homogeneous microstructures
2015-12-01 22:36:19 +00:00
3365edc6cb
simplified and updated abaqus
2015-12-01 15:13:33 +00:00
4eb43908c3
some simplifications as discussed with IBF colleagues
2015-11-25 20:55:17 +00:00
a6cc799375
Removed unnecessary variables
2015-11-17 16:55:23 +00:00
51059abaf0
Martensite volume fraction evolution
2015-11-17 16:30:06 +00:00
a8b157a87c
Did not initialize a variable
2015-11-17 15:20:36 +00:00
4a8477c2c3
Characteristic transformation shear
2015-11-17 14:50:59 +00:00
cf84ceddd8
Added Ndot0 for trans
2015-11-16 10:15:18 +00:00
302178ebb4
some trouble since clean was already in PETSc rules
2015-11-14 09:47:18 +00:00
9d9367ae11
still having trouble, incorporation of PETSC_ARCH not needed
2015-11-13 10:16:35 +00:00
4a52c59854
corrected/removed PETSc variables path
2015-11-12 23:13:25 +00:00
900b072781
removed unused math_qDot (deprecated anyway)
2015-11-10 21:56:06 +00:00
45fb209658
respect PETSC_ARCH
2015-11-10 21:18:53 +00:00
c5b28a22d4
small polishing
2015-11-10 21:18:13 +00:00
aefd8fa390
Added a few more trans dependent variables
2015-11-10 15:23:15 +00:00
2778331d7b
Added trans interaction matrices
2015-11-10 13:30:34 +00:00
4e4136cd58
Consistent variable names
2015-11-10 13:01:03 +00:00
9c91ba4608
Trans interaction matrices
2015-11-09 15:57:20 +00:00
2b10625f0c
Typo
2015-11-09 08:51:05 +00:00
ff7cd17b51
not needed, PETSc is doing the job
2015-11-08 09:43:43 +00:00
4d9a2f8f6b
introduce pointer referencing of plastic state for better readability
2015-11-06 17:00:00 +00:00
7b0c130d6f
started to introduce alias name via pointer for plastic state in phenopowerlaw for more comfortable access
2015-10-30 15:48:30 +00:00
2ae4d24fb3
Use rotated trans elasticity matrices
2015-10-27 12:32:03 +00:00
2ba15e462a
corrected values and typo
2015-10-26 19:27:57 +00:00
c91788d721
new intel compiler was compiling about *** in float output
2015-10-26 17:11:36 +00:00
1f356a6833
added parameters from indentation results
2015-10-22 17:45:28 +00:00
04841cb727
List of dependent trans states
2015-10-21 13:24:26 +00:00
20b0429517
Calculation of equilibrium separation of partials for twin and trans
2015-10-21 09:37:45 +00:00
565d68edd7
Added some TRIP inputs
2015-10-20 15:54:01 +00:00
b9fcd9482b
Renamed some input variables to denote twin
2015-10-20 13:45:06 +00:00
fc8a11f027
corrected typo: alias for first Piola stress should read 1stpiola NOT 1piola
2015-10-19 06:57:23 +00:00
f2b29e4844
correct communication between
...
elements with different
phases
2015-10-14 21:04:40 +00:00
484a34b7f1
added pheno+ module
...
modify crystallite microstructure call
to pass orientations
2015-10-14 18:36:19 +00:00
ea0b16c726
pheno plus for commercial FEM
2015-10-13 21:06:59 +00:00
d6abc00218
add pheno+ module in
2015-10-13 18:52:01 +00:00
8fac635c15
made hybridIA stuff working again
2015-10-13 17:02:07 +00:00
6de8c821f8
improved consistency check on loadcase
2015-10-12 19:02:42 +00:00
9ec7c08bb6
gfortran was complaining abound mismatch
2015-10-07 13:07:21 +00:00
c1a53bc8c0
introduced dummy target to check compilation
2015-10-06 18:00:17 +00:00
0f0b4eaa6b
cleaned
2015-09-28 17:47:00 +00:00
ad16162eb8
simplified
2015-09-24 17:50:11 +00:00
c50a522dd8
some shapes not 100% correct, removed inverse laplace
2015-09-24 17:38:49 +00:00
2e70e6fded
make clear that this isNaN is an intrinsic extension
2015-09-24 08:45:44 +00:00
b528088653
added all 13 slip families for bct(beta)-Sn
...
increased the allowable columns for lattice.f90 to allow for the large interaction matrix.
2015-09-22 18:42:23 +00:00
42ac7902f7
indents wrong in fortran, seeds_fromDistribution bugfix and using ASCII table facilities
2015-09-20 21:03:10 +00:00
41f09dd0f8
should make no difference, but test with ifort was failing. compierler buig?
2015-09-12 18:26:25 +00:00
03896fcd09
redefined required debug variables removed in last commit
2015-09-11 18:20:37 +00:00
840fd0e9fa
report absolute value of accumulated_shear_slip in dot state
2015-09-11 17:44:26 +00:00
8716752bf8
fixed fourier convolution and div/curl calculation for even/odd grid according to Johnston 2011 (MIT, FFTW)
2015-09-11 13:05:46 +00:00
8f32d03a9e
curl calculation overestimated RMS due to factor 2 instead of one for DC component and Nyquist component
2015-09-11 09:19:13 +00:00
e88cedc6ae
clearer naming, debug options for spectral do not work for MPI
2015-09-11 08:52:03 +00:00
751d1d7582
function statement is more clear
2015-09-10 13:01:33 +00:00
580bb76b4c
pInt missing
2015-09-10 10:26:09 +00:00
427b5f4bc1
wrong results for restart
2015-09-10 10:06:14 +00:00
bbe37c842e
name conflict
2015-09-09 21:52:00 +00:00
34980a1d44
more logical structure when reading in
2015-09-08 19:44:32 +00:00
7d996cfa11
fixed new IO_stringPos for Abaqus
2015-09-05 16:26:55 +00:00
caf0ac7e8d
more verbose
2015-09-01 16:53:48 +00:00
7a3a67601f
realloc lhs was not always working, fixed now and added test for new string pos function
2015-08-31 16:30:04 +00:00
197ae53553
lots of memory savings (~ stateside + 250 pReals per integration point) when using analytic jacobian
2015-08-28 10:55:38 +00:00
6848d83d13
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name.
...
functions in the core module are not used, including kdTree
2015-08-28 07:38:48 +00:00
1171dc4344
MPI has a 2GB limit for writing at once, now chunking
2015-08-21 17:51:05 +00:00
8474da26d9
empty line
2015-08-18 16:31:47 +00:00
bae5cfcf89
named if
2015-08-18 16:07:01 +00:00
53c574033a
size of MAXCHUNKS same as maximum allocated size for “heat_time” & “heat_rate”
2015-08-18 03:20:16 +00:00
f64a5e7552
output safe for 64bit integer
2015-08-14 05:03:54 +00:00
31ea4dadf1
IO_stringPos replacement not needing MAXNCHUNKS but making use of new Fortran features for dynamic allocation
2015-08-13 14:54:34 +00:00
563d9e64dd
MAXNCHUNKS capitalized since it is a parameter
2015-08-06 09:24:56 +00:00
3d94ae4074
finite strain correction for respective expansions
2015-08-05 14:12:54 +00:00
09c7affbbc
more efficient rate calculation in updateState when subdt is zero
2015-08-05 11:09:38 +00:00
7928ce3ac3
fixed erroneous rate calculation in updateState when time increment (subdt) was zero during initialization.
2015-08-04 21:27:30 +00:00
65d114e4f0
improved debug statements
2015-08-04 21:26:22 +00:00
905b9f8f42
petsc update (now v3.6.1)
2015-08-04 15:04:53 +00:00
5b57ea7d16
fluctuations proportional to thermal equilibrium concentration
2015-08-03 11:08:59 +00:00
a01d7d8d66
corrections to crystallite_push33ToRef function
2015-08-03 11:07:19 +00:00
1ed276ce2d
quadratic softening for finite dissipation energy
2015-07-30 10:29:25 +00:00
c97b83c9a1
added user defined external heat source
2015-07-27 11:09:37 +00:00
4bb6664de9
tightened default thermal tolerance
2015-07-27 11:08:43 +00:00
c82c4b74f1
updated config files to reflect recent changes in input parameters
2015-07-24 17:14:03 +00:00
5e09954575
remove deprecated spectral load case definition of temperature
2015-07-24 14:57:29 +00:00
7554647c8e
more control over initialisation of field values. specify initial field value in the homogenisation part of the material config file using the appropriate tags
2015-07-24 14:53:50 +00:00
87d42bf447
initialise Fi correctly for initial field values away from equilibrium
2015-07-24 14:47:18 +00:00
0a5ccb3d91
possibility to set damage potential order
2015-07-24 14:43:05 +00:00
7c101cdc31
fixed issues causing compilation errors
2015-07-23 19:22:03 +00:00
7f8a8d5b0f
new function to report initial thermal strain based on current temperature deviation from reference.
...
to be used in constitutive_initialFi, which collects from all kinematics, derives and sums Li terms to report initial Fi.
2015-07-23 15:59:25 +00:00
25a03e128a
added remark to include looping over sources of same type within one phase. (for instance to accommodate different activity periods.)
2015-07-16 17:20:18 +00:00
d2c748e716
TRIP input parameters for fcc to hex
2015-07-10 09:02:48 +00:00
1e9682eef9
corrected upper bound for spectral damage solver as well.
2015-07-09 15:51:34 +00:00
d2bfd26937
1) added the error message of all MPI plans 2) deleted the suffix (MPI) for all MPI plans
2015-07-09 13:38:21 +00:00
7a4415b28a
Sanity checks for SFE
2015-07-08 15:37:16 +00:00
df5cdb6b01
Elasticity matrices for transformed phase
2015-07-08 11:58:52 +00:00
66564c1f75
Read in elastic constants for transformed phase
2015-07-07 13:23:56 +00:00
4435244681
corrected missing initialization. newState = oldState + deltaState. But for irreversibility driving force must be compared with last converged subState.
2015-07-06 09:34:43 +00:00
3d5a85e4af
corrected implementation of deltaState for isoBrittle damage
2015-07-06 08:47:40 +00:00
Martin Diehl