28fc9cc221
put recalculation of elastic stress back in (was accidentally deleted)
...
changed check of residuum improvement for small residuum
2007-08-13 07:12:24 +00:00
85df3b95ee
changed restriction of dTstar, its now based on the residuum, if it does not improve 0.5*dTstar is used instead
2007-08-07 08:02:31 +00:00
a5f49a0578
changed first elastic guess for small Tstar
...
changed residuum to only consider deviatoric part of Tstar
changed scaling when dTstar becomes too large to proportionally scaling all components of dTstar
2007-07-20 13:32:44 +00:00
a1c3df84ab
changed quite a few things, mainly using last result as initial guess now, thinks work a lot better now but still not perfect
2007-06-06 14:38:06 +00:00
0efe7756bc
included fix in constitutive_post_results to handle case of constitutive_Nresults==0
2007-05-22 06:36:52 +00:00
6f3f5ee23c
corrected error in deMandeling of jacoby matrix
2007-05-22 06:34:58 +00:00
09c47f15e3
changed constitutive to open file using IO_open_file and use fileunit 200 instead of 1
...
changed constitutive to skip empty lines when parsing files
activated call quit in IO_error
2007-05-16 14:36:03 +00:00
e1ae35c5e4
corrected typo in constitutive_Sslip_v
...
corrected demandle of jacoby in hypela2
2007-05-15 08:47:44 +00:00
49886a072c
corrected CPFEM_stressIP to use FFN1 at the end of the time increment
...
calculate for ncycle==0 and not only inc==0 in hypela2
2007-05-14 12:13:36 +00:00
6aece52d1c
changed Calculation of the tangent of Lp to always use state(i) and not constitutive_state_new
2007-05-11 13:22:16 +00:00
e5749cbff0
changed scatter check for gauss components to if(noise==0)
2007-04-26 14:46:54 +00:00
09aad5ccb4
changed _Xvalue() return codes on error to huge()
2007-04-26 12:40:06 +00:00
11cc6d1ee1
math_sampleGaussOri did not treat the special case of zero noise (scatter)
2007-04-26 11:46:58 +00:00
6fd473eaa7
scatter was not transformed to rad when read from material file
2007-04-26 11:42:16 +00:00
1fd722e4cf
changed format for reading float from('F') to * in IO_floatValue
2007-04-25 14:38:22 +00:00
033a6a03b3
determine number of grains as texture_Ngrains(mesh_element(4,cp_en)) instead of constitutive_Ngrains (which is not used anymore) in CPFEM_stressIP
2007-04-25 13:58:10 +00:00
1747c0cc98
corrected calculation of constitutive_maxNstatevars
2007-04-25 13:50:28 +00:00
03f244995e
corrected mesh_build_elements (unit) so that state variables are read correctly from input file
2007-04-25 07:33:24 +00:00
3a4ec4c2d1
removed several unused variables
2007-04-24 06:49:13 +00:00
c0ebb4d846
reordered subroutine calls in mesh_init to allocate nodemap and elementmap BEFORE they are used
2007-04-24 06:29:22 +00:00
83d2dea8c6
CPFEM_general call now uses CPFEM_en instead of cp_en -- the required mapping is not in place on first call!
2007-04-23 13:25:09 +00:00
471d53e3ba
_general needs CPFEM_en NOT cp_en..!
...
stress loop counter was not been reset
matmul(A,B) now replaced by stored AB
2007-04-23 13:23:03 +00:00
9e0b9a7096
removed syntax error in _LpAndItsTangent (enforcing symmetry was flawed...)
2007-04-17 08:02:27 +00:00
0f6ff82f30
Checked indices for Jacobi calculation
2007-04-17 07:58:53 +00:00
c76ab95b21
Added in LpandIstTangent 0.5*(transpose+untranspose)
2007-04-17 07:53:06 +00:00
b26704c23f
added missing use statement in _init
2007-04-13 15:37:10 +00:00
4718abfc70
agreed on wrong calc of scaled state residual
2007-04-13 14:22:55 +00:00
09495e852d
added necessary use statements
2007-04-13 14:20:59 +00:00
bbfb86fb8f
added necessary use statements
2007-04-13 14:20:33 +00:00
d90162dc38
Corrected Rstress(i)/state_new(i) into Rstate(i)/state_new(i)
2007-04-13 11:30:49 +00:00
eb292917b9
untied interlinked constitutive and CPFEM modules
...
now constitutive stores EulerAngles and CPFEM generates Fp_old from these
2007-04-11 15:28:46 +00:00
cb4a85319d
corrected small syntax errors
2007-04-11 14:55:06 +00:00
3c5d490148
added _plot_results function to tell about constitutive results (state + current shear per slip system)
2007-04-11 14:52:47 +00:00
f9f3e2bd9b
result processing now in stressIP
...
i.e. no extra polar decompositions performed any more
2007-04-11 14:51:49 +00:00
9704a4e83c
syntax polishing
2007-04-11 13:21:22 +00:00
ff6642ea8c
major edit.
...
combined routines into
CPFEM_general, CPFEM_stressIP (incl cutback scheme),
CPFEM_stressCrystallite, and CPFEM_timeIntegration (Newton scheme)
error management now based on text strings
2007-04-11 10:06:28 +00:00
9da4b4cf24
name change to reflect State and Stress (was outer,inner)
...
added nReg for Jacobi regularization attempts
2007-04-11 10:05:08 +00:00
4743c1cd86
error from inversion routines now boolean type
2007-04-11 10:04:22 +00:00
6095ce0972
minor changes (syntax polishing)
2007-04-11 10:03:43 +00:00
e1f2daf753
_Lp and _dotState now explicitely depends on given state
2007-04-10 12:29:11 +00:00
df3f327ece
added function calls in init()
2007-04-10 11:22:53 +00:00
35266bd513
added nreg to limit number of Jacobi regularizations
2007-04-10 11:22:06 +00:00
84eeddc6ba
small polishing
2007-04-10 11:21:34 +00:00
092bff57d3
Added some modifications in constitutive_LpandItsTangent
2007-04-10 11:20:05 +00:00
db12071577
Finalized constitutive_Assignment (orientation to subips)
2007-04-04 13:57:32 +00:00
6fcf763054
Corrected Tstar calculation by 0.5
2007-04-04 13:56:10 +00:00
ea20fc73c1
Added pure statements
2007-04-04 13:55:21 +00:00
f1653ad88c
debugged quicksort algorithm slightly (now able to do multiple identical keys)
2007-04-04 08:49:48 +00:00
f055a7b037
added IO_continousIntValues
2007-04-04 08:49:00 +00:00
5db61c19d3
fixed buggy build_ipNeighborhood and faceMatch (that sounds so easy...)
2007-04-04 08:47:34 +00:00
Franz Roters