Pratheek Shanthraj
9779127521
added intermediate configuration to crystal plasticity kinematics intended handle intermediate deformation stages consistently in a finite strain framework
...
currently implemented for thermal strains, but will expand on this to add damage strains, phase transformation strains etc.
2014-10-20 15:43:28 +00:00
Pratheek Shanthraj
96f036e34e
updated MPI reporting in line with recent changes
2014-10-10 16:21:10 +00:00
Pratheek Shanthraj
ba65d04bb3
removed redundant crystallite_temperature/heat and materialpoint_heat and using temperature from the thermal module
2014-10-10 12:28:57 +00:00
Pratheek Shanthraj
c8929b47d3
only output from the root processor for parallel runs
2014-10-09 20:23:06 +00:00
Pratheek Shanthraj
c7a746ddc1
constitutive dot state now takes Lp which is needed to the thermal dot state.
...
thermal problems should be working now
2014-09-26 20:49:25 +00:00
Pratheek Shanthraj
5c3f4e3e69
some fixes to get homogenisation none working with recent changes
2014-09-26 15:16:10 +00:00
Pratheek Shanthraj
c2c5c0c521
subsumed constitutive_damage/thermal into constitutive
2014-09-23 10:38:20 +00:00
Pratheek Shanthraj
6ace12be0c
removed unused variables
2014-09-22 20:34:42 +00:00
Luv Sharma
e83a0fb3f7
polished field state related changes and brittle damage
2014-09-22 18:15:19 +00:00
Martin Diehl
7cb7815176
polishing style
2014-09-13 10:04:44 +00:00
Luv Sharma
dc406a01c0
added funtions to get averged properties at integration points.
2014-09-10 08:37:12 +00:00
Franz Roters
b6847ebb01
added use mesh statement needed for output of mesh_element
2014-09-09 15:18:49 +00:00
Philip Eisenlohr
e8a428613c
general polishing; simplification of RKCK45 dotstate assembly for different stages.
2014-09-02 20:11:57 +00:00
Philip Eisenlohr
630d9efffd
Heavy bug fixing: Butcher tableau multiplication in RK45 integration scheme, state updating of adaptive Euler, wrong assignment to evolution rate instead of state in fixed-point iterator.
...
Exchanged possibly non thread-safe matmul in RKCK45, simplified RK45 integration step 1--3 vs 4 handling.
(Thanks to Pratheek, Luv, and Chen for their help!)
2014-09-02 19:46:52 +00:00
Philip Eisenlohr
a360f14581
debug now reporting elFE in addition to elCP
2014-08-27 15:54:11 +00:00
Chen Zhang
5bc49e6441
change the debug level of cutback reporting to non-selective, of central solution from basic to extensive, of stiffness perturbation to selective.
...
fixed few typos and deleted excess whitespace.
2014-08-26 14:44:32 +00:00
Pratheek Shanthraj
7666013758
use subF instead of F for analytic jacobian calculation
2014-08-15 09:04:00 +00:00
Luv Sharma
33b4c0d908
minor typo
2014-08-11 13:33:17 +00:00
Pratheek Shanthraj
8f444b0695
minor changes to multiphysicsOut
2014-08-10 11:14:43 +00:00
Pratheek Shanthraj
ec71d77038
corrections to analytic dPdF calculation
2014-08-10 10:27:35 +00:00
Pratheek Shanthraj
da5f3d8f85
removed some unused variables
2014-08-08 13:54:08 +00:00
Pratheek Shanthraj
2f5d9376dc
call IO_error if inversion error in analytic jacobian calculation
2014-08-07 21:16:37 +00:00
Pratheek Shanthraj
98a297cd9e
removed old inaccurate analytic jacobian with exact jacobian. seems to work as well as perturbed jacobian and better than the old analytic jacobian for some simple tests.
2014-08-07 21:08:34 +00:00
Martin Diehl
a787d66763
updated restart test, deviations seems to be related to tolerance in stress BC
2014-08-04 17:50:01 +00:00
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