Philip Eisenlohr
e507bd0bcd
data structure of UserOutput (negative values of post result indices...)
...
please note the introduction of size identifiers (N, n_1,...n_m) which allow to more easily parse IPs which happen to have different constitutive laws for different constituent grains.
2009-05-07 17:41:56 +00:00
Philip Eisenlohr
4f705f524f
major restructuring of code.
...
homogenization as well as constitutive are now free to choose. the runtime got somewhat longer (25% on simple tests) compared to a hardcoded isostrain homogenization. this might be a point of further optimization at a later stage...
please use homogenization_isostrain.f90 as starting point / example for future developments of homog-schemes.
the homogenization scheme now can additionally output certain results. hence, the userdata structure at each integration point now looks like this:
- sizeHomogPostResults
- block of that size containing homogPostResults
then for each grain:
- sizeGrainPostResults
- block of that size containing crystallitePostResults, which consist of:
+ phaseID
+ volFrac
+ Eulers (3)
+ any constitutive post results requested
2009-05-07 16:27:36 +00:00
Christoph Kords
449b791463
constitutive law for j2 pasticity again corrected, dotState didn't need the factor of sqrt(1.5)
2009-05-06 08:07:18 +00:00
Christoph Kords
4c2eecf280
corrected constitutive law for j2-plasticity, a factor of sqrt(1.5) was missing in the calculation of Lp, dLp/dT and dotState
2009-05-05 14:51:02 +00:00
Christoph Kords
d052a69574
corrected spelling in subroutine "material_parseTexture"
2009-04-28 09:45:52 +00:00
Franz Roters
5063842f1b
added lclass to dimension statement
2009-04-22 14:21:18 +00:00
Christoph Kords
284a339f38
corrected type of 'material_volFrac' to real; 'volFracOfGrain' is not normalized with the mesh_volume of the current microstructure anymore, since this affects the calculation with domain decomposition but is not needed for any other calculations
2009-04-08 06:45:17 +00:00
Christoph Kords
e3b9a35fc3
now with full functionality for element type 6 (6-node pentaeder) and with element data structure partly initialized as function variables and not parameters
2009-04-06 13:25:19 +00:00
Philip Eisenlohr
115ac229d4
now with center of gravity subnode (21)
2009-04-03 15:00:05 +00:00
Franz Roters
d1c1c34faa
corrected code for non parallel calculation
2009-04-03 12:41:01 +00:00
Philip Eisenlohr
09148ef862
unused variable "counter" deleted -- now compiles without error...
2009-04-03 11:07:00 +00:00
Philip Eisenlohr
7ecd0355d0
added array "phase_localConstitution" which is keeping track of phases which emlpoy local or non-local constitutive laws. Non-local constitution is indicated by a key "/nonlocal/" in the phase definition.
2009-04-03 10:34:17 +00:00
Philip Eisenlohr
72f012992d
new function IO_spotTagInPart
2009-04-03 10:30:18 +00:00
Christoph Kords
f8507cfefd
suppressed output of subnode coordinates to *.out file, but still available in 'mesh_tell_statistics'
2009-04-03 09:46:10 +00:00
Christoph Kords
b57a128a88
added subroutine 'IO_skipChunks' that is needed by mesh.f90
2009-04-03 07:04:31 +00:00
Christoph Kords
0e84b07dc7
Simplified 'mesh_get_nodeElemDimensions', mesh_build_elements' and 'mesh_build_sharedElems'.
...
Corrected entries in 'FE_subNodeOnIPFace'
2009-04-03 07:03:25 +00:00
Christoph Kords
b2e8510b6d
corrected entries in 'FE_subNodeOnIPFace'
2009-04-02 14:17:29 +00:00
Christoph Kords
917a8523af
corrected misspelling of function 'IO_continousIntValues'
2009-04-02 13:02:25 +00:00
Christoph Kords
a024826c6c
Element type 21 (20-node hexagon with triquadric interpolation functions) is now included.
...
'mesh_get_nodeElemDimensions', mesh_build_elements', 'mesh_build_sharedElems' now allow for more than 14 nodes per elements which are then written in more than one line in the 'connectivity' field of the input file
2009-04-02 13:00:51 +00:00
Philip Eisenlohr
e5af3874b2
added homogneization_typeInstance and homogenization_Noutput for future use
2009-03-31 09:24:12 +00:00
Philip Eisenlohr
d96c0a6495
added subroutine to throw warnings (instead of terminal errors)
2009-03-31 09:21:57 +00:00
Franz Roters
318fa8008c
corrected check of number of textures
2009-03-31 09:09:24 +00:00
Franz Roters
51cc6714d0
changed pLongInt to 8
2009-03-31 09:07:30 +00:00
Franz Roters
a5a2fdb306
introduced pLongInt for Long Integers
2009-03-31 08:51:14 +00:00
Franz Roters
482b256d84
changed output format for more than 999 microstructures
2009-03-31 08:49:23 +00:00
Philip Eisenlohr
a9824e473d
added function (math_inv3x3) for inverting a 3x3 matrix by Cramer's methods
2009-03-31 07:31:38 +00:00
Christoph Kords
78eed33d6d
changed "tickrate" to type integer(8) for microsecond resolution
2009-03-30 14:54:05 +00:00
Philip Eisenlohr
ee3fdcb2c4
extended linking between nodes and IPs to account for more complex
...
situations in quadratic shape function elements.
2009-03-30 14:50:19 +00:00
Philip Eisenlohr
5c1889e6cd
simplified constitutive_init call
2009-03-26 08:43:31 +00:00
Franz Roters
2fa3f1de4d
renamed makeMe to makeMe.py to make file association of Windows work
2009-03-25 08:39:53 +00:00
Philip Eisenlohr
31d3dae279
FE (Marc) interface is now general. Requires to execute "makeMe" in order to derive the specific interface routines for each Marc version (e.g., 2007r1 and 2008r1) from mpie_cpfem_marc.f90. Pls do not commit explicit mpie_cpfem_marcXXX.f90 anymore and ignore them for SVN control.
2009-03-23 16:08:17 +00:00
Franz Roters
8441ff5372
adopted sequence of including lattice.f90 and material.f90
2009-03-23 09:52:23 +00:00
Philip Eisenlohr
ac42e24252
(re)introduced integer(8) types for clock_tick measurements
2009-03-20 15:19:58 +00:00
Philip Eisenlohr
a57ce74e15
exchanged new lattice_structure ("fcc", "bcc", "hexagonal")
2009-03-20 15:19:21 +00:00
Luc Hantcherli
e5255f4384
Update in order to include the c/a-ratio as additional information about the structure (modules lattice, constitutive_ changed accordingly)
...
Correction of the expression for the velocity of mobile dislocation densities in constitutive_dislobased.f90
Correction of some Fortran syntax, incompatible with Compaq Visual Fortran (e.g. integer(8))
2009-03-20 14:34:24 +00:00
Philip Eisenlohr
121723ffa5
corrected (?) RodrigToR
2009-03-18 14:05:27 +00:00
Philip Eisenlohr
9007fbc73c
added math_tensorproduct
2009-03-17 15:13:17 +00:00
Philip Eisenlohr
eed32fe588
added timing facilities for LpAndItsTangent
2009-03-16 17:38:33 +00:00
Philip Eisenlohr
a8db780af4
remembered the Bal98 (MIT) PhD thesis which might contain an analytical way to get the Jacobian matrix -- as opposed to numerical perturbation...
2009-03-12 07:06:55 +00:00
Luc Hantcherli
919a6cd32a
Change in constitutive_dislobased.f90: now ..._sizeState is 8 times Nslip long, in agreement with the required information about the microstructure.
2009-03-09 22:41:53 +00:00
Philip Eisenlohr
dde0af0b37
sanity check for phase and texture index out of bounds was errorenous...
2009-03-09 13:58:59 +00:00
Christoph Kords
c0899b3097
minor changes to comments
2009-03-06 16:53:31 +00:00
Christoph Kords
a8ab5be60c
no calculation of Lp and dLp_dTstar in first loop
2009-03-06 13:30:22 +00:00
Christoph Kords
f845b9dc87
corrected indices for calculation of dLp_dTstar
2009-03-06 12:29:47 +00:00
Philip Eisenlohr
41d9e2d805
sorry, compilation problem corrected...
2009-03-06 11:47:26 +00:00
Philip Eisenlohr
f84d2abedc
corrected mistake in Lp tangent calculation (Christoph, pls check)
2009-03-06 11:28:48 +00:00
Philip Eisenlohr
b907961032
added examples for j2
2009-03-06 11:15:55 +00:00
Christoph Kords
ae0eafdb90
added calculation for tangent of Lp
2009-03-06 11:07:15 +00:00
Philip Eisenlohr
d12904a412
constitutitions of "j2" and "dislobased" are referenced in constitutive and included to the Marc interface
2009-03-06 10:13:08 +00:00
Christoph Kords
ea8ce09764
changed line endings to unix style
2009-03-06 10:02:36 +00:00