Martin Diehl
1496c62258
updated doxygen config file for FEM solvers, replaced some external statements, removed debugging options from abaqus_v6.env causing Abaqus to crash (didn't detect the bug so far), added new homogenization test to automator, added more debug options + comments to Makefile
2013-02-13 17:54:56 +00:00
Martin Diehl
b78fafc3b5
corrected integer initialized to reals
2013-02-13 15:47:00 +00:00
Martin Diehl
12b81a2b75
moved second file extension (.pes) into DAMASK_abaqus_*.f
2013-02-13 10:56:50 +00:00
Franz Roters
35346b70b5
fixed RGC homogenization in case dt==0
...
removed misplaced parameter statement in IO.f90
2013-02-13 09:36:06 +00:00
Martin Diehl
0856a7aa32
fixed bugs and simplified warning reporting when reading in int, float, and string values
2013-02-12 19:00:41 +00:00
Martin Diehl
59a265ec7f
forgot intent(in) statements needed for pure routines
2013-02-11 10:56:10 +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
c7c81a5ab5
removed unused variables and declared external functions as external
2013-02-11 09:44:17 +00:00
Martin Diehl
a89efaa4a6
added missing bracket and substituted one more forall by do loop
2013-02-09 08:23:47 +00:00
Martin Diehl
f928d0ca0f
some changes on the IKML libs, now using explicit linkin
2013-02-09 08:17:12 +00:00
Martin Diehl
e9f351c69b
forgot to commit
2013-02-08 15:56:58 +00:00
Martin Diehl
147cc8efca
removed debug statement
2013-02-08 15:56:24 +00:00
Martin Diehl
e644c6dbc5
improved reading in of values, now only warnings in case of problematic entries in material.config
...
divergence calculation sqrt scaling optionally introduced for basic scheme spectral solver
2013-02-08 15:55:53 +00:00
Martin Diehl
59e59c90c5
added Option DEBUG to Makefile to turn run time debugging on, compile tests use this option
2013-02-08 15:43:15 +00:00
Philip Eisenlohr
32a416e298
fixed memory leak associated to accumulated shear storage.
2013-02-08 13:33:25 +00:00
Philip Eisenlohr
c7740b41a5
added possible accumulated shear output request to example list of phenopowerlaw
2013-02-07 11:09:14 +00:00
Franz Roters
eacf6f1196
now also fixed the abaqus part
2013-02-07 10:45:10 +00:00
Franz Roters
a036f1eb3a
as some IO functions are no longer pure I changed some forall statements to ordinary do loops
2013-02-07 08:52:47 +00:00
Philip Eisenlohr
722d5574fb
added output capabilities for accumulated shears on slip and twin systems.
2013-02-06 18:09:11 +00:00
Mahesh Balasubramaniam
362cbf3359
Made corresponding changes in the constitutive files and DAMASK_spectral_driver and interface files with respect to the changes made to IO.f90
2013-02-06 16:45:34 +00:00
Mahesh Balasubramaniam
b591dd5f33
made changes to string, int, float interpretations and added warnings.
2013-02-06 16:41:09 +00:00
Franz Roters
c4b877d4f1
replaced missued family index by system index
2013-02-06 08:45:08 +00:00
Mahesh Balasubramaniam
037081831b
last commit did't compile, hope it's working as it should now
2013-02-05 20:07:25 +00:00
Martin Diehl
4997328de1
added use statements
2013-02-05 15:03:36 +00:00
Philip Eisenlohr
865842b0ed
LpandTangent assumed to be called with deviatoric stress input.
...
Now calculates deviator internally...
(Might have led to non volume preserving deformation.)
2013-02-05 13:47:00 +00:00
Philip Eisenlohr
79e7deca55
introduced dedicated calculation of interface normals in IP neighborhood. (2D elements did not work so far..!)
...
fixed definitions of some (2D) interface planes to have all normals pointing outward.
removed "int, private :: i" defined for spectral compilation.
2013-02-05 13:27:37 +00:00
Martin Diehl
cc9eb685fe
fixed wrong setting of coordinates in spectral solvers and removed old python scripts for setup of code and processing
2013-02-05 12:31:44 +00:00
Franz Roters
45344cff6b
added check if nGrains==1 when non-local plasticity is used
2013-02-04 14:34:01 +00:00
Franz Roters
01a7850517
in Abaqus there is a second input file ending (.pes), if this file exists it is to be preferred over the *.inp file (thanks Steffen)
2013-02-04 08:29:58 +00:00
Martin Diehl
7f5e9f70b2
added missing flags
2013-02-03 09:22:25 +00:00
Martin Diehl
6cd6172c0c
fixed bug in dotState causing strange hardening for certain parameters
2013-02-01 15:44:50 +00:00
Martin Diehl
b0263a8aab
added possibility to parse options for setup_code.sh and setup_processing.sh + other improvements in setup functionality
2013-02-01 12:21:56 +00:00
Martin Diehl
c2495d0d4a
changed misplaced private to public statement and renamed qsort to math_qsort
2013-02-01 07:26:21 +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
Mahesh Balasubramaniam
04d48c79b9
symlink_Code.py split seperately from the previous setup_code.py
2013-01-29 15:31:01 +00:00
Mahesh Balasubramaniam
bbd7e72fe8
setup_Makefile split seperately from the previous setup_code.py
2013-01-29 15:30:28 +00:00
Mahesh Balasubramaniam
6dcea3f419
a simple shell script to run setup_Makefile.py,compile_SpectralSolvers.py and symlink_Code.py in the said order
2013-01-29 15:29:43 +00:00
Mahesh Balasubramaniam
c816d5ca74
compiling DAMASK_Spectral split seperately from setup_code.py
2013-01-29 15:28:21 +00:00
Franz Roters
5228af627c
corrected calculation of gmod from c_66 (Voigt approximation for a random polycrystal)
...
approximate twin volume by a disk rather than an elipsoid
2013-01-29 14:59:26 +00:00
Martin Diehl
08a2aa79f7
doxygen comments for homogenization.f90, unified naming ip->i, el->e
2013-01-29 10:28:01 +00:00
Martin Diehl
1594a4bdf8
doxygen comments for isostrain, unified naming ip->i, el->e
2013-01-28 16:36:26 +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
Claudio Zambaldi
ae9af9143d
in python never put import of standard library modules on the same import
...
line as own modules (here:damask)
- und die folgenden Zeilen werden ignoriert --
M setup_code.py
2013-01-24 14:45:23 +00:00
Martin Diehl
db12e6e0c8
removed forgotten debug write statements (thanks Franz!)
2013-01-24 13:20:57 +00:00
Pratheek Shanthraj
7a84209a14
bug fix in microstructure
...
-This line, and those below, will be ignored--
M constitutive_dislotwin.f90
2013-01-24 13:12:51 +00:00
Martin Diehl
6ce78cf806
gfortran now detects "call flush" statements, removed the last of them in mesh.f90
2013-01-23 21:50:01 +00:00
Martin Diehl
d1985e0bef
checked and improved on Franz' one call to CPFEM only improvement
2013-01-23 19:56:45 +00:00
Martin Diehl
96577b18fb
introduced error code for run_test.py
2013-01-23 18:33:46 +00:00
Pratheek Shanthraj
50e874f887
fixed bug in elasticity matrix calculation
2013-01-22 15:48:47 +00:00
Martin Diehl
652a8366fe
splitted lines > 132, added pure statements where suggested by gfortran 4.7
2013-01-22 13:02:23 +00:00
Martin Diehl
df089b94bd
line break in lines > 132
2013-01-22 11:06:39 +00:00
Martin Diehl
2b319f02e6
small corrections on changes related to non-schmid systems
2013-01-22 10:04:15 +00:00
Pratheek Shanthraj
f3bd920c23
added non-schmid structure
2013-01-21 23:50:28 +00:00
Pratheek Shanthraj
60fec0e8ec
added code structure for non-schmid mechanics. work in progress…
2013-01-21 23:11:16 +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
Franz Roters
c195f316e1
corrected calculation of twin stiffness matrix
2013-01-21 08:17:43 +00:00
Martin Diehl
9ee8108b6b
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 11:30:52 +00:00
Martin Diehl
c6a79d2b3d
divergence_correction for basic solver variants has now 3 possibilities:
...
0: uncorrected, slope per sidelength (physical dimension) e = res/dim
1: corrected by sidelength, slope per unitlength e = res/1
2: corrected such that distance between FPs e = 1
alway regarding the medium length of x,y,z direction
2013-01-16 10:40:53 +00:00
Martin Diehl
20bc97b7eb
added doxygen comments and unified naming scheme for ip, element and grain
2013-01-16 10:14:57 +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
Martin Diehl
645b2605e1
corrected reporting to statistics file (*.sta)
2013-01-11 10:40:16 +00:00
Martin Diehl
4c145376a0
forgot some integers
2013-01-10 18:50:14 +00:00
Martin Diehl
6e3e06bed7
corrected PETSc macro expansion
2013-01-10 15:36:55 +00:00
Martin Diehl
563b1f5e4b
added some warning and explicit size of arrays
2013-01-10 13:33:43 +00:00
Martin Diehl
0d5e91ac87
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers
...
improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-09 22:19:32 +00:00
Martin Diehl
55b88e47b7
fixed wrong temperature when using spectral solver
2013-01-09 18:08:08 +00:00
Martin Diehl
00246ade4e
missing range in array fixed
2013-01-09 14:43:27 +00:00
Martin Diehl
093cf92338
remove omp statements during initialization, "line" for reading in from material.config now has a default (empty) value to prevent conditional jump depend on this value
2013-01-08 22:11:59 +00:00
Martin Diehl
4eaa97b33c
minor changes on buggy latest commit (wrong use statement)
2013-01-08 21:54:25 +00:00
Martin Diehl
19d86ca06a
removed unnecsessary omp statements, removed one transpose in constitutitve
2013-01-08 11:09:20 +00:00
Martin Diehl
ce7a0571fd
fixed bug in forwarding fields for AL solver
2013-01-08 10:12:03 +00:00
Martin Diehl
d90efef535
fixed bug in cutback scheme
2013-01-07 21:42:00 +00:00
Martin Diehl
738f363263
corrected small issues with PETSc debugging
2013-01-03 16:17:23 +00:00
Martin Diehl
7a43d1b6ad
added funtionality to specify working directory to spectral solver
2013-01-02 17:02:12 +00:00
Christoph Kords
766202bae9
corrected debug output of aged state: first reports for selected debug_el,debug_ip
2012-12-28 12:15:59 +00:00
Christoph Kords
f14fd45ce1
new material parameter "egdejog": fraction of annihilated screw dipoles that forms edge jogs
2012-12-23 12:56:15 +00:00
Martin Diehl
64d167fa90
loadcase rotation now working for AL solver
2012-12-17 10:18:39 +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
Martin Diehl
1baf8dea5d
load case rotation no working for Basic PETSc solver
2012-12-15 23:52:06 +00:00
Martin Diehl
5c0c7121e7
made PETSc silent by removing -snes_view from default options.
...
introduced PETSc option for debugging that introduces some debugging options into the petsc options and move PETSc initialization from numerics to DAMASK_spectral_utilities
2012-12-15 18:07:49 +00:00
Martin Diehl
566b680319
made new spectral solver the default (no need for SOLVER=NEW anymore) and added additional warnings switche to be used with gfortran >= 4.6
...
added comments
2012-12-15 16:21:10 +00:00
Martin Diehl
b1f28ba5e7
made restart working for all solvers and added corresponding test
...
corrected error handling for FFT based function in case of odd numbers etc.
2012-12-14 17:30:22 +00:00
Martin Diehl
7b87987751
introduced option for regridding to numerics.f90, working as follows:
...
first cut back is tried if material point model oder BVP solver does not converge.
If no regridding is enabled after max cut back, in case of non-converged material point the simulation stops and in case of non-converged BVP solver it continues.
set regridMode to 2 to enable regridding if BVP solver OR materialPoint model do not converge,
set regridMode to 1 to enable regridding if materialPoint model do not converge, non-converged BVP solver will be ignored as in the standard case.
For regridding, the load case need to have a restart freq set.
enabled restarting for Basic PETSc variant
2012-12-14 15:18:04 +00:00
Christoph Kords
64d9c70dc7
corrected nonconforming use of openmp parallelization
2012-12-14 14:30:08 +00:00
Martin Diehl
d559653722
small improvements on tests
2012-12-13 20:20:04 +00:00
Martin Diehl
2b4eb2a257
changed shell variable containing name of Workstation from $HOST to $HOSTNAME as this is also working on Ubuntu
2012-12-13 13:53:54 +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
03a88d4798
fixed timesyncing for special case of ip that synchronizes its time step, but does not have a neighbor
2012-12-11 13:38:36 +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
37028455e9
always use the current state for leaving flux calculation when no neighbor is present
2012-12-09 12:24:32 +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
149c747688
use non-corrected densities for rate check in dotState
2012-12-06 17:14:35 +00:00
Christoph Kords
1f5a30d6f0
corrected typo in random multiplication term
2012-12-06 14:04:18 +00:00
Christoph Kords
20f7053ad3
correction of interaction coefficients due to line tension effect only applies to coefficients 3, 4 and 5
2012-12-04 18:23:46 +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