Pratheek Shanthraj
3d83ac64fd
added code for local damage
2014-08-01 14:54:57 +00:00
Martin Diehl
a8c0ca66fc
corrected linking: use PETSC MPI (when available) rather then system wide version.
...
little polishing (e.g. debug output for new state)
2014-07-23 13:26:05 +00:00
Martin Diehl
0252fea3d7
fixed small bug in newstate
2014-07-10 08:47:00 +00:00
Martin Diehl
146d151357
disabled multiphysics output because homogenization output was not accessible anymore
2014-07-09 13:34:33 +00:00
Luv Sharma
7dfda9d8af
remved bug in phenopower law state initialization and a minor change in Makefile
2014-07-07 14:21:58 +00:00
Martin Diehl
27d861decc
fixed calculation of post results size introduced with new state
2014-07-03 13:17:29 +00:00
Luv Sharma
d5952138e3
corrected wrong dimensions of aTolstate and inclusion order of files in abaqus and marc interfaces
2014-07-03 12:52:33 +00:00
Martin Diehl
8fa2dcffbd
changed to new state, please report bugs to Luv or Martin
2014-07-02 12:27:39 +00:00
Pratheek Shanthraj
ac36190f10
no need to allocate a large array for constitutive(thermal,damage)_sizePostResults (changes apply only to new state)
2014-06-30 14:47:30 +00:00
Pratheek Shanthraj
676cba13a8
modified state integrators to perform staggered integration of plasticity damage and thermal internal variables
2014-06-24 23:21:25 +00:00
Pratheek Shanthraj
6347adcb65
fixed error in previous commit
2014-06-23 22:09:15 +00:00
Pratheek Shanthraj
d3356bde4e
moved mappings for new state layout to material
2014-06-22 18:58:29 +00:00
Martin Diehl
2160b37b1e
fixed warning message in constitutive
2014-06-17 15:24:44 +00:00
Martin Diehl
3533138936
removed delta state for J2 and phenopowerlaw when using new state
2014-06-17 06:54:49 +00:00
Martin Diehl
5bd80123b4
fixed two severe bugs in new state, now a little bit faster (399 vs. 378 sec) and more memory efficicient (251312 vs. 236832 kbytes) for Phenopowerlaw 16x16x16 example (short load only)
2014-06-16 19:19:38 +00:00
Luv Sharma
5755ed547b
applied new State for nonlocal model and a trivial changes in crystallite
2014-06-13 20:53:17 +00:00
Luv Sharma
8130f4b21f
changes to make newstate aware of damage variables
2014-06-11 16:32:09 +00:00
Luv Sharma
401b31c951
new state related changes
2014-06-11 12:27:41 +00:00
Pratheek Shanthraj
639ca89133
DAMASK (except nonlocal) now sees and interacts with FEM solver
2014-06-06 00:38:29 +00:00
Martin Diehl
164252213b
more work on the new state
2014-05-27 14:46:03 +00:00
Philip Eisenlohr
5d21ef5a84
fixed selective debugging logic flaw in crystallite_stressAndItsTangent
2014-05-19 16:43:32 +00:00
Philip Eisenlohr
88f94cbbec
increased debug threshold for some outputs (decreased verbosity level at "basic")
2014-05-15 13:01:54 +00:00
Martin Diehl
4bfced1a70
fixed small allocation flaws for new state
2014-05-12 13:00:37 +00:00
Martin Diehl
838a8be1b9
removed forgotten debug statements from last commit
...
fixed material.config for abaqus example
2014-05-12 09:28:32 +00:00
Luv Sharma
d820a5aaa6
fixed some errors in changes related to newstate
2014-05-12 00:44:44 +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
7c79b31f6c
remove not needed use statement
2014-03-13 05:50:56 +00:00
Martin Diehl
9afc1e3a15
fixed wrong omp statement
2014-03-12 15:29:14 +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
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
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
d45aea4467
moved reading in of lattice stru
2014-02-10 14:31:19 +00:00
Christoph Kords
b9b87a785c
added crystallite_F to debugging output
2014-02-03 13:19: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
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
Martin Diehl
831e905b16
added missing lp output to crystallite
2013-12-20 11:13:12 +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
Pratheek Shanthraj
6f6bd86078
fixed bug in heat generation and analytic jacobian calculations
2013-12-18 10:19:48 +00:00
Martin Diehl
d5bdcf5a6d
some small improvements, still something wrong, probably in crystallite.f90
2013-12-13 14:14:17 +00:00
Martin Diehl
84d51d743c
fixed small issues, ATTENTION: Still a bug introduced with rev. 2783 not fixed
2013-12-13 13:19:17 +00:00
Martin Diehl
65ae979920
indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input
2013-12-12 17:09:59 +00:00
Martin Diehl
0bc99a9622
change latticeName to latticeID
2013-11-27 16:20:27 +00:00
Martin Diehl
2f7efa2055
introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file
2013-11-27 08:04:05 +00:00
Christoph Kords
883669bd77
in preguess of FPintegrator: resetting of previous dotstates can be done in same loop as collectdotstate
...
polishing: indentation level and capitalization
2013-11-21 10:58:41 +00:00
Pratheek Shanthraj
724380c624
Coupled phase field implementation for the basic petsc solver. needs clearing up… work in progress
2013-11-12 00:44:23 +00:00
Martin Diehl
3ecc8103f0
introduced dummy heat calculation, available as an crystallite output
...
temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-18 18:57:28 +00:00
Martin Diehl
03e26af37b
small bufix, didn't compile
2013-10-16 20:44:04 +00:00
Martin Diehl
dc95c82d4a
removed temperature integration and corresponding data structures and debugging options
...
temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver.
introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver.
went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 13:04:59 +00:00
Philip Eisenlohr
42cd4eb38f
fixed some e12.6 string formats to avoid compile time warnings
2013-10-08 13:51:36 +00:00
Martin Diehl
3bf0ed84c8
some comments and minor improvements.
...
removed tol_gravityNodePos from prec because it's not used
tol_math_check is now used by spectral driver to check rotation
2013-09-14 10:59:35 +00:00
Martin Diehl
c93373a4ac
removed long lines
2013-08-08 19:33:21 +00:00
Franz Roters
4f9dbfa193
last commit also solves problem with infinite loop in case of a cutback in homogenization
...
corrected indentation of a comment
2013-08-02 07:59:55 +00:00
Franz Roters
a08d263865
remnoved some unused fuctions from use statements
2013-08-02 06:18:41 +00:00
Philip Eisenlohr
98528f9a89
added possibility for multi-level inclusion of files in *.config and loadcase files.
...
include subfiles by stating
{path/to/include}
2013-06-26 19:19:00 +00:00
Martin Diehl
6f7740a243
made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables
2013-06-11 16:35:04 +00:00
Christoph Kords
55735b9580
grainrotation output now uses corrected rotation math
2013-06-07 09:35:00 +00:00
Martin Diehl
c7ba8a2a9b
introduced plain mode (no ping pong) again and added test for MSC.MArc 2012 as prove that its working
2013-05-17 17:52:46 +00:00
Martin Diehl
cfac73043d
some lines in crystallite.f90 exceeded limit of 132 characters
...
fixed bug in Lib_VTK_IO.f90 (informed author as well)
added zip file containing source of Lib_VTK_IO.f90
2013-05-08 17:48:00 +00:00
Christoph Kords
3b788720e5
added crystallite output "neighboringelement" and "neighboringip"
...
some polishing in crystallite_postResults
2013-05-08 12:02:30 +00:00
Martin Diehl
b98ce3d22d
fixed out of range bug introduced with singleRun changes in rev. 2356
2013-04-30 12:14:07 +00:00
Martin Diehl
533fcec96f
named some loops and made BUTCHER TABLEAUs parameters in crystallite_integrateStateRKCK45
2013-04-29 11:17:30 +00:00
Martin Diehl
9fb6d8167b
simplified logic for single run, now not checking for extra arguments but simply checking if start and end elem/IP are the same
2013-04-26 13:23:36 +00:00
Christoph Kords
a7706090d2
neighborhood relations were screwed up for some elementtypes as of version 2318, might have caused segmentation faults with nonlocal model, now fixed
2013-04-22 13:35:35 +00:00
Martin Diehl
8b2d9d8155
further improvements on implementing the plain mode
2013-04-16 17:07:27 +00:00
Christoph Kords
44eb771583
added crystallite output of current ip coordinates
2013-03-27 10:12:28 +00:00
Franz Roters
02f092492e
updated copyright header
2013-03-22 17:35:05 +00:00
Martin Diehl
7e0ea5dd04
hickups in doxygen documentation fixed
2013-03-05 14:35:26 +00:00
Martin Diehl
ada2beb8b8
reorganized calculation modes for CPFEM, now having better readable and cleaner structure
2013-03-01 11:48:29 +00:00
Christoph Kords
875abcd06f
replaced some array assignments, which might have caused trouble when using the "timeSyncing" option, by parallelized loops of scalar value assignments
2013-02-27 16:36:12 +00:00
Christoph Kords
3736ab1d97
added shortcuts for two do loops in the "timesyncing" procedure
...
math_pDecomposition enclosed by omp critical statement; somehow, this subroutine is not thread safe, not clear what exactly the problem is
2013-02-27 10:32:37 +00:00
Pratheek Shanthraj
07ff1f7c93
reporting current time for most module inits
2013-02-25 16:34:59 +00:00
Martin Diehl
393c60565d
one line was beyond 132 characters
2013-02-22 05:38:02 +00:00
Martin Diehl
42b96354db
doxygen comments
2013-02-21 23:08:36 +00:00
Martin Diehl
f0b4281400
removed unused variables, added pure statements, and declared external functions as external (using gfortrans debug option)
2013-02-11 10:43:45 +00:00
Martin Diehl
e74b5da19a
changed all remaining routines to fortran-fast arrays (geometry reconstruction etc.)
...
changed all remaining routines in f2py to more clever determination of array size (requires f2py >= 2.0)
enabled 3D visualize to work with odd resolution by switching to linear reconstruction
PLEASE NOTE: Redefinition of routines for f2py might cause trouble -> DELETE DAMASK_ROOT/lib/damask/core.so in this case
further changes: added pure statement where possible, polished, unified use of "Q" for "Quaternion" and reordered math to have similar routines together
2013-01-31 16:28:08 +00:00
Martin Diehl
3ffd6499b1
change bash-only "let" to "`expr`", abaqus_v6.env now suppresses warnings about long lines (there are only comments anyway)
2013-01-28 15:30:51 +00:00
Pratheek Shanthraj
fd94c786f0
moved stiffness tensor calculation to lattice
...
introduced 'isotropic' and 'orthorhombic' lattice types to use corresponding symmetries in stiffness tensor. intended to be used with non-crystal plasticity models (j2, constitutive_none with isotropic, cubic or orthotropic elasticity).
2013-01-21 21:57:26 +00:00
Christoph Kords
8a45a90775
new crystallite output "grainrotationX", "grainrotationY", "grainrotationZ" that gives the deviation of the initial orientation in degrees around the sample reference x/y/z axis, while "grainrotation" still gives the same as axis angle pair in crystal coordinates
2013-01-16 08:45:41 +00:00
Christoph Kords
7dd1130e92
reset syncSubFrac flag when turning terminallyIll
2012-12-16 16:18:36 +00:00
Christoph Kords
92d2c83334
more debugging info for timesyncing
2012-12-16 11:30:19 +00:00
Christoph Kords
4f12073dcb
once again corrected parallelization statements
...
added some debugging writes for timesyncing
2012-12-16 10:54:13 +00:00
Christoph Kords
f744fc1d47
fixed bug in parallelization: flushing of crystallite_todo was missing at a couple of places
2012-12-11 15:05:08 +00:00
Christoph Kords
6a9b19e4c0
fixed bug in timesyncing procedure: if any ip that synchronizes its time step did not converge, all nonlocals become terminally ill; before, this led to an infinite loop in crystallite
2012-12-11 13:42:22 +00:00
Christoph Kords
35a0483507
fixed bug in last revision: FE_geomType was not in namespace
2012-12-10 17:22:16 +00:00
Philip Eisenlohr
10ed848258
fixed 2 bugs: missing mapping with FE_geomtype(mesh_element(2,...))
2012-12-10 17:03:28 +00:00
Christoph Kords
98d02dcffb
removed call to microstructure in stateJump function, since this is already done in the state integrator subroutine
2012-12-07 07:50:41 +00:00
Christoph Kords
0986a66449
accidentally also included some changes in crystallite_integrateStress in last commit; reverted these changes in this revision
2012-12-03 16:34:41 +00:00
Christoph Kords
f2dacb5f40
polished debugging output
2012-12-03 16:29:54 +00:00
Christoph Kords
f92d34dd22
removed some unused variables
2012-11-29 18:44:00 +00:00
Christoph Kords
e10000a338
Introduced an alternative cutback scheme for the nonlocal model, that allows to keep the results of most of the ips that immediately converged, and only do cutbacks in regions where some ips are in big trouble. Only works for nonlocal model and state integrator 2.
2012-11-27 18:36:55 +00:00
Christoph Kords
e584c0e6e8
deleted space in front of preprocessor line
2012-11-22 20:04:33 +00:00
Christoph Kords
38ec462f4a
order of integration now same for all state integrators: first calculate the dotState and then the deltaState
2012-11-22 13:04:19 +00:00
Philip Eisenlohr
7d196fbb25
polishing
2012-11-22 09:58:36 +00:00
Philip Eisenlohr
d189a6a7a0
polishing
2012-11-21 16:57:57 +00:00