9e2396ffc8Documentation now with all files
Claudio Zambaldi
2010-04-29 09:12:20 +0000
c744c01ec9Here 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
aaf1bf0edffixed (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
c34c07a6ffstatedamper 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
77dc16d15fprint aged state of 1 1 1, which always exists, to avoid segmentation fault
Christoph Kords
2010-04-29 07:34:20 +0000
a36951a6acchanged common block files to release version
Franz Roters
2010-04-28 07:28:12 +0000
e8ab495a6cIntroduced 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
cab2261693New 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
249042c2d3New 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
1730a90e12corrected 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
9b5047c828corrected two OMP CRITICAL terminations
Philip Eisenlohr
2010-04-06 06:47:15 +0000
1bf740eb9dan allocation statement was deleted from rev 518 to rev 520!!!
Christoph Kords
2010-03-30 13:15:18 +0000
67c523e4b1enclosed output statements in an "if (verboseDebugger)"
Christoph Kords
2010-03-25 08:27:42 +0000
9026cc4016again corrections. but this version must work.
Denny Tjahjanto
2010-03-24 16:23:21 +0000
414050303bforgotten: the change in the input config
Denny Tjahjanto
2010-03-24 13:28:46 +0000
40b1478dacthe latest RGC model + corrections for "element homogeneous" feature
Denny Tjahjanto
2010-03-24 13:20:12 +0000
3aa2dd5fefremoved 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
2c6c8fa04eensure positive eulerangles in math_QuaternionToEuler
Christoph Kords
2010-03-19 16:11:53 +0000
f4e5d32cdeadded flag verboseDebugger to debug in order to have more control about debugging statements
Christoph Kords
2010-03-19 14:14:08 +0000
418bb9ecbesmall improvement in readability
Claudio Zambaldi
2010-03-19 12:56:00 +0000
948c119ee9changed 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
e6fdfdfc36iJacoStiffness belongs to numerics, not debug sorry for not having checked that before
Christoph Kords
2010-03-17 17:00:54 +0000
1a992e1aefnonlocal 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
e3dd7551dcadded output for dislocation evolution rate associated with flux: "rho_dot_flux"
Christoph Kords
2010-03-10 09:49:40 +0000
4752bbe42edeleted 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
947af80a2eremoved old debugging statements
Christoph Kords
2010-02-26 12:02:49 +0000
8c8ed34356reworked crystallite part to allow for flexible user output
Philip Eisenlohr
2010-02-25 17:39:11 +0000
538faecf45new output "dislocationvelocity"
Christoph Kords
2010-02-23 17:23:07 +0000
c1717329aconly 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
7526a24a1bnow 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
aab5598d96detabbed file for better readability
Christoph Kords
2010-02-19 13:44:38 +0000
8c5852dedfcorrected 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
7e632e8132corrcected definition of forceLocalStiffnessCalculation
Franz Roters
2010-02-18 14:43:40 +0000
e1a25e0f73renaming of abaqus to _std and adding of VUMAT interface for explicit.
Philip Eisenlohr
2010-02-18 10:23:42 +0000
af06e8e58ffixed some dangerous constructs in computation mode assignments
Philip Eisenlohr
2010-02-18 10:23:02 +0000
78808a7d47corrected corrupted header $id
Philip Eisenlohr
2010-02-18 10:18:15 +0000
4cb7254a21added 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
c3cd75c2c2small polishing
Philip Eisenlohr
2010-02-18 10:12:45 +0000
9c63005c8fallow for longer pathes when openening files corrected error messages concerning file openening
Franz Roters
2010-02-18 08:29:57 +0000
97db70cf23constitutive_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
a2037db9c6corrected sign of screw part of dislocation stress
Christoph Kords
2010-02-12 12:54:00 +0000
f1daf9f2e7added comment, my first 'test' SVN commit
Claudio Zambaldi
2010-02-10 16:40:50 +0000
f186c83068introduced racing condition breaker for the case of multithreaded CPFEM init.
Philip Eisenlohr
2010-02-02 15:00:08 +0000
101d4c6bb4updated list
Franz Roters
2010-02-02 13:05:17 +0000
02b98100d1added / 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
2098c106admisorientation is set to 400 degrees if our neighbor has a different phase
Christoph Kords
2010-01-05 16:03:45 +0000
aee0721ab2corrected math_misorientation
Christoph Kords
2010-01-05 16:02:00 +0000
ac080e0b52corrected 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
043356e8a9homogenization_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
1f7aebfa4dnew 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
3ab5cdc770Additional 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
39d1932e2bcorrected line parsing now counts all parts
Franz Roters
2009-12-15 16:03:53 +0000
01b7e88e5fadapted the material.config file to last changes in constitutive_nonlocal
Christoph Kords
2009-12-15 10:22:17 +0000
26e9056f90added stiffness state loop counter to nonlocal stiffness calculation
Christoph Kords
2009-12-15 09:03:04 +0000
d784153e0cconstitutive_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
6e1c731175crystallite 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
cad278c3bcIntroduce a randomness in the odd stress, using pre-factor rnd.
Denny Tjahjanto
2009-12-09 10:33:00 +0000
d847f63e0dAdded 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
9b1b612f67added some examples
Franz Roters
2009-12-08 15:31:50 +0000
1bf3b978ebmodifying the criteria for the sum of microstructure fraction.
Denny Tjahjanto
2009-12-02 13:08:14 +0000
59d22d47b2introduction 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
6702cb3baaIntroduction 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
cb88019aa6introduced 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
67f87486b1in constitutive_dislotwin.f90: in constitutive_nonlocal.f90:
Luc Hantcherli
2009-11-06 17:46:12 +0000
b03f229613climb equation corrected
Philip Eisenlohr
2009-11-02 14:20:21 +0000
f96f9332cfadded 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
8f40c516bbadded to the previous one
Denny Tjahjanto
2009-10-30 10:00:17 +0000
2068346f6djust small changes in the 'default' value of RGC numerical parameters.
Denny Tjahjanto
2009-10-30 09:54:52 +0000
1e5a0318f3don't do any cutbacks in homogenization when my materialpoint consists of only one grain.
Christoph Kords
2009-10-27 08:31:45 +0000
c766ba2e3aminimum 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
f50d3291f9sorry, 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
9eb1a6c3fbnew, more beautiful version of the UserOutput.pdf
Denny Tjahjanto
2009-10-22 15:46:25 +0000
d344aaf3a0added 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
1dcb80753ctidying up code and reversing order in debug timing output
Philip Eisenlohr
2009-10-22 09:14:17 +0000
67f4bdfaa0volume 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
5d01e30f77deleted debugging statements that print out interaction/hardening matrixes
Philip Eisenlohr
2009-10-21 13:41:49 +0000