9d14dbde2bstarted 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
7b076185b2This 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
e8cd7cef94fixed a small bug: text about finding symmetric solver was always written to out-file
Franz Roters
2008-11-27 15:11:44 +0000
824cec5ca1correct error in case of outdated FFN1: element number is cp_en not e
Franz Roters
2008-11-11 14:57:20 +0000
db245e8935a small text file to remember ideas for code improvement
Philip Eisenlohr
2008-11-10 08:41:32 +0000
3ec8ec6b15New 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
6e13ed0566I 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
9692f2f3d3added mpie_cpfem_marc2008r1_sequential.f90
Franz Roters
2008-10-30 11:18:10 +0000
1a5a68c87aadded 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
fa7cf61fd9FEsolfing.f90: - added flag for symmetric solver - added subroutine to detect symmetric solver
Franz Roters
2008-09-22 08:36:51 +0000
3216976d62Question: do we have to symmetrize H_bar?
Luc Hantcherli
2008-08-19 14:02:14 +0000
b3814020f3added slip and twin systems for hcp lattice structure added straightforward 48x48 interaction matrix
Yun Jo Ro
2008-08-06 14:48:11 +0000
7faaa8532cthis 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
db69217ce1removed 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
cf1c7ce82areplaced almost all calls of matmul by our own code, gives better performance in parallel computations
Franz Roters
2008-07-23 12:49:40 +0000
6abff4b3bachanged hybridIA_population from being wrong integer type to correct real type.
Philip Eisenlohr
2008-07-16 16:35:13 +0000
bdcf76cd74fixed stupidity in hybridIA evaluation. added informative output during ODF parsing.
Philip Eisenlohr
2008-07-16 14:49:17 +0000
7a7a5d2ae3messed up the parameter list of TimeIntegration with the last submit
Franz Roters
2008-07-15 08:32:44 +0000
4773cf7d09corrected 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
63d81b92b7added math_mul33x33, 66x66, and 99x99 to avoid OMP_Critical "evilmatmul"
Philip Eisenlohr
2008-07-08 19:38:22 +0000
de9e35cdb8modified cutback counting restore Lpguess_old in case of iInner > nInner
Franz Roters
2008-06-25 06:26:41 +0000
434ac7f06eIn 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
2472de77c2OpenMP 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
6921c57c7drev 225 had a small syntax error...
Philip Eisenlohr
2008-06-17 11:54:34 +0000
3676b453caadded two more element types (157 and 136)
Philip Eisenlohr
2008-06-17 09:11:54 +0000
5834a2965dfixed 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
60f3571266added code lines for OpenMP, however it does not work yet!
Franz Roters
2008-05-26 13:11:25 +0000
51e596c81dsmall polishing to be compatible with NAG compiler (basis for automatic differentiation)
Philip Eisenlohr
2008-04-29 12:18:17 +0000
9626c25bfbcrystallite: restoration of last known good Lp_guess CPFEM_Taylor: exit whenever cutback limit exceeded
Luc Hantcherli
2008-04-28 16:04:08 +0000
cdb2dd8808crystallite: 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
fa2d6b9b6dMake 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
664a29c183Subscript Crystal_ changed to Lattice_ according to the module Lattice
Luc Hantcherli
2008-04-10 12:13:37 +0000
9d2ce61698A revision for CPFEM_Taylor.f90.
Denny Tjahjanto
2008-04-10 11:22:17 +0000
d2312e81ffdeleted 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
c063ce5bc1add a kronecker delta function
Denny Tjahjanto
2008-03-27 11:54:34 +0000
5462d2c754just a correction for a little typo :P
Denny Tjahjanto
2008-03-27 11:53:53 +0000
f3eda7cdaainclude 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
92ee2cbf7dto 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
cc7e13d351New 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
7faf4093a5reflect 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
3640a48095changed 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
6d721dc16cModifications of cpfem_marc2005r3 according the ones in cpfem_marc2007r1
Luc Hantcherli
2008-03-14 21:56:46 +0000
1267ce78b6CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.)
Philip Eisenlohr
2008-03-14 21:32:57 +0000
7975d62586CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.)
Philip Eisenlohr
2008-03-14 21:29:31 +0000
0a08d9eff0Corrected errors from the last release of constitutive_dislo Changed parameter designation in mattex file
Luc Hantcherli
2008-03-14 16:23:08 +0000
6950eee59baltered stress and strain "language" to first PK and deformation gradient.
Denny Tjahjanto
2008-03-14 14:19:10 +0000
6f06133b90removed some unused variables
Franz Roters
2008-03-14 08:23:43 +0000
778cbcd3cachanged file read formats from explicit types to more general (and not harmful) '*'
Denny Tjahjanto
2008-03-12 13:53:00 +0000
5641cfd477changed precision assignment to intrinsic (and standardized) "selected_*_kind"
Philip Eisenlohr
2008-03-07 15:21:46 +0000
c74e071f5eupdated constitutive_xxx to reflect new 9x9 dLp_dTstar - used to be 3x3x3x3
Philip Eisenlohr
2008-02-21 12:50:37 +0000
9c1e8a7944added acceleration capability after time-step cut backing
Philip Eisenlohr
2008-02-21 10:33:34 +0000
9707fd0f8ffull update, i.e. my development snapshot
Philip Eisenlohr
2008-02-19 12:58:46 +0000
f4edf4bd0c2007 shortened stress_recovery = ... added (crucial) deMandeling for return values of CPFEM_general
Philip Eisenlohr
2008-02-18 18:49:06 +0000
2b567ad20dCPFEM 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
199c65cdbachanged pert_e to new perturbation value pert_Fg
Denny Tjahjanto
2008-02-18 10:46:18 +0000
e07d38fdecooops. missed nInner and nOuter as loop counters...
Denny Tjahjanto
2008-02-18 10:41:37 +0000
9c9cdeb58bcorrected 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
cfaa0e696dexchanged TINY (from prec) with intrinsic "tiny" function (Fortran90)
Philip Eisenlohr
2008-02-15 12:42:27 +0000
a41a4a75efCPFEM: 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
f39c68ddedCPFEM_stressCrystallite must be called with t not dt in CPFEM_stressIP
Franz Roters
2008-02-08 08:25:19 +0000
6a3a9e8f98Added more detailed reporting
Hongmei Li
2008-01-11 18:32:02 +0000
1233d01cd5mesh_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
3ef451824cjust a test
Luc Hantcherli
2007-12-14 13:28:41 +0000
ab0bde576eCommit of IO
Luc Hantcherli
2007-12-11 14:23:33 +0000
32fc32865fCommit of mesh
Luc Hantcherli
2007-12-11 14:23:21 +0000
08a1c38c73Finalizing the implementation of deformation twinning in constitutive_dislo.f90 Assume to be ready to use
Luc Hantcherli
2007-12-10 12:55:43 +0000
d0f6c81d66constitutive_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
44d6235777Constitutive_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
b69aa2b112Construction of the elasticity matrices for rotated/twin regions. Definition of constitutive_Ntwin and constitutive_MaxNtwin
Luc Hantcherli
2007-12-07 13:59:15 +0000
e18f304c5aCrystal.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
f63f6949ccCorrected undefined "state" in HomogenisedC Added "forall" loops in C construction
Luc Hantcherli
2007-12-07 10:54:19 +0000
d615406722Corrected 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
975c113ae8Homogenised C per grain can now depend on microstructure (important for deformation twinning)
Luc Hantcherli
2007-12-07 10:34:21 +0000
3c0d51542eChanged to SI conformity. Mattex file changed bg->burgers
Luc Hantcherli
2007-11-15 17:04:30 +0000
ad3ccf22c2clarified tolerances in prec (abs vs rel) all units conform to SI now !! (especially MPa --> Pa)
Philip Eisenlohr
2007-11-15 11:21:37 +0000
e86a288a21CPFEM_dummy_stress now user-given constant in _init (SI units !!)
Philip Eisenlohr
2007-11-15 11:12:58 +0000
896c37ede2changed structure of CPFEM_general to allow for non local algorithm and parrallelisation
Franz Roters
2007-11-08 08:26:02 +0000
8619a92009made sumVolfrac real also in constitutive_dislo.f90
Franz Roters
2007-11-07 15:02:54 +0000
56e82e0fe2changed sumVolfrac to real
Franz Roters
2007-11-06 11:21:46 +0000
d94f6c80dedeleted unused variable constitutive_results constitutive_maxNresults gets now calculated as maxval instead of hardwired
Philip Eisenlohr
2007-11-06 10:21:03 +0000
4dbf6ac3faadded volume fraction of orientation to output array (now CPFEM_Nresults = 4)
Philip Eisenlohr
2007-11-06 08:30:04 +0000
d2471ab479corrected formating error in output from mesh
Franz Roters
2007-10-26 13:39:54 +0000
909ebf2cd0removed unused variables corrected allocation error in constitutive_dislo.f90
Franz Roters
2007-10-24 15:24:49 +0000
20693b9d38hybrid IA sampling now respects sample symmetry
Philip Eisenlohr
2007-10-24 15:13:10 +0000
bf35fa4868former 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
f0617b74f6added output of mesh statistics at end of input file parsing
Philip Eisenlohr
2007-10-24 09:00:42 +0000
67a5711ed1added awareness of set definitions (new function build_elemSetMapping) added new style tables for initial conditions
William Counts
2007-10-23 13:09:46 +0000
ca0e2e55dffixed 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
9043b4374eCorrected the treatment of interaction coefficients in constitutive_dislo
Luc Hantcherli
2007-10-22 15:32:45 +0000
57ed908129removed unused variables made hardening matrix finally (?) work
Franz Roters
2007-10-22 14:53:10 +0000