Commit Graph

828 Commits

Author SHA1 Message Date
Martin Diehl 0434ef3daa some comments to be seen in documentation 2012-08-15 13:38:38 +00:00
Pratheek Shanthraj c2c31a0a59 added SolverBasicPETSC module to the new FFT solver (solves the basic scheme using PETSc when mysolver = basicpetsc in numerics.config)
some bug fixes
2012-08-14 16:58:23 +00:00
Christoph Kords f5e7d576f1 introduced annihilation of screw dipoles; so far no thermal activation, but immediate annihilation; each annihilated screw dipole leaves behind an edge jog of same length as the dipole height on the collinear system 2012-08-14 12:26:20 +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
Christoph Kords eae9a34c6b corrected debugging output: extreme values of returned stress and jacobian belong to cpfem debugging, not homogenization 2012-08-14 11:48:43 +00:00
Pratheek Shanthraj 259eb7850e more work on array ordering in new solver. (now working but not tested) 2012-08-10 17:01:58 +00:00
Martin Diehl cd0da03ebc changed comments to doxygen style 2012-08-10 15:58:17 +00:00
Krishna Komerla 08c7be7d15 adaptive regrid now working, updates F and F_lastInc based on mapping. 2012-08-10 15:48:27 +00:00
Martin Diehl b324a98014 changed order of most arrays to fortran-fast, whole new solver still work in progress 2012-08-09 13:04:56 +00:00
Martin Diehl fe4d4d9525 changed comments in prec and IO to be read by doxygen, fine tuned output of doxygen and added precompilation for documentation 2012-08-09 11:01:53 +00:00
Pratheek Shanthraj d0933dad7b added filter option to the spectral solver to mitigate spurious oscillations due to gibb's phenomenon. activate by setting myfilter in config file appropriately (currently only 'none' and 'cosine' options coded). more cleaning up of AL code 2012-08-07 17:23:13 +00:00
Pratheek Shanthraj 925915e30d some cleaning up 2012-08-07 12:02:57 +00:00
Pratheek Shanthraj 32f8c9c6ea moved petsc option string to numerics and added new numerics parameters to config file 2012-08-06 17:27:53 +00:00
Martin Diehl 93658c72db set auto values back to auto 2012-08-06 12:44:07 +00:00
Martin Diehl 7e3a837640 moved calculation and output of geometry related data to mesh.f90, moved output of command line related information to DAMASK_spectral_interface.f90.
improved new solver structure (use with make SOLVER=NEW)
2012-08-06 12:43:05 +00:00
Christoph Kords a48d180f8c allow prefactor for self-diffusion to be equal to zero 2012-08-06 09:03:33 +00:00
Pratheek Shanthraj e2a2621235 updates for AL version of spectral solver 2012-08-06 08:53:12 +00:00
Christoph Kords e8df3678aa added possibility to assign the transmissivity of grain boundaries a fixed value by use of the keyword "grainboundaryTransmissivity" in the material.config; grain boundaries are identified as interfaces between areas of different material_texture 2012-08-03 14:32:49 +00:00
Martin Diehl 948d9c03d0 new version of modular solver structure 2012-08-03 09:25:48 +00:00
Mahesh Balasubramaniam 3ffba037d4 Reading of euler angles from the loadcase corrected. 2012-08-02 17:52:07 +00:00
Mahesh Balasubramaniam 740357172a Reading of Euler angles from the loadcase corrected. 2012-08-02 17:46:50 +00:00
Martin Diehl ee1bde0cd7 enabling regridding more than once by introducing deallocation of arrays
added J2 test (stub from Taymour)
2012-07-31 15:37:49 +00:00
Martin Diehl 5a89c783d4 deallocating already allocated arrays in case of regridding/restarting (i.e. when calling init more than once) 2012-07-31 14:00:15 +00:00
Martin Diehl 991a7bbd2d reordered initialization for spectral method, corrected bug of deformation BC parsing when prescribing velocity gradient resulting in wrong average deformation 2012-07-30 15:51:48 +00:00
Martin Diehl d4163dd16f further work on modularization, basic scheme is now calculating the same results as DAMASK_spectral.f90 2012-07-30 14:06:22 +00:00
Pratheek Shanthraj c341ddd855 more modularization…. 2012-07-26 13:58:47 +00:00
Martin Diehl 4ed68bb4ae next round in modularization 2012-07-25 14:01:39 +00:00
Martin Diehl e1faad8bcf next round in modularizzation of solvers 2012-07-24 17:07:10 +00:00
Philip Eisenlohr 882c44a61d changed the majority of debugging output to level "extensive" 2012-07-24 15:29:23 +00:00
Christoph Kords 15c5131874 corrected formulation of maximum height of screw dislocation dipoles 2012-07-24 14:50:11 +00:00
Martin Diehl 9c0a161ec0 in case of regridding, undeformed configuration is written out correctly. Restart improve by reading in additional data from converged step. 2012-07-24 12:21:41 +00:00
Christoph Kords c4c00f04f4 no dipole formation by stress decrease;
dipole dissociation rate proportional to "deltaDipoleLimit / oldDipoleLimit", not "deltaDipoleLimit / newDipoleLimit"
2012-07-24 06:57:37 +00:00
Martin Diehl 3d59fec305 further modularisation of spectral solver (only for basic scheme so far) 2012-07-23 10:12:31 +00:00
Martin Diehl bd9e81fbec intermediate, not working state of new solver structure 2012-07-20 15:33:13 +00:00
Martin Diehl ad3f9d8050 removed old standalone AL solver and introduced new structure for solver zoo 2012-07-19 17:24:56 +00:00
Martin Diehl 7d4ff278cb worked on the restart capabilities 2012-07-19 14:16:59 +00:00
Martin Diehl 13caf2d389 made linear shape reconstruction working again, pretty similar results like corrected FFT reconstruction 2012-07-18 18:39:59 +00:00
Martin Diehl a22163c2ae reverted unwanted commit of Makefile 2012-07-18 12:34:46 +00:00
Martin Diehl 95d292b286 2012-07-18 12:31:07 +00:00
Philip Eisenlohr 7c6fc121fd condensed error reporting for constitutive_XYZ_init
removed erroneous check for structure>3
2012-07-17 17:36:24 +00:00
Philip Eisenlohr 0fdd43b54b minor clarification of file content and its use 2012-07-05 10:02:20 +00:00
Philip Eisenlohr 4d09ef0648 changed variable name "debug_what" to "debug_level" 2012-07-05 09:54:50 +00:00
Philip Eisenlohr b69f446caa added constitutive_none
added wrongly deleted "endif" statements
removed tmglib and refblas from linker command for lapack (Koen??)
xxxROOT is changed by setup_code unless commented-out... So had to delete '#' again.
2012-07-03 13:59:38 +00:00
Philip Eisenlohr aeb1e18268 passing of multiple make options possible now 2012-07-03 13:57:05 +00:00
Philip Eisenlohr 7f10ab0996 corrected to consistent naming: "IMKL"
restored xxxROOT to commented-out auto entries.
commented-out currently unused PETSC parts.
added "refblas" and "tmglib" to LAPACK library inclusion (seems to be necessary since 3.2.1).
2012-07-03 12:36:54 +00:00
Franz Roters 823013d485 new plasticity: none
see material.config for necessary parameters (lattice_structure, elastic constants)
2012-07-03 11:16:38 +00:00
Philip Eisenlohr b4103b205d parts of the material.config can now be set to echo their content back into the log-file (i.e. STDOUT / unit 6)
use keyword "/echo/" before first section in the respective parts that you want to see echoed back into the log file.
2012-06-26 10:24:54 +00:00
Martin Diehl 767842acbf now real time monitoring of DAMASK_spectral output to screen
working with GNU Fortran as well (new checking for return code)
2012-06-25 09:03:19 +00:00
Martin Diehl 0606845f15 reordered calculation of current stiffness to enable restarting with different resolution 2012-06-22 09:20:23 +00:00
Martin Diehl bd462d12c6 fixed some bugs in the regridding routine 2012-06-20 18:33:20 +00:00
Martin Diehl d650972904 corrected some bugs concerning the regridding
prevented FEsolving from potentially write to a none existing file
started to introduce petsc into the make chain (nothing happens if PETSC_DIR is not set)
2012-06-20 12:49:46 +00:00
Nan Jia f5eb8fcb25 fixed bug caused by the inconsistent sequence of slip and twin systems 2012-06-19 16:06:25 +00:00
Martin Diehl 481268c8c2 fixed closing of file (unit 222) statement in FEsolving.f90
added missing keyword in mesh.f90 with ifndef Spectral statement (needed for Marc and Abaqus)
2012-06-19 14:33:24 +00:00
Krishna Komerla 06be437bc9 added minRes to regridding function and writing out of new geometry file
updated f2py wrapper to enable the use of init functions.
added 2 new error messages to io
2012-06-19 13:31:15 +00:00
Martin Diehl 8537e87b7e added shell scripts to run tests from crontab, improved numbering in spectral compile tests and added some first files for testing restart capabilities of spectral solver.
added configuration file for generation of code documentation using doxygen
restart of spectral solver is fixed and seem to work now also for restart at significant deformation.
spectral solver now gives exit code 2 if some increments did not converge
2012-06-18 15:27:01 +00:00
Philip Eisenlohr d9c6f9e134 switch on restartRead logical when requested by "--restart" cmdline switch 2012-06-18 09:27:36 +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
Martin Diehl b2fd3e1180 introduced preprocessor identifiers Marc, Abaqus, and Spectral to enable conditional compilation. This allows deleted dummy functions that are used by one solver only.
Mainly affected modules are IO and mesh. Most of the changes in mesh result from reordering the functions when grouping them depending on their solver.
Further advantage is that FE solver do not need FFTW and kdtree2 anymore. The include files for these two libraries moved to DAMASKROO/lib now as I figured out how to use a include path in the Makefile.
Put all the files I got when testing compilation with abaqus in a folder which to become the abaqus compilation test.
2012-06-15 16:10:21 +00:00
Martin Diehl eb39d332bd ULTRA OPTIMIZATION wasnt working anymore, seems to be related to Bertholds update of ws12 and ws13. had to remove -static from optimization string to be able to link again. Added some support for IKLM, but this didn't solve the problem 2012-06-13 17:19:16 +00:00
Martin Diehl 3ef17c208b set default lapack to lapack instead of acml, added definition of makro for spectral solver for later use in preprocessor 2012-06-13 13:45:58 +00:00
Martin Diehl 052b526bde made compilation of spectral code a command line switch and removed interactive question 2012-06-13 13:42:00 +00:00
Franz Roters 7754416b35 now both getSolverJobName and getSolverWorkingDirectoryName use unit 5! 2012-06-13 09:58:06 +00:00
Franz Roters d3d0cfdab1 common blocks from 2012
updated setup to also create DAMASK_marc2012.f90
2012-06-13 07:30:27 +00:00
Franz Roters a5c27f22ca now uses unit 5 for file name inquiery (unit 6 points to stdout in 2012, might be a bug though) 2012-06-13 07:28:43 +00:00
Franz Roters 5b7f2e122e fixed usage of OpenMP function library 2012-06-12 09:44:05 +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
Martin Diehl fa7f9866df removed phase contrast/preconditioning
added information on itmin in output, impoved output
set exit code to 0 on successful termination (seems to be unix standard)
exit codes:
0: successful termination
1: error (using IO_error)
2: require regrid 
updated the AL solver, still VERY experimental
2012-06-05 16:34:20 +00:00
Martin Diehl 4b6800b89a removed some unwanted open statements when writing out the sizeState 2012-06-05 13:13:33 +00:00
Franz Roters 01522cf6f2 check for valid elasticity model
improved error message in case of unknown elaticity/plasticity model
2012-06-02 14:23:28 +00:00
Nan Jia 3af9165b08 useless debugging inforomation (asking for elastic moduli output) is removed 2012-05-31 09:25:04 +00:00
Christoph Kords 145d971d7e added some debug output 2012-05-30 08:05:36 +00:00
Christoph Kords 5cc73a15a9 small correction in debug output 2012-05-30 07:41:22 +00:00
Martin Diehl 0127828f70 corrected initialization order 2012-05-29 15:08:18 +00:00
Taymor El Achkar e1c5671e1f relaxed error message triggering such that purely elastic calculations and more than one constitutive instance using hex lattice are possible. no twinning possible if no slip system exists 2012-05-24 14:08:48 +00:00
Martin Diehl 1003cfa7f8 added missing space in makefile, made removing of old files in Makefile silent and updated files for tests to current versions 2012-05-23 18:08:34 +00:00
Philip Eisenlohr c607441717 (likely) fixed a bug in the FFT-based geometry reconstruction. For (hopefully) correct math see Appendix B in paper. 2012-05-22 18:35:15 +00:00
Nan Jia bbf790a0a5 Input parameter 'QedgePerSbSystem' which denotes activation energy for each shear band system is added 2012-05-22 16:10:28 +00:00
Franz Roters 9201133913 activated keyword expansion 2012-05-22 12:45:37 +00:00
Martin Diehl cec3357e8f corrected bug causing wrong element construction when using spectral solver. Also improved error messages for invalid resolution, dimension, and homogenization. 2012-05-21 09:06:02 +00:00
Christoph Kords 20f203f799 allow minimum stable dipole height equal to zero; ensure that maximum stable dipole height is not smaller than the minimum value 2012-05-20 13:57:35 +00:00
Christoph Kords 5b02d4e8eb dipole dissociation and formation by change in stress as new mechanism in deltaState; had to add dipole height to state variables, which is however updated by deltaState instead of microstructure; alternatively microstructure had to know the current stress state 2012-05-18 14:35:52 +00:00
Christoph Kords 4da6907648 moved remobilization of dead dislocations from dotState to deltaState 2012-05-18 13:35:44 +00:00
Christoph Kords a5c1624648 check for invalid mesh type once in constitutive_init and not every time constitutive_nonlocal_dotState is called 2012-05-18 13:15:23 +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 351c2c6e65 deltaState depends on the current state, no the state at the beginning of the timestep 2012-05-17 12:18:30 +00:00
Christoph Kords f6d5efeed3 removed arguments Fe and Fp from collectDeltaState function call 2012-05-17 11:04:22 +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
Christoph Kords 84d4652a07 all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState). 2012-05-16 14:43:26 +00:00
Franz Roters df2d520630 stripping abaqus comments is now aware of damask comments 2012-05-15 14:59:26 +00:00
Franz Roters 701cf2e71b replaced duplicated math entry by missing mesh entry 2012-05-15 14:58:47 +00:00
Martin Diehl 0d745adfa0 corrected NaN in prec_single.f90 (was already silent but with typo), polished output and added switch in order to check range on all files except of prec.f90 (in which NaN is defined) 2012-05-11 12:46:17 +00:00
Franz Roters 7ee368f324 check whether there is any constitutive output 2012-05-11 07:38:20 +00:00
Franz Roters f8f9c089d6 check for NaN in math_QuaternionInSST 2012-05-11 07:37:32 +00:00
Franz Roters babf2647fb use non signaling version of NaN 2012-05-11 07:35:52 +00:00
Martin Diehl ce83026965 corrected bug in fft based geometry reconstruction causing mirrored geometries 2012-05-09 14:56:25 +00:00
Taymor El Achkar d11aa15a67 changed indices counting (fortran to python style, easier division to get rid of periodic copies)
corrected f2py pyf file
2012-05-09 10:09:56 +00:00
Krishna Komerla fbd075dc6e corrected automatic determination of resolution 2012-05-08 15:12:43 +00:00
Krishna Komerla c752dd5474 regridding is now working, changed the subroutine into a function
changed order of arrays in nearest neighbor search to make it fortran fast
constitutive.f90 and homogenization.f90 write state size out during initialization
setup/setup_processing.py is using byterecl to be compatible with binary files written out by solver
2012-05-08 14:57:06 +00:00
Taymor El Achkar a9a72cee97 fixed bug in nearest neighbor search, corrected error message for kdtree2.f90 2012-05-08 13:16:59 +00:00