Commit Graph

221 Commits

Author SHA1 Message Date
Christoph Kords daf89806b1 fixed bug in deltaState; used to lead to unpredictable values in accumulated shear, since deltaState was not defined at this state position 2013-06-11 21:03:23 +00:00
Martin Diehl 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
Christoph Kords a6361738fb (again!) fixed bug that was introduced in rev 2085: messed up indices of stiffness matrix in calculation of Poisson's ratio 2013-05-26 14:31:31 +00:00
Christoph Kords ab997f19ac fixed bug in nonlocal state init that was introduced in rev 2438 and used to reset states from other constitutions 2013-05-24 13:43:44 +00:00
Christoph Kords c0539d2383 replaced all remaining occurrences of state indices 2013-05-24 11:48:34 +00:00
Christoph Kords e2d970ce57 combined some state indices to an array with a more generic name 2013-05-24 09:02:30 +00:00
Christoph Kords ad4a66cb4e use named state variable indices in "deltaState" 2013-05-23 21:46:15 +00:00
Christoph Kords 332b207270 use named state variable indices in "kinetics" and "microstructure" 2013-05-23 21:40:00 +00:00
Christoph Kords 5900076d72 rename state variable index for critical resolved shear stress from "iTau" to "iTauF" 2013-05-23 21:32:36 +00:00
Christoph Kords ad2a245cdf use named state variable indices in "microstructure" 2013-05-23 21:28:45 +00:00
Christoph Kords c02645c996 use named state variable indices in "aTolState" 2013-05-23 21:13:56 +00:00
Christoph Kords f04a531e9b use named state variable indices in "stateInit" 2013-05-23 21:10:31 +00:00
Christoph Kords 65524b7f84 corrected named indices for state variables; however, not yet used 2013-05-23 20:30:06 +00:00
Christoph Kords e249366ded introduced named indices for state variables; however, not yet used 2013-05-23 20:15:23 +00:00
Christoph Kords 72129cae00 removed prefix "constitutive_nonlocal_" from all private variables 2013-05-23 19:56:36 +00:00
Christoph Kords b559666436 replaced FE_maxNipNeighbors with mesh_maxNipNeighbors, which is sufficient for memory allocation, since this gives the maximum of neighbors in all cp elements 2013-05-23 17:46:21 +00:00
Christoph Kords 77e9b80a16 bug fix for calculation of fluxdensity and slipnormal/direction in constitutive_nonlocal_postResults: indices for Fe were missing 2013-05-23 14:52:57 +00:00
Christoph Kords 41faa60323 fixed bug in backstress calculation for elements with periodic neighbors 2013-05-23 12:36:48 +00:00
Christoph Kords 3aaf60cffe polishing 2013-05-23 12:25:56 +00:00
Christoph Kords ec2503253e added output of slip direction and slip normal
changed name of fluxdensity outputs
2013-05-23 08:19:36 +00:00
Christoph Kords e758ee32be improved accuracy of accumulated shear output: take shear rate from LpAndItsTangent for time integration 2013-05-21 10:04:52 +00:00
Christoph Kords 3759024b42 forgot to set neighboring_rho (also reverting changes made in rev 2384) 2013-05-19 19:23:31 +00:00
Christoph Kords 0504fb3b21 reverted changes made in rev2384 2013-05-17 12:54:47 +00:00
Christoph Kords 6fac7c347c fixed bug in constitutive_nonlocal_kinetics that was introduced in rev2088 along with the non-Schmid behavior: nonSchmid matrix of last slip system was used for all slip systems; led to extremely bad convergence due to flawed dLp_dT 2013-05-10 22:29:12 +00:00
Christoph Kords d9cc42daef slight change in flux formulation: no flux to neighbor, not only if there is a sign change in the velocity, but also if the neighboring velocity is exactly zero 2013-05-10 08:42:42 +00:00
Martin Diehl b2a3c9235b fixed some warnings issued by gfortran (type conversion, not needed use statements) 2013-05-08 09:23:47 +00:00
Christoph Kords 6e64b5c476 better not allow for deads to contribute to flux, since not clear at which interface they got stuck 2013-05-07 18:47:17 +00:00
Christoph Kords 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
Christoph Kords ad879ebcf9 forgot to delete all occurrences of obsolete "constitutive_nonlocal_ accumulatedShear" in last commit 2013-04-08 08:51:12 +00:00
Christoph Kords 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
Christoph Kords 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
Christoph Kords 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
Christoph Kords 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
Christoph Kords 19655c2d92 significantRho check was flawed for dipole density 2013-03-25 17:36:04 +00:00
Franz Roters 02f092492e updated copyright header 2013-03-22 17:35:05 +00:00
Christoph Kords 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
Pratheek Shanthraj 07ff1f7c93 reporting current time for most module inits 2013-02-25 16:34:59 +00:00
Philip Eisenlohr 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
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
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
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
Pratheek Shanthraj 50e874f887 fixed bug in elasticity matrix calculation 2013-01-22 15:48:47 +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
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
Christoph Kords f14fd45ce1 new material parameter "egdejog": fraction of annihilated screw dipoles that forms edge jogs 2012-12-23 12:56:15 +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