Commit Graph

16393 Commits

Author SHA1 Message Date
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
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
Christoph Kords cb82ebedef added error message for non-positive taylor factor 2009-03-05 17:07:32 +00:00
Philip Eisenlohr 1420b773f0 needs renaming 2009-03-05 16:49:13 +00:00
Philip Eisenlohr ffa6af1853 update of sample material.config parts 2009-03-05 16:42:05 +00:00
Philip Eisenlohr 5a39497742 added sample material.config sections at module start 2009-03-05 16:36:31 +00:00
Luc Hantcherli 34aad5c945 Same with constitutive_dislobased.f90 2009-03-05 16:15:24 +00:00
Philip Eisenlohr ca1b9787b6 removed now obsolete files from former structure 2009-03-05 16:15:09 +00:00
Christoph Kords 122760601c not yet operational, though necessary to compile the rest 2009-03-05 16:12:13 +00:00
Luc Hantcherli 31f15ee146 Add constitutive_dislobased.f90, containing the current dislocation based model for dislocation glide.
Add corresponding changes in constitutive.f90 and IO.f90

Compilation errors in math.f90 corrected
2009-03-05 16:06:01 +00:00
Christoph Kords 599ba2453c added math_mul6x6 and math_mul3x3 to replacce more general calls to "dot_product" 2009-03-05 14:37:59 +00:00
Christoph Kords c9d43cccec now with j2 2009-03-05 14:36:01 +00:00
Franz Roters a7135ed83b removed some unused variables
correct some spelling errors in material.config
adopted marc2008 interface routine
removed sequential interface routines as they are not needed any more
2009-03-05 14:18:50 +00:00
Philip Eisenlohr 9aa89bd9b5 forgot to add those two (rather crucial) files .. ;-) 2009-03-04 14:01:36 +00:00
Christoph Kords e10964cc5d corrected error in <microstructure> template: "material" --> "phase" 2009-03-04 13:44:14 +00:00
Philip Eisenlohr 46de79d1f9 some tasks for Denny 2009-03-04 11:50:36 +00:00
Philip Eisenlohr 9351508ae1 overhaul to allow for multiphase/homogenization setup. 2009-03-04 11:48:54 +00:00
Eralp Demir 5a09badb71 corrected another small mistake, sorry for the inconvenience 2009-01-26 13:14:25 +00:00
Eralp Demir f3e5be4b54 small mistake corrected! 2009-01-26 13:07:10 +00:00
Eralp Demir 3625afff3b add symmetry operators and a subroutine for misorientation calculation 2009-01-26 12:58:58 +00:00
Philip Eisenlohr 5fc9816ff7 corrected Abaqus name to C3D8 2009-01-20 08:28:07 +00:00
Philip Eisenlohr bb32110cb0 further documentation of element type 5 (4 IP tetrahedron) and type 6 (6 IP pentahedron) 2009-01-20 08:10:56 +00:00
Philip Eisenlohr cbd7c279d4 changed Line Ending to Unix style (LF) 2009-01-19 19:10:58 +00:00
Philip Eisenlohr abb2e3ef30 introduction of logical "wantsConstitutiveResults" 2009-01-19 19:10:02 +00:00
Philip Eisenlohr 82e65648a2 major debugging and streamlining of the calculation of ipVolumes and
ipAreas
2009-01-19 18:42:31 +00:00
Philip Eisenlohr e3e42fdf44 shell script taking two arguments and setting CPFEM_xxx and
constitutive_yyy symbolic links from those
2009-01-19 13:37:37 +00:00
Philip Eisenlohr 4b67529aaf added full description of element type 2 2009-01-16 17:36:37 +00:00
Philip Eisenlohr 8788cc4f8a description of element type 2 2009-01-16 17:35:14 +00:00
Eralp Demir 483896218f added functions
-- mesh_build_subNodeCoords
-- mesh_build_ipVolimes
-- mesh_build_ipAreas

new globals mesh_ipVolume and mesh_ipArea / mesh_ipAreaNormal contain the volume per IP, area of the interface to aach neighboring IP and the outward pointing unit normal vector of these interfaces
2009-01-16 15:29:57 +00:00
Eralp Demir de04c4cccd added functions
-- math_volTetrahedron: calc the volume of a tetrahedron from four vertices
-- math_vectorproduct: return vector product of two dim(3) vectors

rendered most of the math lib to be "pure" functions
2009-01-16 15:27:13 +00:00
Eralp Demir 1529d6115d rendered most of the functions to be "pure" 2009-01-16 15:25:37 +00:00
Philip Eisenlohr 9d14dbde2b started to document things with the description of the "all-purpose" element 7 (or 3DC8) as interpreted by mesh.f90 2009-01-16 08:29:03 +00:00
Franz Roters 7b076185b2 This is a major update.
I restructured the subroutine in view of non local constitutive models.
ATTENTION OpenMP parallelization does not work with this version, needs some more work.
ATTENTION CPFEM_GIA8.f90 is not addopted to the new scheme and will NOT work, should actually be renamed CPFEM_RGC.f90 ;-)
I removed crystallite.f90, the routines are now in CPFEM_*.f90.
I removed the marc2005 interface routine since it is outdated.
2008-11-28 07:39:39 +00:00
Franz Roters e8cd7cef94 fixed a small bug: text about finding symmetric solver was always written to out-file 2008-11-27 15:11:44 +00:00
Franz Roters 824cec5ca1 correct error in case of outdated FFN1: element number is cp_en not e 2008-11-11 14:57:20 +00:00
Philip Eisenlohr db245e8935 a small text file to remember ideas for code improvement 2008-11-10 08:41:32 +00:00
Luc Hantcherli 3ec8ec6b15 New version of constitutive_dislo.f90, where some dangerous spots were reformulated.
PROBLEM: [001] in tension does not have a symmetrical change of shape, which is definitely not good at all. I checked my implementation and I really do not see any errors. Could someone check it out? Can it comes from the numerical integration?
2008-11-07 16:34:29 +00:00
Luc Hantcherli 6e13ed0566 I will proceed step by step.
constitutive_dislo.f90 contains a modified version of Anxin's dislocation based model. As far as tested it, I consider it now error free. Mechanical twinning is NOT implemented YET.
mattex.mpie file is modified in accordance with the changes in constitutive_dislo.f90
It would be a great help if someone else can check the implementation. I may have overseen something.
2008-11-06 13:09:12 +00:00
Franz Roters 9692f2f3d3 added mpie_cpfem_marc2008r1_sequential.f90 2008-10-30 11:18:10 +00:00