Commit Graph

6723 Commits

Author SHA1 Message Date
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
Franz Roters 1a5a68c87a added interface routine for Marc2008r1
this includes two files: concom_f90 and creeps_f90
so you either create links of that name to the new files concom2008r1 and creeps2008r1 or rename copies of those files accordingly
2008-10-29 15:28:10 +00:00
Franz Roters fa7cf61fd9 FEsolfing.f90: - added flag for symmetric solver
- added subroutine to detect symmetric solver

mesh.f90: - added subroutine call in mesh_init to detect symmetric solver during input file parsing

mpie_cpfem_marc2005r3.f90
mpie_cpfem_marc2007r1.f90
mpie_cpfem_marc2007r1_sequential.f90: - resorted include order of other source files
                                      - symmetrize d in case a symmetric solver is used

constitutive_pheno.f90: - included code to output shear rates and shear activity as post results

CPFEM_GIA8.f90
CPFEM_Taylor.f90
CPFEM_Taylor_sequential.f90: - symmetrize H_bar
                             - generalized reference to CPFEM_results in call of SingleCrystallite
2008-09-22 08:36:51 +00:00
Luc Hantcherli 3216976d62 Question: do we have to symmetrize H_bar? 2008-08-19 14:02:14 +00:00
Yun Jo Ro b3814020f3 added slip and twin systems for hcp lattice structure
added straightforward 48x48 interaction matrix
2008-08-06 14:48:11 +00:00
Philip Eisenlohr 7faaa8532c this provides adapted functionality to go with the sequential wrapper.
moved recording of temperature, ffn, and ffn1 into case of computation
mode 1,2 to be directly used in the subsequent (and now FE sequential)
call to MaterialPoint.
2008-08-01 07:58:59 +00:00
Philip Eisenlohr db69217ce1 removed computation mode 3 (collect) such that now every cycle will
compute with either mode=1 (new inc) or mode=2
2008-08-01 07:54:32 +00:00
Philip Eisenlohr 10358fc1f0 (hopefully) improved the cutbacking scheme somewhat. introduced a bunch
of debugging to the .out file
2008-08-01 07:47:01 +00:00
Franz Roters cf1c7ce82a replaced almost all calls of matmul by our own code, gives better performance in parallel computations 2008-07-23 12:49:40 +00:00
Philip Eisenlohr 6abff4b3ba changed hybridIA_population from being wrong integer type to correct real
type.
2008-07-16 16:35:13 +00:00
Philip Eisenlohr bdcf76cd74 fixed stupidity in hybridIA evaluation. added informative output during
ODF parsing.
2008-07-16 14:49:17 +00:00
Franz Roters 7a7a5d2ae3 messed up the parameter list of TimeIntegration with the last submit 2008-07-15 08:32:44 +00:00
Franz Roters 4773cf7d09 corrected determination and evaluation of updateJaco
introduced check whether FFN1 changed and force recollect in that case
2008-07-14 14:38:19 +00:00
Philip Eisenlohr 63d81b92b7 added math_mul33x33, 66x66, and 99x99 to avoid OMP_Critical "evilmatmul" 2008-07-08 19:38:22 +00:00
Franz Roters de9e35cdb8 modified cutback counting
restore Lpguess_old in case of iInner > nInner
2008-06-25 06:26:41 +00:00
Luc Hantcherli 7453a18b42 in lattice:
- add slip-twin interaction type matrix
2008-06-17 14:42:45 +00:00
Luc Hantcherli 434ac7f06e In dislo:
- model for stacking fault energy computation w.r.t temperature and orientation
- classical thermodynamical approach to twin nucleation (delta_E)
- rename some parameters

In crystallite:
- add msg='ok', even if components in consistent tangent are not re-computed
- allow for pseudo-convergence in the outer-loop (case of slow convergence of dislocation densities)

in mattex:
- add new parameters
2008-06-17 14:40:46 +00:00
Franz Roters 2472de77c2 OpenMP seems to work now
watch out for critical sections
write and matmul statements need to be marked accordingly!!!
2008-06-17 13:15:36 +00:00
Philip Eisenlohr 6921c57c7d rev 225 had a small syntax error... 2008-06-17 11:54:34 +00:00
Philip Eisenlohr 3676b453ca added two more element types (157 and 136) 2008-06-17 09:11:54 +00:00
Philip Eisenlohr 5834a2965d fixed bug in mesh_build_sharedElems which so far counted nodes in
distorted elements twice -- part of this problem has already been fixed
with revision 182 but then only in getNodeElemDimensions.
M    mesh.f90
2008-06-16 20:49:48 +00:00