Commit Graph

351 Commits

Author SHA1 Message Date
Philip Eisenlohr 5d21ef5a84 fixed selective debugging logic flaw in crystallite_stressAndItsTangent 2014-05-19 16:43:32 +00:00
Philip Eisenlohr 88f94cbbec increased debug threshold for some outputs (decreased verbosity level at "basic") 2014-05-15 13:01:54 +00:00
Martin Diehl 4bfced1a70 fixed small allocation flaws for new state 2014-05-12 13:00:37 +00:00
Martin Diehl 838a8be1b9 removed forgotten debug statements from last commit
fixed material.config for abaqus example
2014-05-12 09:28:32 +00:00
Luv Sharma d820a5aaa6 fixed some errors in changes related to newstate 2014-05-12 00:44:44 +00:00
Martin Diehl 4bb5439994 removed copyright from *.f90 and *.f files, will be added using the deployMe script for release/nightly build. Now based on presence of svn:keywords property set (usually "Id"). Similar thing will be done for .py, .sh, and .config files but more involved because they might be executable 2014-03-29 08:20:36 +00:00
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
Franz Roters 7ee368f324 check whether there is any constitutive output 2012-05-11 07:38:20 +00:00
Christoph Kords 91a70b0fb3 before turning terminally ill report at which g,i,e crystallite did not converge 2012-03-30 13:16:45 +00:00
Pratheek Shanthraj 131c9ac93e dPdF calculations made consistent with constitutive_TandItsTangent 2012-03-21 15:00:36 +00:00
Franz Roters 8a2f2c5a95 stress iteration loop now uses generalized elasticity by calling TandItsTangent
the anlalytical tangent calculation should now be adopted to also use TandItsTangent
2012-03-21 10:57:27 +00:00
Pratheek Shanthraj c2f5cebacb simplified analytic jacobian calculation. removed Lpfrac, time_sensitive. introduced rate_sensitivity flag when calling crystallite_stressAndItsTangent that is currently set to .false. and is to be set according to which dPdF the FE solver is asking for 2012-03-14 13:56:50 +00:00
Franz Roters c9ed989dad replaced keyword constitution with plasticity in preparation of upcomming generalized elasticity 2012-03-12 14:09:37 +00:00
Martin Diehl bd9667bd4b added new, flexible debugging scheme.
now all modules have their own debug specification.
compiles and runs, I hope nothing is broken
did a lot of polishing
2012-03-08 20:25:28 +00:00
Pratheek Shanthraj 17e9698659 added time_sensitive flag to ignore or add time sensitive component to analytic jacobian (should be close to perturbed jacobian when ignored)
now uses sub_dt for analytic jacobian calculations
bug fixes
2012-03-01 19:42:43 +00:00
Nan Jia 211f66a728 added elastic stiffness matrix as crystallite output 2012-03-01 13:06:09 +00:00
Pratheek Shanthraj 040e244993 Added stress-order terms to analytic stiffness 2012-02-27 17:40:28 +00:00
Pratheek Shanthraj 99d6dcecb5 fixed bug in calculation of analytic jacobian (should work much better now). parallelized analytic jacobian calculation loop 2012-02-24 14:07:46 +00:00
Pratheek Shanthraj 28d654d61e forgot to remove some print statements in previous commit 2012-02-22 20:22:54 +00:00
Pratheek Shanthraj b9a82156c9 crystallite.f90: can now calculate analytic jacobian by setting analyticJaco = 1 in numerics.config
math.f90: added math_mul3333xx3333…

