60ef9d5ab0
updated installation instructions
2011-09-13 07:58:54 +00:00
7be2701989
* new output variable "internalstress", which gives the internal stress tensor
...
* use "math_invert3x3" instead of "math_inv3x3" for inversion of Fe
* for dislocation stress calculation: first regular case, then special case of dead dislocations in central ip
* "dv_dtau" now given for each dislocation type, so is a (ns,4) array
* deleted unused variables in "_LpAndItsTangent"
* corrected contribution of deads in "_LpAndItsTangent"
* the NaN variables defined in math did not give a proper NaN value, so use 0.0/0.0 again
* neighbors with nonlocal constitution but local properties (i.e. /nonlocal/ flag not set) are also considered for incoming fluxes
2011-09-07 11:30:28 +00:00
b412239d9f
reflects new output format of DAMASK_spectral to find
...
*** Loadcase...
lines.
2011-09-05 13:39:44 +00:00
48ac3956ce
added '\n' escaping of ID string in VTK header
2011-09-05 13:38:10 +00:00
f458f026db
forgot to include string module, too
2011-09-05 08:54:29 +00:00
0ac5109c94
forgot to include string module
2011-09-05 08:50:36 +00:00
e8628996ff
added IO error in case mixed BC that allow rotation are given
2011-09-02 13:50:05 +00:00
320ec4d0b4
added some information regarding spectral method
...
added the information about the new glide system for hex
2011-09-02 13:49:02 +00:00
fe9754a41c
previous commit was wrong file
2011-09-02 10:46:37 +00:00
43e2684cfb
added new glide system for hex lattice
2011-09-02 10:43:49 +00:00
0469d37fc3
cleaned and restructured output on screen, changed L_inf norm to L_2 for calculation of error in Fourier space.
...
removed unnecessary function from math.f90
added some documents/manuals regarding spectral method
2011-08-31 14:37:01 +00:00
f527bdcc78
changed standard values according to numerics.f90
2011-08-30 11:21:39 +00:00
8f90f1f1f9
corrected description of algorithm to invert matrix
2011-08-30 07:29:13 +00:00
4fb1cb8f87
merged calcmodes, i.e. equilibrium and fulfilling of stress BC is done in one step,
...
made convergence independent of size and resolution,
polishing output in DAMASK_spectral.f90
added function to compute eigenvalues without eigenvectors and function to convert a 3x3 logical to a 9 vector in math.f90
removed obsolete variable in numerics.f90
2011-08-26 14:06:37 +00:00
380a536b45
homogenization_*, constitutive_*, DAMASK_spectral_interface: did some polishing regarding init output
...
makefile: corrected spelling mistake
2011-08-26 13:57:29 +00:00
4eb8e25b6d
no fitting recent file renames
2011-08-26 12:25:11 +00:00
d4e7d6e02f
added $Id$ and some further comments
2011-08-26 07:29:55 +00:00
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
eb0ce3992b
debugging of addDebugInformation (now working, but not tested in detail)
...
Added Id to postprocessingMath
2011-08-25 18:18:38 +00:00
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
282d2d7b97
corrected wrong Id tag
2011-08-25 16:47:12 +00:00
bbee6c4658
forget to check in addDebugInformation
2011-08-25 16:44:36 +00:00
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
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
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
39b97f7b67
added new script "spectral_scaleGeom"
2011-08-19 13:32:27 +00:00
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
2dd8a353bb
now also works for single stepped loadcases, too
2011-08-19 10:03:42 +00:00
483234047c
renamed to be consistent with other file names
2011-08-19 07:53:35 +00:00
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
589666d19e
added python script to simply (without antialising) scale a given *.geom file
2011-08-19 07:29:50 +00:00
def056f30f
added $Id$ escaping
2011-08-18 08:00:19 +00:00
10493ea8b0
added $Id$ login name escaping
2011-08-18 07:53:07 +00:00
164da22730
added some warnings if Mentat or PIL cabability missing
2011-08-18 07:52:07 +00:00
ff3bac2496
does not double size of 2D data for spectral method any more
2011-08-17 16:06:50 +00:00
3298e6821b
added forgotten "averageDown" script to list of shortcuts
2011-08-16 18:20:52 +00:00
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
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
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
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
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
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
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
4f1cb3a8af
exchanged forall with nested do loops... Bye, bye smatrix warning, will miss you!!
2011-08-03 17:58:16 +00:00
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
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
64435b8a97
loadcase takes 'temperature' (or 'temp') as input. Previously it was hard-wired.
2011-08-02 13:58:28 +00:00
422d6d9c6c
slight correction of debugging output
2011-08-02 12:36:08 +00:00
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
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
Franz Roters