Commit Graph

  • 5d5e7eb8f8 corrected uncritical mix-up of Ntwin and Nslip in reshaping fcc interaction (slip--twin) Philip Eisenlohr 2009-10-21 13:41:09 +0000
  • 4e98935287 changed 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
  • 9892fd9a10 sorry: initialization of periodicBC was erroneous Christoph Kords 2009-10-20 16:38:05 +0000
  • e1586d5059 fixed (probably severe) bug in hexagonal twin system transformation into ortho-hexagonal coordinates. Philip Eisenlohr 2009-10-20 16:13:25 +0000
  • c1ee34d235 constitutive_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
  • 48e642bb15 twin 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
  • 7914f80931 stiffness 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
  • a94e551c5f absolute 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
  • ba27947142 material 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
  • b5fc0a9769 small rearrangement of cutback dealings. may give minute speed improvement... Philip Eisenlohr 2009-10-19 12:46:02 +0000
  • 5a5d461a30 changed 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
  • 992077e88e fixed 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
  • 974116808b fixed a potential memory leak for hexagonal structures. added some status output to constitutive_xx Philip Eisenlohr 2009-10-15 20:02:52 +0000
  • 5adcb5e36c now 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
  • 7f84961c93 less debug clutter... Philip Eisenlohr 2009-10-14 13:23:52 +0000
  • c3cb72d356 put some comments to hybridIA Philip Eisenlohr 2009-10-14 13:21:03 +0000
  • 00450b2c16 Added 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
  • 3bf3d7ddbd New version of constitutive_dislotwin after first-round checks Luc Hantcherli 2009-10-13 15:05:21 +0000
  • 2c1a610b14 removed 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
  • ca1c61bf98 Changed 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
  • 2c06cad4d0 now $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
  • 02602144d2 can now be executed from different directories...but will operate on "its" one Philip Eisenlohr 2009-10-12 14:18:54 +0000
  • 3f41b271fb The 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
  • 4ff497e5ba mesh_build_sharedElems does not need an input argument anymore Christoph Kords 2009-10-08 09:23:15 +0000
  • e7247551f0 changed "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
  • 8cb156dc35 in 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
  • bb8a71557b Added comment on FCC different interaction types 1 --> self 2 --> coplanar 3 --> collinear 4 --> Hirth locks 5 --> glissile junctions 6 --> Lomer locks Luc Hantcherli 2009-10-07 17:00:35 +0000
  • 9b3a59646a constitutive_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
  • 7d6729f1c5 Added summary of dislocation model as implemented in constitutive_dislotwin.f90 Luc Hantcherli 2009-10-07 14:28:55 +0000
  • cc6aa1aeba Updated constitutive_dislotwin Deleted constitutive_dislobased Luc Hantcherli 2009-10-06 15:19:19 +0000
  • bc5a408c1a Changed 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
  • c5865bb7ec corrected calculation of dTdLp Franz Roters 2009-10-02 06:40:12 +0000
  • c67412c36b added missing constitutive outputs for nonlocal law Christoph Kords 2009-10-01 16:31:06 +0000
  • b09b2b17f3 convergence 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
  • ec524a1dd2 homogenization_RGC.f90 Denny Tjahjanto 2009-09-07 15:28:32 +0000
  • 5f7e5b8409 homogenization_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
  • a16dbe4d32 material_config: added w0_slip in [Aluminum_phenopowerlaw] Franz Roters 2009-08-31 15:38:09 +0000
  • 6d3c620e46 forgot to "use" the new init-routines Franz Roters 2009-08-31 15:21:15 +0000
  • 2029b23f98 added version information to all files Franz Roters 2009-08-31 15:09:15 +0000
  • bb0840ecb2 constitutive_phenopowerlaw.f90: check for non-zero w0_slip Franz Roters 2009-08-31 14:13:10 +0000
  • 4357e29e07 yet another missing "then" in homogenization Franz Roters 2009-08-28 15:26:18 +0000
  • c0d1140234 added a missing "then" Franz Roters 2009-08-28 15:14:27 +0000
  • 3ea17dd2fd CPFEM.f90 >>> A very small, not important thing. Denny Tjahjanto 2009-08-28 15:05:14 +0000
  • f27d66a9ae constitutive_nonlocal: - dipole dislocations with evolution Christoph Kords 2009-08-28 13:50:47 +0000
  • 644265381e mpie_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
  • b25396374a homogenization_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
  • f70013ee9b lattice.f90: correction of typo from the last update (Philip) Denny Tjahjanto 2009-08-26 09:02:22 +0000
  • 7720680f3a makeMe: some polishment lattice: found memory leak (lattice_interactionTwinSlip had wrong size) and added further commenting Philip Eisenlohr 2009-08-26 07:28:43 +0000
  • 1dbd0865db constitutive_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
  • 387195e036 that was tooo quick, now it works hopefully Franz Roters 2009-08-13 13:32:17 +0000
  • f400cdac5b corrected paramter check for phenopowerlaw Franz Roters 2009-08-13 13:27:14 +0000
  • cafc39b183 replaces Nast with Nreps in hybridIA Denny Tjahjanto 2009-08-13 13:21:22 +0000
  • 0ef4ec27b3 there was a lonesome endif that needed a comment Franz Roters 2009-08-13 12:56:07 +0000
  • 93543c21e9 In crystallite_stressAndItsTangent: state is now correctly collected during perturbation method Luc Hantcherli 2009-08-13 10:04:14 +0000
  • 86211bf0ce homogenization_RGC.f90: adding comments to improve clarity Denny Tjahjanto 2009-08-12 14:13:20 +0000
  • 2aae7b7574 added dislocation multiplication to dotState in constitutive_nonlocal.f90 Christoph Kords 2009-08-12 11:22:02 +0000
  • 6171361c7e terminallyIll was missing in the use statement Franz Roters 2009-08-12 08:10:28 +0000
  • 8ed3ddc03b now with first draft of nonlocal constitutive law Christoph Kords 2009-08-11 16:31:57 +0000
  • 4689966c79 correction 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
  • 1a7eb3158b wrapped lines in order to ensure maximum number of allowed characters per line = 132 Christoph Kords 2009-08-03 06:37:37 +0000
  • 360fb069ba List 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
  • f0729a2e52 Cosmetic changes in constitutive_dislobased.f90 - more precise output labels - do changes into forall Luc Hantcherli 2009-07-28 10:08:13 +0000
  • 1cd0b74b08 divide by 3.0 instead of 3 as it is _pReal Franz Roters 2009-07-27 09:23:39 +0000
  • 53cff7458c polishing of constitutive output size determination in _init() Philip Eisenlohr 2009-07-24 15:23:45 +0000
  • 9bbde133a1 Updated 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
  • 97d8a321e0 some cosmetics... Philip Eisenlohr 2009-07-23 13:33:53 +0000
  • f337847f35 quite some changes: Philip Eisenlohr 2009-07-22 16:07:19 +0000
  • 290410b3fc Example of material.config, delivering the required information for the new version of constitutive_dislobased Luc Hantcherli 2009-07-14 11:31:06 +0000
  • ce515beb39 THIS 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
  • d3343ef795 Corrected 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
  • 1c35103524 Corrected 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
  • 4aed2ade80 Major Update: all modules are now correctly submitted Luc Hantcherli 2009-07-01 10:55:31 +0000
  • a16b8a619d Major update: corrected treatment of temperature Luc Hantcherli 2009-07-01 10:29:35 +0000
  • 2e783df5ed corrected typo in the prime number function Christoph Kords 2009-06-29 15:29:07 +0000
  • a6ccfe2e44 restructured 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
  • 06ed8c8177 corrected the constitutive law for j2-plasticity in Lp and its tangent and dotState. Christoph Kords 2009-06-23 10:39:29 +0000
  • 0873231e78 corrected NaN check in crystallite_updateState return false in case NaN occurs in crystallite_updateState Franz Roters 2009-06-23 06:52:51 +0000
  • 660dce0f09 now 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
  • b5b4b32961 move operation to "phenoPowerLaw" Philip Eisenlohr 2009-06-19 07:03:43 +0000
  • 8e70963008 move operation into "Code" part of documentation Philip Eisenlohr 2009-06-19 07:02:57 +0000
  • ef8b672a00 numerics.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
  • 904ea78ac5 Tstar_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
  • ada92a9b74 added 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
  • 204e296ecd modify document for ConstitutivePhenoPowerLaw Yun Jo Ro 2009-06-15 08:09:49 +0000
  • caf568eb89 adapted 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
  • a8ff024d97 modify document for ConstitutivePhenoPowerLaw Yun Jo Ro 2009-06-10 14:43:30 +0000
  • 118e68e881 modify document for ConstitutivePhenoPowerLaw Yun Jo Ro 2009-06-09 13:48:46 +0000
  • c17249c953 Add 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
  • 3196049496 in 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
  • 306cd95992 sort 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
  • db731a4be3 Taylor 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
  • 9d4521b7ec modify document for ConstitutivePhenoPowerLaw Yun Jo Ro 2009-06-08 09:29:54 +0000
  • 5a32d644fd modify document for ConstitutivePhenoPowerLaw Yun Jo Ro 2009-06-05 12:50:09 +0000