Commit Graph

1498 Commits

Author SHA1 Message Date
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
Pratheek Shanthraj f3ea92f50c fixed error in previous commit 2013-12-18 09:38:44 +00:00
Pratheek Shanthraj 61981617d7 separated forwarding and solution subroutines for better control at the load step looping level 2013-12-18 09:35:05 +00:00
Pratheek Shanthraj ff6211b78c rolled back phase field changes 2013-12-18 09:09:32 +00:00
Martin Diehl 6fa9ed8f48 homogenization: added enums and sourced allocation for RGC, some higher verbosity for isostrain 2013-12-18 07:28:01 +00:00
Martin Diehl c4a592867f updated to PETSc 3.4.3, will not compile on all workstations except for maws0X 2013-12-17 15:37:14 +00:00
Martin Diehl 676e621af2 new version of Lib_VTK_IO 2013-12-17 13:24:34 +00:00
Martin Diehl 7b9d76cdcc added missing allocation of outputID caused crashing 2013-12-17 08:02:21 +00:00
Martin Diehl ba89df7784 continued with sourced allocation and EOL improvements 2013-12-16 11:58:03 +00:00
Martin Diehl 118dff9f23 fixed wrong initialized symmetry value 2013-12-16 10:56:56 +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 7ba3cc0ab7 removed symlink to marc 2010 2013-12-13 10:04:25 +00:00
Martin Diehl acad3d4a33 using sourced allocation 2013-12-12 23:03:37 +00:00
Martin Diehl da337962af fixed bug for Abaqus input parsing 2013-12-12 22:29:40 +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
Christoph Kords 9b9f4dd624 adapted changes from rev 2776 for constitutive_nonlocal: improved usage of "enum", fixed bug in file reading, usage of "source" specifier for "allocate" 2013-12-12 10:16:50 +00:00
Christoph Kords dd77690a68 fixed bug introduced in rev2777: renamed all remaining occurrences of variable "ncomponents" to "nconstituents" 2013-12-12 10:06:39 +00:00
Martin Diehl 102712d91f added enum for dislotwin output, fixed bug when using recursive file input function 2013-12-11 23:42:33 +00:00
Martin Diehl 95d6430b09 improved on enums (introduced them for output in j2), fixed wrong recursion depth bug in IO, fixed a bug in reading in for none, j2, phenopowerlaw 2013-12-11 22:03:09 +00:00
Martin Diehl 7885ebaf8f added 'reset' flag for recursive function IO_read, need to reset stack when stopping exectution (will be used by constitutive and homogenization) 2013-12-11 16:49:20 +00:00
Franz Roters efa7a6dae2 added example material.config again 2013-12-05 15:11:30 +00:00
Franz Roters bf6e749abd adopted deployMe.sh to new svn server
editorial changes in LICENSE and CPFEM.f90
2013-12-05 14:31:28 +00:00
Christoph Kords 61ebfb4c85 removed obsolete output 'boundarylayer' from material.config 2013-11-29 12:30:06 +00:00
Christoph Kords b0695c0904 enumerators for case switch of outputs as already used in phenopowerlaw
removed obsolete output 'boundarylayer'
2013-11-29 12:29:47 +00:00
Christoph Kords f56cdafccb polishing indentation levels 2013-11-29 12:25:36 +00:00
Martin Diehl 2d624d022f fixed small bug in detecting the lattice structure. only first 3 letters count (e.g. 'ort' and 'orthography' will give the orthorombic lattice type) 2013-11-28 08:56:02 +00:00
Martin Diehl 345e5cd47c split up material.config and added a small helper script to create wiki content from all material.configS
(Work done by Luv, just his TortoiseSVN is not yet installed)
2013-11-27 17:49:34 +00:00
Martin Diehl 0bc99a9622 change latticeName to latticeID 2013-11-27 16:20:27 +00:00
Martin Diehl 6c2ab1f1f8 added ID (integer, enum) case switching for phenopowerlaw and isostrain 2013-11-27 11:39:28 +00:00
Martin Diehl 080edffd87 file was missing in last commit 2013-11-27 08:05:23 +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
Martin Diehl c9ce8d2c52 some improvements on HDF5 test, polishing math 2013-11-22 10:31:52 +00:00
Christoph Kords e8de8f69a8 renaming of "maxNinstance" to be consistent with naming scheme introduced in rev2662 2013-11-21 13:40:14 +00:00
Christoph Kords 3cc6781eaf renaming of "instance" and thelike to be consistent with naming scheme introduced in rev2662 2013-11-21 13:35:43 +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
Martin Diehl 9083aa53fb fixed long lines due to makro expansion 2013-11-19 16:06:53 +00:00
Pratheek Shanthraj 9c34c1b871 error reporting for phase fields only if there are active phase fields 2013-11-14 10:36:49 +00:00
Pratheek Shanthraj f51c600189 convergence check for phase field only if there are active phase fields 2013-11-14 10:32:41 +00:00
Pratheek Shanthraj 3b113f6a62 more errors :) 2013-11-13 19:38:25 +00:00
Pratheek Shanthraj 8444b4116d corrected error in previous commit 2013-11-13 19:26:36 +00:00
Pratheek Shanthraj 50db944c0c reworked phase field interface to damask spectral solvers. now specify 'thermal a b c' or 'fracture a b c' to activate either phase field where a b c are the initial value, diffusion coefficient and mobility respectively.
Right now only thermal and fracture phase fields implemented and only in the basic petsc solver
2013-11-13 19:21:35 +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 9ac8846e41 compiling with hdf5 is now working for DAMASK (if HDF5=ON) and for playground example in private 2013-11-07 12:52:36 +00:00
Philip Eisenlohr c857381881 added "-d", "--directory" as possible cmd switches for working directory.
fixed (sometimes?) non-working check for valid working directory path.
2013-10-24 18:56:56 +00:00
Martin Diehl ea980f428e fixed long line coming from PETSc macro expansion 2013-10-23 17:52:58 +00:00
Pratheek Shanthraj 57cf472982 Improved convergence checking when using newton solver (does not oversolve the problem like before) 2013-10-23 15:22:12 +00:00
Philip Eisenlohr eb0d2b7e24 fixed bug that (sometimes, always?) prevented reading in of debug.config, when variable 'line' was still '#EOF#' from former IO_read... 2013-10-23 11:23:41 +00:00
Philip Eisenlohr bdafc6f1b8 condensed error msg numbering,
introduced new error for "microstructure without any constituent"
2013-10-23 11:21:48 +00:00
Martin Diehl 205eb66ba8 bugfix for last commit, too many variables declared as private 2013-10-23 10:29:38 +00:00
Martin Diehl 96b03841e8 redefined macro for intel compiler, should work with windows (but Nicolo had some problems)
made configure verbose
2013-10-23 08:38:14 +00:00
Martin Diehl 7b28dedfe3 removed temperature calculation and debugging statements 2013-10-22 15:52:58 +00:00
Martin Diehl ea664688f8 introduced dummy temperature calculation. 2013-10-18 20:56:10 +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 1c7ed337d7 corrected linebreak in comment 2013-10-17 12:52:46 +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
Martin Diehl 6a1c40d540 introduced error when element/IP/component(grain) selected for debugging does not exist 2013-10-16 12:38:00 +00:00
Philip Eisenlohr 0a292d8d8b fixed final bug in hex twintwin matrix 2013-10-15 22:02:39 +00:00
Franz Roters 4d184b1ee4 improved comments 2013-10-14 14:35:41 +00:00
Martin Diehl 295d0cd28e removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC 2013-10-14 10:54:45 +00:00
Martin Diehl 029502c12e fixed typo causing compilation error 2013-10-14 06:11:50 +00:00
Martin Diehl b651f334fe set useful default values to remove dummy functions in individual homogenization schemes (stateInit, averageTemperature, and updateState (isostrain only)) 2013-10-11 16:01:53 +00:00
Franz Roters 80cb78c1fd corrected activation of slipbands; both velocity and resistance have to be non zero 2013-10-11 13:52:00 +00:00
Martin Diehl 3b5a5ab812 changed error in case of number of found interaction types is less than expected to warning 2013-10-11 09:17:03 +00:00
Martin Diehl bb0b35e466 made compilation_info.f90 windows compatible and added abaqus_v6_windows.env as a template for using DAMASK with windows.
Deleted old fftw and ACML installation instruction (can be found on damask.mpie.de)
2013-10-10 09:38:31 +00:00
Martin Diehl ed47d25adc remove dotTemperature (returned always 0 anyway) to make it easier to replace it with a heating rate funtction
remove potentially dangerous short circuiting in homogenization files
added doxygen comment for nonlocal (giving the functions a name)
fixed wrong definition of maxnchungs introcuded into dislotwin with last commit, causes trouble with hex
2013-10-09 06:12:16 +00:00
Martin Diehl 82faf74363 added doxygen documentation and unified variable names and some common parts of the code 2013-10-08 16:27:26 +00:00
Philip Eisenlohr 0a7052a7da fixed two errors in hexagonal slipslip and twintwin interaction matrix 2013-10-08 13:52:19 +00:00
Philip Eisenlohr 42cd4eb38f fixed some e12.6 string formats to avoid compile time warnings 2013-10-08 13:51:36 +00:00
Philip Eisenlohr a2169e338f fixed type-casting problems (mainly with Lib_VTK_IO) 2013-09-26 17:21:46 +00:00
Martin Diehl ca2a8d0a03 fixed typos in Abaqus related functions 2013-09-23 09:01:46 +00:00
Martin Diehl 595c8860a2 renamed myInstance -> matID, myStructure -> structID to be consistent with other constitutive models 2013-09-20 17:03:11 +00:00
Martin Diehl ad32aa3ced fixed calculation of CurlRMS when grid(1)==1 2013-09-20 16:17:25 +00:00
Martin Diehl a0f304021f fixed restarting for Polarisation and AL, updated corresponding test and polished output in python (some statements are screwed up when making them Python3 compatible) 2013-09-20 14:22:37 +00:00
Martin Diehl 40c399bf2e fixed bug in last commit and remove one more goto marker 2013-09-19 15:28:55 +00:00
David Mercier ce926aaf1f Minor contribution :
- add SVN keyword to the material.config file for cp-Ti
- add comments in lattice file
2013-09-19 09:17:09 +00:00
Martin Diehl ca0b3b7b81 changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw 2013-09-19 07:46:01 +00:00
Martin Diehl 70c888f6ea fixed a typo in documentation 2013-09-18 19:54:39 +00:00
Martin Diehl 56d21286a6 fixed comment on values changed automatically 2013-09-18 14:33:40 +00:00
Martin Diehl 3eec0ecdcd introduced #EOF# as module wide parameter in IO.f90
renamed read/write binary file to read/write real/int
removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 14:07:55 +00:00
Martin Diehl 19a353edb4 made error in case of invalid precision more verbose 2013-09-18 13:59:42 +00:00
David Mercier a8eb6a985e - Addition of an example of a material.config file for hcp material (cp-Ti).
- New Error Message #214 when elastic constants are not correctly defined.
- Error message when the c/a ratio for hcp is higher than 2 (no physical sense).
2013-09-17 12:20:50 +00:00
Martin Diehl cb16b30409 set output order of spectral parameters to more reasonable one 2013-09-17 09:19:16 +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
Christoph Kords 5f973a21c0 require all hardening coefficients/ non schmid coefficients to be specified in material.config, do not silently ignore any
warning message now more meaningful (?) when using less than maximum number of slip families
2013-09-12 14:47:09 +00:00
Martin Diehl ddfc40afe6 improve restart test, deletion of results was messed up. fixed restart for BasicPETSc, AL and Polarization currently don't run when restarting 2013-09-05 12:07:58 +00:00