Philip Eisenlohr
61f8a5fcbe
corrected UpperCase "SolidSolutionStrength" in config file parsing.
2010-10-15 12:59:00 +00:00
Martin Diehl
3c502561ee
added matrix multiplication 3333x33 to math.f90
...
added some parameters for spectral method to numerics.f90 (tolerance)
changed error message concerning spectral method in IO.f90
corrected calculation of stress BC in mpie_spectral.f90
2010-10-13 16:04:44 +00:00
Christoph Kords
724960686f
* in nonlocal model: dislocation flux now with valid description also for large angle grain boundaries; transfer of dislocations from one slip system to the other according to new variable "constitutive_nonlocal_compatibility" which depends on the angle between slip plane normals and slip directions and is updated once per cycle in function "crystallite_orientations"
...
* reactivated debugging functionality for "non-standard" integration methods
2010-10-12 13:08:54 +00:00
Christoph Kords
efd92d9b51
allocation of special state variables for "non-standard integration" was erroneous. used to produce segmentation fault for integration method 4 and 5.
2010-10-12 12:11:24 +00:00
Christoph Kords
fce7590c17
* now able to choose method for state integration (integrator and integratorStiffness in numerics.config); standard is "Fixed Point Iteration", which is basically the same as before; others available are "Explicit Euler", "AdaptiveEuler", "Classical Runge-Kutta" and "Runge-Kutta Cash-Karp"
...
* now remembering stiffness similar to how we do it for Lp etc.; avoids undefined stiffness values for nonconverged stiffness calculation
* non-local stuff:
* changed non-local kinetics (Gilman2002)
* enforce zero shearrate for overall carrrier density below relevant density
* enforce zero density for those states that become negative and were below relevant density before
* dislocation velocity is not limited by V^(1/3) / dt anymore
2010-10-01 12:18:49 +00:00
Martin Diehl
9927cd7adb
changed format of resolution. instead of reading exponent of two, now the number given equals directly to the number of FPs. Allows us to use arbitrary number of FPs.
...
example:
OLD: reading 5 results in 2^5=32 FPs
NEW: reading 33 results in 33 FPs
2010-10-01 10:42:15 +00:00
Alankar Alankar
2cc7bc38d7
Fixed a missing ','. Assigned plane and direction names.
2010-10-01 06:32:44 +00:00
Philip Eisenlohr
1ae252c1ce
altered hexagonal slip systems
...
(needs comparison to former set, which is potentially wrong)
2010-09-30 17:31:19 +00:00
Claudio Zambaldi
b85f9b0770
material.config: changed texture of microstructure 1 from 1 to 3, because 1 needs a ODF file which is not in the repository and makes the one_element_model.proc complain.
2010-09-30 12:55:38 +00:00
Philip Eisenlohr
a8dfdc6487
added _mul33xx33 double contraction: A_ij * B_ij = scalar
2010-09-30 09:32:49 +00:00
Christoph Kords
cb4b4afc88
corrected header for math_qRot
2010-09-30 08:46:58 +00:00
Philip Eisenlohr
b259fbd9c6
1) added distribution of leapfrog breaks
...
2) lattice_symmetryType is now a function (former lookup array was buggy)
3) stricter check of state var values (>0!) and memory deallocation done
2010-09-30 07:31:53 +00:00
Alankar Alankar
e562df35a9
Added more state dependent variables in post results.
2010-09-29 06:35:08 +00:00
Philip Eisenlohr
3a4743f1d3
fixed bug in grain naming scheme. occurred for gaps in consecutive numberings...
2010-09-28 16:56:49 +00:00
Martin Diehl
5a51082977
polishing, changed calculation of stress BC.
...
prevent updating of gamma hat in case of perfect plasticity
2010-09-24 13:27:53 +00:00
Alankar Alankar
d94236b2d6
Declaration of filename etc. was missing. Fixed now.
2010-09-23 09:13:46 +00:00
Philip Eisenlohr
83a921eec2
crystallite: polishing
...
rest: added "$ID" and line ending settings...
2010-09-23 08:05:50 +00:00
Philip Eisenlohr
a467fb4ae3
forgot "$ID" and setting line endings...
2010-09-23 08:03:00 +00:00
Philip Eisenlohr
89136afade
debug now has config file, too
2010-09-23 08:01:41 +00:00
Philip Eisenlohr
23746e7e78
chopped line beginnings
...
introduced check for PIL, if this fails, no image will be output
2010-09-23 07:59:52 +00:00
Martin Diehl
6ea8623f65
added functions math_mul33x3_complex math_Plain99to3333 to math.f90.
...
mpie_spectral.f90 uses both functions
math_Plain99to3333 is used for inversion of c0
math_mul33x3_complex is used for equilibrium check in fourier space
also did some cleaning up on mpie_spectral.f90
2010-09-22 12:04:43 +00:00
Alankar Alankar
80618c9814
Cleaned-up, rearranged state variables
2010-09-22 09:00:40 +00:00
Alankar Alankar
7548858ffa
Some useful changes in damper for quicker convergence.
2010-09-22 08:54:36 +00:00
Martin Diehl
604992a9e1
fixes bugs in screen output, changed gmsh output, now the deformed mesh is written out
2010-09-22 08:51:34 +00:00
Martin Diehl
6920d8d939
cleaning up, added new output to gmsh and to screen, prepared calculation of deformed mesh
...
new convergence criteria: divergence of stressfield (in fourier space)
2010-09-21 15:38:25 +00:00
Alankar Alankar
8488cafe3a
now consistent with numerics.f90 values
2010-09-21 14:37:52 +00:00
Alankar Alankar
18c86ca0b1
Updated corresponding files for titanmod constitutive
2010-09-13 09:29:03 +00:00
Philip Eisenlohr
9a20f742ea
leaner memory requirement to store states and their rates when calculating the crystallite stiffness
2010-09-13 09:13:25 +00:00
Philip Eisenlohr
d8a4cd1739
added optional argument --baseindex to specify what colorset should be influenced. identify colors, for instance, start at idx 214...
2010-09-13 07:31:27 +00:00
Leyun Wang
e604e0946e
2010-09-12 20:12:53 +00:00
Alankar Alankar
006dce7d8a
printing euler angles (and defgrads and orientations) for individual ips now works
2010-09-08 11:55:57 +00:00
Martin Diehl
1693bfca47
changed calculation of defgrad from cauchy-stress to pk-stress, now working in large-strain-formulation.
...
output of msh-file changed to deformed configuration, removed output of defgrad as a field
2010-09-07 16:37:55 +00:00
Alankar Alankar
3c5f38643d
Brushed up accountability of twinning to Lp
2010-09-07 14:44:37 +00:00
Christoph Kords
7d4c7f7fa7
corrected bug in debug: stressloop info was erroneous
...
mesh init now showing much more ip statistics in verbose mode
2010-09-07 09:06:02 +00:00
Philip Eisenlohr
4a31b175da
1) introduced default colors
...
2) complaints are more verbose
2010-09-06 16:51:11 +00:00
Philip Eisenlohr
a3c92061f9
1) terminallyIll was reset before FE did cutback --> useless extra calculations of same problem over and over...
...
2) local stiffness calculation is now standard for non-local grains
3) stressLoopDistribution discriminates between (a) central solution and (b) stiffness perturbation
4) debugger is switched on as standard... (but verboseDebugger not!)
2010-09-06 16:06:41 +00:00
Martin Diehl
00922705eb
now using current average stiffness for reference material stiffness
2010-09-06 10:00:59 +00:00
Philip Eisenlohr
c19524f264
1) fixed terminallyIll bug occurring with Ngrains=1 homogenization
...
2) brushed up the output to be more easily readable/understandable
2010-09-01 21:04:02 +00:00
Alankar Alankar
4f29e8c2fe
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 10:07:52 +00:00
Martin Diehl
a5c228fd02
changed back transform to complex-to-real, removed redundant variables, reduced size of arrays where possible
2010-09-01 08:05:11 +00:00
Martin Diehl
85febf0803
added fftw3 as fft(library will not versioned, should be in a linkable folder) , did some corrections on the code, splitted main file up (allows use of makefile), added makefile
...
changes on mpie_spectral.f90:
new structure, changed variable names, now using defgrad instead of disgrad, cleaned up, removed augmented Lagrange.
ToDo: Implement Augmented Lagrange again (but then a working version), implement Large strain, think about complex-to real-transform backwards, try to implement MP-support
2010-08-27 16:39:38 +00:00
Christoph Kords
05568ef73b
calculation of state damper was erroneous for use with nonlocal and other constitution at the same time
2010-08-26 12:54:31 +00:00
Christoph Kords
197d7ea8d9
rearranged output statements in crystallite init for better readability in output file
2010-08-26 09:52:40 +00:00
Alankar Alankar
ec4269bdd8
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 11:36:51 +00:00
Philip Eisenlohr
da53656c9a
added C3D20R (type 9) to list of possible non-local elements (Christoph, should we invent a reasonable error number for this? "-1" may not be optimal...)
2010-08-19 23:00:26 +00:00
Philip Eisenlohr
c352146fe6
exchanged forall (warning causers) to nested do loops
...
deleted (currently unused) symmetrization of H tensor in CPFEM.f90 --> Denny pls advise?
2010-08-19 21:35:38 +00:00
Philip Eisenlohr
af2437c0af
removed hardwired installation directory from submitX scripts and added corresponding flexibility into make script
2010-08-17 21:22:24 +00:00
Philip Eisenlohr
3b671faae7
added some flexibility (and bug-fixing) to the pre/post processing. the file 'MSCpath' contains the LOCAL path to the MSC installation (/msc in our case)
2010-08-17 18:21:22 +00:00
Franz Roters
f78b07448f
some modifications in line with Davids analytical model
...
input some additional variables via material.config
2010-08-17 14:23:55 +00:00
Philip Eisenlohr
bb9899e7de
1) added element type 57 (C3D20R) reduced quadratic hexahedral
...
2) need discussion about the present omission of the FE_nodesAtIP for type 8 --> added comment into source code
2010-08-16 22:53:24 +00:00