cdb2dd8808
crystallite:
...
added full check whether residuum in inner loop is NaN
SingleCrystallite now takes/stores Lp guess
CPFEM_Taylor (missing in GIA!!):
former fully plastic guess for Lp has been exchanged by keeping the last converged Lp (global array CPFEM_Lp) to serve as new best guess for the next time step. This speeds up the inner loop of TimeIntegration.
2008-04-16 17:00:28 +00:00
fa2d6b9b6d
Make use of average consistent tangent dPdF_bar to fill eventual holes in the individual dPdF of one orientation (during consistent tangent calculation)
2008-04-11 13:04:25 +00:00
664a29c183
Subscript Crystal_ changed to Lattice_ according to the module Lattice
2008-04-10 12:13:37 +00:00
9d2ce61698
A revision for CPFEM_Taylor.f90.
2008-04-10 11:22:17 +00:00
d2312e81ff
deleted CFPFEM.f90 and crystal.f90 from repository
...
adopted mattex.mpie to use lattice_structure instead of crystal_structure
2008-04-09 15:35:16 +00:00
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
# two symbolic links are required: constitutive.f90 and CPFEM.f90
2008-04-07 14:54:29 +00:00
c063ce5bc1
add a kronecker delta function
2008-03-27 11:54:34 +00:00
5462d2c754
just a correction for a little typo :P
2008-03-27 11:53:53 +00:00
f3eda7cdaa
include a new function epsilon_ijk, i.e., the permutation function (used to compute the cross product of two tensors)
2008-03-26 13:35:01 +00:00
92ee2cbf7d
to include 2 material parameters, i.e., grain size and the length of burgers vector used in GIA model
2008-03-26 09:33:12 +00:00
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
2008-03-25 14:19:22 +00:00
7faf4093a5
reflect change of mesh_element(2,e) in mesh.f90
...
-change of line, and those below, will be ignored--
M constitutive_pheno.f90
2008-03-25 13:59:22 +00:00
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!
2008-03-25 12:52:27 +00:00
6d721dc16c
Modifications of cpfem_marc2005r3 according the ones in cpfem_marc2007r1
2008-03-14 21:56:46 +00:00
1267ce78b6
CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.)
...
renamed CPFEM_stressIP to CPFEM_MaterialPoint
renamed CPFEM_stressCrystallite to CPFEM_Crystallite
introduced new global variables to keep track of FE state within module
FEsolving
2008-03-14 21:32:57 +00:00
7975d62586
CPFEM_general now independed of the wrapper code (UMAT, hypela2, etc.)
...
renamed CPFEM_stressIP to CPFEM_MaterialPoint
renamed CPFEM_stressCrystallite to CPFEM_Crystallite
introduced new global variables to keep track of FE state within module
FEsolving
2008-03-14 21:29:31 +00:00
0a08d9eff0
Corrected errors from the last release of constitutive_dislo
...
Changed parameter designation in mattex file
2008-03-14 16:23:08 +00:00
6950eee59b
altered stress and strain "language" to first PK and deformation gradient.
...
plot_results are now part of the crystallite result and not explicitly called from the homogenization layer.
2008-03-14 14:19:10 +00:00
6f06133b90
removed some unused variables
2008-03-14 08:23:43 +00:00
778cbcd3ca
changed file read formats from explicit types to more general (and not
...
harmful) '*'
2008-03-12 13:53:00 +00:00
5641cfd477
changed precision assignment to intrinsic (and standardized)
...
"selected_*_kind"
2008-03-07 15:21:46 +00:00
c74e071f5e
updated constitutive_xxx to reflect new 9x9 dLp_dTstar - used to be 3x3x3x3
2008-02-21 12:50:37 +00:00
9c1e8a7944
added acceleration capability after time-step cut backing
2008-02-21 10:33:34 +00:00
9707fd0f8f
full update, i.e. my development snapshot
2008-02-19 12:58:46 +00:00
f4edf4bd0c
2007
...
shortened stress_recovery = ...
added (crucial) deMandeling for return values of CPFEM_general
2005
matched to 2007 version
2008-02-18 18:49:06 +00:00
2b567ad20d
CPFEM
...
changed consistent tangent scheme to D.Tjahjanto Diss version
fixed nasty bug in decision when to calculate new values (CPFEM_general)
prec
precision values should now be on safe side.
2008-02-18 18:11:24 +00:00
199c65cdba
changed pert_e to new perturbation value pert_Fg
2008-02-18 10:46:18 +00:00
1106bb27eb
introduced logical "debugger"
2008-02-18 10:45:46 +00:00
e07d38fdec
ooops. missed nInner and nOuter as loop counters...
2008-02-18 10:41:37 +00:00
9c9cdeb58b
corrected nStress and nState to reflect more general nInner and nOuter (now vompatible with changes in prec.f90)
2008-02-18 10:39:39 +00:00
cfaa0e696d
exchanged TINY (from prec) with intrinsic "tiny" function (Fortran90)
2008-02-15 12:42:27 +00:00
a41a4a75ef
CPFEM:
...
time integration scheme changed to start from guessing of Lp.
Newton-Raphson algorithm now accelerates as long as residuum improves.
IO:
added error 650 for polar decomposition problems
prec:
some cleanup. Introduced "relevantStrain" for convergence checking in CPFEM
2008-02-15 10:04:49 +00:00
f39c68dded
CPFEM_stressCrystallite must be called with t not dt in CPFEM_stressIP
2008-02-08 08:25:19 +00:00
6a3a9e8f98
Added more detailed reporting
2008-01-11 18:32:02 +00:00
720d8da82b
Corrected syntax
2008-01-10 19:00:45 +00:00
b8b171c95b
Introduce debug module, contains distributions of nCutback, nStressLoop, and nStateLoop
2008-01-10 18:53:57 +00:00
3a7ec38a3d
Delete unused physical constants from pheno
...
Rename Boltzmann Constant Kb->kB
2008-01-10 18:51:58 +00:00
1233d01cd5
mesh_maxNsharedElems is now correct. Added check for multiple identical nodes in an element, e.g. triangular quad-element and the like.
2008-01-10 17:12:33 +00:00
fe70a82d6d
!!!!! IMPORTANT !!!!!
...
All subroutines were committed at once:
- constitutive_pheno works
- constitutive_dislo without twinning also works
This release should serve as reference
2007-12-14 13:36:04 +00:00
3ef451824c
just a test
2007-12-14 13:28:41 +00:00
90195fe1c8
Commit dislo
2007-12-11 14:30:21 +00:00
ab0bde576e
Commit of IO
2007-12-11 14:23:33 +00:00
32fc32865f
Commit of mesh
2007-12-11 14:23:21 +00:00
08a1c38c73
Finalizing the implementation of deformation twinning in constitutive_dislo.f90
...
Assume to be ready to use
2007-12-10 12:55:43 +00:00
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
2007-12-07 15:12:43 +00:00
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
2007-12-07 14:17:42 +00:00
b69aa2b112
Construction of the elasticity matrices for rotated/twin regions.
...
Definition of constitutive_Ntwin and constitutive_MaxNtwin
2007-12-07 13:59:15 +00:00
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" ...
2007-12-07 12:36:17 +00:00
f63f6949cc
Corrected undefined "state" in HomogenisedC
...
Added "forall" loops in C construction
2007-12-07 10:54:19 +00:00
d615406722
Corrected HomogenisedC:
...
1- No dependence on material volume fractions!!!
2- Can now depend on state of microstructure
changed some loops from Nstatevars to Nslip
2007-12-07 10:37:06 +00:00
Philip Eisenlohr