10358fc1f0
(hopefully) improved the cutbacking scheme somewhat. introduced a bunch
...
of debugging to the .out file
2008-08-01 07:47:01 +00:00
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
6abff4b3ba
changed hybridIA_population from being wrong integer type to correct real
...
type.
2008-07-16 16:35:13 +00:00
bdcf76cd74
fixed stupidity in hybridIA evaluation. added informative output during
...
ODF parsing.
2008-07-16 14:49:17 +00:00
7a7a5d2ae3
messed up the parameter list of TimeIntegration with the last submit
2008-07-15 08:32:44 +00:00
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
63d81b92b7
added math_mul33x33, 66x66, and 99x99 to avoid OMP_Critical "evilmatmul"
2008-07-08 19:38:22 +00:00
de9e35cdb8
modified cutback counting
...
restore Lpguess_old in case of iInner > nInner
2008-06-25 06:26:41 +00:00
7453a18b42
in lattice:
...
- add slip-twin interaction type matrix
2008-06-17 14:42:45 +00:00
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
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
6921c57c7d
rev 225 had a small syntax error...
2008-06-17 11:54:34 +00:00
3676b453ca
added two more element types (157 and 136)
2008-06-17 09:11:54 +00:00
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
d2335e4a0b
2008-06-10 10:23:32 +00:00
60f3571266
added code lines for OpenMP, however it does not work yet!
2008-05-26 13:11:25 +00:00
51e596c81d
small polishing to be compatible with NAG compiler (basis for automatic differentiation)
2008-04-29 12:18:17 +00:00
9626c25bfb
crystallite: restoration of last known good Lp_guess
...
CPFEM_Taylor: exit whenever cutback limit exceeded
2008-04-28 16:04:08 +00:00
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
Philip Eisenlohr