Commit Graph

  • ba48e5d3ab fixed syntax errors in last commit -- hopefully no further ones... Philip Eisenlohr 2010-05-04 13:03:35 +0000
  • cecffffe2b disorientation with symmetryType == 0 returns plain misorientation Philip Eisenlohr 2010-05-04 12:54:13 +0000
  • adaeb0e1ef rotation axis is undefined for small angles and therefore gets a NaN Claudio Zambaldi 2010-05-03 18:27:19 +0000
  • 9a257fb68b moved "UserOutput" into figures of Manual Philip Eisenlohr 2010-04-29 15:31:40 +0000
  • 2494da8f51 moved Manual to its own folder. Philip Eisenlohr 2010-04-29 15:28:29 +0000
  • ea718573ef return 180° disorientation about 100 for neighbors with different phases Christoph Kords 2010-04-29 12:28:51 +0000
  • 8ab101da50 safe use of acos function in math_QuaternionToAxisAngle Christoph Kords 2010-04-29 10:33:39 +0000
  • 15e5dcf8f2 math_QuaternionToAxisAngle now safe for small or close to 2pi rotations. Philip Eisenlohr 2010-04-29 10:01:09 +0000
  • e8719cb6b8 Table driven input works wirh 2010! That leaves DDM to do for tomorrow ;-) Franz Roters 2010-04-29 09:43:31 +0000
  • b7b2b4e9c3 small bug fix. Claudio Zambaldi 2010-04-29 09:42:50 +0000
  • 9e2396ffc8 Documentation now with all files Claudio Zambaldi 2010-04-29 09:12:20 +0000
  • c744c01ec9 Here is a suggestion for a manual to the CPFEM code. Please remember to add any solved problems to the troubleshooting section... Claudio Zambaldi 2010-04-29 09:10:47 +0000
  • aaf1bf0edf fixed (first) problem with Marc2010 single job with old style input and extended precision should work now! Franz Roters 2010-04-29 08:15:29 +0000
  • c34c07a6ff statedamper has to be local (specific for each e,i,g); with a global damping we may produce spurious convergence Christoph Kords 2010-04-29 07:41:29 +0000
  • 77dc16d15f print aged state of 1 1 1, which always exists, to avoid segmentation fault Christoph Kords 2010-04-29 07:34:20 +0000
  • 79f286f24c small typo Claudio Zambaldi 2010-04-28 17:45:30 +0000
  • 0b83fa0fb2 corrected (?) disorientation calc and introduced some new assisting functions Claudio Zambaldi 2010-04-28 17:19:58 +0000
  • 35cebfb132 extended output string length for init msg Claudio Zambaldi 2010-04-28 17:19:06 +0000
  • a36951a6ac changed common block files to release version Franz Roters 2010-04-28 07:28:12 +0000
  • e8ab495a6c Introduced new parameter NaN, not a number. It should be assigned for invalid results, such as currently for misorientation of unknown crystal symmetry. Claudio Zambaldi 2010-04-19 13:40:22 +0000
  • cab2261693 New array crystallite_symmetryID(i,g,e) is now filled during initialization run with 1 for bcc and fcc phase and 2 for hexagonal phase. The values are needed for misorientation calculations to apply the correct symmetry operators for cubic and hexagonal phases. Claudio Zambaldi 2010-04-19 10:03:34 +0000
  • 249042c2d3 New output can be requested from crystallite: (output) grainrotation it gives the deviation from the initial grain orientation in axis-angle representation with the angle in degrees. Claudio Zambaldi 2010-04-12 11:14:36 +0000
  • 653837046e new: math_QuaternionToAxisAngle Claudio Zambaldi 2010-04-12 11:07:25 +0000
  • adf6206b71 corrected symmetry operations for hex. Christoph Kords 2010-04-12 08:04:26 +0000
  • fc3f35cf3b debugger = .false. Claudio Zambaldi 2010-04-08 09:30:25 +0000
  • 1730a90e12 corrected ip numbers in FE_ipNeighbor for element type 21; used to get wrong neighborhood for this element type Christoph Kords 2010-04-06 11:45:23 +0000
  • 9b5047c828 corrected two OMP CRITICAL terminations Philip Eisenlohr 2010-04-06 06:47:15 +0000
  • 1bf740eb9d an allocation statement was deleted from rev 518 to rev 520!!! Christoph Kords 2010-03-30 13:15:18 +0000
  • 67c523e4b1 enclosed output statements in an "if (verboseDebugger)" Christoph Kords 2010-03-25 08:27:42 +0000
  • 9026cc4016 again corrections. but this version must work. Denny Tjahjanto 2010-03-24 16:23:21 +0000
  • 414050303b forgotten: the change in the input config Denny Tjahjanto 2010-03-24 13:28:46 +0000
  • 40b1478dac the latest RGC model + corrections for "element homogeneous" feature Denny Tjahjanto 2010-03-24 13:20:12 +0000
  • 3aa2dd5fef removed special characters from twip steel section in material.config changed crirtical twin stress in constitutive_dislotwin.f90 Franz Roters 2010-03-24 08:17:27 +0000
  • 2c6c8fa04e ensure positive eulerangles in math_QuaternionToEuler Christoph Kords 2010-03-19 16:11:53 +0000
  • f4e5d32cde added flag verboseDebugger to debug in order to have more control about debugging statements Christoph Kords 2010-03-19 14:14:08 +0000
  • 418bb9ecbe small improvement in readability Claudio Zambaldi 2010-03-19 12:56:00 +0000
  • 948c119ee9 changed internal representation of orientation and misorientation from euler angles to quaternions (this should also fix some bugs in the math_misorientation subroutine). includes a couple of new functions in math.f90 and some changes in crystallite.f90. Christoph Kords 2010-03-18 12:23:17 +0000
  • e6fdfdfc36 iJacoStiffness belongs to numerics, not debug sorry for not having checked that before Christoph Kords 2010-03-17 17:00:54 +0000
  • 1a992e1aef nonlocal stiffness calculation produced segmentation fault for cycle number > 16. Corrected calculation of perturbation indices k and l. Christoph Kords 2010-03-10 09:53:41 +0000
  • e3dd7551dc added output for dislocation evolution rate associated with flux: "rho_dot_flux" Christoph Kords 2010-03-10 09:49:40 +0000
  • 4752bbe42e deleted unused <global> part from file Philip Eisenlohr 2010-03-09 15:36:07 +0000
  • 61bd0224c1 - corrected an if statement in the state loop - nonlocal stiffness calculation: we perturb all material points at the same time, so instead of N^2 loops we just need N - set "forceLocalStiffnessCalculation" to false as standard Christoph Kords 2010-03-04 17:27:39 +0000
  • a0d28ebc18 - flux calculation is now also compatible to neighborhood of local constitution - flux density interpolation now depends on the position of the interface between ttwo neighboring material points - simplified flux calculation scheme - introduced sanity check for dislocation velocity to ensure v*dt< cellsize Christoph Kords 2010-03-04 17:14:47 +0000
  • 947af80a2e removed old debugging statements Christoph Kords 2010-02-26 12:02:49 +0000
  • 8c8ed34356 reworked crystallite part to allow for flexible user output Philip Eisenlohr 2010-02-25 17:39:11 +0000
  • 538faecf45 new output "dislocationvelocity" Christoph Kords 2010-02-23 17:23:07 +0000
  • c1717329ac only need 2 logicals in crystallite_stressAndItsTangent to reflect all possible integration states, so crystallite_onTrack was obsolete and replaced by crystallite_todo; broken state update now directly produces a cutback Christoph Kords 2010-02-23 09:46:39 +0000
  • 7526a24a1b now with functionality to output a "JOB.outputHomogenization". This contains size descriptions of all requested output variables to be parsed by a further script to generate meaningful names for user defined vars... (watch out!) Claudio Zambaldi 2010-02-19 18:03:16 +0000
  • aab5598d96 detabbed file for better readability Christoph Kords 2010-02-19 13:44:38 +0000
  • 8c5852dedf corrected messed up output of grain orientation in t16-file. Philip Eisenlohr 2010-02-18 16:09:45 +0000
  • e9aac0c4af (1) change phaseOfGrain from real to integer and (2) introducing a condition to switch off the grain reshuffling when element homogeneous is requested. Denny Tjahjanto 2010-02-18 15:54:10 +0000
  • d60d8a4103 small readability improvement Claudio Zambaldi 2010-02-18 15:06:08 +0000
  • 7e632e8132 corrcected definition of forceLocalStiffnessCalculation Franz Roters 2010-02-18 14:43:40 +0000
  • e1a25e0f73 renaming of abaqus to _std and adding of VUMAT interface for explicit. Philip Eisenlohr 2010-02-18 10:23:42 +0000
  • af06e8e58f fixed some dangerous constructs in computation mode assignments Philip Eisenlohr 2010-02-18 10:23:02 +0000
  • 78808a7d47 corrected corrupted header $id Philip Eisenlohr 2010-02-18 10:18:15 +0000
  • 4cb7254a21 added new computation modes 8 and 9 for Abaqus_exp. These correspond to mode 1 and 2 but do collection and calculation in one step. Philip Eisenlohr 2010-02-18 10:15:08 +0000
  • c3cd75c2c2 small polishing Philip Eisenlohr 2010-02-18 10:12:45 +0000
  • 9c63005c8f allow for longer pathes when openening files corrected error messages concerning file openening Franz Roters 2010-02-18 08:29:57 +0000
  • 97db70cf23 constitutive_nonlocal: - reworked contribution of immobile dislocation density for rate equations - flux is now calculated on the basis of interpolated velocities and densities at the interface; both incoming and outgoing fluxes are considered, so every material point only changes his own dotState - dislocation velocity is now globally defined and calculated by subroutine constitutive_nonlocal_kinetics; the subroutine is called inside _LpAndItsTangent as well as _microstructure; therefore, microstructure now needs Tstar_v as additional input; in the future one should perhaps create a subroutine constitutive_kinetics that calls constitutive_nonlocal_kinetics separately, to clearly distinguish between microstructural and kinetic variables - better use flux density vector as output variable instead of scalar flux values for each interface - added output variables internal and external resolved stress Christoph Kords 2010-02-17 13:21:36 +0000
  • a2037db9c6 corrected sign of screw part of dislocation stress Christoph Kords 2010-02-12 12:54:00 +0000
  • f1daf9f2e7 added comment, my first 'test' SVN commit Claudio Zambaldi 2010-02-10 16:40:50 +0000
  • f186c83068 introduced racing condition breaker for the case of multithreaded CPFEM init. Philip Eisenlohr 2010-02-02 15:00:08 +0000
  • 101d4c6bb4 updated list Franz Roters 2010-02-02 13:05:17 +0000
  • 02b98100d1 added / updated files for use with marc2010r1 Franz Roters 2010-01-11 08:35:28 +0000
  • ed8446e0f3 - corrected contribution of immobile singles to dipole formation and annihilation rates - remobilization of immobile singles immediately increases the mobile single dislocation density and therefore directly affects all other mechanisms in dotState Christoph Kords 2010-01-06 09:54:00 +0000
  • dce840c0a5 - introduction of immobile (used) single dislocation densities as new state variables - changed nomenclature (rho -> rhoSgl) to distinguish precisely between single dislocation density and total dislocation density - changed material.config accordingly Christoph Kords 2010-01-05 16:07:24 +0000
  • 2098c106ad misorientation is set to 400 degrees if our neighbor has a different phase Christoph Kords 2010-01-05 16:03:45 +0000
  • aee0721ab2 corrected math_misorientation Christoph Kords 2010-01-05 16:02:00 +0000
  • ac080e0b52 corrected small mistake in calculation for ip center of gravity / ip coordinates; used to ignore subnodes that had coordinates (0.0, 0.0, 0.0) Christoph Kords 2009-12-22 13:23:20 +0000
  • 043356e8a9 homogenization_RGC.f90 >> just switching-off one of the debugging statement. crystallite.f90 >> optimizing the perturbation algorithm for local tangent. Denny Tjahjanto 2009-12-22 12:28:02 +0000
  • 1f7aebfa4d new subroutine crystallite_orientations: calculates the orientation matrix and the euler angles and - in case of single grain ips - the misorientation for each neighboring ip. Christoph Kords 2009-12-18 15:46:33 +0000
  • 3ab5cdc770 Additional feature to the RGC scheme: The volumetric consistency constraint. This is implemented by penalizing any volumetric discrepancy which occurs due to severe relaxation. The penalty for volume inconsistency is described using a power-law model. The setting for this function can be done in the numerics.config via three new parameters, i.e., maxvoldiscrepancy_RGC, voldiscrepancymod_RGC and discrepancypower_RGC. For monintoring, an extra output variable in prescribed in the material.config in the RGC section, namely "(output) volumediscrepancy". Denny Tjahjanto 2009-12-16 16:20:53 +0000
  • 39d1932e2b corrected line parsing now counts all parts Franz Roters 2009-12-15 16:03:53 +0000
  • 01b7e88e5f adapted the material.config file to last changes in constitutive_nonlocal Christoph Kords 2009-12-15 10:22:17 +0000
  • 26e9056f90 added stiffness state loop counter to nonlocal stiffness calculation Christoph Kords 2009-12-15 09:03:04 +0000
  • d784153e0c constitutive_nonlocal: - corrected flux term - multiplication is now aware of dislocation type - corrected change rate for "dipole size" dupper - corrected term for dipole dissociation by stress change - added transmissivity term in fluxes which accounts for misorientation between two neighboring grains (yet hardcoded transmissivity according to misorientation angle) - added more output variables Christoph Kords 2009-12-15 08:20:31 +0000
  • 07303d9506 - corrected math_misorientation subroutine Christoph Kords 2009-12-15 07:34:53 +0000
  • 596b7f3727 - rearranged array with symmetry operations - new subroutine "math_misorientation" which calculates misorientations (based on old subroutine "misorientation", which was not used previously) - rendered some functions pure Christoph Kords 2009-12-14 11:02:10 +0000
  • 6e1c731175 crystallite does not report volume fraction but real volume (i.e. sum of all volumes gives volume of your geometry) Philip Eisenlohr 2009-12-09 15:20:46 +0000
  • cad278c3bc Introduce a randomness in the odd stress, using pre-factor rnd. Denny Tjahjanto 2009-12-09 10:33:00 +0000
  • d847f63e0d Added summary of the constitutive model for mechanical twinning as implemented in DisloTwin.f90 (both pdf and corresponding tex source) Luc Hantcherli 2009-12-08 15:58:39 +0000
  • 9b1b612f67 added some examples Franz Roters 2009-12-08 15:31:50 +0000
  • 1bf3b978eb modifying the criteria for the sum of microstructure fraction. Denny Tjahjanto 2009-12-02 13:08:14 +0000
  • 59d22d47b2 introduction of a possibility to have homogeneous element (all ips in one element are identical, sort of reduced integration) and bugs-fixing in crystallite.f90, homogenization_RGC.f90, numerics.f90. Denny Tjahjanto 2009-11-24 15:00:25 +0000
  • 6702cb3baa Introduction of a numerical viscosity into the RGC scheme to improve the convergent behavior and the numerical stability of the scheme (see: homogenization_RGC.f90). These changes have been incorporated in all related subroutines (homogenization.f90, IO.f90, numerics.config and numerics.f90). Denny Tjahjanto 2009-11-17 13:42:38 +0000
  • cb88019aa6 introduced a flexibility in cut-backing scheme in homogenization.f90 and in crystallite.f90: (1) subStepSizeHomog and subStepSizeCryst := size of substep when cut-back is applied (initially was hard-coded). (2) stepIncreaseHomog and stepIncreaseCryst := step increase when calculation for substep converge (was also hardcoded). Denny Tjahjanto 2009-11-10 13:36:27 +0000
  • 67f87486b1 in constitutive_dislotwin.f90: in constitutive_nonlocal.f90: Luc Hantcherli 2009-11-06 17:46:12 +0000
  • b03f229613 climb equation corrected Philip Eisenlohr 2009-11-02 14:20:21 +0000
  • f96f9332cf added some comments and made sure that texture distribution among a large number of IPs is statistically sound. Philip Eisenlohr 2009-11-02 08:03:14 +0000
  • 8f40c516bb added to the previous one Denny Tjahjanto 2009-10-30 10:00:17 +0000
  • 2068346f6d just small changes in the 'default' value of RGC numerical parameters. Denny Tjahjanto 2009-10-30 09:54:52 +0000
  • 1e5a0318f3 don't do any cutbacks in homogenization when my materialpoint consists of only one grain. Christoph Kords 2009-10-27 08:31:45 +0000
  • c766ba2e3a minimum substepping size in homogenization and crystallite is now controlled by two independent parameters subStepMinHomog and subStepMinCryst respectively. Christoph Kords 2009-10-26 16:43:43 +0000
  • f50d3291f9 sorry, back to the previous version. grain deformation output is part of crystallite.f90, no longer done by homogenization. material.config reflects the changes. Denny Tjahjanto 2009-10-22 16:59:24 +0000
  • 9eb1a6c3fb new, more beautiful version of the UserOutput.pdf Denny Tjahjanto 2009-10-22 15:46:25 +0000
  • d344aaf3a0 added possibility to output individual grain deformation gradient in homogenization_RGC.f90 and the homogenization.f90 has been modified accordingly. see material.config for the key-words to output the individual grain deformation gradient in the homogenization block. Denny Tjahjanto 2009-10-22 09:24:05 +0000
  • 1dcb80753c tidying up code and reversing order in debug timing output Philip Eisenlohr 2009-10-22 09:14:17 +0000
  • 67f4bdfaa0 volume fraction where new twins may form is now limited to (1-F) instead of 1... "relevantResistance" is set to 1 Pa as default Philip Eisenlohr 2009-10-22 08:58:14 +0000
  • 5d01e30f77 deleted debugging statements that print out interaction/hardening matrixes Philip Eisenlohr 2009-10-21 13:41:49 +0000