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
Christoph Kords
63451f6edf
added "absavg" (=average of the absolute values) to the predefined mappings
2011-07-15 08:07:33 +00:00
Martin Diehl
09ba92c26e
added and rearranged error messages, polished output and simplified calculation of f depending on fdot or L
...
guessing along former trajectory is now on per default, 'guessreset' and 'dropguessing' switch it off.
2011-07-14 09:37:31 +00:00
Martin Diehl
56340fd487
changed input of loadcase. Now fdot (time derivative) can be used instead of velocity gradient. Velocity gradient needs to have each line fully or not at all defined, as for other loadcases the stress BCs are not known in advance. Also added the possibility to keep guessing along trajectory of former loadcase.
...
changed back to use the compliance of initial linear material behavior.
added counter of non-converged steps
renamed compiler flags in makefile
2011-07-13 16:33:12 +00:00
Martin Diehl
8153cd50b4
added spectral numeric parameters to numerics.config, changed values for spectral numeric parameters in numerics.f90 to more suitable values for new method for determination of divergence in Fourier space
2011-07-11 14:28:56 +00:00
Martin Diehl
a08eaecf4e
polishing (names of variables, spaces, comments etc.)
...
renamed mpie to DAMASK in IO.f90 error message
2011-07-07 15:27:35 +00:00