4a0c1c1717
comments
2015-03-15 15:30:14 +00:00
12cb5e56fe
fixed interpretation 'spectralOut' header. Seem to have been buggy all the time, but manifested itself only after keyword change. Also added ':' after keywords to be more robust
2015-03-15 15:29:38 +00:00
ac3d82493e
changed output keywords form resolution to grid and from dimension to size
2015-03-13 14:45:33 +00:00
e4d70e971d
cp is commercially pure ;)
2015-03-13 13:49:09 +00:00
d8debc3396
updated reference geom files (used old names), updated numerics.config to reflect use of analytic jaco
2015-03-13 06:47:26 +00:00
aa125bac5b
aliases for grid and size not supported any more
2015-03-12 22:30:24 +00:00
00cba25a44
improved update of ip coordinates for spectral solver.
...
do not need to create a new fftw plan all the time, using data already defined for the convolution
2015-03-12 22:28:33 +00:00
263f997cf2
seems to be an old output
2015-03-12 08:51:56 +00:00
b338f1aeb8
added some outputs
2015-03-09 23:10:36 +00:00
6b0a2ed5a4
Changed initialization of projection matrix for transformation based on the number of slip systems
2015-03-09 12:05:38 +00:00
a5ee02ff01
added some aliases for output to finally converge from 'shear_rate' to 'shearrate' and from 'accumulatedshear' to 'accumulated_shear'
2015-03-09 11:29:44 +00:00
2ab90c28a5
replaced example ODF by reasonable data in fully correct format
2015-03-06 15:56:58 +00:00
35f9e91e73
some corrections to the damage models
2015-03-06 13:12:32 +00:00
bbb5ff6ae9
changes related to intermediate configuration kinematics:
...
- switched Fi and Li from state variables to crystallite variables
- Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration
- T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi
- Updated analytic jacobian to take into account tangents wrt Fi
- Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 13:09:00 +00:00
4fd7e0165e
more meaningful name
2015-03-05 11:07:55 +00:00
36f14a2d1d
Removed nr = ns condition
2015-02-26 14:16:29 +00:00
2baac6fc10
Fix for 3881 because gfortran was complaining
2015-02-24 10:15:07 +00:00
3c2a95efd8
if matrix inversion error encountered during analytic tangent calculation, give a warning and return elastic tangent
2015-01-30 14:00:32 +00:00
f1f4f939a3
analytic jacobian default for all solvers
2015-01-29 13:59:49 +00:00
dd8458a775
updated analytic jacobian calculation to correctly take into account intermediate configuration Fi. improved convergence of Li loop in stress integration
2015-01-29 13:58:25 +00:00
0b59519a2a
updated damage models:
...
- coupling to plasticity handled within damage module instead of plasticity module
- anisotropic models more stable
2015-01-29 13:56:09 +00:00
12c586afb9
removed unused variables
2015-01-23 13:44:16 +00:00
1f666dffc9
fixed bug: output size was not set
2015-01-23 13:08:25 +00:00
e9c3137a3e
disabled a faulty ntrans loop (matmul was complaining)
2015-01-23 11:23:35 +00:00
43095b1f1b
added threshold_stress_slip as output (got lost once upon a time)
2015-01-22 09:34:11 +00:00
1c2dd37e28
deleted unused variables in dislo UCLA
2015-01-21 15:14:00 +00:00
e57cb2b4f4
Wrong number of slip families in BCC material config files
2015-01-20 17:22:52 +00:00
19dc945065
Inroduced Fd in kinematic decomposition.
2015-01-19 13:39:47 +00:00
54dc4943de
adaptive penalty to bound Cv in (0,1)
2015-01-16 18:31:14 +00:00
f16244d240
corrected state integration
2015-01-16 17:51:43 +00:00
eecb78665c
removing some unused code and cleaning up
2015-01-16 17:32:47 +00:00
6411d36633
updated vacancy-damage model
2015-01-16 17:04:01 +00:00
93e50366bb
fixed automatic determination of compiler
2015-01-16 11:19:30 +00:00
701f35c326
reverting the changes related to undamaged_Fi during rev. 3850/3851.
2015-01-16 00:09:22 +00:00
85465333ec
2015-01-15 22:42:56 +00:00
b1ef2f9dd9
2015-01-15 22:40:17 +00:00
62ff5d1d45
gfortran/ifort is now selected depending on the compiler provided by PETSc
2015-01-15 11:35:53 +00:00
ced00954fe
added new model by david (LLNL/UCLA)
2015-01-15 10:56:15 +00:00
125276283b
fix in analytically integrated expression.
2015-01-15 09:22:24 +00:00
19611660c2
passing undamaged Fe for calculation of stress.(Pratheek please check.)
2015-01-13 09:54:44 +00:00
9c8873f16c
added functions to get degraded stiffness and undamaged intermediate deformation gradient.
2015-01-13 09:46:25 +00:00
80bb9a2e66
function to calculate mises stress.
2015-01-13 09:43:05 +00:00
c2b8627981
Currently, there is no updating for the dot state of vacancy generation, set sizeDotState = 0
2015-01-11 20:01:43 +00:00
a2552d6241
Polishing.
2015-01-11 19:10:16 +00:00
e20271b89a
fix for temperature in spectral solver
2015-01-08 13:26:00 +00:00
9635ae9f3d
fix in integration.
2015-01-07 14:51:06 +00:00
8610531aae
removed unused variable.
2015-01-04 20:13:16 +00:00
460c0d73fa
Making use of extended plastic state with pointers, removed helper functions.
2015-01-04 19:26:33 +00:00
e6d1b21b57
introduced extended plastic state type (including pointer/aliases)
2014-12-19 14:01:03 +00:00
e4cb387571
polishing comments
2014-12-18 18:44:47 +00:00
72cb9072e1
added some debug statements
2014-12-18 18:41:02 +00:00
7fe20fe805
intent(out) specifier for getXX functions
2014-12-18 18:26:16 +00:00
fdbaf00fa7
cleaned up naming and added elastic material with volume preserving stiffness
2014-12-18 18:20:11 +00:00
299518adfa
cleaned tungsten parameter set, renamed to be consistent
2014-12-18 16:55:43 +00:00
3f14ebe43d
helper routine to get heat generation rate needed for MARC/Abaqus
2014-12-17 13:37:13 +00:00
e6e72ab3a7
In mesh_marc_build_elements use mesh_NcpElems instead of mesh_Nelems (line 2906).
...
Fixes error when not all elements are CP elements.
Some formating and comment improvements.
2014-12-15 11:51:32 +00:00
7c8f2ede31
fixed missing initialization causing a NaN in gdot and crashes
2014-12-11 15:04:16 +00:00
5bc93efd63
Fix bug about vacancy generation in subroutine constitutive_microstructure
2014-12-11 13:07:15 +00:00
bdb6bf8090
recasting vacancy dot state as microstructure state
2014-12-10 13:40:30 +00:00
39a3a5b5a0
cleaned up FEM examples
2014-12-09 18:59:37 +00:00
3064c14ff6
fig a bug of energy barrier of vacancy generation.
2014-12-09 18:56:19 +00:00
3f8678c7c4
plastic dot states evaluated at effective stress not damaged stress
2014-12-09 18:23:48 +00:00
55e2de6ffd
made dotstate and state pointers (targets are not allowed in derived types) and introduced plastic state having aliases to parts of the state.
2014-12-09 12:12:53 +00:00
d1e1191df4
corrected bug: renamed file list for commercial FEM solvers to start with small letter
2014-12-08 16:36:22 +00:00
2458aa9c3e
BBar stabilisation now optional (default is off)
2014-12-08 16:27:23 +00:00
c34252c71c
changed constitutive_* prefix for plasticity laws to plastic_* to be consistent with novel thermal_, damage_, vacancy_, etc. laws.
2014-12-08 15:55:30 +00:00
0bd4ab48f8
numerics with respect to vacancy diffusion.
2014-12-08 15:37:44 +00:00
db0e02b94d
changed to lower canse convention
2014-12-08 15:37:02 +00:00
67c3eb3a2e
bounding the analytically integrated damage variable
2014-12-08 14:05:38 +00:00
347dac74c6
added vacancy state integrators
2014-12-08 10:45:12 +00:00
9301b56aa3
initialize restoration points of vacancy concentration at material points, postResults of VC at material points.
2014-12-08 09:27:32 +00:00
62f254c9de
age results and write restart data for vacancy concentration
2014-12-08 09:09:54 +00:00
55c2b729a0
while node coordinates are always 3D displacements are passed into hypela2 as 2D array (dispt) for 2D simulations
2014-12-05 07:35:28 +00:00
334ade5260
Every line in a recipe must begin with a tab character.
...
http://stackoverflow.com/questions/14576237/errors-in-make-file-missing-separator-stop
2014-12-03 13:37:04 +00:00
9d12e026e8
renumbered some errors, fixed non-existing errors
2014-12-03 00:42:35 +00:00
01a8c500ca
detabbing…
2014-12-03 00:41:12 +00:00
fa6f784102
deleted unused stuff
2014-12-02 21:49:33 +00:00
82942d1618
bug fix for previous commit
2014-12-02 17:23:47 +00:00
b796e46f98
residualStiffness is not just an FEM parameter
2014-12-02 17:17:35 +00:00
2fc264d90a
forgot file in previous commit
2014-11-28 20:07:22 +00:00
37dd574790
recast power law and viscous dot state variables as derived quantities using analytical solutions of the evolution equation
2014-11-28 20:06:24 +00:00
8d64e513da
improved commenting and variable names.
...
values for unused slip/twin families are ignored in material.config (do not need to be numeric anymore…)
2014-11-26 21:23:02 +00:00
f8b29bdaf0
fixed bug in type definition for ‘drivingForce’
2014-11-26 18:40:31 +00:00
54de033c45
using analytic solution for local/nonlocal coupled damage (as opposed to having crystallite iterate to find a consistent solution)
2014-11-26 17:29:45 +00:00
61bb9a373e
all damage models consistently semi-explicit
2014-11-26 11:32:50 +00:00
6e90260817
updated local damage post results for newer models
2014-11-26 11:32:19 +00:00
a91ebafa0e
semi-explicit damage evolution
2014-11-26 10:36:14 +00:00
646f1c48dc
typo, IO_warning and not IO_error
2014-11-26 10:16:30 +00:00
9714a172d1
updated to get local damage from microstructure state
2014-11-26 00:42:20 +00:00
c89ef9de61
relaxed some sanity checks
2014-11-26 00:29:54 +00:00
26c01f64e5
updated file list
2014-11-26 00:06:16 +00:00
32ee3a4a66
corrected constitutive_getVacancyPotentialDrivingForce function
2014-11-26 00:04:06 +00:00
36dc59b09f
reverted back to rate independent damage evolution. updated example config files to recent changes
2014-11-25 23:56:52 +00:00
a4584fecd2
potential driving force for conservative cahn hilliard vacancy diffusion
2014-11-25 20:13:33 +00:00
20437b2ae0
added brittle(elastic energy release rate)/ductile(vacancy condensation) phase field damage model coupled to vacancy concentration
...
simplified vacancy_generation
2014-11-25 17:23:37 +00:00
65b22f9a40
calculate analytic jacobian from xxx0 and not xxxsub0
2014-11-25 17:20:52 +00:00
e6edc3d8cd
Some sanity checks
2014-11-25 13:20:12 +00:00
296146490f
changed some default parameters for FEM
2014-11-21 14:20:54 +00:00
3d4b2d2f52
tolerance for damage can be any positive value. Set to 0.001 if found negative.
2014-11-21 13:12:54 +00:00
ee1d409bad
critical shear can be set to zero for family not in use.
2014-11-21 09:03:54 +00:00
Martin Diehl