Martin Diehl
7c79b31f6c
remove not needed use statement
2014-03-13 05:50:56 +00:00
Martin Diehl
9afc1e3a15
fixed wrong omp statement
2014-03-12 15:29:14 +00:00
Martin Diehl
ff1b1c1a50
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code
...
forgot to commit dislotwin last time, now seems to work
2014-03-11 23:55:40 +00:00
Martin Diehl
2b589c3d71
moved reading in of lattice type and elastic constants to lattice module
...
removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config
constitutive modules will only be initialized if needed
homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-08 20:50:31 +00:00
Martin Diehl
a0d75ee05e
moved quaternion disorientation to lattice because it requires the knowledge of the lattice structure.
2014-02-28 13:28:27 +00:00
Martin Diehl
d45aea4467
moved reading in of lattice stru
2014-02-10 14:31:19 +00:00
Christoph Kords
b9b87a785c
added crystallite_F to debugging output
2014-02-03 13:19:49 +00:00
Christoph Kords
391d83fa2b
moved rate sensitivity block to where it belongs to, namely the analytic jacobian calculation
2014-01-22 15:34:10 +00:00
Pratheek Shanthraj
4b7890f133
fixed wrong indices in analytic jacobian calculation following christoph's changes. convergence is restored at least for the spectral solver
2014-01-22 10:16:55 +00:00
Christoph Kords
3dfdbaff5b
Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged.
...
In contrast, the analytical jacobian will probably be affected by this change!
@Pratheek: Can you check with me how this can be fixed?
2014-01-22 08:38:13 +00:00
Christoph Kords
c10d0d15cb
crystallite_integrateStateFPI: convergence check always ran about all e,i,g; now checking only those specified in FEsolving_execElem, FEsolving_execIP
2014-01-21 18:45:41 +00:00
Christoph Kords
ff648492c1
perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode
...
(missed some part in last commit)
2014-01-21 18:28:21 +00:00
Christoph Kords
029826d84f
perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode
2014-01-21 16:35:12 +00:00
Christoph Kords
351aa6ba29
setting todo and converged flag only for those e,i,g that are to be executed (FE_elemExec etc)
...
used to spoil convergence checking in non ping pong mode and lead to slow performance
2014-01-21 16:09:00 +00:00
Christoph Kords
77cba87bf3
in crystallite_integrateStress: pass full 2nd PK stress to LpAndItsTangent instead of only the deviatoric part; otherwise constitutive laws that depend on the hydrostatic pressure do not work; constitutive laws that require a deviatoric stress need to calculate it internally (as is e.g. already done in the J2 model)
2014-01-16 10:36:40 +00:00
Martin Diehl
831e905b16
added missing lp output to crystallite
2013-12-20 11:13:12 +00:00
Martin Diehl
af075aaa9d
found bug (array index out of bounds), Philip/Franz please check if correct (Gfortran was complaining)
2013-12-20 08:33:46 +00:00
Pratheek Shanthraj
6f6bd86078
fixed bug in heat generation and analytic jacobian calculations
2013-12-18 10:19:48 +00:00
Martin Diehl
d5bdcf5a6d
some small improvements, still something wrong, probably in crystallite.f90
2013-12-13 14:14:17 +00:00
Martin Diehl
84d51d743c
fixed small issues, ATTENTION: Still a bug introduced with rev. 2783 not fixed
2013-12-13 13:19:17 +00:00
Martin Diehl
65ae979920
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 17:09:59 +00:00
Martin Diehl
0bc99a9622
change latticeName to latticeID
2013-11-27 16:20:27 +00:00
Martin Diehl
2f7efa2055
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 08:04:05 +00:00
Christoph Kords
883669bd77
in preguess of FPintegrator: resetting of previous dotstates can be done in same loop as collectdotstate
...
polishing: indentation level and capitalization
2013-11-21 10:58:41 +00:00
Pratheek Shanthraj
724380c624
Coupled phase field implementation for the basic petsc solver. needs clearing up… work in progress
2013-11-12 00:44:23 +00:00
Martin Diehl
3ecc8103f0
introduced dummy heat calculation, available as an crystallite output
...
temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-18 18:57:28 +00:00
Martin Diehl
03e26af37b
small bufix, didn't compile
2013-10-16 20:44:04 +00:00
Martin Diehl
dc95c82d4a
removed temperature integration and corresponding data structures and debugging options
...
temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver.
introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver.
went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 13:04:59 +00:00
Philip Eisenlohr
42cd4eb38f
fixed some e12.6 string formats to avoid compile time warnings
2013-10-08 13:51:36 +00:00
Martin Diehl
3bf0ed84c8
some comments and minor improvements.
...
removed tol_gravityNodePos from prec because it's not used
tol_math_check is now used by spectral driver to check rotation
2013-09-14 10:59:35 +00:00
Martin Diehl
c93373a4ac
removed long lines
2013-08-08 19:33:21 +00:00
Franz Roters
4f9dbfa193
last commit also solves problem with infinite loop in case of a cutback in homogenization
...
corrected indentation of a comment
2013-08-02 07:59:55 +00:00
Franz Roters
a08d263865
remnoved some unused fuctions from use statements
2013-08-02 06:18:41 +00:00
Philip Eisenlohr
98528f9a89
added possibility for multi-level inclusion of files in *.config and loadcase files.
...
include subfiles by stating
{path/to/include}
2013-06-26 19:19:00 +00:00
Martin Diehl
6f7740a243
made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables
2013-06-11 16:35:04 +00:00
Christoph Kords
55735b9580
grainrotation output now uses corrected rotation math
2013-06-07 09:35:00 +00:00
Martin Diehl
c7ba8a2a9b
introduced plain mode (no ping pong) again and added test for MSC.MArc 2012 as prove that its working
2013-05-17 17:52:46 +00:00
Martin Diehl
cfac73043d
some lines in crystallite.f90 exceeded limit of 132 characters
...
fixed bug in Lib_VTK_IO.f90 (informed author as well)
added zip file containing source of Lib_VTK_IO.f90
2013-05-08 17:48:00 +00:00
Christoph Kords
3b788720e5
added crystallite output "neighboringelement" and "neighboringip"
...
some polishing in crystallite_postResults
2013-05-08 12:02:30 +00:00
Martin Diehl
b98ce3d22d
fixed out of range bug introduced with singleRun changes in rev. 2356
2013-04-30 12:14:07 +00:00
Martin Diehl
533fcec96f
named some loops and made BUTCHER TABLEAUs parameters in crystallite_integrateStateRKCK45
2013-04-29 11:17:30 +00:00
Martin Diehl
9fb6d8167b
simplified logic for single run, now not checking for extra arguments but simply checking if start and end elem/IP are the same
2013-04-26 13:23:36 +00:00
Christoph Kords
a7706090d2
neighborhood relations were screwed up for some elementtypes as of version 2318, might have caused segmentation faults with nonlocal model, now fixed
2013-04-22 13:35:35 +00:00
Martin Diehl
8b2d9d8155
further improvements on implementing the plain mode
2013-04-16 17:07:27 +00:00
Christoph Kords
44eb771583
added crystallite output of current ip coordinates
2013-03-27 10:12:28 +00:00
Franz Roters
02f092492e
updated copyright header
2013-03-22 17:35:05 +00:00
Martin Diehl
7e0ea5dd04
hickups in doxygen documentation fixed
2013-03-05 14:35:26 +00:00
Martin Diehl
ada2beb8b8
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 11:48:29 +00:00
Christoph Kords
875abcd06f
replaced some array assignments, which might have caused trouble when using the "timeSyncing" option, by parallelized loops of scalar value assignments
2013-02-27 16:36:12 +00:00
Christoph Kords
3736ab1d97
added shortcuts for two do loops in the "timesyncing" procedure
...
math_pDecomposition enclosed by omp critical statement; somehow, this subroutine is not thread safe, not clear what exactly the problem is
2013-02-27 10:32:37 +00:00
Pratheek Shanthraj
07ff1f7c93
reporting current time for most module inits
2013-02-25 16:34:59 +00:00
Martin Diehl
393c60565d
one line was beyond 132 characters
2013-02-22 05:38:02 +00:00
Martin Diehl
42b96354db
doxygen comments
2013-02-21 23:08:36 +00:00
Martin Diehl
f0b4281400
removed unused variables, added pure statements, and declared external functions as external (using gfortrans debug option)
2013-02-11 10:43:45 +00:00
Martin Diehl
e74b5da19a
changed all remaining routines to fortran-fast arrays (geometry reconstruction etc.)
...
changed all remaining routines in f2py to more clever determination of array size (requires f2py >= 2.0)
enabled 3D visualize to work with odd resolution by switching to linear reconstruction
PLEASE NOTE: Redefinition of routines for f2py might cause trouble -> DELETE DAMASK_ROOT/lib/damask/core.so in this case
further changes: added pure statement where possible, polished, unified use of "Q" for "Quaternion" and reordered math to have similar routines together
2013-01-31 16:28:08 +00:00
Martin Diehl
3ffd6499b1
change bash-only "let" to "`expr`", abaqus_v6.env now suppresses warnings about long lines (there are only comments anyway)
2013-01-28 15:30:51 +00:00
Pratheek Shanthraj
fd94c786f0
moved stiffness tensor calculation to lattice
...
introduced 'isotropic' and 'orthorhombic' lattice types to use corresponding symmetries in stiffness tensor. intended to be used with non-crystal plasticity models (j2, constitutive_none with isotropic, cubic or orthotropic elasticity).
2013-01-21 21:57:26 +00:00
Christoph Kords
8a45a90775
new crystallite output "grainrotationX", "grainrotationY", "grainrotationZ" that gives the deviation of the initial orientation in degrees around the sample reference x/y/z axis, while "grainrotation" still gives the same as axis angle pair in crystal coordinates
2013-01-16 08:45:41 +00:00
Christoph Kords
7dd1130e92
reset syncSubFrac flag when turning terminallyIll
2012-12-16 16:18:36 +00:00
Christoph Kords
92d2c83334
more debugging info for timesyncing
2012-12-16 11:30:19 +00:00
Christoph Kords
4f12073dcb
once again corrected parallelization statements
...
added some debugging writes for timesyncing
2012-12-16 10:54:13 +00:00
Christoph Kords
f744fc1d47
fixed bug in parallelization: flushing of crystallite_todo was missing at a couple of places
2012-12-11 15:05:08 +00:00
Christoph Kords
6a9b19e4c0
fixed bug in timesyncing procedure: if any ip that synchronizes its time step did not converge, all nonlocals become terminally ill; before, this led to an infinite loop in crystallite
2012-12-11 13:42:22 +00:00
Christoph Kords
35a0483507
fixed bug in last revision: FE_geomType was not in namespace
2012-12-10 17:22:16 +00:00
Philip Eisenlohr
10ed848258
fixed 2 bugs: missing mapping with FE_geomtype(mesh_element(2,...))
2012-12-10 17:03:28 +00:00
Christoph Kords
98d02dcffb
removed call to microstructure in stateJump function, since this is already done in the state integrator subroutine
2012-12-07 07:50:41 +00:00
Christoph Kords
0986a66449
accidentally also included some changes in crystallite_integrateStress in last commit; reverted these changes in this revision
2012-12-03 16:34:41 +00:00
Christoph Kords
f2dacb5f40
polished debugging output
2012-12-03 16:29:54 +00:00
Christoph Kords
f92d34dd22
removed some unused variables
2012-11-29 18:44:00 +00:00
Christoph Kords
e10000a338
Introduced an alternative cutback scheme for the nonlocal model, that allows to keep the results of most of the ips that immediately converged, and only do cutbacks in regions where some ips are in big trouble. Only works for nonlocal model and state integrator 2.
2012-11-27 18:36:55 +00:00
Christoph Kords
e584c0e6e8
deleted space in front of preprocessor line
2012-11-22 20:04:33 +00:00
Christoph Kords
38ec462f4a
order of integration now same for all state integrators: first calculate the dotState and then the deltaState
2012-11-22 13:04:19 +00:00
Philip Eisenlohr
7d196fbb25
polishing
2012-11-22 09:58:36 +00:00
Philip Eisenlohr
d189a6a7a0
polishing
2012-11-21 16:57:57 +00:00
Christoph Kords
760af1fbbf
crystallite_requested now masked by execIP, now elemhomogeneous works again
2012-11-19 15:13:48 +00:00
Christoph Kords
9fb94cd707
breaking crystallite loop in case of terminally ill (this functionality was probably lost since rev1886)
2012-11-19 08:39:48 +00:00
Philip Eisenlohr
d9a98417ca
switched element library to geomType based.
...
saves to copy same geometry description for different elements that are essentially similar regarding the IP number but differ in total node count.
introduced quadratic tetrahedron (Marc element 127 -- element 157 might also work, but did not perform well in fully elastic calc so far)
2012-11-15 22:45:20 +00:00
Christoph Kords
11264f3fe4
forgot to check for convergence flag in FPI integrator
2012-11-08 20:15:19 +00:00
Christoph Kords
60438cdbcc
criterion for crystallite cutback loop changed from crystallite_subStep to crystallite_todo; this avoids problems when the last substep of the crystallite increment is below subStepMinCryst (which is allowed during last wind forward)
...
fixed bug in the adaptive Euler integrator and the RK4 integrator, which could not reach the stateJump function, because of erroneous check for converged flag; also corrected the state counter for the RK4 integrator
2012-11-08 20:00:29 +00:00
Christoph Kords
c775edaa6d
better avoid initialization of residuum_old with huge, since it posed problems in Abaqus
2012-11-08 13:26:22 +00:00
Christoph Kords
dad9922f54
fixed bug in crystallite_FPI: stateDamper always has to be defined for each grain
...
added some OMP FLUSH statements were necessary
replaced openmp do by forall construct where possible; this is much safer and perhaps even as fast for small loops
2012-11-07 15:43:29 +00:00
Christoph Kords
bb033c5fe7
Simplified algorithm of crystallite_integrateStress while preserving (almost) same functionality.
...
Removed "leapfrogging" (increase of step for next guess, when last guess was ok); Replaced Armijo rule testing for step size by simple check if the residuum got better, since the former virtually did not have any effect; consistently using the 2-norm of the residuum rather than infinity-norm for the convergence check throughout the function
2012-11-06 12:35:45 +00:00
Christoph Kords
639d6e0655
inserted missing flush(crystallite_todo) in parallel section; might have caused trouble in parallel calculations before
2012-11-05 14:45:52 +00:00
Martin Diehl
086fe138b1
substituted matrix inversion to solve equation by direct solution routine from LAPACK
2012-10-31 09:56:26 +00:00
Christoph Kords
2be331b74d
In crystallite for state integrators 2 and 3: now that "stateJump" is called before integrating dotState, we have to do state = state + dotState *dt instead of state = subState0 + dotState * dt; otherwise the deltaState is lost
2012-10-22 13:34:15 +00:00
Christoph Kords
0980920b66
In crystallite: for state integrators 2 and 3 the deltaState always has to be done at the beginning of the integration, not at the end; otherwise the stiffness calculation can be completely wrong, when the state is changed after the stress integration for the central solution.
...
Also changed some debugging levels in crystallite and mainly in constitutive_nonlocal
2012-10-22 07:59:35 +00:00
Christoph Kords
36676029e1
for state integrator 2 (explicit euler) and 3 (adaptive 2nd order euler): don't do any state evolution during stiffness perturbation, but only stress integration, since for those cases the state is not restored after every perturbation step
2012-10-19 11:20:31 +00:00
Philip Eisenlohr
dd5f453994
mainly cosmetics. added some comments. changed debug levels for some outputs. corrected (probably non-essential) check for NaN in FPI stressIntegrator.
2012-10-18 09:53:26 +00:00
Christoph Kords
5f0440fde3
crystallite_integrateStress now raises an error if Fp becomes NaN
2012-10-17 16:19:42 +00:00
Martin Diehl
0bcb8f59db
added doxygen comments, some polishing, added "protected" statements where applicable
2012-10-02 12:53:25 +00:00
Martin Diehl
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
Martin Diehl
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
Christoph Kords
ad7381cfdd
maximum relative step length in stress integration cannot decrease below 1
...
intermediate tensor "A" is calculated only once before the Lp loop in "integrateStress"
2012-08-21 08:53:36 +00:00
Christoph Kords
c1c5941733
stressAndItsTangent:
...
when winding forward to last substep in increment, also do the calculation if this substep is below "subStepMinCryst";
output of "wind forward" and "cutback" info only for crystallites that are still "todo"
integrateStress:
maximum component in Lp determines absolute convergence criterium:
aTol = max(maxLpComp * rTol_fromConfig , aTol_fromConfig)
this ensures that small values in Lp do not spoil the convergence, similar to the effect of the relevantStrain parameter before (not needed anymore in the stress integration)
2012-08-14 12:07:59 +00:00
Philip Eisenlohr
4d09ef0648
changed variable name "debug_what" to "debug_level"
2012-07-05 09:54:50 +00:00
Franz Roters
380f40561b
first calculation of Tstar still used hard coded Hooke, now usese constitutive_TandItsTangent(...)
2012-06-18 07:26:34 +00:00
Christoph Kords
fc7b4d6471
new function stateJump, which takes care of immediate changes in the state (deltaState); this introduces discontinuities in the state evolution; therefore this is always and only once done after each integration step, so no evaluation of deltaState for intermediate steps of e.g. the Runge-Kutta integrator; otherwise integration becomes a pain without significant gain in accuracy
...
deltaState now successfully tested for nonlocal model; was not correct for integrators 1,4,5 before
2012-06-06 15:11:30 +00:00
Christoph Kords
8c849d62ee
small corrections for how to treat the new deltaState in FPI and euler integration; all integrators successfully tested with j2
2012-05-18 13:00:15 +00:00
Christoph Kords
abbae76c51
instantaneous jumps in the state by constitutive_deltaState are now incorporated for all state integrators. still they (should) not influence the result, since all constitutive laws simply return zero for the deltaState
2012-05-17 15:25:21 +00:00
Christoph Kords
a54439e3b5
dotState does not have to be reset to zero. this is a remnant from older versions when the dotState for the nonlocal model used to be updated by the neighboring integration point
2012-05-16 15:35:14 +00:00