Franz Roters
7b27606000
modified constitutive description in line with other dislocation density based models
2014-03-30 15:04:06 +00:00
Martin Diehl
4bb5439994
removed copyright from *.f90 and *.f files, will be added using the deployMe script for release/nightly build. Now based on presence of svn:keywords property set (usually "Id"). Similar thing will be done for .py, .sh, and .config files but more involved because they might be executable
2014-03-29 08:20:36 +00:00
Martin Diehl
aae3f95f76
added script to add Copyright information (complete header) to f90 files having the svn property "MPIE" set.
...
Development version don't need an header anymore
2014-03-26 14:57:33 +00:00
Martin Diehl
62497b648a
removed forgotten debug statement
2014-03-26 09:09:28 +00:00
Martin Diehl
58ad750b2a
small changes regarding restart and HDF5
2014-03-26 08:41:45 +00:00
Martin Diehl
99fe2dfccd
small changes to HDF5, put link to homepage into README not to have double effort
2014-03-25 17:21:47 +00:00
Martin Diehl
2659ee51d4
restart is working for all solvers, but slight deviations still occur. Reason is most probably, that for CPFEM init, a 0s timestep is required which seems to alter the state a tiny bit, leading to small deviations. For AL and Polarization, this is even more severe since the regular call to CPFEM has F-residual as the current deformation gradient, however the init step uses only F.
2014-03-25 15:44:16 +00:00
Martin Diehl
c7994db5ef
writing out f_aim and f_aim_last inc for restart, was calculated as average before
2014-03-24 09:59:30 +00:00
Franz Roters
04b2011195
reintroduced dislocation hardening for bcc
2014-03-17 15:42:46 +00:00
Martin Diehl
bac0e24d21
added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given
2014-03-13 23:50:55 +00:00
Martin Diehl
5d4017bbfc
introduced homogenization_none to substitute isostrain with ngrains 1.
...
cleaned up homogenization isostrain (has no state)
2014-03-13 23:20:50 +00:00
Duancheng Ma
670c0caabc
2014-03-13 12:12:19 +00:00
Martin Diehl
3f7a389ff7
changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used
2014-03-13 06:43:49 +00:00
Martin Diehl
7c79b31f6c
remove not needed use statement
2014-03-13 05:50:56 +00:00
Martin Diehl
e1a31457ba
fixed potential array mismatch in assignment
2014-03-13 05:49:07 +00:00
Martin Diehl
fced0168f0
added a little bit more HDF5 functionality
...
made abbreviations in configure options consistently capitals
2014-03-12 16:51:01 +00:00
Martin Diehl
9afc1e3a15
fixed wrong omp statement
2014-03-12 15:29:14 +00:00
Martin Diehl
3aea8b39e9
added some HDF5 functionality (needs to be activated with preprocessor makro)
2014-03-12 07:33:51 +00:00
Martin Diehl
ff1b1c1a50
fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code
...
forgot to commit dislotwin last time, now seems to work
2014-03-11 23:55:40 +00:00
Martin Diehl
ef8fbf4dda
introduced case in dislotwin for bcc (peierls stress as critical stress)
...
renamed some parameters, now exponents (p,q) for slip are per family, shear band got own ones. exponent for twin (r) now per family
fixed bugs in lattice ("empty" interactions should be 1 not zero)
see both Phase examples to check what is needed
2014-03-11 17:40:59 +00:00
Franz Roters
6781d7d8a5
last line got lost
...
as it is fcc it should actually be only one slip system family
2014-03-11 10:49:48 +00:00
Martin Diehl
c2e81414e3
updated TWIP steel data set to have only 2 slip system families, added data set for tungsten
2014-03-11 10:34:36 +00:00
Martin Diehl
721183e702
removed some forgotten debug statements
2014-03-08 22:51:32 +00:00
Martin Diehl
2b589c3d71
moved reading in of lattice type and elastic constants to lattice module
...
removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config
constitutive modules will only be initialized if needed
homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-08 20:50:31 +00:00
Martin Diehl
32493675d6
fixed bug in constitutive introduced in rev 2988 causing wrong names in *.outputConstitutive
2014-03-05 08:06:21 +00:00
Christoph Kords
0265978941
renamed "maxNmatIDs" accordingly to "maxNinstances"
2014-03-04 13:47:04 +00:00
Martin Diehl
a0d75ee05e
moved quaternion disorientation to lattice because it requires the knowledge of the lattice structure.
2014-02-28 13:28:27 +00:00
Martin Diehl
92bf1565fc
formatting
2014-02-28 13:05:11 +00:00
Martin Diehl
b0f191c88c
also rename matID/i to instance like in the other constitutive models
2014-02-28 13:03:21 +00:00
Martin Diehl
b9722b866c
improved error handling when getting number of threads
2014-02-28 10:30:07 +00:00
Martin Diehl
650b71ffa9
renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling
2014-02-28 10:18:40 +00:00
Martin Diehl
fe13bb16fe
added some consistency check and remove unnamed constants
2014-02-27 18:57:46 +00:00
Martin Diehl
ad203d32c1
circular dependencie (copy and paste error) removed
2014-02-26 12:04:52 +00:00
Martin Diehl
35be468b24
fixed detection of working directory for ifort
2014-02-25 21:10:55 +00:00
Martin Diehl
a554bc5267
now ignoring c/a ratio keyword
2014-02-19 07:40:21 +00:00
Martin Diehl
69992638a8
added missing covera_ration in hex test. was not necessary before for constitutive_none, now it is
2014-02-17 09:45:54 +00:00
Martin Diehl
8bbba55506
function version seems to be quasei-standard, not subroutine for getcwd
2014-02-14 19:59:37 +00:00
Martin Diehl
4efb1e5617
fixed bug: wrong initialization of lattice structure, had effect on j2 and phenopowerlaw
2014-02-14 18:05:59 +00:00
Franz Roters
d3d4931061
corrected checking for ~/.damask/damask.conf
2014-02-14 15:58:01 +00:00
Franz Roters
4399fb9535
improved comment on Euler angles, now includes unit
2014-02-10 15:48:00 +00:00
Martin Diehl
d8311d46b3
.. last commit message was missing because I wanted to abort, so:
...
move reading in of lattice structure to lattice (elastic constants, c/a ratio, lattice structure). Still needs a little bit of documentation, the lattice structure names are not easy to understand.
tested it for j2 and phenopowerlaw and seems ok, but before continuation with the other constitutives models might simplifiy a little bit.
2014-02-10 14:38:05 +00:00
Martin Diehl
d45aea4467
moved reading in of lattice stru
2014-02-10 14:31:19 +00:00
Martin Diehl
20bea80f0c
changed order to have PETSc variables before others, fixed problem when MPI and FFTW were installed in standard location but PETSc build own MPI
2014-02-09 18:30:31 +00:00
Martin Diehl
33ab0ff0fe
introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway
2014-02-09 17:08:13 +00:00
Martin Diehl
a07c738778
polished a little bit
2014-02-08 10:48:09 +00:00
Martin Diehl
8741f2368d
first (not too complicated) workaround for the shortcomings of the new intel compiler. still crashes, but later ;)
2014-02-06 17:48:01 +00:00
Martin Diehl
bc4cc20c55
slightly polished makefile and numerics, recompile fftw with icc; still not running with ifort 14.01
2014-02-06 10:41:34 +00:00
Martin Diehl
1ea48bb8ff
added accumulatedshear as output for twin and slip
2014-02-05 15:26:22 +00:00
Martin Diehl
83e4bc3247
updated test for restart, now also checking if number of increments is the same
...
fixed output of test class
fixed small bug? in asciitable.py, @philip: please check
2014-02-03 19:30:28 +00:00
Martin Diehl
1444739d60
prevent ifort from complaining too loud about old syntax in openmpi files
2014-02-03 16:47:16 +00:00
Christoph Kords
caca65148d
fixed bug in restart: stiffness values that were read from file used to be overwritten by "utilities_constitutiveResponse"
...
Still not fixed in polarization and AL solver !!!
2014-02-03 15:57:04 +00:00
Christoph Kords
b9b87a785c
added crystallite_F to debugging output
2014-02-03 13:19:49 +00:00
Franz Roters
c5b31b373c
no longer needed
2014-02-03 10:38:59 +00:00
Christoph Kords
aa4eee2b50
fixed bug in math_rotate_forward3333
2014-01-31 14:28:20 +00:00
Martin Diehl
c7e59821d1
fixed compilation bug when using IMKL single core
2014-01-30 15:38:10 +00:00
Martin Diehl
986591682c
remove equals sign in rpath definition to work with mac, changed linker options in compile_CoreModule to work with linux again, please let me know if mac doesn't work anymore (it should)
2014-01-30 13:07:44 +00:00
Martin Diehl
9e5a2d8e10
fixed runpath
2014-01-29 19:52:53 +00:00
Martin Diehl
f82e5ec26e
some small changes on the installation scripts, IMKL now works with gfortran, more verbosity and more sanity checks
2014-01-29 12:20:56 +00:00
Martin Diehl
ee99216edc
polishing compile_CoreModule, bug fix for Makefile
2014-01-28 07:34:13 +00:00
Martin Diehl
02441bef0f
introduced rpath to tell executable where libraries are located
2014-01-27 14:10:49 +00:00
Pratheek Shanthraj
b9e55b4495
moved setting of CPFEM forwarding flag to the forwarding subroutine
2014-01-27 11:05:36 +00:00
Martin Diehl
656e3cc6ca
some more small changes, including bug fix in configure + relative linking in Makefile for spectral solver
2014-01-27 10:12:29 +00:00
Martin Diehl
4850e5987d
remove LD_LIBRARY_PATH ( http://linuxmafia.com/faq/Admin/ld-lib-path.html ) and fixed wrong named xxx_ROOT
2014-01-25 00:15:04 +00:00
Martin Diehl
d66b41c686
fixed installation, will test now
2014-01-24 12:30:27 +00:00
Martin Diehl
c99ca8b0ee
fixed magnesium parameters
2014-01-23 13:33:09 +00:00
Martin Diehl
e6d88e4458
prevent twin volume fraction from going above 1.0
2014-01-22 15:47:49 +00:00
Christoph Kords
391d83fa2b
moved rate sensitivity block to where it belongs to, namely the analytic jacobian calculation
2014-01-22 15:34:10 +00:00
Pratheek Shanthraj
4b7890f133
fixed wrong indices in analytic jacobian calculation following christoph's changes. convergence is restored at least for the spectral solver
2014-01-22 10:16:55 +00:00
Christoph Kords
3dfdbaff5b
Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged.
...
In contrast, the analytical jacobian will probably be affected by this change!
@Pratheek: Can you check with me how this can be fixed?
2014-01-22 08:38:13 +00:00
Christoph Kords
c10d0d15cb
crystallite_integrateStateFPI: convergence check always ran about all e,i,g; now checking only those specified in FEsolving_execElem, FEsolving_execIP
2014-01-21 18:45:41 +00:00
Christoph Kords
ff648492c1
perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode
...
(missed some part in last commit)
2014-01-21 18:28:21 +00:00
Christoph Kords
029826d84f
perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode
2014-01-21 16:35:12 +00:00
Christoph Kords
351aa6ba29
setting todo and converged flag only for those e,i,g that are to be executed (FE_elemExec etc)
...
used to spoil convergence checking in non ping pong mode and lead to slow performance
2014-01-21 16:09:00 +00:00
Philip Eisenlohr
74cbfb6d48
added few comments
2014-01-17 01:38:35 +00:00
Christoph Kords
77cba87bf3
in crystallite_integrateStress: pass full 2nd PK stress to LpAndItsTangent instead of only the deviatoric part; otherwise constitutive laws that depend on the hydrostatic pressure do not work; constitutive laws that require a deviatoric stress need to calculate it internally (as is e.g. already done in the J2 model)
2014-01-16 10:36:40 +00:00
Christoph Kords
a011b0a4a2
corrected tangent of plastic velocity gradient dLp_dT for slip (twin volume fraction was missing)
2014-01-16 09:53:45 +00:00
Christoph Kords
9223cf52b4
corrected non-Schmid effect: used to calculate resolvedstress based on Mandel notation which can be incorrect for nonsymmetric non-Schmid tensors
2014-01-16 09:50:40 +00:00
Christoph Kords
b24dd62431
corrected tangent of plastic velocity gradient dLp_dT for slip (twin volume fraction was missing)
2014-01-16 09:29:21 +00:00
Christoph Kords
563c4ac94b
corrected tangent of plastic velocity gradient dLp_dT, which improves convergence a lot
2014-01-16 09:18:26 +00:00
Martin Diehl
278f854e41
introduced functionality to get output results by point from the HDF5 file
2014-01-15 22:17:32 +00:00
Martin Diehl
3fffb2ae3e
added parameter set for Magensium
2014-01-15 12:29:44 +00:00
Martin Diehl
9f3a17b1b6
fixed bug in hexagonal slip system definition
2014-01-15 12:08:46 +00:00
Franz Roters
98f677916c
fixed crash in RtoEuler for PHI=180°
2014-01-15 08:32:43 +00:00
Martin Diehl
449c66c70f
did some testing...
2014-01-14 19:32:55 +00:00
Martin Diehl
01c8c434f7
finished work on the new installation scheme (if no bugs are found ;))
2014-01-14 19:03:41 +00:00
Martin Diehl
aab0208c65
more modifications, only configure needs modifications
2014-01-12 17:29:00 +00:00
Martin Diehl
00eb9b2cf6
forgot Makefile
2014-01-11 17:18:59 +00:00
Christoph Kords
a1927b60a7
allow flux between neighboring material points if one of both got a velocity of zero
2014-01-08 16:35:10 +00:00
Christoph Kords
b41593f279
blocked density does no longer contribute to dislocation flux; used to generate source term for mobile density if burgers vector had any non-zero component perpendicular to blocking boundary.
2014-01-08 16:31:50 +00:00
Christoph Kords
629295242f
added two new outputs: "rho_dot_flux_mobile" and "rho_dot_sgl_mobile"
...
corrected outputs of dislocation density rates: take into account that blocked dislocation density is a signed quantity
2014-01-07 15:55:28 +00:00
Martin Diehl
79307de6a7
added some code parts for HDF5 in Fortran
2013-12-30 10:06:01 +00:00
Martin Diehl
5b4b132176
introduced preprocessor statements in mesh to hide unneeded variables as a preparation for HDF5
2013-12-27 20:03:28 +00:00
Martin Diehl
831e905b16
added missing lp output to crystallite
2013-12-20 11:13:12 +00:00
Pratheek Shanthraj
de3fc70a02
moved restart writing to the forwarding routine
2013-12-20 10:49:14 +00:00
Martin Diehl
98d5a99a43
reverted unwanted changes to Makefile (configure)
2013-12-20 10:31:33 +00:00
Martin Diehl
d595f59f0e
improved sanity checks, mainly for RGC
2013-12-20 08:36:15 +00:00
Martin Diehl
af075aaa9d
found bug (array index out of bounds), Philip/Franz please check if correct (Gfortran was complaining)
2013-12-20 08:33:46 +00:00
Franz Roters
543639bb8f
added support for Marc/Mentat2013.1
2013-12-19 13:37:35 +00:00
Martin Diehl
61c6839723
fixed over-sensitive error in homogeniztion and ifort option in makefile
2013-12-19 08:49:47 +00:00
Pratheek Shanthraj
6f6bd86078
fixed bug in heat generation and analytic jacobian calculations
2013-12-18 10:19:48 +00:00