386ca7ebc3
added some basic I/O functions to mpie_spectral.f90
2010-06-10 08:50:04 +00:00
740db98090
- use upwind differences scheme for fluxes instead of central differences
...
- introduced possibility to "reflect" fluxes at free surfaces in order to have quasi-periodic fluxes
2010-06-09 08:56:00 +00:00
11b98fbfa6
added ID functionality and correct line endings
2010-06-08 10:10:57 +00:00
82687d31c8
first shot at mpie_spectral wrapper
2010-06-08 10:08:15 +00:00
7e78932976
renamed "mpie_cpfem_init" subroutine within mpie_interface module to more general name "mpie_interface_init"
2010-06-08 09:35:13 +00:00
71c6179ee5
dislocation stress had wrong sign in some components
2010-06-08 09:34:23 +00:00
3ce1882d29
when calculating dislocation stress at grain boundaries, densities are now extrapolated similarly to like it was already done at free surfaces
2010-06-07 16:01:37 +00:00
3989cc2688
dislocation stress now considers dislocation segments of length V^(1/3) and ensures non-singular solution for points that are collinear to the dislocation line by means of a core spreading parameter "a". stress is scaled by a reference cutoff radius "r" to avoid mesh-dependent stress values.
...
corrected formula for conversion of anisotropic elastic constants to isotropic Youngs modulus and poissons ratio
2010-06-07 14:32:23 +00:00
5283ba5546
if it is pReal it should be 1.0 not just 1
2010-05-31 07:12:45 +00:00
c1fed9d197
polishing of verbosity
2010-05-28 16:05:16 +00:00
76b392e12f
dislocation stress calculation in constitutive_nonlocal_microstructure: corrected mapping from lattice to slip system frame and the neighboring dislocation density formulation
2010-05-28 13:12:36 +00:00
0d69623b74
grpmf...
2010-05-27 14:32:16 +00:00
89f09ed784
use testMe in this subdir to check for compilation errors.
2010-05-27 14:31:04 +00:00
920bccc556
single IP, reduced integration element for 3D
2010-05-26 17:30:16 +00:00
bb73a74325
there was a line break missing so that it did not compile
2010-05-26 17:10:26 +00:00
9932c90c73
restoring the crystallite.f90, IO.f90 and material.f90 to the version 575. the version 576 of these files contain some uncleared debugging statements.
2010-05-26 16:34:44 +00:00
d114a600c3
fixed some errors in QuaternionToEuler, RToQuaternion to get the correct texture representation in the output.
2010-05-26 15:52:54 +00:00
e5f0af638e
fixed calculation of dislocation stress for elements with a single ip and no neighbors; used to get segmentation fault
2010-05-26 09:17:10 +00:00
0c5bc83469
in nonlocal_microstructure: at the model surface the excess density is now extrapolated from opposite neighbor instead of assuming zero density; this results in better modeling of the dislocation stress at the surface
...
restructured nonlocal_dotstate, to be able to easily switch on and off particular effects in the microstructure evolution
nonlocal_dotstate now enforces a cutback when single density runs the risk of becoming negative; in the case of a state already below the relevantState dotState is set to zero
introduced two new output variables: rho_dot_edge and rho_dot_screw
2010-05-21 08:51:15 +00:00
dba4ae7ef1
new numerics parameter "defgradTolerance" replaces "relevantStrain" in cpfem, since "relevantStrain" is also used in crystallite but in a different context.
...
"defgradTolerance" gets the same standard value as "relevantStrain" so calculations should be consistent with older ones
2010-05-20 14:55:11 +00:00
3a23dc4bbd
first call to hypela2 in inc 0 is done with computation mode 2 (=plain calc), but since rev. 539 the second call was done with computation mode 6 (=only return values) instead of 2. This doesnt work, because computation mode 2 might calculate stress and tangent for all materialpoints, but the mapping from lattice to current configuration is only done for the element/IP that triggers the call.
...
So switched back to the former order of calls which also works for marc2010:
el 1 IP 1 | el 1 IP 2
|
inc mode lovl | inc mode lovl
------------------ | ------------------
0 2 6 | 0 2 6
1 6 4 | 1 6 4
1 4 6 | 1 3 6
1 6 4 | 1 6 4
1 1 6 | 1 2 6
Still to be checked:
according to Franz the lovl order could be different in contact analysis (e.g. starting with lovl 4 instead of 6)
2010-05-20 09:14:59 +00:00
51da3b6242
sorry, in rev 568 some variables were not declared in crystallite_postResults
2010-05-18 12:36:09 +00:00
472684c30e
make element type 8 (= hex8 reduced integration) available for nonlocal constitutive law
2010-05-18 08:13:23 +00:00
19b6fae88b
added a missing exclamation mark in front of an $OMP directive
2010-05-18 08:07:55 +00:00
d96780967b
added some output variables in crystallite: F, Fe, Fp, Ee, P, S
2010-05-18 07:57:13 +00:00
15328e51cf
Like this Euler angles AND grain rotation should be OK.
2010-05-11 15:06:21 +00:00
483584ceb3
improved version of general I/O functions
...
renamed module cpfem_interface to mpie_interface
2010-05-11 06:57:15 +00:00
e810e5cfa1
changed fileopen routines in IO.f90 to use solver dependant functions defined in the interface routines
...
DID NOT test ABAQUS versions yet, not shure whether (V)GETOUTDIR returns path with or without terminating slash
2010-05-10 15:02:59 +00:00
0416c9a616
reworked build_ipNeighborhood
...
added C3D8R reduced integration hexahedral element
2010-05-10 14:54:59 +00:00
2405a51042
crystallite:
...
- grainrotation calculation now is done with symmetryID 0, i.e. without symmetry reduction since we want the absolute misorientation.
- While math has everything in radians, post results eulerangles and axisangle are given in degrees.
And: grain rotation seems OK after the previous changes in math module.
2010-05-07 12:01:46 +00:00
97f206d0bc
new errors regarding Spectral method input file parsing
...
mesh setup from Spectral input file
$Id$ in FEsolving activated
2010-05-06 16:40:47 +00:00
5a0963fb81
added variable working dir (if apply_script is called from somewhere else) and made the script executable
...
why using "comp_user?*" instead of "comp_user*" ?? (@Franz...)
2010-05-06 14:17:53 +00:00
e4067f2d26
debugged zoo of rotation operations and transformations
...
all angles now in radians
introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
2010-05-06 14:07:21 +00:00
4954228b3a
forgot to change this file accordingly (main.msb revmoved)
2010-05-06 13:19:45 +00:00
f1950cacc9
corrected apply script
...
removed binary menu file, it is now created by the apply script
2010-05-06 13:17:28 +00:00
f23b4d16d1
updated section with dislotwin paramters
2010-05-06 09:01:54 +00:00
383a36b226
added comment on reformating to common block files
...
added modifications necessary for running the subroutine with Marc/Mentat2010
2010-05-06 08:50:54 +00:00
393a67e587
now stating instance when sanity checks run into trouble
2010-05-05 13:36:59 +00:00
1ee5df41fd
Sorry, that was the wrong file, now the correct one...
2010-05-04 17:34:51 +00:00
48c554811b
A simple one-element model that you can use to try out changed code before commits. Hopefully, the start of a testing suite.
2010-05-04 17:32:17 +00:00
0a36527f61
fixed bugs in math_quaternionDisorientation
...
though running, there still seem to be some problems with the quaternion functions, since crystallite_grainRotation gives some dubious results!
2010-05-04 16:02:05 +00:00
ba48e5d3ab
fixed syntax errors in last commit -- hopefully no further ones...
2010-05-04 13:03:35 +00:00
cecffffe2b
disorientation with symmetryType == 0 returns plain misorientation
...
complain about unknown symmetryType
2010-05-04 12:54:13 +00:00
adaeb0e1ef
rotation axis is undefined for small angles and therefore gets a NaN
2010-05-03 18:27:19 +00:00
9a257fb68b
moved "UserOutput" into figures of Manual
2010-04-29 15:31:40 +00:00
2494da8f51
moved Manual to its own folder.
...
set up include sub-structure.
based on KOMA-script book class and heavily borrowing from Wiley book recently written.
2010-04-29 15:28:29 +00:00
ea718573ef
return 180° disorientation about 100 for neighbors with different phases
2010-04-29 12:28:51 +00:00
8ab101da50
safe use of acos function in math_QuaternionToAxisAngle
2010-04-29 10:33:39 +00:00
15e5dcf8f2
math_QuaternionToAxisAngle now safe for small or close to 2pi rotations.
2010-04-29 10:01:09 +00:00
e8719cb6b8
Table driven input works wirh 2010!
...
That leaves DDM to do for tomorrow ;-)
2010-04-29 09:43:31 +00:00
Martin Diehl