numerics.f90: to read in analyticJaco and Lp_frac
2012-02-22 20:11:09 +00:00
Martin Diehl d68a3d1869 added LF=UNIX to a bunch of files 2012-02-22 13:30:00 +00:00
Martin Diehl d8ffc29236 again: adding pInt, removing redundant use statments, chang in dble to real(,pReal) 2012-02-21 16:31:37 +00:00
Martin Diehl 6c0f9d163b polishing: removed variable names like 'unit' and 'data' that are keywords of fortran and ensured that integer and real precision matches independent of machine standard.
removed cut_off parameter for damask_spectral
removed outpot of derived divergence measures and added RMS output in brackets
added comments and options to the makefile
2012-02-15 18:58:38 +00:00
Christoph Kords c786336af3 reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore.
all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 17:41:27 +00:00
Krishna Komerla 94100e8d8e missing line break 2012-02-13 15:21:03 +00:00
Martin Diehl 156ec4582a polishing, adding _pInt etc. where applicable
post_results now handels zero increment different (like FEM, it is always there even if it is not counted)
2012-02-13 14:18:07 +00:00
Philip Eisenlohr fa32a00d56 changed symmetryID initialization to integer (instead of float) value 2012-02-10 17:46:43 +00:00
Christoph Kords e7d407c2f6 crystallite_Fp has to be initialized in crystallite_init, cause it is used by constitutive_microstructure 2012-02-03 14:52:09 +00:00
Martin Diehl c5128e94eb some f2008 standard related correction i didn't check in last time 2012-02-01 20:20:05 +00:00
Martin Diehl 800e291240 made code standard conform to Fortran 2008 (ignoring warning concerning comments beyond character 132). Basically, changing "x" format specifier to "Nx" ("1x") plus removing $ format specifier
added compiler switches for gfortran and ifort to check for standard conformity
old gnu compilers <4.4 are not longer supported because they don't provide the c binding for fftw
2012-01-31 19:18:55 +00:00
Martin Diehl 80583fefb5 added compilation_info.f90, to store compilation date,time and compiler info in each file
added #include statement at each init() routine
2012-01-31 14:54:49 +00:00
Christoph Kords 1330576a01 added new math function "math_deviatoric33" which returns the deviatoric part of a 3x3 tensor
renamed some math functions, so that we have a universal naming scheme: for matrix multiplications use an "x" (e.g. math_mul33x3); don't use the "x" to describe the shape of the tensor that the function is applied to (e.g. math_invert33 instead of math_invert3x3)
2012-01-26 13:50:00 +00:00
Christoph Kords 6fec57bbd7 replaced variable fraction by "timefraction" 2012-01-26 12:48:59 +00:00
Christoph Kords 4980f57fd5 added total Green-Lagrange strain to crystallite outputs 2012-01-20 10:25:35 +00:00
Christoph Kords 8739b9522e dislocation stress does not add to stress field anymore, but is only available as constitutive output for purposes of postprocessing. instead, we now use a backstress term (which is added to the resolved shear stress) that depends on the gradient of excess density on the same slip system. 2012-01-17 10:26:57 +00:00
Philip Eisenlohr a7d9d711f6 edited lines exceeding 132 chars (before any comment started) 2012-01-11 16:56:35 +00:00
Martin Diehl 5ebeb96e85 made DAMASK to work with gfortran:
-removed to long lines
-restructured f2py modules and merged make_DAMASK2Python into setup processing
-setup_code.py now sets library path in makefile and asks for compile switches for spectral code
-substituted \ in format strings with $

restructured DAMASK_spectral:
-more logical output and structure of code
-better input for spectral debug parameters
2011-12-06 16:58:17 +00:00
Christoph Kords 9d1bc584d0 corrected stress integration scheme: now use norm(Lpresiduum) as a target function for Armijo's rule instead of whole tensor Lp; also corrected the guess for the improvement in Lp 2011-11-04 12:57:12 +00:00
Christoph Kords ca3d21a3b6 Depending on which state integrator one uses for the stiffness calculation, the initial state has to be chosen accordingly: e.g. for FPI choose last converged state, but for explicit RK choose converged state from start of increment (in case of explicit euler no state integration at all, but only stress integration). For this purpose we also need to remember "Fe" which now follows the cutbacking procedure as it is used for "Fp". 2011-11-04 12:44:50 +00:00
Martin Diehl c35ea33f8e did a lot of polishing:
- removed unnecessary "return" before end of subroutine or function:
- changed undetermined array length (:) to (1:3)

To prevent problems with some code analysing tools:
- "3D oneliner loops" (with ";) only for "do" and "enddo" at the same time
- removed line continuation in OMP statements

made the makefile more flexible, removed heap-arrays switch
2011-09-13 15:54:06 +00:00
Christoph Kords b300cc7faa ip volume is now based on the determinant of F. "mesh_ipVolume" represents the initial volume and is multiplied with det(F) wherever the current volume is needed. Since this works for all solver types, the "volume" output in crystallite is now also correct for spectral method and abaqus. 2011-08-19 11:18:56 +00:00
Christoph Kords 1ffb59a96a Calculation of current ip volumes now working. Crystallite output also reflects current grain volume, not reference volume. However, this is only available for Marc. Abaqus and spectral method still return the reference ip volume. The ip coordinates though are correctly updated for all solver types. 2011-08-10 16:37:17 +00:00
Christoph Kords ef7405fe21 Stress integration now uses Armijo rule to find an appropriate correction of Lp: decreases step in case that residuum does not improve significantly, accelerates as usual in case of good convergence. This turned out to improve convergence behavior. 2011-08-02 11:29:08 +00:00
Martin Diehl 75c67f53f7 introduced alternative location for material configuration.
If present, the >>>JOB.materialConfig<<< file takes precedence over the standard >>>material.config<<<
2011-08-02 10:14:16 +00:00
Martin Diehl 564eb5009f shape mismatch in crystallite, invalid names for variables(pos, size) in homogenization corrected. polishing of DAMASK_spectral_interface and makefile 2011-08-01 18:10:55 +00:00
Philip Eisenlohr a26d9c844b max size of crystallite result is now based on actual use within the mesh, not just picking the largest chunk from the material.config file... 2011-06-08 09:16:48 +00:00
Philip Eisenlohr 8041587a72 changed internal debug verbosity in accord with debug.config listing. 2011-06-06 15:27:35 +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
Philip Eisenlohr affd383ef8 simplified disorientation calculation. performed FEM check that "grainrotation" reflects the *active* rotation from the starting orientation to the current one with respect to the lab coord system (as expected). 2011-04-14 19:09:44 +00:00
Christoph Kords 5814262f55 moved some write statements within stateIntegrationEuler 2011-04-14 09:35:56 +00:00
Philip Eisenlohr 349f022100 reinstantiated calculation of grainrotation (still needs check for correctness, i.e. active or passive rotation output) 2011-04-13 16:08:52 +00:00