Commit Graph

980 Commits

Author SHA1 Message Date
Philip Eisenlohr f3c2547fba name change to better reflect what happens: comparison between local deformation field and a (necessarily compatbile) geometry reconstructed from the whole deformation field. 2011-08-26 07:17:35 +00:00
Martin Diehl eb0ce3992b debugging of addDebugInformation (now working, but not tested in detail)
Added Id to postprocessingMath
2011-08-25 18:18:38 +00:00
Philip Eisenlohr 7aba05ed9f reworked to match former script layouts and logics. (not yet tested, use at your own risk) 2011-08-25 17:55:36 +00:00
Martin Diehl 282d2d7b97 corrected wrong Id tag 2011-08-25 16:47:12 +00:00
Martin Diehl bbee6c4658 forget to check in addDebugInformation 2011-08-25 16:44:36 +00:00
Martin Diehl a5f176cf18 put divergence calculation using FDM to postprocessingMath, added divergence calculation using differentation in Fourier space to postprocessingMath
added functions to calculate the shape and volume mismatch of compatible and non-compatible geometry reconstruction to postprocessingMath. They can be accessed using addDebugInformation
2011-08-25 15:31:05 +00:00
Martin Diehl cc3bc4a216 made handling of unknown dimension (i.e. resolution ==1) the same as for preprocessing (last commit) 2011-08-25 15:28:33 +00:00
Martin Diehl 272f7cc655 added auto scaling of dimension to voronoi_tessellation.f90. If dimension(i) == 0, it will get the smallest step of the remaining dimensions in that direction. Same is now done by patchFromReconstructedBoundaries
setup_processing.py now gives executables the extension *.exe and removes the old executables before compiling
2011-08-25 12:16:06 +00:00
Martin Diehl 39b97f7b67 added new script "spectral_scaleGeom" 2011-08-19 13:32:27 +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
Philip Eisenlohr 2dd8a353bb now also works for single stepped loadcases, too 2011-08-19 10:03:42 +00:00
Martin Diehl 483234047c renamed to be consistent with other file names 2011-08-19 07:53:35 +00:00
Martin Diehl 878ea443d6 added options to control written output (mesh(default), box, points) and to control reconstruction. --undeformed does not calculate average deformation of ve, --scaling controls the scaling of the fluctuation. Default is 1.0, with 0.0 no fluctuations are shown 2011-08-19 07:33:22 +00:00
Martin Diehl 589666d19e added python script to simply (without antialising) scale a given *.geom file 2011-08-19 07:29:50 +00:00
Philip Eisenlohr def056f30f added $Id$ escaping 2011-08-18 08:00:19 +00:00
Philip Eisenlohr 10493ea8b0 added $Id$ login name escaping 2011-08-18 07:53:07 +00:00
Philip Eisenlohr 164da22730 added some warnings if Mentat or PIL cabability missing 2011-08-18 07:52:07 +00:00
Martin Diehl ff3bac2496 does not double size of 2D data for spectral method any more 2011-08-17 16:06:50 +00:00
Nan Jia 3298e6821b added forgotten "averageDown" script to list of shortcuts 2011-08-16 18:20:52 +00:00
Philip Eisenlohr bedce496a8 deleted unused option "dimension" since this is figured out by parsing the respective spectral or geom input files anyway. 2011-08-15 08:34:55 +00:00
Philip Eisenlohr e851c934ee interrogate output from DAMASK_spectral for "error divergence:" evolution.
results in an ASCII datafile containing the increment counter with corresponding divergence value. the evolution within an increment i is stretched between i and i+1 (and likely jumps then for the next increment i+1 at the value i+1)
2011-08-14 10:44:15 +00:00
Martin Diehl f99bf63397 removed functions added for debugging of divergence calculation to math.f90
corrected calculation of stress BC condition. Depending on given BC, the stiffness matrix is reduced and than inversed. Then it is filled with zeros and used for the calculation of the correct change of deformation gradient. All calculation is done using dP/dF
2011-08-10 17:45:37 +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
Martin Diehl 34de2e301b DEBUGGING VERSION for calculation of Divergence (various methods implemented). DO NOT USE. Wait for next update coming in a few minutes 2011-08-10 16:02:13 +00:00
Martin Diehl bb1ca2330b added option to enable use of AMD math core library, added comments on the various parameters 2011-08-10 16:00:02 +00:00
Onur Guevenc 31b60b8814 changed debug intensity so that the artificially high stress and stiffness values in collection cycles do not get printed out anymore... 2011-08-03 17:59:31 +00:00
Onur Guevenc 4f1cb3a8af exchanged forall with nested do loops... Bye, bye smatrix warning, will miss you!! 2011-08-03 17:58:16 +00:00
Onur Guevenc 907965272b added option to build DAMASK_spectral without openMP (make OPENMP=OFF)
fixed bug in bc_temperature assignment that was hitting memory.
Temperature is taken from the first loadcase and evolves from there in an adiabatic fashion for the moment. I.e. T-specifications from later loadcases are ignored...
2011-08-03 17:57:28 +00:00
Martin Diehl b5abd6f203 adding new option OPTIMIZATION with values OFF, DEFENSIVE and AGGRESSIVE to makefile. Polishing of makefile and DAMASK_spectral.f90 2011-08-03 07:54:32 +00:00
Onur Guevenc 64435b8a97 loadcase takes 'temperature' (or 'temp') as input. Previously it was hard-wired. 2011-08-02 13:58:28 +00:00
Christoph Kords 422d6d9c6c slight correction of debugging output 2011-08-02 12:36:08 +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
Christoph Kords 6f859e99de * internal stress calculation now considers dead dislocations always at the interface, not at the center of the IP volume; used to merge them together with "normal" dislocations for stress calculation.
* dislocation flux is blocked if we encounter a sign change in the resolved shear stress from the central ip to the neighbor
* do not set density to zero if below certain threshold; this creates an artificial sink term
2011-08-02 11:17:45 +00:00
Christoph Kords 0373fa64e4 have to initialize constitutive_state during constitutive_init, since the first call to constitutive_microstructure (where we already need the state) is done before the cutback scheme starts. 2011-08-02 11:10:41 +00:00
Christoph Kords 9116a1edcb tried to find and execute mpie_cpfem_marc.f90. changed to new name DAMASK_... 2011-08-02 11:02:28 +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
Martin Diehl 3adb7ab382 corrected makefile, now working again without giving standard values explicitly.
did some small modification in order to make it possible to compile with gfortran. Changed NaN=0.0/0.0 to bitwise representation (3 different ways)
2011-08-01 10:11:32 +00:00
Philip Eisenlohr 06a4ac2565 fixed bug (segmentation fault) occurring for non-renumbered meshes: ipNeighborhood building did use FE IDs of twin nodes to address memory -- dangerous..! 2011-07-31 15:42:59 +00:00
Christoph Kords fb121b1435 * boundary condition masking changed
* damper initialized with one 
* inversion of Mandelized stiffness tensor does not work, have to use plain tensor
* new functions in math that allow for conversion between Mandel and Plain tensors
2011-07-29 15:57:39 +00:00
Franz Roters aa714a3d84 some textconstants were wider than specified in the format string, this might have caused Abaqus to crash due to an illegal memory reference (signal 11) 2011-07-29 14:54:26 +00:00
Franz Roters 92d7dcb3f4 added modfiles for Marc/Mentat2010.2
Marc/Mentat2010.2 is installed on msuws1 only right now, please test!
simply start metat2010.2 instead of mentat2010
2011-07-28 07:27:39 +00:00
Martin Diehl 72d20875de added some switches and variables to the makefile to make it more flexible
DAMASK_spectral.f90 is a "debug version" with a number of different criteria to determine divergence. will be removed later on.
2011-07-25 16:30:21 +00:00
Philip Eisenlohr a58158fc7c bug fix related to hash key error when requesting --increments range 2011-07-25 12:53:54 +00:00
Onur Guevenc 8cb8637d41 fixed running with default --range & using --separation options. 2011-07-22 11:59:03 +00:00
Onur Guevenc 105a09f3df only consider increments actually present in (spectral) result
added switch to change from range of positions to range of increments

polished help output
2011-07-21 15:45:41 +00:00
Philip Eisenlohr 2e9c605571 somehow forgot to move the numerics.config into the config subdir... 2011-07-21 08:40:45 +00:00
Franz Roters e92e5cae53 restarting now works with Abaqus (standard, cannot test explicit due to lack of license)
you have to specify the job you are restarting from in the job description (cae), if you prepare your input file by hand this is the first line after *Heading
example: if the first job was using Oldjob.inp the first entry in the job description needs to be Oldjob (without the .inp)
as for Marc restart works only from last converged increment, i.e. ther restart writing should be specified like this:
*retsart, write, frequency=1, overlay
Overlay is not essential but saves a lot of disk space and as stated before you can only restart from the last converged increment anyway
2011-07-18 09:15:20 +00:00
Philip Eisenlohr 8c65cb9ea6 small output polishing 2011-07-15 12:25:38 +00:00
Christoph Kords cb96559d99 mappings should rather be called "sumabs" and "avgabs" than "abssum" and "absavg" 2011-07-15 11:27:47 +00:00