61f8a5fcbe
corrected UpperCase "SolidSolutionStrength" in config file parsing.
2010-10-15 12:59:00 +00:00
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
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
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
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
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
2cc7bc38d7
Fixed a missing ','. Assigned plane and direction names.
2010-10-01 06:32:44 +00:00
1ae252c1ce
altered hexagonal slip systems
...
(needs comparison to former set, which is potentially wrong)
2010-09-30 17:31:19 +00:00
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
a8dfdc6487
added _mul33xx33 double contraction: A_ij * B_ij = scalar
2010-09-30 09:32:49 +00:00
cb4b4afc88
corrected header for math_qRot
2010-09-30 08:46:58 +00:00
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
e562df35a9
Added more state dependent variables in post results.
2010-09-29 06:35:08 +00:00
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
d94236b2d6
Declaration of filename etc. was missing. Fixed now.
2010-09-23 09:13:46 +00:00
83a921eec2
crystallite: polishing
...
rest: added "$ID" and line ending settings...
2010-09-23 08:05:50 +00:00
a467fb4ae3
forgot "$ID" and setting line endings...
2010-09-23 08:03:00 +00:00
89136afade
debug now has config file, too
2010-09-23 08:01:41 +00:00
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
80618c9814
Cleaned-up, rearranged state variables
2010-09-22 09:00:40 +00:00
7548858ffa
Some useful changes in damper for quicker convergence.
2010-09-22 08:54:36 +00:00
604992a9e1
fixes bugs in screen output, changed gmsh output, now the deformed mesh is written out
2010-09-22 08:51:34 +00:00
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
8488cafe3a
now consistent with numerics.f90 values
2010-09-21 14:37:52 +00:00
18c86ca0b1
Updated corresponding files for titanmod constitutive
2010-09-13 09:29:03 +00:00
9a20f742ea
leaner memory requirement to store states and their rates when calculating the crystallite stiffness
2010-09-13 09:13:25 +00:00
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
3c5f38643d
Brushed up accountability of twinning to Lp
2010-09-07 14:44:37 +00:00
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
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
00922705eb
now using current average stiffness for reference material stiffness
2010-09-06 10:00:59 +00:00
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
4f29e8c2fe
Introduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
2010-09-01 10:07:52 +00:00
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
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
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
197d7ea8d9
rearranged output statements in crystallite init for better readability in output file
2010-08-26 09:52:40 +00:00
ec4269bdd8
Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
2010-08-23 11:36:51 +00:00
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
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
f78b07448f
some modifications in line with Davids analytical model
...
input some additional variables via material.config
2010-08-17 14:23:55 +00:00
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
9f071c3e01
Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
2010-08-06 15:53:45 +00:00
298cecbfec
corrected syntax errors (long lines, line continuation by \) and logical mistake in mpie_cpfem_marc / abq_std which went unnoticed in ifort. Using SunStudio f90 surfaced those...
2010-08-03 23:47:00 +00:00
60decea760
Calibration (Prismatic slip) of the new dislocation density based CPFE model for single crystal alpha-Ti (Titanmod). Formulated edge and screw dislocations mobilities in presence of kink pairs. Updated comments. The model shows the plateau in stress evolution in prismatic orientation which is observed in experiments.
2010-07-15 07:16:15 +00:00
49926d5d66
F77 --> F90 polishing/condensation as far as possible.
...
next step is change of Re/Im FFT to Re only...
2010-07-13 15:29:26 +00:00
4d110126da
adopted improvements done by Arun Prakash.
...
mesh:
elemType identification based on lower case
Abaqus now reports more errors
IO:
new function to inquire whether inputfile contains "parts"
new function to assemble multiply included inputfile into a flat one
awareness of range generation in element numbers
error reporting
2010-07-13 10:26:07 +00:00
eb7830dc8f
New Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
2010-07-13 08:19:25 +00:00
397db06fb3
cpfe_general now returns also Kirchhoff-stress P and dPdF
...
added dummy variables P and dPdF to mpie_cpfe_* to be able to call cpfem_general correctly
2010-07-07 09:58:18 +00:00
80016f8429
added P and dPdF for the call to cpfem_general, made some comments in the code, aligned some loops
2010-07-07 09:10:54 +00:00
Philip Eisenlohr