Martin Diehl
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
Martin Diehl
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
Pratheek Shanthraj
72bc1df832
forgot in previous commit
2013-02-28 17:37:26 +00:00
Martin Diehl
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
Martin Diehl
77c0b32af9
removed possible (but not-working) output of schmid_factor_shearband from dislotwin
2013-02-14 08:53:17 +00:00
Philip Eisenlohr
c7740b41a5
added possible accumulated shear output request to example list of phenopowerlaw
2013-02-07 11:09:14 +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
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
738f363263
corrected small issues with PETSc debugging
2013-01-03 16:17:23 +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
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
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
Christoph Kords
9eb89b0050
made random multiplication switchable via material.config
2012-12-03 12:59:38 +00:00
Christoph Kords
62ddbf26df
slightly changed the multiplication term in the nonlocal model for the starvation case
2012-11-29 18:50:25 +00:00
Christoph Kords
f666f8dcf3
added possibility to scale the plastic shearrate by the ratio of mobile to total density; keyword in material.config is "deadZone"
2012-10-29 12:49:28 +00:00
Christoph Kords
3ecaa9e371
modified multiplication term in order to mimic a discrete distribution of dislocation sources; namely, at high resolution, when there is (on average) less than one source per integration point, a random generator decides whether there is a source in the finite volume or not; if there is one, then the multiplication turns independent of the density and depends only on the velocity/volume times a scalar factor; recovers density*velocity/lambda for large finite volumes
2012-10-19 11:40:17 +00:00
Martin Diehl
324dfda5a2
added comments (doxygen conform) to phenopowerlaw, added warning on specifying h0_sliptwin as it has no effect
2012-10-11 14:49:12 +00:00
Christoph Kords
40a2e652e7
introduced material parameter "linetension", which indicates how strong the effect of the line tension on the hardening coefficients is
2012-10-05 14:42:41 +00:00
Christoph Kords
1502a71f0c
added possibility to start with an initial random distribution of dislocation segments with specific overall density
2012-10-02 12:57:24 +00:00
Franz Roters
9fc46b4950
added /echo/ statements
2012-10-02 08:12:53 +00:00
Christoph Kords
8336d0705c
new numerics parameter "unitlength" that scales the model size in DAMASK; this enables to use very small models in marc without getting any numerical problems for example with contact which does not seem to work when using very small numbers
2012-09-24 06:13:26 +00:00
Martin Diehl
1cfc467705
added new cut back parameter for spectral solver and instrumented variables for use with doxygen
2012-09-13 09:48:38 +00:00
Christoph Kords
f153866030
introduced factor to control edge contribution to multiplication;
...
dislocation density below a single dislocation per IP considered not significant
2012-09-04 16:56:37 +00:00
Christoph Kords
d80b45c1d9
added relativistic correction of velocity (maximum velocity given by parameter in material.config); probability of thermal activation in velocity law not allowed to be greater than one.
...
output of creation rate for edge jogs now for own system, not collinear.
check for negative densities in dotState now uses absoluteToleranceRho instead of significantRho.
2012-08-30 07:33:13 +00:00
Christoph Kords
aa380dcb08
added separate output of dipole formation rate for edge and screw
2012-08-27 15:57:22 +00:00
Christoph Kords
4f67d04c69
new parameter "significantRho" for nonlocal constitutive law ; density below this value will hardly contribute to dislocation glide or any dislocation reaction ; old parameter "absoluteToleranceRho" is now used only for its initial purpose, namely the absolute tolerance in the state integration
2012-08-23 05:48:21 +00:00
Christoph Kords
b2aacf9ca4
rates in constitutive_results are taken directly from last converged step and not calculated again; added output of rates for annihilation of screws and deposition of edge jogs on collinear system
2012-08-16 11:03:22 +00:00
Christoph Kords
bd754dbefd
corrected debug output of evolution rates; added possibility to change the safety factor for the CFL flux condition via the material.config file; updated the example material.config due to recent changes in nonlocal constitutive
2012-08-16 09:13:38 +00:00
Pratheek Shanthraj
d0933dad7b
added filter option to the spectral solver to mitigate spurious oscillations due to gibb's phenomenon. activate by setting myfilter in config file appropriately (currently only 'none' and 'cosine' options coded). more cleaning up of AL code
2012-08-07 17:23:13 +00:00
Pratheek Shanthraj
32f8c9c6ea
moved petsc option string to numerics and added new numerics parameters to config file
2012-08-06 17:27:53 +00:00
Philip Eisenlohr
882c44a61d
changed the majority of debugging output to level "extensive"
2012-07-24 15:29:23 +00:00
Philip Eisenlohr
0fdd43b54b
minor clarification of file content and its use
2012-07-05 10:02:20 +00:00
Franz Roters
823013d485
new plasticity: none
...
see material.config for necessary parameters (lattice_structure, elastic constants)
2012-07-03 11:16:38 +00:00
Martin Diehl
1003cfa7f8
added missing space in makefile, made removing of old files in Makefile silent and updated files for tests to current versions
2012-05-23 18:08:34 +00:00
Nan Jia
bbf790a0a5
Input parameter 'QedgePerSbSystem' which denotes activation energy for each shear band system is added
2012-05-22 16:10:28 +00:00
Martin Diehl
cec3357e8f
corrected bug causing wrong element construction when using spectral solver. Also improved error messages for invalid resolution, dimension, and homogenization.
2012-05-21 09:06:02 +00:00
Franz Roters
701cf2e71b
replaced duplicated math entry by missing mesh entry
2012-05-15 14:58:47 +00:00
Martin Diehl
9a92d343dc
added Onurs modified description for tausat in j2 law.
...
Polishing phenopowerlaw
2012-04-11 14:01:02 +00:00
Martin Diehl
f2da887899
added missing ":" in mesh.f90, introduced absolute stress tolerance for spectral solver in numerics.config/numerics.f90
2012-04-11 12:57:25 +00:00
Franz Roters
7947ae36b3
changed default value of err_div_tol
2012-03-22 10:33:54 +00:00
Martin Diehl
dbdc7ddfa2
debug.config, debug.f90, DAMASK_abaqus_exp.f, DAMASK_abaqus_std.f: changed to new debug scheme (wasn't working)
...
lattice.f90, FEsolving.f90: explicitly defined public functions and variables, all others are now private
numerics.f90: changed output format of real numbers, now instead of 0.1eX 1.0e(X-1) is printed to screen
Makefile: now using correct Optimization flags for OPTIMIZATION=AGGRESSIVE
DAMASK_spectral_AL.f90: improved, but still testing. Stress BCs now seem to be handled correctly
2012-03-20 18:01:31 +00:00
Christoph Kords
5263366615
adding switch in material.config for short range stress correction
2012-03-15 14:58:12 +00:00
Christoph Kords
503f6ed040
i think it's elasticity not elaticity, right?
2012-03-15 08:46:51 +00:00
Franz Roters
36edb4dd9f
introduced keyword elasticity in material.config
...
added constitutive_TandItsTangent and special version constitutive_hooke_TandItsTangent
those are however not yet used and need to be tested
2012-03-14 16:16:11 +00:00
Pratheek Shanthraj
c2f5cebacb
simplified analytic jacobian calculation. removed Lpfrac, time_sensitive. introduced rate_sensitivity flag when calling crystallite_stressAndItsTangent that is currently set to .false. and is to be set according to which dPdF the FE solver is asking for
2012-03-14 13:56:50 +00:00
Franz Roters
c9ed989dad
replaced keyword constitution with plasticity in preparation of upcomming generalized elasticity
2012-03-12 14:09:37 +00:00
Martin Diehl
bd9667bd4b
added new, flexible debugging scheme.
...
now all modules have their own debug specification.
compiles and runs, I hope nothing is broken
did a lot of polishing
2012-03-08 20:25:28 +00:00
Pratheek Shanthraj
17e9698659
added time_sensitive flag to ignore or add time sensitive component to analytic jacobian (should be close to perturbed jacobian when ignored)
...
now uses sub_dt for analytic jacobian calculations
bug fixes
2012-03-01 19:42:43 +00:00
Nan Jia
211f66a728
added elastic stiffness matrix as crystallite output
2012-03-01 13:06:09 +00:00
Franz Roters
1042ac018f
output values used for Lp_frac and analyticJaco
...
improved comment for Lp_frac on numerics.config
2012-02-24 08:21:18 +00:00
Martin Diehl
a98832100f
introduced itmin parameter for spectral code for defining minimum number of cycles
...
removed simplified_algorthim flag because the basic scheme using the polarization field will not be implemented
introduced divergence_correction flag for making divergence criterion resolution-independent (still experimental and not set by default)
corrected output and restart frequency (now modulo on incs of current load case)
2012-02-23 16:43:17 +00:00
Franz Roters
989a2b2d39
added default entries for Lp_frac and analyticJaco
2012-02-23 10:53:15 +00:00
Christoph Kords
8dc6e256c0
choose textures and phases of example materials according to what the labels suggest
2012-02-21 13:16:59 +00:00
Christoph Kords
5ddade27fc
changed Peierls stress to more realistic values for fcc
2012-02-03 13:12:45 +00:00
Christoph Kords
d62eddc0cd
corrected kinetics law and changed parameters. For solid solution hardening there are 3 parameters: the activation energy, the concentration of obstacles that determines the activation length and meanfreepath, and the obstacle size that determines the activation volume. For the Peierls mechanism there is: the width of doublekinks that determines the activation volume and the Peierls stress for edge and screw.
...
Still testing needed to check whether the current formulation makes sense or not.
2012-02-03 12:50:54 +00:00
Christoph Kords
df931890e0
corrected typos
2012-01-26 12:50:04 +00:00
Christoph Kords
5a658d1b82
dislocation velocity output for each dislocation type (edge+,edge-,screw+,screw-)
...
for fcc constitutive_nonlocal_kinetics is called only once and the result is used for each dislocation type, for all other lattice structures we have one call for each dislocation type
2012-01-26 07:43:36 +00:00
Christoph Kords
e5407894b5
Reworked kinetics. Now there is a universal law which can be used for both fcc and bcc (and probably even hcp). It mainly consists of an athermal forest cutting part , solid solution and peierls as thermally activated processes, and viscous glide. Therefore, lots of new parameters in the material.config. Not sure yet if the values make sense though.
...
Also renamed some of the old parameters.
Unknown constitutive output raises error.
2012-01-25 17:04:37 +00:00
Christoph Kords
4980f57fd5
added total Green-Lagrange strain to crystallite outputs
2012-01-20 10:25:35 +00:00
Christoph Kords
8739b9522e
dislocation stress does not add to stress field anymore, but is only available as constitutive output for purposes of postprocessing. instead, we now use a backstress term (which is added to the resolved shear stress) that depends on the gradient of excess density on the same slip system.
2012-01-17 10:26:57 +00:00
Franz Roters
515986c4d6
renamed w0 and w0_slip to a and a_slip
2011-11-23 14:48:39 +00:00
Martin Diehl
5ef73e164a
restructured algorithm: moved into loop to reallocate fields and replan FFTW in case resolution changes during runtime ==> regridding
...
introduced parameters for selective debugging of spectral code and partly introduced the advanced divergence calculation again which is controlled by debug.config
added switch in numerics to control divergence behavior (uncorrected and corrected by phenomenological factor)
added precision directive to all values I found
2011-11-15 17:54:18 +00:00
Christoph Kords
7033e84fc7
new outputs "accumulatedshear", "velocity_edge", and "velocity_screw" (the latter two replace the output "dislocationvelocity")
2011-11-09 09:22:52 +00:00
Franz Roters
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
Franz Roters
9d091bdc3f
added parameters for shear band model
2011-09-14 13:35:28 +00:00
Martin Diehl
f527bdcc78
changed standard values according to numerics.f90
2011-08-30 11:21:39 +00:00
Philip Eisenlohr
2e9c605571
somehow forgot to move the numerics.config into the config subdir...
2011-07-21 08:40:45 +00:00
Christoph Kords
9a64f3dd85
* comments are now recognized by IO even when they do not appear at the beginning of the line
...
* no need for the zero entries in the material.config anymore
2011-05-30 09:09:19 +00:00
Christoph Kords
bbf75f3fb5
cannot omit entries for second to fourth slip system family, though they are zero
2011-05-30 07:53:41 +00:00
Franz Roters
08d39342e4
reworked restarting for compatibility with abaqus (not yet fully working)
...
added new orientation feature for direct simulation:
component type (random) asigns random orientation to an entire grain
2011-05-28 09:44:43 +00:00
Christoph Kords
a910c469df
update of nonlocal material parameters
2011-05-20 08:11:22 +00:00
Philip Eisenlohr
1baae024cb
added two more <crystallite> examples
2011-05-13 16:54:08 +00:00
Philip Eisenlohr
d366a53f4e
probably installation is a better place for the abq.env file
2011-05-11 17:10:45 +00:00
Philip Eisenlohr
e5a2d829b0
new substructure with 'include' and 'config' directories
...
renaming of mpie_... to DAMASK for main usersubroutines
extension of element outputs from 5 to 8 digits (FFT!!!)
2011-05-11 16:38:45 +00:00