295d0cd28e
removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC
2013-10-14 10:54:45 +00:00
80cb78c1fd
corrected activation of slipbands; both velocity and resistance have to be non zero
2013-10-11 13:52:00 +00:00
3b5a5ab812
changed error in case of number of found interaction types is less than expected to warning
2013-10-11 09:17:03 +00:00
ed47d25adc
remove dotTemperature (returned always 0 anyway) to make it easier to replace it with a heating rate funtction
...
remove potentially dangerous short circuiting in homogenization files
added doxygen comment for nonlocal (giving the functions a name)
fixed wrong definition of maxnchungs introcuded into dislotwin with last commit, causes trouble with hex
2013-10-09 06:12:16 +00:00
82faf74363
added doxygen documentation and unified variable names and some common parts of the code
2013-10-08 16:27:26 +00:00
3eec0ecdcd
introduced #EOF# as module wide parameter in IO.f90
...
renamed read/write binary file to read/write real/int
removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 14:07:55 +00:00
5f973a21c0
require all hardening coefficients/ non schmid coefficients to be specified in material.config, do not silently ignore any
...
warning message now more meaningful (?) when using less than maximum number of slip families
2013-09-12 14:47:09 +00:00
c93373a4ac
removed long lines
2013-08-08 19:33:21 +00:00
813ed130b0
nonSchmid behavior now ready to use for bcc with phenopowerlaw and nonlocal constitutive model.
...
nonSchmid tensors according to Koester,Ma,Hartmaier,2012.
extended lattice_Sslip with non-Schmid tensors (the full non-symmetric tensors are required for the tangent and cannot be retrieved from the symmetrized Mandel notation)
2013-08-05 09:23:21 +00:00
4a291dc372
implementation of twin nucleation criteria according to Davids PhD thesis
...
applies for fcc crystal structure only
2013-07-01 13:06:01 +00:00
fd8d85896a
error msg in case of unknown keyword was not properly trimmed --> 64kB of spaces...
2013-06-27 16:41:00 +00:00
98528f9a89
added possibility for multi-level inclusion of files in *.config and loadcase files.
...
include subfiles by stating
{path/to/include}
2013-06-26 19:19:00 +00:00
296a5040b9
fixed bug in initialization when reading in material.config.
...
Fortran does not short-circuit, e.g. (if a>0 .and. b(a) ==c) might cause an out-of-bounds error when a=-1 it'll try to access b(-1)!
2013-06-11 20:16:40 +00:00
a21dd816c7
init reporting of constitutive_*.f90 had less/to much spaces, renamed label to LABEL because it is a parameter.
...
removed debug output of geom_fromEuclideanDistance.py
2013-05-28 17:31:55 +00:00
399a0218c7
bug fix on commit 2424 that caused compilation errors
2013-05-22 09:40:23 +00:00
a104fda3b9
Adding accumulated shear due to "slip" and "twin" as two outputs
2013-05-19 20:32:06 +00:00
871df5427d
Correction for shear_rate_twin in postresults: characteristic twin shear was missing
2013-05-19 15:45:34 +00:00
0d95d33646
use disk shape for twin also in initialization
2013-04-30 07:34:05 +00:00
b5b6bf56d4
reintroduced Calculation of forest projections for edge dislocations, got lost after rev.1908
2013-04-22 15:17:24 +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
02f092492e
updated copyright header
2013-03-22 17:35:05 +00:00
07ff1f7c93
reporting current time for most module inits
2013-02-25 16:34:59 +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
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
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
c4b877d4f1
replaced missued family index by system index
2013-02-06 08:45:08 +00:00
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
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
50e874f887
fixed bug in elasticity matrix calculation
2013-01-22 15:48:47 +00:00
652a8366fe
splitted lines > 132, added pure statements where suggested by gfortran 4.7
2013-01-22 13:02:23 +00:00
60fec0e8ec
added code structure for non-schmid mechanics. work in progress…
2013-01-21 23:11:16 +00:00
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
c195f316e1
corrected calculation of twin stiffness matrix
2013-01-21 08:17:43 +00:00
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
1896c07c70
corrected some bugs in the last commits
2012-11-14 14:33:41 +00:00
e4ce63d689
switched matrix indices in SlipSlip, SlipTwin, TwinSlip, and TwinTwin such that first index refers to first name and second to second... (used to be transposed for memory efficiency)
2012-11-14 10:22:34 +00:00
7c6fc121fd
condensed error reporting for constitutive_XYZ_init
...
removed erroneous check for structure>3
2012-07-17 17:36:24 +00:00
f5eb8fcb25
fixed bug caused by the inconsistent sequence of slip and twin systems
2012-06-19 16:06:25 +00:00
3af9165b08
useless debugging inforomation (asking for elastic moduli output) is removed
2012-05-31 09:25:04 +00:00
e1c5671e1f
relaxed error message triggering such that purely elastic calculations and more than one constitutive instance using hex lattice are possible. no twinning possible if no slip system exists
2012-05-24 14:08:48 +00:00
bbf790a0a5
Input parameter 'QedgePerSbSystem' which denotes activation energy for each shear band system is added
2012-05-22 16:10:28 +00:00
84d4652a07
all constitutive modules now contain a new function "deltaState", which in the future allows to have an instantaneous and incremental change of the state (additional to the rate based evolution with dotState).
2012-05-16 14:43:26 +00:00
70a44fbbdd
constitutive models need to skip new keyword elasticity when reading their arameters
2012-03-15 09:22:24 +00:00
433d8a338c
missed several occurences of constitution
2012-03-12 14:43:19 +00:00
c9ed989dad
replaced keyword constitution with plasticity in preparation of upcomming generalized elasticity
2012-03-12 14:09:37 +00:00
79663a7f76
polishing: adding pInt, removing unused use-statements etc
2012-02-21 16:00:00 +00:00
db06c797cd
store number of outputs for each instance of the constitution in constitutive_[j2/titanmod/...]_Nouput and use this value to loop over all valid outputs when determining the size of the outputs array
2012-02-14 15:19:59 +00:00
d9522bf588
minor changes, substituting implicit type conversions with explicit ones
2012-02-14 12:17:47 +00:00
30c55c3142
added check for unknown keywords again (in agreement with Philip ;-)
...
indeed 'constitution' and '/nonlocal/' are the only keywords which have to be ignored, since they are parsed elsewhere
2012-02-14 09:22:37 +00:00
9b73cb3c65
removed parameter check complaining about unknown keyword as default case. there are just too many keywords that do not trigger a variable assignment to cope with them (e.g. /nonlocal/, constitution, etc.)
...
I suggest to kill similar logic from other files if present somewhere...
2012-02-13 23:30:59 +00:00
c786336af3
reordered (and partly redistributed) error message identifiers, deleted those which are not in use anymore.
...
all constitutive as well numerics now raises an error, if an unknown keyword is found in the respective config file
2012-02-13 17:41:27 +00:00
156ec4582a
polishing, adding _pInt etc. where applicable
...
post_results now handels zero increment different (like FEM, it is always there even if it is not counted)
2012-02-13 14:18:07 +00:00
37ac7bf1b4
fixing inconsistencies in variable assignments
2012-02-10 11:56:05 +00:00
c5128e94eb
some f2008 standard related correction i didn't check in last time
2012-02-01 20:20:05 +00:00
800e291240
made code standard conform to Fortran 2008 (ignoring warning concerning comments beyond character 132). Basically, changing "x" format specifier to "Nx" ("1x") plus removing $ format specifier
...
added compiler switches for gfortran and ifort to check for standard conformity
old gnu compilers <4.4 are not longer supported because they don't provide the c binding for fftw
2012-01-31 19:18:55 +00:00
80583fefb5
added compilation_info.f90, to store compilation date,time and compiler info in each file
...
added #include statement at each init() routine
2012-01-31 14:54:49 +00:00
1330576a01
added new math function "math_deviatoric33" which returns the deviatoric part of a 3x3 tensor
...
renamed some math functions, so that we have a universal naming scheme: for matrix multiplications use an "x" (e.g. math_mul33x3); don't use the "x" to describe the shape of the tensor that the function is applied to (e.g. math_invert33 instead of math_invert3x3)
2012-01-26 13:50:00 +00:00
a7d9d711f6
edited lines exceeding 132 chars (before any comment started)
2012-01-11 16:56:35 +00:00
5ebeb96e85
made DAMASK to work with gfortran:
...
-removed to long lines
-restructured f2py modules and merged make_DAMASK2Python into setup processing
-setup_code.py now sets library path in makefile and asks for compile switches for spectral code
-substituted \ in format strings with $
restructured DAMASK_spectral:
-more logical output and structure of code
-better input for spectral debug parameters
2011-12-06 16:58:17 +00:00
58917b5abe
corrected error messages
...
you can now switch of shear banding by seting sb_velocity to zero
stacking fault energy parameters are now read from material.config
2011-09-26 09:55:08 +00:00
2e263dee90
merged in the shear band code, it compiles but is not otherwise tested yet
2011-09-16 15:55:18 +00:00
1dd2591982
corrected wrong check of twin systems
2011-09-13 08:18:43 +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
b4678112cd
sorry, last version did not compile. now with correct encoding.
2011-04-13 11:51:46 +00:00
dd2057bac5
changed header to unicode utf-8 encoding
2011-04-13 11:24:36 +00:00
c1b8391110
changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
2011-04-07 07:20:28 +00:00
fcdb805225
added copyright text to all f90 (free) format files
2011-04-04 14:09:54 +00:00
11138e3ee2
debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output)
...
0 : only version infos and all from "hypela2"/"umat"
1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info
2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines
3 : basic outputs from "homogenization.f90"
4 : extensive outputs from "homogenization.f90"
5 : basic outputs from "crystallite.f90"
6 : extensive outputs from "crystallite.f90"
7 : basic outputs from the constitutive files
8 : extensive outputs from the constitutive files
If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 10:31:17 +00:00
366d52bd71
* former "relevantRho" and "relevantResistance" is renamed to "atol_rho" and "atol_resistance" and now used as an absolute tolerance for the state residuum. before it was used rather as significant state, so whenever the state dropped below that value it was considered converged. (In dislotwin and titanmod constitutive law there is only one value atol_rho which is used for all the states, though the state array consists not only of densities. We need further parameters here!)
...
* somewhat simplified convergence check for adaptive euler and runge-kutta
2010-10-26 13:16:37 +00:00
61f8a5fcbe
corrected UpperCase "SolidSolutionStrength" in config file parsing.
2010-10-15 12:59:00 +00:00
f78b07448f
some modifications in line with Davids analytical model
...
input some additional variables via material.config
2010-08-17 14:23:55 +00:00
298cecbfec
corrected syntax errors (long lines, line continuation by \) and logical mistake in mpie_cpfem_marc / abq_std which went unnoticed in ifort. Using SunStudio f90 surfaced those...
2010-08-03 23:47:00 +00:00
3aa2dd5fef
removed special characters from twip steel section in material.config
...
changed crirtical twin stress in constitutive_dislotwin.f90
2010-03-24 08:17:27 +00:00
67f87486b1
in constitutive_dislotwin.f90:
...
in constitutive_nonlocal.f90:
Derivatives of shear rates w.r.t. resolved shear stress HAVE to be positive.
Computation of dgdot_dtauslip is now correct.
2009-11-06 17:46:12 +00:00
00450b2c16
Added constraint on the dipole distance
...
if dipole distance > dislocation mean free distance then
dipole distance = dislocation mean free distance
2009-10-13 18:14:44 +00:00
3bf3d7ddbd
New version of constitutive_dislotwin after first-round checks
2009-10-13 15:05:21 +00:00
7d05845d5d
# now $ID aware
...
# LF as line separation and some tidy up
2009-10-12 16:26:16 +00:00
3f41b271fb
The new dislocation model is now implemented in constitutive_dislotwin.f90
...
(seems to run fine, but still under investigation)
2009-10-09 20:28:23 +00:00
cc6aa1aeba
Updated constitutive_dislotwin
...
Deleted constitutive_dislobased
2009-10-06 15:19:19 +00:00
Martin Diehl