Commit Graph

  • 9d14dbde2b started to document things with the description of the "all-purpose" element 7 (or 3DC8) as interpreted by mesh.f90 Philip Eisenlohr 2009-01-16 08:29:03 +0000
  • 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. Franz Roters 2008-11-28 07:39:39 +0000
  • e8cd7cef94 fixed a small bug: text about finding symmetric solver was always written to out-file Franz Roters 2008-11-27 15:11:44 +0000
  • 824cec5ca1 correct error in case of outdated FFN1: element number is cp_en not e Franz Roters 2008-11-11 14:57:20 +0000
  • db245e8935 a small text file to remember ideas for code improvement Philip Eisenlohr 2008-11-10 08:41:32 +0000
  • 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? Luc Hantcherli 2008-11-07 16:34:29 +0000
  • 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. Luc Hantcherli 2008-11-06 13:09:12 +0000
  • 9692f2f3d3 added mpie_cpfem_marc2008r1_sequential.f90 Franz Roters 2008-10-30 11:18:10 +0000
  • 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 Franz Roters 2008-10-29 15:28:10 +0000
  • fa7cf61fd9 FEsolfing.f90: - added flag for symmetric solver - added subroutine to detect symmetric solver Franz Roters 2008-09-22 08:36:51 +0000
  • 3216976d62 Question: do we have to symmetrize H_bar? Luc Hantcherli 2008-08-19 14:02:14 +0000
  • b3814020f3 added slip and twin systems for hcp lattice structure added straightforward 48x48 interaction matrix Yun Jo Ro 2008-08-06 14:48:11 +0000
  • 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. Philip Eisenlohr 2008-08-01 07:58:59 +0000
  • db69217ce1 removed computation mode 3 (collect) such that now every cycle will compute with either mode=1 (new inc) or mode=2 Philip Eisenlohr 2008-08-01 07:54:32 +0000
  • 10358fc1f0 (hopefully) improved the cutbacking scheme somewhat. introduced a bunch of debugging to the .out file Philip Eisenlohr 2008-08-01 07:47:01 +0000
  • cf1c7ce82a replaced almost all calls of matmul by our own code, gives better performance in parallel computations Franz Roters 2008-07-23 12:49:40 +0000
  • 6abff4b3ba changed hybridIA_population from being wrong integer type to correct real type. Philip Eisenlohr 2008-07-16 16:35:13 +0000
  • bdcf76cd74 fixed stupidity in hybridIA evaluation. added informative output during ODF parsing. Philip Eisenlohr 2008-07-16 14:49:17 +0000
  • 7a7a5d2ae3 messed up the parameter list of TimeIntegration with the last submit Franz Roters 2008-07-15 08:32:44 +0000
  • 4773cf7d09 corrected determination and evaluation of updateJaco introduced check whether FFN1 changed and force recollect in that case Franz Roters 2008-07-14 14:38:19 +0000
  • 63d81b92b7 added math_mul33x33, 66x66, and 99x99 to avoid OMP_Critical "evilmatmul" Philip Eisenlohr 2008-07-08 19:38:22 +0000
  • de9e35cdb8 modified cutback counting restore Lpguess_old in case of iInner > nInner Franz Roters 2008-06-25 06:26:41 +0000
  • 7453a18b42 in lattice: - add slip-twin interaction type matrix Luc Hantcherli 2008-06-17 14:42:45 +0000
  • 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 Luc Hantcherli 2008-06-17 14:40:46 +0000
  • 2472de77c2 OpenMP seems to work now watch out for critical sections write and matmul statements need to be marked accordingly!!! Franz Roters 2008-06-17 13:15:36 +0000
  • 6921c57c7d rev 225 had a small syntax error... Philip Eisenlohr 2008-06-17 11:54:34 +0000
  • 3676b453ca added two more element types (157 and 136) Philip Eisenlohr 2008-06-17 09:11:54 +0000
  • 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 Philip Eisenlohr 2008-06-16 20:49:48 +0000
  • d2335e4a0b Denny Tjahjanto 2008-06-10 10:23:32 +0000
  • 60f3571266 added code lines for OpenMP, however it does not work yet! Franz Roters 2008-05-26 13:11:25 +0000
  • 51e596c81d small polishing to be compatible with NAG compiler (basis for automatic differentiation) Philip Eisenlohr 2008-04-29 12:18:17 +0000
  • 9626c25bfb crystallite: restoration of last known good Lp_guess CPFEM_Taylor: exit whenever cutback limit exceeded Luc Hantcherli 2008-04-28 16:04:08 +0000
  • cdb2dd8808 crystallite: added full check whether residuum in inner loop is NaN SingleCrystallite now takes/stores Lp guess Philip Eisenlohr 2008-04-16 17:00:28 +0000
  • fa2d6b9b6d Make use of average consistent tangent dPdF_bar to fill eventual holes in the individual dPdF of one orientation (during consistent tangent calculation) Luc Hantcherli 2008-04-11 13:04:25 +0000
  • 664a29c183 Subscript Crystal_ changed to Lattice_ according to the module Lattice Luc Hantcherli 2008-04-10 12:13:37 +0000
  • 9d2ce61698 A revision for CPFEM_Taylor.f90. Denny Tjahjanto 2008-04-10 11:22:17 +0000
  • d2312e81ff deleted CFPFEM.f90 and crystal.f90 from repository adopted mattex.mpie to use lattice_structure instead of crystal_structure Franz Roters 2008-04-09 15:35:16 +0000
  • 12dfbaf6b4 # separated CPFEM_Crystallite from CPFEM. Now in separate module crystallite.f90 as "SingleCrystallite" # improved SingleCrystallite to advance by true cutbacking (instead of improving guess and integrating always from t_0) # module "crystal" renamed to "lattice" together with its prefix for variables # extension of "computationMode" to deal with cutbacks (CPFEM_general). # cutback and new inc detection for MARC is based on common block variable cptim (and inc), not incsub anymore! # generalized GrainInterAction as new homogenization scheme Denny Tjahjanto 2008-04-07 14:54:29 +0000
  • c063ce5bc1 add a kronecker delta function Denny Tjahjanto 2008-03-27 11:54:34 +0000
  • 5462d2c754 just a correction for a little typo :P Denny Tjahjanto 2008-03-27 11:53:53 +0000
  • f3eda7cdaa include a new function epsilon_ijk, i.e., the permutation function (used to compute the cross product of two tensors) Denny Tjahjanto 2008-03-26 13:35:01 +0000
  • 92ee2cbf7d to include 2 material parameters, i.e., grain size and the length of burgers vector used in GIA model Denny Tjahjanto 2008-03-26 09:33:12 +0000
  • cc7e13d351 New update of constitutive_dislo.f60 considering the changes in mesh.f90 (see map_ElemType) Parameters assignation, hardening matrices and constitutive equations are PROPERLY implemented (tested) Comparison to my own old subroutine: on going Luc Hantcherli 2008-03-25 14:19:22 +0000
  • 7faf4093a5 reflect change of mesh_element(2,e) in mesh.f90 -change of line, and those below, will be ignored-- Philip Eisenlohr 2008-03-25 13:59:22 +0000
  • 3640a48095 changed element type mapping to be performed right after input file parsing. mesh_element(2,e) now contains INTERNAL type index -- no subsequent mapping necessary anymore! Philip Eisenlohr 2008-03-25 12:52:27 +0000
  • 6d721dc16c Modifications of cpfem_marc2005r3 according the ones in cpfem_marc2007r1 Luc Hantcherli 2008-03-14 21:56:46 +0000
  • 1267ce78b6 CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.) Philip Eisenlohr 2008-03-14 21:32:57 +0000
  • 7975d62586 CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.) Philip Eisenlohr 2008-03-14 21:29:31 +0000
  • 0a08d9eff0 Corrected errors from the last release of constitutive_dislo Changed parameter designation in mattex file Luc Hantcherli 2008-03-14 16:23:08 +0000
  • 6950eee59b altered stress and strain "language" to first PK and deformation gradient. Denny Tjahjanto 2008-03-14 14:19:10 +0000
  • 6f06133b90 removed some unused variables Franz Roters 2008-03-14 08:23:43 +0000
  • 778cbcd3ca changed file read formats from explicit types to more general (and not harmful) '*' Denny Tjahjanto 2008-03-12 13:53:00 +0000
  • 5641cfd477 changed precision assignment to intrinsic (and standardized) "selected_*_kind" Philip Eisenlohr 2008-03-07 15:21:46 +0000
  • c74e071f5e updated constitutive_xxx to reflect new 9x9 dLp_dTstar - used to be 3x3x3x3 Philip Eisenlohr 2008-02-21 12:50:37 +0000
  • 9c1e8a7944 added acceleration capability after time-step cut backing Philip Eisenlohr 2008-02-21 10:33:34 +0000
  • 9707fd0f8f full update, i.e. my development snapshot Philip Eisenlohr 2008-02-19 12:58:46 +0000
  • f4edf4bd0c 2007 shortened stress_recovery = ... added (crucial) deMandeling for return values of CPFEM_general Philip Eisenlohr 2008-02-18 18:49:06 +0000
  • 2b567ad20d CPFEM changed consistent tangent scheme to D.Tjahjanto Diss version fixed nasty bug in decision when to calculate new values (CPFEM_general) Philip Eisenlohr 2008-02-18 18:11:24 +0000
  • 199c65cdba changed pert_e to new perturbation value pert_Fg Denny Tjahjanto 2008-02-18 10:46:18 +0000
  • 1106bb27eb introduced logical "debugger" Denny Tjahjanto 2008-02-18 10:45:46 +0000
  • e07d38fdec ooops. missed nInner and nOuter as loop counters... Denny Tjahjanto 2008-02-18 10:41:37 +0000
  • 9c9cdeb58b corrected nStress and nState to reflect more general nInner and nOuter (now vompatible with changes in prec.f90) Denny Tjahjanto 2008-02-18 10:39:39 +0000
  • cfaa0e696d exchanged TINY (from prec) with intrinsic "tiny" function (Fortran90) Philip Eisenlohr 2008-02-15 12:42:27 +0000
  • a41a4a75ef CPFEM: time integration scheme changed to start from guessing of Lp. Newton-Raphson algorithm now accelerates as long as residuum improves. Philip Eisenlohr 2008-02-15 10:04:49 +0000
  • f39c68dded CPFEM_stressCrystallite must be called with t not dt in CPFEM_stressIP Franz Roters 2008-02-08 08:25:19 +0000
  • 6a3a9e8f98 Added more detailed reporting Hongmei Li 2008-01-11 18:32:02 +0000
  • 720d8da82b Corrected syntax Luc Hantcherli 2008-01-10 19:00:45 +0000
  • b8b171c95b Introduce debug module, contains distributions of nCutback, nStressLoop, and nStateLoop Luc Hantcherli 2008-01-10 18:53:57 +0000
  • 3a7ec38a3d Delete unused physical constants from pheno Rename Boltzmann Constant Kb->kB Luc Hantcherli 2008-01-10 18:51:58 +0000
  • 1233d01cd5 mesh_maxNsharedElems is now correct. Added check for multiple identical nodes in an element, e.g. triangular quad-element and the like. Philip Eisenlohr 2008-01-10 17:12:33 +0000
  • fe70a82d6d !!!!! IMPORTANT !!!!! All subroutines were committed at once: - constitutive_pheno works - constitutive_dislo without twinning also works Luc Hantcherli 2007-12-14 13:36:04 +0000
  • 3ef451824c just a test Luc Hantcherli 2007-12-14 13:28:41 +0000
  • 90195fe1c8 Commit dislo Luc Hantcherli 2007-12-11 14:30:21 +0000
  • ab0bde576e Commit of IO Luc Hantcherli 2007-12-11 14:23:33 +0000
  • 32fc32865f Commit of mesh Luc Hantcherli 2007-12-11 14:23:21 +0000
  • 08a1c38c73 Finalizing the implementation of deformation twinning in constitutive_dislo.f90 Assume to be ready to use Luc Hantcherli 2007-12-10 12:55:43 +0000
  • d0f6c81d66 constitutive_Microstructure() computes now twin "plop" volumes from the given microstructure. Mattex.dat has 2 parameters for deformation twinning: - average grain size - average twin stack thickness Luc Hantcherli 2007-12-07 15:12:43 +0000
  • 44d6235777 Constitutive_dislo.f90 contains now a modified version of HomogenisedC that accounts for twin reoriented regions. ONLY constitutive_dislo was changed, NOT constitutive.pheno Luc Hantcherli 2007-12-07 14:17:42 +0000
  • b69aa2b112 Construction of the elasticity matrices for rotated/twin regions. Definition of constitutive_Ntwin and constitutive_MaxNtwin Luc Hantcherli 2007-12-07 13:59:15 +0000
  • e18f304c5a Crystal.f90 contains now the twin systems (at least for fcc, not implemented for bcc and hcp). Schmid and transformations matrices for twin systemss are defined and constructed. Twin systems are specified using the type of variables as slip systems: "MaxNTwin", "MaxMaxNtwin" ... Luc Hantcherli 2007-12-07 12:36:17 +0000
  • f63f6949cc Corrected undefined "state" in HomogenisedC Added "forall" loops in C construction Luc Hantcherli 2007-12-07 10:54:19 +0000
  • d615406722 Corrected HomogenisedC: 1- No dependence on material volume fractions!!! 2- Can now depend on state of microstructure Luc Hantcherli 2007-12-07 10:37:06 +0000
  • 975c113ae8 Homogenised C per grain can now depend on microstructure (important for deformation twinning) Luc Hantcherli 2007-12-07 10:34:21 +0000
  • 3c0d51542e Changed to SI conformity. Mattex file changed bg->burgers Luc Hantcherli 2007-11-15 17:04:30 +0000
  • ad3ccf22c2 clarified tolerances in prec (abs vs rel) all units conform to SI now !! (especially MPa --> Pa) Philip Eisenlohr 2007-11-15 11:21:37 +0000
  • e86a288a21 CPFEM_dummy_stress now user-given constant in _init (SI units !!) Philip Eisenlohr 2007-11-15 11:12:58 +0000
  • 896c37ede2 changed structure of CPFEM_general to allow for non local algorithm and parrallelisation Franz Roters 2007-11-08 08:26:02 +0000
  • 8619a92009 made sumVolfrac real also in constitutive_dislo.f90 Franz Roters 2007-11-07 15:02:54 +0000
  • 56e82e0fe2 changed sumVolfrac to real Franz Roters 2007-11-06 11:21:46 +0000
  • d94f6c80de deleted unused variable constitutive_results constitutive_maxNresults gets now calculated as maxval instead of hardwired Philip Eisenlohr 2007-11-06 10:21:03 +0000
  • 4dbf6ac3fa added volume fraction of orientation to output array (now CPFEM_Nresults = 4) Philip Eisenlohr 2007-11-06 08:30:04 +0000
  • d2471ab479 corrected formating error in output from mesh Franz Roters 2007-10-26 13:39:54 +0000
  • 909ebf2cd0 removed unused variables corrected allocation error in constitutive_dislo.f90 Franz Roters 2007-10-24 15:24:49 +0000
  • 20693b9d38 hybrid IA sampling now respects sample symmetry Philip Eisenlohr 2007-10-24 15:13:10 +0000
  • bf35fa4868 former revision of mesh.f90 caused segmentation fault. probably due to string assignment across memory. fixed now... William Counts 2007-10-24 11:25:58 +0000
  • f0617b74f6 added output of mesh statistics at end of input file parsing Philip Eisenlohr 2007-10-24 09:00:42 +0000
  • 67a5711ed1 added awareness of set definitions (new function build_elemSetMapping) added new style tables for initial conditions William Counts 2007-10-23 13:09:46 +0000
  • ca0e2e55df fixed broken IO_continousIntValues to be aware of element (or other) sets -This line, and those below, will be ignored-- William Counts 2007-10-23 13:08:27 +0000
  • 9043b4374e Corrected the treatment of interaction coefficients in constitutive_dislo Luc Hantcherli 2007-10-22 15:32:45 +0000
  • 57ed908129 removed unused variables made hardening matrix finally (?) work Franz Roters 2007-10-22 14:53:10 +0000