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
Denny Tjahjanto
d2335e4a0b
2008-06-10 10:23:32 +00:00
Franz Roters
60f3571266
added code lines for OpenMP, however it does not work yet!
2008-05-26 13:11:25 +00:00
Philip Eisenlohr
51e596c81d
small polishing to be compatible with NAG compiler (basis for automatic differentiation)
2008-04-29 12:18:17 +00:00
Luc Hantcherli
9626c25bfb
crystallite: restoration of last known good Lp_guess
...
CPFEM_Taylor: exit whenever cutback limit exceeded
2008-04-28 16:04:08 +00:00
Philip Eisenlohr
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
Luc Hantcherli
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
Luc Hantcherli
664a29c183
Subscript Crystal_ changed to Lattice_ according to the module Lattice
2008-04-10 12:13:37 +00:00
Denny Tjahjanto
9d2ce61698
A revision for CPFEM_Taylor.f90.
2008-04-10 11:22:17 +00:00
Franz Roters
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
Denny Tjahjanto
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
Denny Tjahjanto
c063ce5bc1
add a kronecker delta function
2008-03-27 11:54:34 +00:00
Denny Tjahjanto
5462d2c754
just a correction for a little typo :P
2008-03-27 11:53:53 +00:00
Denny Tjahjanto
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
Denny Tjahjanto
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
Luc Hantcherli
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
Philip Eisenlohr
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
Philip Eisenlohr
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
Luc Hantcherli
6d721dc16c
Modifications of cpfem_marc2005r3 according the ones in cpfem_marc2007r1
2008-03-14 21:56:46 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
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
Luc Hantcherli
0a08d9eff0
Corrected errors from the last release of constitutive_dislo
...
Changed parameter designation in mattex file
2008-03-14 16:23:08 +00:00
Denny Tjahjanto
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
Franz Roters
6f06133b90
removed some unused variables
2008-03-14 08:23:43 +00:00
Denny Tjahjanto
778cbcd3ca
changed file read formats from explicit types to more general (and not
...
harmful) '*'
2008-03-12 13:53:00 +00:00
Philip Eisenlohr
5641cfd477
changed precision assignment to intrinsic (and standardized)
...
"selected_*_kind"
2008-03-07 15:21:46 +00:00
Philip Eisenlohr
c74e071f5e
updated constitutive_xxx to reflect new 9x9 dLp_dTstar - used to be 3x3x3x3
2008-02-21 12:50:37 +00:00
Philip Eisenlohr
9c1e8a7944
added acceleration capability after time-step cut backing
2008-02-21 10:33:34 +00:00
Philip Eisenlohr
9707fd0f8f
full update, i.e. my development snapshot
2008-02-19 12:58:46 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
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
Denny Tjahjanto
199c65cdba
changed pert_e to new perturbation value pert_Fg
2008-02-18 10:46:18 +00:00
Denny Tjahjanto
1106bb27eb
introduced logical "debugger"
2008-02-18 10:45:46 +00:00
Denny Tjahjanto
e07d38fdec
ooops. missed nInner and nOuter as loop counters...
2008-02-18 10:41:37 +00:00
Denny Tjahjanto
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
Philip Eisenlohr
cfaa0e696d
exchanged TINY (from prec) with intrinsic "tiny" function (Fortran90)
2008-02-15 12:42:27 +00:00
Philip Eisenlohr
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
Franz Roters
f39c68dded
CPFEM_stressCrystallite must be called with t not dt in CPFEM_stressIP
2008-02-08 08:25:19 +00:00
Hongmei Li
6a3a9e8f98
Added more detailed reporting
2008-01-11 18:32:02 +00:00
Luc Hantcherli
720d8da82b
Corrected syntax
2008-01-10 19:00:45 +00:00
Luc Hantcherli
b8b171c95b
Introduce debug module, contains distributions of nCutback, nStressLoop, and nStateLoop
2008-01-10 18:53:57 +00:00
Luc Hantcherli
3a7ec38a3d
Delete unused physical constants from pheno
...
Rename Boltzmann Constant Kb->kB
2008-01-10 18:51:58 +00:00
Philip Eisenlohr
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
Luc Hantcherli
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
Luc Hantcherli
3ef451824c
just a test
2007-12-14 13:28:41 +00:00
Luc Hantcherli
90195fe1c8
Commit dislo
2007-12-11 14:30:21 +00:00
Luc Hantcherli
ab0bde576e
Commit of IO
2007-12-11 14:23:33 +00:00
Luc Hantcherli
32fc32865f
Commit of mesh
2007-12-11 14:23:21 +00:00
Luc Hantcherli
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
Luc Hantcherli
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