Commit Graph

4082 Commits

Author SHA1 Message Date
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
Yun Jo Ro b4092966fe 2009-03-26 07:00:10 +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 03bd466aec changed line ending 2009-03-16 17:39:52 +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 b19b8081a6 changed line ending to unix style 2009-03-06 11:16:09 +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