51552b027d
fixed bug: initialize computation mode to zero, used to be unset for lovl=4 and not(timinc < theDelta .and. theInc == inc), which unpredictably caused trouble with marc
2013-04-23 14:42:22 +00:00
5837b25e27
accidently commited constitutive.f90 wasn't working, fixed now
2013-04-23 07:37:14 +00:00
d5e4a354e6
misplaced space corrected
2013-04-22 15:35:15 +00:00
b5b6bf56d4
reintroduced Calculation of forest projections for edge dislocations, got lost after rev.1908
2013-04-22 15:17:24 +00:00
8da7544978
marc doesn't need error prone setup_code any more, moved Marc includes from code/include to lib/includeMarc
...
DAMASK_marcXXXX.f90 now is not a copy any more but contains one definition and an include statement
setup_code of spectral solver is now replaced by the standard 1) configure 2) make 3) make install
2013-04-22 14:41:33 +00:00
a7706090d2
neighborhood relations were screwed up for some elementtypes as of version 2318, might have caused segmentation faults with nonlocal model, now fixed
2013-04-22 13:35:35 +00:00
49ef6d7193
spectral_build_neighborhood now "hidden" in spectral compile part. (should not interfere with Marc/Abaqus compile anymore)
2013-04-22 09:01:58 +00:00
ca99274a5d
fluxes now always periodic for spectral solver, no need to explicitly state this in the geom file by use of "periodic" keyword
2013-04-22 08:32:51 +00:00
60fa91be9a
calculation of cellnode positions now available in python via damask.core
...
first call damask.core.mesh.mesh_init_postprocessing(meshfilename) to initialize all necessary mesh variables
then damask.core.mesh.mesh_build_cellnodes(nodes) calculates the cellnode positions for a given list of node positions
the meshfile that is needed for the init is created automatically by mesh_init in DAMASK
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-21 18:48:59 +00:00
d203e8dece
removed extra comma (might have caused issues?)
2013-04-19 16:05:13 +00:00
e8f238a413
added specialized ipNeighborhood calculation for spectral solver case
2013-04-19 16:04:20 +00:00
12760026a4
last version still did not compile: IO_write_jobFile was not used by mesh_init
2013-04-19 13:02:37 +00:00
55bc741d1d
last version did not compile: fixed small typo in mesh_init
2013-04-19 12:59:45 +00:00
4bc5e6717b
fixed small bug when calculating mesh_Ncells
2013-04-19 12:41:06 +00:00
103c770ee6
fixed bug in cell generation: "FE_Ncellnodes" was not correct for elem type 17 and 136
...
too much memory allocated for mesh_cell
2013-04-18 17:56:14 +00:00
6b46a9c338
introduced output of initial geometry (vtk file) to mesh
2013-04-18 16:40:49 +00:00
8b2d9d8155
further improvements on implementing the plain mode
2013-04-16 17:07:27 +00:00
2ba986be2a
introduced new dummy module "libs" for inclusion of further library sources. contains kdtree2 and Lib_VTK_IO at and is enabled for the spectral solver only at the moment
2013-04-16 12:45:12 +00:00
2a6ac44496
Major restructuring of mesh.f90 due to introduction of "cells". A cell is the volume (area in 2D elements) that surrounds an ip. It is described by up to 8 cell nodes and is used to calculate ip volume, area, and face normals. It can also be used to generate a mesh for e.g. paraview that contains ip output as cell centered data.
...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-15 08:13:20 +00:00
eaccf88e12
fixed bug that was introduced in the element type mapping in version 2305
2013-04-10 13:22:02 +00:00
9b1fd579d4
corrected phase numbering in example material.config
2013-04-10 12:09:09 +00:00
ddcc795461
allowed arbitrary grid (formerly know as resolution), i.e. any integer>1 in each dimension is possible.
...
removed square root correction for divergence (doesn't make sense)
fixed FFT debug of spectral solver
2013-04-10 10:19:16 +00:00
4338cd13bc
- added element type 54 for macro (2d 8node 4ip reduced integration)
...
- automatically raising error for unknown element type inside mapping function
2013-04-10 09:38:40 +00:00
338e338c9a
- fixed bug in volume calculation of 2d elements
...
- subroutine "mesh_build_ipVolumes" uses openmp parallelization
2013-04-10 08:24:53 +00:00
d80d416c32
2d elements now have a volume corresponding to a thickness of 1; used to have zero volume
...
now also 2d elements can be used with nonlocal model
2013-04-09 18:07:30 +00:00
5b508b5ee4
no need to exclude any geometry type from usage of nonlocal model
2013-04-09 13:10:31 +00:00
924d943edc
simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional)
...
commented and cleaned up the marc interface.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 10:08:00 +00:00
7b5a202e8c
enabled alternative (and soon standard) keywords grid (resolution) and size (dimension).
...
removed errors for odd resolution
2013-04-08 14:22:32 +00:00
ad879ebcf9
forgot to delete all occurrences of obsolete "constitutive_nonlocal_ accumulatedShear" in last commit
2013-04-08 08:51:12 +00:00
10e50bf41d
improved accuracy of "accumulated shear": now added as a state variable, which facilitates integration (dotstate=shearrate); needs absolute tolerance value for state integration in material.config
2013-04-04 13:37:14 +00:00
3c25723139
- make sure that formerly blocked density is not able to move to neighbor until it is remobilized
...
- added some explicit indices instead of ":"
2013-04-03 16:29:48 +00:00
6d3f05b387
small changes in nonlocal constitutive model:
...
- make sure neighboring density is always greater equal zero
- no extra Bauschinger in postResults
2013-04-03 16:22:55 +00:00
5b8257a7f9
added debugging possibility for MSC.Marc, rename parameters to CAPITALS
2013-03-31 13:06:49 +00:00
2f76365ac8
fixed wrong string concatenation
2013-03-30 19:43:06 +00:00
965ab38217
corrected error bounds for spectral parameter "maxCutback"
2013-03-28 16:04:47 +00:00
d59e072209
general polishing
2013-03-28 13:50:20 +00:00
7f37f55904
added missing line continuation in numerics and explicit private statements in the other files
2013-03-28 10:37:00 +00:00
d2ae06015a
moved "Random is random" to output of random seed, cleaned up a little
2013-03-28 10:02:11 +00:00
745ca5a6a5
added "public" (and where possible "protected") statements to functions/subroutines/variables/parameters where needed
...
added -fmodule-private to compile options for gfortran, i.e. gfortran won't compile anymore if an entity is not explicitly declared as public but used in another module to ensure stronger encapsulation
2013-03-28 07:40:30 +00:00
b8f8d66f82
total density used in backstress calculation now based on average of neighborhood
...
fixed small bug in state_init: random distribution of density was probably not working correctly, as some variables were not properly initialized
2013-03-27 13:04:01 +00:00
39a70e8a19
fixed bug in coordinate calculation when no average F was given (scaled wrong by ncp_elems**2)
...
other changes: just polishing + some more comments
2013-03-27 12:28:55 +00:00
b6aecdac17
corrected debugging output of node twins
2013-03-27 10:25:37 +00:00
e25af5db04
added function for finding the (non-periodic) nearest neighbor
2013-03-27 10:14:44 +00:00
44eb771583
added crystallite output of current ip coordinates
2013-03-27 10:12:28 +00:00
15b6242585
set more appropriate default values for spectral solver tolerances (abs stress BC was at huge())
2013-03-26 16:36:52 +00:00
33dc8251c0
saved some memory not copying data within one function
2013-03-25 19:02:12 +00:00
97fe4400ff
corrected comment explaining structure of nBlock
2013-03-25 18:36:21 +00:00
e6f6c22a30
polished Abaqus Standard, changed inclusion of OMP functionality from "spectral style" (include) to "marc style" (use) to prevent explicit from crashing
2013-03-25 18:25:15 +00:00
966ad2826b
merged two-stepped subroutine call into one call, added comments and cleaned up.
...
working for small example (with and without openMP), don't know if the results make any sense
2013-03-25 17:45:58 +00:00
19655c2d92
significantRho check was flawed for dipole density
2013-03-25 17:36:04 +00:00
d55ba3e3e6
now using correct exit routine for Abaqus/Explicit
2013-03-25 13:52:22 +00:00
02f092492e
updated copyright header
2013-03-22 17:35:05 +00:00
5b96c1d62a
scaling of deformation gradient fluctuations now possible in x,y,z independently (give array).
2013-03-22 15:09:55 +00:00
be655ae536
Changed AL solver errors to absolute. Stress boundary error now max of relative and absolute error
2013-03-22 14:46:55 +00:00
2a6f66d504
Abaqus explicit now with correct interface, seems to work now
2013-03-21 17:24:58 +00:00
cee90bbdb3
simplified dislocation kinetics: dropped relativistic correction of velocity, since the infinite attack frequency already sets a limit for the velocity; dropped second; also dropped the correction term for backward jumps, since it was probably not meaningful
2013-03-21 12:52:29 +00:00
0b10c52b77
restore Jacoby after cutback, was lost during restructering
2013-03-15 06:27:04 +00:00
dd06d07e1c
set default number of iteration for spectral solver to 250 instead of 40, abaqus doxygen now does the correct inclusion (everything appeared twice before)
2013-03-14 23:10:02 +00:00
1779510dda
remove some unused variables, set default value for divergence calculation to 1.e-5 and made calculation independent of dim and res
2013-03-12 18:48:28 +00:00
3e2ecbe0a3
added time stamp to init of lattice
2013-03-07 22:45:00 +00:00
b9e0326240
added function to calculate the gradient, addGradient.py will follow
2013-03-07 07:58:42 +00:00
1ce6028ad3
simplified CurlFFT and DivergenceFFT functions, and did the last changes to Utilities related to the new structure of the spectral solver
2013-03-06 19:34:30 +00:00
60633ffd98
some doxygen corrections
2013-03-06 14:41:15 +00:00
1f4d7c2ca4
changed reference stiffness to min max avg from volume avg. using terminally ill information in convergence check
2013-03-06 14:31:13 +00:00
7e0ea5dd04
hickups in doxygen documentation fixed
2013-03-05 14:35:26 +00:00
5c1185a5d2
renamed F_lambda to F_tau for clarity. using stress tangent to precondition equations
2013-03-04 09:49:40 +00:00
ada2beb8b8
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 11:48:29 +00:00
33e75972cc
unified indices used for better understanding
2013-03-01 09:36:45 +00:00
72bc1df832
forgot in previous commit
2013-02-28 17:37:26 +00:00
0209a1d5a7
restructured AL solver to resemble polarization scheme and introduced alpha and beta parameters in numerics to determine what kind of solution you want. alpha = beta = 1 results in the AL scheme. alpha = beta = 2 results in the accelerated scheme. safe values are anything in between. changed the solver to use stress calculations on the campatible field which seems to result in faster convergence.
2013-02-28 17:35:02 +00:00
59da88460e
preparation for introduction of plain mode and restructuring of calc modes for CPFEM_general
2013-02-27 20:41:14 +00:00
875abcd06f
replaced some array assignments, which might have caused trouble when using the "timeSyncing" option, by parallelized loops of scalar value assignments
2013-02-27 16:36:12 +00:00
62e55483f7
fixed bug in lattice: line continuation character was missing
2013-02-27 10:49:23 +00:00
3736ab1d97
added shortcuts for two do loops in the "timesyncing" procedure
...
math_pDecomposition enclosed by omp critical statement; somehow, this subroutine is not thread safe, not clear what exactly the problem is
2013-02-27 10:32:37 +00:00
bac186c5b7
removed 2003 standard initialization expression not supported by ifort so far
2013-02-27 07:35:53 +00:00
0be6706483
some more detailed doxygen comments
2013-02-26 19:01:31 +00:00
04c2b22766
removed time information from interface routines, should work now
2013-02-26 12:49:07 +00:00
9ec87e8d1f
moved include statement of IO.f90 before corresponding use statement, still untested
2013-02-26 06:28:28 +00:00
72418f7b60
removed -c from options to compile_spectralSolver and added functionality to remove "-c" or "--compile" from command line arguments befor calling compile_spectralSolver
2013-02-25 18:46:36 +00:00
07ff1f7c93
reporting current time for most module inits
2013-02-25 16:34:59 +00:00
5b16f57727
added int vec pointer for more efficient memory allocation of homog and micro counting (elemsOfHomogMicro) in material.f90
2013-02-25 13:13:52 +00:00
8f2e164fd2
some more documentation, mainly on the exit codes
2013-02-23 00:24:30 +00:00
393c60565d
one line was beyond 132 characters
2013-02-22 05:38:02 +00:00
42b96354db
doxygen comments
2013-02-21 23:08:36 +00:00
1f86e098ba
doxygen documentation, enabled uncommented debug statement to be turned on by switches, minor polishing
2013-02-20 22:06:15 +00:00
a98be4e0cf
IO_error and IO_warning crashed for string of zero length (or only containing whitespaces), fixed now.
...
Also added BEGIN and END statement to numerics in case of non-presen DAMASK_NUM_THREADS, would prevent this error and is better for debugging
2013-02-20 21:56:59 +00:00
dfe34fd162
specified quit as external as it is not part of a module
2013-02-20 14:37:12 +00:00
ecbc15abd0
doxygen related changes: corrected arrows in material.f90 and replaced $REVISION$ by substitution key to be changed after test
2013-02-19 22:12:05 +00:00
d8da2f60d8
added core module function math_periodicNearestNeighborDistances
2013-02-19 14:56:26 +00:00
7461e13c1e
relocated spectral_quit() into driver code.
...
added (dummy) core_quit() for Python damask module.
2013-02-19 14:54:34 +00:00
b57d001e71
added hints on offending positions when parsing material.config
2013-02-15 08:26:38 +00:00
c1cf446774
parsing of material.config file now expects adequate number of (slip/twin) family entries according to given lattice type.
...
used to read up to maxN, but that caused unnecessary IO_warnings...
2013-02-14 22:24:55 +00:00
cd0325baf1
introduced atol_twinfrac for dislotwin
...
set values in material.config example
atol_rho 1.0,
atol_twinFrac 1.0e-7
2013-02-14 09:57:26 +00:00
77c0b32af9
removed possible (but not-working) output of schmid_factor_shearband from dislotwin
2013-02-14 08:53:17 +00:00
d715ec769b
removed debug statements from symlink_Code.py, added autodetection for fortran compiler and disabled compiling on default
2013-02-13 18:45:49 +00:00
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
b78fafc3b5
corrected integer initialized to reals
2013-02-13 15:47:00 +00:00
12b81a2b75
moved second file extension (.pes) into DAMASK_abaqus_*.f
2013-02-13 10:56:50 +00:00
35346b70b5
fixed RGC homogenization in case dt==0
...
removed misplaced parameter statement in IO.f90
2013-02-13 09:36:06 +00:00
0856a7aa32
fixed bugs and simplified warning reporting when reading in int, float, and string values
2013-02-12 19:00:41 +00:00
Christoph Kords