Franz Roters
f400cdac5b
corrected paramter check for phenopowerlaw
2009-08-13 13:27:14 +00:00
Denny Tjahjanto
cafc39b183
replaces Nast with Nreps in hybridIA
...
formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
2009-08-13 13:21:22 +00:00
Franz Roters
0ef4ec27b3
there was a lonesome endif that needed a comment
2009-08-13 12:56:07 +00:00
Luc Hantcherli
93543c21e9
In crystallite_stressAndItsTangent:
...
state is now correctly collected during perturbation method
2009-08-13 10:04:14 +00:00
Denny Tjahjanto
86211bf0ce
homogenization_RGC.f90: adding comments to improve clarity
2009-08-12 14:13:20 +00:00
Christoph Kords
2aae7b7574
added dislocation multiplication to dotState in constitutive_nonlocal.f90
...
cleaned up debugging output statements to *.out file
2009-08-12 11:22:02 +00:00
Franz Roters
6171361c7e
terminallyIll was missing in the use statement
2009-08-12 08:10:28 +00:00
Christoph Kords
8ed3ddc03b
now with first draft of nonlocal constitutive law
...
debugging memory leak closed
debugging counters corrected
center of gravity stored in mesh
state updated is now split into a collecting loop and an execution
updateState and updateTemperature fill sequentially separate logicals and evaluate afterwards to converged
added 3x3 transposition function, norm for 3x1 matrix and 33x3 matrix multiplication in math
non-converged crystallite triggers materialpoint cutback (used to respond elastically)
non-converged materialpoint raises terminal illness which in turn renders whole FE increment useless by means of odd stress/stiffness and thus waits for FE cutback
2009-08-11 16:31:57 +00:00
Denny Tjahjanto
4689966c79
correction in homogenization_RGC.f90: in the use of material module
...
homogenization_typeInstance was included two times in previous version.
2009-08-06 12:47:14 +00:00
Christoph Kords
1a7eb3158b
wrapped lines in order to ensure maximum number of allowed characters per line = 132
2009-08-03 06:37:37 +00:00
Denny Tjahjanto
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
2009-07-31 12:02:20 +00:00
Luc Hantcherli
f0729a2e52
Cosmetic changes in constitutive_dislobased.f90
...
- more precise output labels
- do changes into forall
2009-07-28 10:08:13 +00:00
Franz Roters
1cd0b74b08
divide by 3.0 instead of 3 as it is _pReal
2009-07-27 09:23:39 +00:00
Philip Eisenlohr
53cff7458c
polishing of constitutive output size determination in _init()
2009-07-24 15:23:45 +00:00
Luc Hantcherli
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
2009-07-24 12:02:20 +00:00
Philip Eisenlohr
97d8a321e0
some cosmetics...
2009-07-23 13:33:53 +00:00
Philip Eisenlohr
f337847f35
quite some changes:
...
# non-greedy memory allocation
# generation of outputConstitutive to allow for script-based T16 extraction
# exchange of phenomenological by more general phenopowerlaw
# lattice is based on slip and twin families which can be treated as individual entities (switched on/off, separate hardening, etc.)
# nicer debugging output
# changed some error/warning codes
# plus potentially some minor additional brushes here and there
2009-07-22 16:07:19 +00:00
Luc Hantcherli
290410b3fc
Example of material.config, delivering the required information for the new version of constitutive_dislobased
2009-07-14 11:31:06 +00:00
Luc Hantcherli
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.
2009-07-14 11:26:52 +00:00
Luc Hantcherli
d3343ef795
Corrected definition of vector n and d per deformation system; n is the system plane normal and d the corresponding direction
2009-07-14 11:20:35 +00:00
Christoph Kords
fe2f3e4d36
- added sanity check for temperature tolerance
...
- added temperature tolerance in numerics.config
2009-07-02 16:37:48 +00:00
Luc Hantcherli
1c35103524
Corrected CPFEM: assignment of constitutive_state and _state0 consistent with homogenization_Ngrains (what is actually needed), instead of homogenization_maxMgrains
2009-07-02 08:57:35 +00:00
Luc Hantcherli
4aed2ade80
Major Update: all modules are now correctly submitted
2009-07-01 10:55:31 +00:00
Luc Hantcherli
a16b8a619d
Major update: corrected treatment of temperature
2009-07-01 10:29:35 +00:00
Christoph Kords
2e783df5ed
corrected typo in the prime number function
2009-06-29 15:29:07 +00:00
Franz Roters
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
2009-06-25 07:17:59 +00:00
Christoph Kords
06ed8c8177
corrected the constitutive law for j2-plasticity in Lp and its tangent and dotState.
2009-06-23 10:39:29 +00:00
Franz Roters
0873231e78
corrected NaN check in crystallite_updateState
...
return false in case NaN occurs in crystallite_updateState
2009-06-23 06:52:51 +00:00
Philip Eisenlohr
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
2009-06-19 07:09:39 +00:00
Philip Eisenlohr
b5b4b32961
move operation to "phenoPowerLaw"
2009-06-19 07:03:43 +00:00
Philip Eisenlohr
8e70963008
move operation into "Code" part of documentation
2009-06-19 07:02:57 +00:00
Christoph Kords
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.
...
added some output to the init subroutines of debug.f90, FEsolving.f90, lattice.f90, material.f90 and mesh.f90.
2009-06-18 14:28:02 +00:00
Christoph Kords
904ea78ac5
Tstar_v is now restored at different stages of cutbacking; we need this because preguess of state relies on consistent Tstar_v
2009-06-16 09:03:30 +00:00
Christoph Kords
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.
2009-06-15 13:11:21 +00:00
Yun Jo Ro
204e296ecd
modify document for ConstitutivePhenoPowerLaw
2009-06-15 08:09:49 +00:00
Christoph Kords
caf568eb89
adapted crystallite_stressAndItsTangent to do a pre-guess for the state before the actual state loop (with order stress integration, state update)
2009-06-10 15:08:33 +00:00
Yun Jo Ro
a8ff024d97
modify document for ConstitutivePhenoPowerLaw
2009-06-10 14:43:30 +00:00
Yun Jo Ro
118e68e881
modify document for ConstitutivePhenoPowerLaw
2009-06-09 13:48:46 +00:00
Yun Jo Ro
c17249c953
Add document for explaining how we sort the order of the slip/twin system for HCP in lattice.f90
2009-06-09 13:13:11 +00:00
Christoph Kords
3196049496
in crystallite_stressAndItsTangent we now update the state first (explicit integration) and then integrate the stress (implicit integration)
...
also inside the stress integration, it is now possible to define the frequency of the Jacobian update oin the LpLoop through iJacoLpResiduum. The frequency of the Jacobian update for the stiffness in the crystallite loop is controlled by the parameter iJacoStiffness.
also updated the corresponding stuctograms in the documentation
2009-06-09 11:05:29 +00:00
Yun Jo Ro
306cd95992
sort compressive and tensile twin-2 group in Hex structure and also correct the name for each twin type
2009-06-09 08:42:02 +00:00
Franz Roters
db731a4be3
Taylor homogenization is now called isostrain, changed material.config accordingly.
...
corrected error message in case of unknown homogenization
2009-06-08 13:28:00 +00:00
Yun Jo Ro
9d4521b7ec
modify document for ConstitutivePhenoPowerLaw
2009-06-08 09:29:54 +00:00
Yun Jo Ro
5a32d644fd
modify document for ConstitutivePhenoPowerLaw
2009-06-05 12:50:09 +00:00
Yun Jo Ro
4c6bb8a6fb
2009-06-05 12:48:30 +00:00
Christoph Kords
a9b38a7b40
adapted structogram for integrateStress and renamed the read.me file to readme.txt
2009-06-04 15:11:46 +00:00
Philip Eisenlohr
2c74dfc5c2
place to put any documentation regarding constitutive laws.
...
phenoPowerLaw as simple template...
please add your own subdirs accordingly
2009-06-04 13:47:40 +00:00
Yun Jo Ro
5adf84e8d8
modify document for ConstitutivePhenoPowerLaw
2009-06-04 12:37:16 +00:00
Christoph Kords
f4095e4127
added Nassi-Schneidermann-Diagrams for two more functions
2009-06-02 15:48:21 +00:00
Christoph Kords
7b638d7b94
sorry, forgot to declare all variables, promise to be more careful to always commit a running version of the subroutine ;-)
2009-06-02 09:55:45 +00:00