5d5e7eb8f8corrected uncritical mix-up of Ntwin and Nslip in reshaping fcc interaction (slip--twin)
Philip Eisenlohr
2009-10-21 13:41:09 +0000
4e98935287changed lattice_interaction matrices to (other:me) notation fixed small mistake in hexagonal twin--slip interaction matrix
Philip Eisenlohr
2009-10-21 13:10:12 +0000
9892fd9a10sorry: initialization of periodicBC was erroneous
Christoph Kords
2009-10-20 16:38:05 +0000
e1586d5059fixed (probably severe) bug in hexagonal twin system transformation into ortho-hexagonal coordinates.
Philip Eisenlohr
2009-10-20 16:13:25 +0000
c1ee34d235constitutive_nonlocal: - read in activation energy for dislocation glide from material.config - changed naming of dDipMin/Max to dLower/dUpper - added new outputs: rho_dot, rho_dot_dip, rho_dot_gen, rho_dot_sgl2dip, rho_dot_dip2sgl, rho_dot_ann_ath, rho_dot_ann_the, rho_dot_flux, d_upper_edge, d_upper_screw, d_upper_dot_edge, d_upper_dot_screw - poisson's ratio is now calculated from elastic constants - microstrucutre has state as first argument, since this is our output variable - periodic boundary conditions are taken into account for fluxes and internal stresses. for the moment, flag has to be set in constitutive_nonlocal. - corrected calculation for dipole formation by glide - added terms for dipole formation/annihilation by stress decrease/increase
Christoph Kords
2009-10-20 14:36:03 +0000
48e642bb15twin resistance evolution did not use magnitude of shear rate (dot gamma) but the signed quantity... leads to non-physical oscillations of twin resistance.
Philip Eisenlohr
2009-10-20 12:15:37 +0000
7914f80931stiffness calculation for structure > 3 was missing: "case(3)" --> now "case(3:)" as in all other constitutive_xyz..!
Philip Eisenlohr
2009-10-19 14:47:52 +0000
a94e551c5fabsolute tolerance was just fine -- residuum is in Lp not in stress. Hence, unit is strain and not Pa...
Philip Eisenlohr
2009-10-19 12:53:56 +0000
ba27947142material point now cuts back by 1/8 crystallite does not accelerate anymore, since, typically, longer step immediately fails and uses resources in vain. (future: remember number of successful steps to increase step size after x of those...)
Philip Eisenlohr
2009-10-19 12:50:59 +0000
b5fc0a9769small rearrangement of cutback dealings. may give minute speed improvement...
Philip Eisenlohr
2009-10-19 12:46:02 +0000
5a5d461a30changed abs tolerance in stress loop to 1 Pa (from 1e-8 Pa ..?)
Philip Eisenlohr
2009-10-19 12:44:20 +0000
6e613303a3
Christoph Kords
2009-10-19 08:01:02 +0000
ba02dfca1e
Christoph Kords
2009-10-19 07:56:47 +0000
992077e88efixed cut back issue with Marc. Necessity to introduce compMode 7 to capture tangent restoration at lovl=4 call.
Philip Eisenlohr
2009-10-17 09:25:36 +0000
974116808bfixed a potential memory leak for hexagonal structures. added some status output to constitutive_xx
Philip Eisenlohr
2009-10-15 20:02:52 +0000
5adcb5e36cnow contains the number of grains "m" as first entry (thanks to Koen for appending the document)
Philip Eisenlohr
2009-10-15 11:53:07 +0000
7f84961c93less debug clutter...
Philip Eisenlohr
2009-10-14 13:23:52 +0000
c3cb72d356put some comments to hybridIA
Philip Eisenlohr
2009-10-14 13:21:03 +0000
00450b2c16Added constraint on the dipole distance if dipole distance > dislocation mean free distance then dipole distance = dislocation mean free distance
Luc Hantcherli
2009-10-13 18:14:44 +0000
3bf3d7ddbdNew version of constitutive_dislotwin after first-round checks
Luc Hantcherli
2009-10-13 15:05:21 +0000
2c1a610b14removed Koen's loading subroutine from abaqus interface
Franz Roters
2009-10-13 06:55:15 +0000
0a34da292e# no need for materialpoint_results in hypela2 # homogenization updated to be compatible with abaqus interface
Philip Eisenlohr
2009-10-12 17:01:42 +0000
ca1c61bf98Changed constitutive_dislobased into constitutive_dislotwin
Luc Hantcherli
2009-10-12 16:44:17 +0000
8b66b8def2# small remark to set "-free" switch for compile # rename to Abaqus-compatible ".f" extension
Philip Eisenlohr
2009-10-12 16:32:57 +0000
7d05845d5d# now $ID aware # LF as line separation and some tidy up
Philip Eisenlohr
2009-10-12 16:26:16 +0000
2c06cad4d0now $ID aware
Philip Eisenlohr
2009-10-12 16:25:48 +0000
21cad32137# new interface for Abaqus # IO has some additional functionality for Abaqus parsing # ping pong scheme in FE interface now similar (and more human understandable) in both versions # mesh has better splitting of different tasks, plus operation on database whenever possible # FEsolver as new global var to indicate FEM solver type # computation mode reshuffling: 6 is now Marc special case of recycling...
Philip Eisenlohr
2009-10-12 16:01:49 +0000
02602144d2can now be executed from different directories...but will operate on "its" one
Philip Eisenlohr
2009-10-12 14:18:54 +0000
3f41b271fbThe new dislocation model is now implemented in constitutive_dislotwin.f90 (seems to run fine, but still under investigation)
Luc Hantcherli
2009-10-09 20:28:23 +0000
4ff497e5bamesh_build_sharedElems does not need an input argument anymore
Christoph Kords
2009-10-08 09:23:15 +0000
e7247551f0changed "build_sharedElements" to use (then existing) internal element data structure instead of reading again from the input file. The data structure is already established by "build_elements"...
Philip Eisenlohr
2009-10-08 09:01:59 +0000
8cb156dc35in Documentation\ConstitutiveLaw\DisloTwinLaw Added tex source that generates PDF-file describing dislocation glide as implemented in constitutive_dislotwin.f90
Luc Hantcherli
2009-10-08 07:33:01 +0000
9b3a59646aconstitutive_nonlocal - take orientation gradients into account when calculating dislocation stress and dislocation fluxes - hard coded value for nu - changed kinetics (parameter G0 is currently defined as a parameter, needs to be read from material.config) - added some output statements
Christoph Kords
2009-10-07 15:31:52 +0000
7d6729f1c5Added summary of dislocation model as implemented in constitutive_dislotwin.f90
Luc Hantcherli
2009-10-07 14:28:55 +0000
bc5a408c1aChanged constitutive_dislobased into constitutive_dislotwin: in constitutive.90 in constitutive_dislotwin.f90 in mpie_cpfem_marc.f90
Luc Hantcherli
2009-10-06 15:16:03 +0000
c5865bb7eccorrected calculation of dTdLp
Franz Roters
2009-10-02 06:40:12 +0000
c67412c36badded missing constitutive outputs for nonlocal law
Christoph Kords
2009-10-01 16:31:06 +0000
b09b2b17f3convergence of state in crystallite is now tested as follows: (state < relevant state) or (residuum < relative tolerance * state) since the relevant value for the state variables depend on their nature and can vary by large scales (e.g. volume fraction: 1e-10, dislocation density: 1e5) it is not possible to set a unique value. instead the constitutive law has to decide what is relevant. therefore, all constitutive laws now read in parameters from the material.config that determine the values for relevantState [@luc: in dislobased law relevant State is for the moment generally set to 1e-200, so no additional parameters necessary in material.config. if you also want this feature, we can still implement it, no big deal]
Christoph Kords
2009-09-18 15:37:14 +0000
5f7e5b8409homogenization_RGC.f90 << correction for a very stupid mistake. replacing stress P with the absolute value of stress abs(P) during convergence check. it should work upon deep-drawing.
Denny Tjahjanto
2009-09-03 20:26:20 +0000
a16dbe4d32material_config: added w0_slip in [Aluminum_phenopowerlaw]
Franz Roters
2009-08-31 15:38:09 +0000
6d3c620e46forgot to "use" the new init-routines
Franz Roters
2009-08-31 15:21:15 +0000
2029b23f98added version information to all files
Franz Roters
2009-08-31 15:09:15 +0000
bb0840ecb2constitutive_phenopowerlaw.f90: check for non-zero w0_slip
Franz Roters
2009-08-31 14:13:10 +0000
4357e29e07yet another missing "then" in homogenization
Franz Roters
2009-08-28 15:26:18 +0000
c0d1140234added a missing "then"
Franz Roters
2009-08-28 15:14:27 +0000
3ea17dd2fdCPFEM.f90 >>> A very small, not important thing.
Denny Tjahjanto
2009-08-28 15:05:14 +0000
f27d66a9aeconstitutive_nonlocal: - dipole dislocations with evolution
Christoph Kords
2009-08-28 13:50:47 +0000
644265381empie_cpfem_marc.f90 >>> changing the compilation sequence of modules: numerics.f90 and IO.f90 move up, while math.f90 down.
Denny Tjahjanto
2009-08-27 15:30:40 +0000
b25396374ahomogenization_RGC.f90 >>> adding some lines, mostly for debugging purpose. no critical change. constitutive_phenopowerlaw.f90 >>> adding new parameter: constitutive_phenopowerlaw_w0_slip, i.e., the hardening rate exponent. homogenization.f90 >>> most important change is to add an if-else statement (line 379-380) to switch crystallite_requeted = .false. for already converged material point iteration (el/ip). the rest of the changes are cosmetics and debugging stuffs. crystallite.f90 >>> similar to homogenization.f90, the most important change is to add additional if-else statement (line 574) in the jacobian (perturbation) loop. now the jacobian calculation will only be performed when crystallite_requested = .true.. the rest is only cosmetic.
Denny Tjahjanto
2009-08-27 12:10:06 +0000
f70013ee9blattice.f90: correction of typo from the last update (Philip)
Denny Tjahjanto
2009-08-26 09:02:22 +0000
7720680f3amakeMe: some polishment lattice: found memory leak (lattice_interactionTwinSlip had wrong size) and added further commenting
Philip Eisenlohr
2009-08-26 07:28:43 +0000
1dbd0865dbconstitutive_nonlocal.f90 - completed postResults output function - connecting vector of neighboring material points is mapped to intermediate configuration of my neighbor
Christoph Kords
2009-08-24 08:16:01 +0000
387195e036that was tooo quick, now it works hopefully
Franz Roters
2009-08-13 13:32:17 +0000
f400cdac5bcorrected paramter check for phenopowerlaw
Franz Roters
2009-08-13 13:27:14 +0000
cafc39b183replaces Nast with Nreps in hybridIA
Denny Tjahjanto
2009-08-13 13:21:22 +0000
0ef4ec27b3there was a lonesome endif that needed a comment
Franz Roters
2009-08-13 12:56:07 +0000
93543c21e9In crystallite_stressAndItsTangent: state is now correctly collected during perturbation method
Luc Hantcherli
2009-08-13 10:04:14 +0000
2aae7b7574added dislocation multiplication to dotState in constitutive_nonlocal.f90
Christoph Kords
2009-08-12 11:22:02 +0000
6171361c7eterminallyIll was missing in the use statement
Franz Roters
2009-08-12 08:10:28 +0000
8ed3ddc03bnow with first draft of nonlocal constitutive law
Christoph Kords
2009-08-11 16:31:57 +0000
4689966c79correction in homogenization_RGC.f90: in the use of material module homogenization_typeInstance was included two times in previous version.
Denny Tjahjanto
2009-08-06 12:47:14 +0000
1a7eb3158bwrapped lines in order to ensure maximum number of allowed characters per line = 132
Christoph Kords
2009-08-03 06:37:37 +0000
360fb069baList of changes/modifications: * IO.f90 :: Adding error messages for RGC homogenization * crystallite.f90 :: Modifying convergent criteria in crystalline_updateState and crystalline_updateTemperature * material.f90 :: Adding IO_lc in homogenization_type(Name) * debug.f90 :: Adding debbugging statement and counter for material point loop * homogenization.f90 :: Adding homogenization_RGC blocks * homogenization_isostrain.f90 :: Modifying argument of homogenization_isostrain_stateInit * mpie_cpfem_marc.f90 :: Adding homogenization_RGC include * numerics.f90 :: Adding numerical parameters for RGC scheme * math.f90 :: Changing function name: math_permut to math_civita
Denny Tjahjanto
2009-07-31 12:02:20 +0000
f0729a2e52Cosmetic changes in constitutive_dislobased.f90 - more precise output labels - do changes into forall
Luc Hantcherli
2009-07-28 10:08:13 +0000
1cd0b74b08divide by 3.0 instead of 3 as it is _pReal
Franz Roters
2009-07-27 09:23:39 +0000
53cff7458cpolishing of constitutive output size determination in _init()
Philip Eisenlohr
2009-07-24 15:23:45 +0000
9bbde133a1Updated constitutive_dislobased.90 and lattice.f90 according to the changes done in the previous release - introduction of system families - analogous structure to constitutive_phenopowerlaw.f90 - interaction types 0-1 changed into 1-2
Luc Hantcherli
2009-07-24 12:02:20 +0000
97d8a321e0some cosmetics...
Philip Eisenlohr
2009-07-23 13:33:53 +0000
f337847f35quite some changes:
Philip Eisenlohr
2009-07-22 16:07:19 +0000
290410b3fcExample of material.config, delivering the required information for the new version of constitutive_dislobased
Luc Hantcherli
2009-07-14 11:31:06 +0000
ce515beb39THIS IS A MAJOR UPDATE The different blocks required for the twinning model are now implemented (I guess correctly...) Keywords in the material.config are changed. Since the flow rule for twin systems remains under investigation, this part is susceptible to further changes.
Luc Hantcherli
2009-07-14 11:26:52 +0000
d3343ef795Corrected definition of vector n and d per deformation system; n is the system plane normal and d the corresponding direction
Luc Hantcherli
2009-07-14 11:20:35 +0000
fe2f3e4d36- added sanity check for temperature tolerance - added temperature tolerance in numerics.config
Christoph Kords
2009-07-02 16:37:48 +0000
1c35103524Corrected CPFEM: assignment of constitutive_state and _state0 consistent with homogenization_Ngrains (what is actually needed), instead of homogenization_maxMgrains
Luc Hantcherli
2009-07-02 08:57:35 +0000
4aed2ade80Major Update: all modules are now correctly submitted
Luc Hantcherli
2009-07-01 10:55:31 +0000
a16b8a619dMajor update: corrected treatment of temperature
Luc Hantcherli
2009-07-01 10:29:35 +0000
2e783df5edcorrected typo in the prime number function
Christoph Kords
2009-06-29 15:29:07 +0000
a6ccfe2e44restructured the repository: renamed trunk into code, moved documentation one directory up, therefore you can now checkout either code or documentation only if you like
Franz Roters
2009-06-25 07:17:59 +0000
06ed8c8177corrected the constitutive law for j2-plasticity in Lp and its tangent and dotState.
Christoph Kords
2009-06-23 10:39:29 +0000
0873231e78corrected NaN check in crystallite_updateState return false in case NaN occurs in crystallite_updateState
Franz Roters
2009-06-23 06:52:51 +0000
660dce0f09now safe to request values from config-parsing which are not present. 3rd and up value of key "xyz" would return "0" or "0.0" in this example: xyz 1.234 4.321
Philip Eisenlohr
2009-06-19 07:09:39 +0000
b5b4b32961move operation to "phenoPowerLaw"
Philip Eisenlohr
2009-06-19 07:03:43 +0000
8e70963008move operation into "Code" part of documentation
Philip Eisenlohr
2009-06-19 07:02:57 +0000
ef8b672a00numerics.config is now optional; if this file does not exist in the working directory, standard values from numerics.f90 are used; either way the output file tells which values were taken.
Christoph Kords
2009-06-18 14:28:02 +0000
904ea78ac5Tstar_v is now restored at different stages of cutbacking; we need this because preguess of state relies on consistent Tstar_v
Christoph Kords
2009-06-16 09:03:30 +0000
ada92a9b74added a new module called numerics.f90 which reads in all numerical "parameters" from the file numerics.config (also being added). From now on this file has to be located in the working directory of the FEM-model and has to contain all necessary parameters.
Christoph Kords
2009-06-15 13:11:21 +0000
204e296ecdmodify document for ConstitutivePhenoPowerLaw
Yun Jo Ro
2009-06-15 08:09:49 +0000
caf568eb89adapted crystallite_stressAndItsTangent to do a pre-guess for the state before the actual state loop (with order stress integration, state update)
Christoph Kords
2009-06-10 15:08:33 +0000
a8ff024d97modify document for ConstitutivePhenoPowerLaw
Yun Jo Ro
2009-06-10 14:43:30 +0000
118e68e881modify document for ConstitutivePhenoPowerLaw
Yun Jo Ro
2009-06-09 13:48:46 +0000
c17249c953Add document for explaining how we sort the order of the slip/twin system for HCP in lattice.f90
Yun Jo Ro
2009-06-09 13:13:11 +0000
3196049496in crystallite_stressAndItsTangent we now update the state first (explicit integration) and then integrate the stress (implicit integration)
Christoph Kords
2009-06-09 11:05:29 +0000
306cd95992sort compressive and tensile twin-2 group in Hex structure and also correct the name for each twin type
Yun Jo Ro
2009-06-09 08:42:02 +0000
db731a4be3Taylor homogenization is now called isostrain, changed material.config accordingly. corrected error message in case of unknown homogenization
Franz Roters
2009-06-08 13:28:00 +0000
9d4521b7ecmodify document for ConstitutivePhenoPowerLaw
Yun Jo Ro
2009-06-08 09:29:54 +0000
5a32d644fdmodify document for ConstitutivePhenoPowerLaw
Yun Jo Ro
2009-06-05 12:50:09 +0000