Commit Graph

324 Commits

Author SHA1 Message Date
Martin Diehl 2fbe60b949 anticipate (proper) change in 23_BasticPETSc_2_PETSc as PRIVATE repository is ahead 2018-05-26 09:54:50 +02:00
Martin Diehl c7c39922f0 pointer assignment was done twice
pointer is re-assigned automatically, but I found it confusing.
Also using automatic left hand side reallocation to simplify handling of outputID
2018-05-26 09:44:14 +02:00
Martin Diehl 8184d51a99 following style of more complex constitutive laws
offset for different states needs to be computed, so it makes sense to define global and local aliases together.
no need to introduce variables for state0 and aTolstate, they are only used once
2018-05-26 09:31:36 +02:00
Eureka Pai c09a7fd157 replaced param(instance) with p => pointer, corrected errors 2018-05-24 18:31:32 -04:00
Eureka Pai 982c0fb90a replaced param(instance) with p => pointer, corrected error 2018-05-24 18:24:58 -04:00
Eureka Pai f97800658f replaced param(instance) with p => pointer 2018-05-24 17:56:09 -04:00
Jaeyong Jung c79a0d39c0 Merge branch 'development' into 23_BasticPETSc_2_PETSc 2018-05-24 15:53:00 +02:00
Martin Diehl 5bf6ede661 polishing 2018-05-24 07:57:10 +02:00
Martin Diehl 1bad719abe undefined variables can be removed, style as in other plastic modules 2018-05-24 07:51:43 +02:00
Martin Diehl f0309ac4dd Bessel functions not used 2018-05-24 07:51:16 +02:00
Martin Diehl 1c75a2e9cd using syntax with better error handling 2018-05-24 07:13:12 +02:00
Martin Diehl 0ee34d608c fixing all appaerances 2018-05-24 00:06:03 +02:00
Mingxuan Lin bf7efad250 Fix errors in CMake script 2018-05-23 21:46:03 +02:00
Martin Diehl 30fe87a88a post-release cleaning 2018-05-22 19:16:12 +02:00
Philip Eisenlohr 324df2eeba Merge branch 'PETSc3.9' into 'development'
PETSc 3.9 working

See merge request damask/DAMASK!25
2018-05-19 17:39:26 +02:00
Martin Diehl 5187ee25b6 destroying unused object causes problem with PETSc 3.9
not clear, why destroying at the end of the program is anyway needed, so removed it
2018-05-19 11:20:04 +02:00
Philip Eisenlohr d79d366355 Merge branch '20_modify_RGC' into 'development'
20 modify rgc

See merge request damask/DAMASK!24
2018-05-18 00:15:47 +02:00
Jaeyong Jung 6800e779cb second try! 2018-05-17 16:27:36 +02:00
Martin Diehl 537bb8df48 employing use statements (almost no external functions anymore)
adopted function calls to newer PETSc versions
2018-05-17 12:04:21 +02:00
Martin Diehl 8193e27e8d line break looked strange 2018-05-16 18:34:49 +02:00
Martin Diehl 31311d63db gfortran 8 was complaining, asigned pointer is not contiguous 2018-05-16 00:09:33 +02:00
Jaeyong Jung a65e4bdad4 first modification. compile succeeded 2018-05-14 15:45:12 +02:00
Martin Diehl 6b27460a3b Merge branch 'development' into 12-fixOrientationSampling 2018-05-05 19:03:51 +02:00
Martin Diehl 9072179c61 Merge branch 'development' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2018-05-04 16:07:15 +02:00
Jaeyong Jung 3f71288e0a Now, homogenization_RGC is fixed! 2018-05-04 14:09:49 +02:00
Zhuowen Zhao 76e6fa24fb added check to the "Plasticity_DetectChanges" test and removed commented-out parts for kinematic hardening 2018-05-03 22:39:49 -04:00
Jaeyong Jung 4c42510ad6 compiles without AL solver 2018-04-26 14:11:45 +02:00
Zhuowen Zhao 65f17f4166 merge development into kinematic-hardening branch 2018-04-25 14:14:17 -04:00
Zhuowen Zhao 7bc736aad7 killed modification not relevant to kinematic hardening 2018-04-25 13:26:38 -04:00
Yi-Chin Yang 35a4fdc358 "et al." instead of "et. al" abbreviates "et alii", meaning "and others" 2018-04-22 10:07:49 +02:00
Yi-Chin Yang e90861956d "et al." instead of "et. al" abbreviates "et alii", meaning "and others" 2018-04-22 09:32:59 +02:00
Sharan Roongta 488ff6bb2d Rectifying a typo 2018-04-18 15:30:54 +02:00
Sharan Roongta 71d23274e7 A citation on twinning was missing 2018-04-18 15:22:19 +02:00
Sharan Roongta 451e9f245a Removing some typos in citations 2018-04-18 15:01:03 +02:00
Sharan Roongta 90a1af5096 Cited HybridIA algorithm 2018-04-18 14:18:48 +02:00
Sharan Roongta e77e0b6f20 Merge branch '17-citations-in-source-code' of magit1.mpie.de:damask/DAMASK into 17-citations-in-source-code 2018-04-17 15:54:20 +02:00
Sharan Roongta 3fca8da6aa Citations in plastic_twin 2018-04-17 15:51:59 +02:00
Yi-Chin Yang 4f5b5d9d4a Merge branch '17-citations-in-source-code' of magit1.mpie.de:damask/DAMASK into 17-citations-in-source-code 2018-04-17 15:42:45 +02:00
Yi-Chin Yang b7d622c32d No tab allowed 2018-04-17 15:39:08 +02:00
Yi-Chin Yang 258be943c7 No tab allowed 2018-04-17 15:25:27 +02:00
Yi-Chin Yang fc3abd4f39 References for spectral solvers 2018-04-17 15:22:18 +02:00
Yi-Chin Yang 41967da684 edit 2018-04-17 15:09:12 +02:00
Yi-Chin Yang d6d9d6e1bb edit 2018-04-17 15:07:16 +02:00
Yi-Chin Yang 86683674c0 mod 2018-04-17 14:52:50 +02:00
Sharan Roongta 3bd09fe837 fixed typo 2018-04-17 14:43:38 +02:00
Yi-Chin Yang afcaafa798 Reference for homogenization_RGC added 2018-04-17 08:20:41 +02:00
Yi-Chin Yang 0701b535fc DAMASK overview reference independently of selected solver 2018-04-16 16:16:10 +02:00
Martin Diehl ea9434432e added reference for the tungsten model 2018-04-09 15:04:58 +02:00
Martin Diehl 8509b13f16 avoid numerical noise 2018-03-13 09:21:41 +01:00
Martin Diehl 302a5017d4 Merge branch 'development' into 12-fixOrientationSampling 2018-03-13 04:12:03 +01:00
Martin Diehl 2ceb3aeea0 not all variables were consistently renamed 2018-03-13 03:56:02 +01:00
Martin Diehl 44e5644e78 fixed random Gaussian sampling
sampling needs to be performed from unfiform misorientation, NOT uniformly distributed rotations
for Fiber, compute uniform tilt of Fiber axis
2018-03-13 03:36:24 +01:00
Zhuowen Zhao 361cf680a5 merged more recent development 2018-03-09 18:19:41 -05:00
Zhuowen Zhao f26fd1d1dc Fixed a physics issue: sense change with respect to difference between resolved stress and backstress 2018-03-09 18:01:09 -05:00
Martin Diehl 0ccce7facc fixed missing angle initialization, simplified and commented 2018-03-08 17:33:54 +01:00
Philip Eisenlohr 88c46d500e Merge branch '13-higher-order-thermal-expansion' into 'development'
Resolve "higher-order thermal expansion"

Closes #13

See merge request damask/DAMASK!19
2018-02-28 21:08:01 +01:00
Martin Diehl e716441787 indexing now follows description 2018-02-26 20:18:38 +01:00
Martin Diehl cd4115bada now using the correct rotation formula 2018-02-26 11:14:02 +01:00
Martin Diehl e1c4c7e6e8 wrong indices in matrix rotation
fixing this resolves issue 14: "Poor performance after load case rotation"
now inline with http://www.continuummechanics.org/rotationmatrix.html
C'_mnop = R_mi R_nj R_ok R_pl C_ijkl
2018-02-25 21:01:23 +01:00
Philip Eisenlohr bc998846ab thermal expansion coefficient alpha = a + b(T-T_0) + c(T-T_0)^2 2018-02-25 14:14:03 -05:00
Martin Diehl c2ced6a274 not needed anymore 2018-02-25 14:45:35 +01:00
Martin Diehl d5d71d2321 math_transpose33 has no benefit here 2018-02-25 13:42:02 +01:00
Martin Diehl 42500eabff simplified (math_transpose33 does not provide any benefit here) 2018-02-25 12:48:58 +01:00
Martin Diehl eedec74a6a correct debug reporting 2018-02-25 12:31:23 +01:00
Martin Diehl 86768c38d9 source code should not be executable 2018-02-25 10:16:27 +01:00
Martin Diehl 4f4fa5daf8 simplified 2018-02-22 14:54:49 +01:00
Martin Diehl e51de7ffd8 explicitly select halton bases in call 2018-02-22 14:16:36 +01:00
Martin Diehl c6c66bb653 using proper Gauss sampling also for the fiber components
function most probably still contains a bu
2018-02-22 00:42:06 +01:00
Martin Diehl 5c908e44ec correct scaling for FWHM
FWHM was wrong by a factor of approx 2.
Analytic expression can be found in https://en.wikipedia.org/wiki/Gaussian_function
2018-02-22 00:10:38 +01:00
Martin Diehl 9173d12d14 correct algorithm for sampling of uniform orientations and fix for Halton series
Halton series gives strange results for large prime numbers, now always starting with 2 for first dimension, 3 for second etc.
Consecutive Halton numbers for rejection sampling seem to cause problems (i.e. introduce patterns).
Algorithm for uniformly distributed orientations with FWHM specified is taken from https://math.stackexchange.com/questions/131336.
WIP: Gauss filtering is currently not implemented!
2018-02-21 20:32:52 +01:00
Martin Diehl ae27660e86 simplified halto procedure (still needs testing) 2018-02-21 18:47:39 +01:00
Martin Diehl 98df2d1427 better description and names 2018-02-21 16:49:29 +01:00
Martin Diehl a52f54a9a0 fixed typo in prime number list and extended to 1600 values
source https://people.sc.fsu.edu/~jburkardt/f_src/halton/halton.f90
2018-02-21 14:37:02 +01:00
Martin Diehl 6b1c1af609 no early return required here 2018-02-21 08:55:33 +01:00
Martin Diehl 051f3592ea Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2018-02-21 07:53:00 +01:00
Philip Eisenlohr ad1179117d modified meaning of restart increment to "restart from" instead of "restart at" 2018-02-16 16:41:07 -05:00
Philip Eisenlohr 0a763ff116 included changes to correct cutback issue of spectral solver 2018-02-16 09:36:18 -05:00
Martin Diehl 9b1e72e7bc more sensible cutoff for scatter 2018-02-15 11:26:19 +01:00
Martin Diehl 230ffe8f35 Bessel functions for Gaussian scatter 2018-02-14 13:03:50 +01:00
Martin Diehl c544fa1b34 typo, cause wrong finalization sequence 2018-02-12 15:03:35 +01:00
Martin Diehl fe07a34c69 source code should not be executable 2018-02-12 12:42:18 +01:00
Martin Diehl 4e62d7793a more detailed comments 2018-02-12 11:46:01 +01:00
Martin Diehl b6abbfca9d early return not needed here 2018-02-12 11:45:42 +01:00
Martin Diehl 2254746177 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2018-02-12 11:23:05 +01:00
Zhuowen Zhao 7f487bb77b merge development into kinematic hardening branch 2018-02-07 11:30:57 -05:00
Fengbo Han 208c4affa4 using math_equivStrain33 and math_equivStress33 instead of Mises 2018-02-07 17:11:43 +01:00
Fengbo Han b834b2a00d removed unnecessary lines 2018-02-07 13:37:26 +01:00
Fengbo Han 9249e7db2b Merge branch 'development' into HMS 2018-02-07 11:44:17 +01:00
Fengbo Han 190a2baf9f when using yield stop criteria, if rotation of the load frame is specified, the output results in .yield and .stressstrain files are also rotated 2018-02-07 11:35:16 +01:00
Martin Diehl b455fcd0a3 acos has the same limits/stability as indirect asin calculation, enforing limits directly instead of scaling 2018-02-03 09:21:19 +01:00
Martin Diehl 9908734710 wrong random number, simplified 2018-02-03 09:04:49 +01:00
Martin Diehl ba6fd70654 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2018-02-03 09:02:20 +01:00
Franz Roters d80a255736 new Marc2017 file format finally working! 2018-02-02 15:06:13 +01:00
Martin Diehl d80e15bd76 compiler_options and compiler_version supported by new Intel compiler 2018-02-02 12:36:09 +01:00
Zhuowen Zhao 700edc313f Merge remote-tracking branch 'origin/spectralSolver-cutbackfix' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2018-01-18 18:18:34 -05:00
Philip Eisenlohr 307aa7e7d1 fixed forgotten file writing of F_aimDot 2018-01-18 18:16:16 -05:00
Zhuowen Zhao 7c755a0241 Merge remote-tracking branch 'origin/spectralSolver-cutbackfix' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2018-01-18 13:59:42 -05:00
Philip Eisenlohr b36151cc32 fixing spectral cutback hiccup and multiple cleanups
flush(6) at better places, added dedicated CPFEM_age subroutine, cleaned up cutback logic, fixed broken assignment of old timeinc, continueCalculation is now a logical, rearrnaged interfaces for utilities_constitutiveResponse and utilities_calculateRate, handling of stressBC more understandable, added more comments and explanations
2018-01-18 11:14:06 -05:00
Philip Eisenlohr 93073ed661 summarized multiple logicals into one 2018-01-18 10:47:52 -05:00
Franz Roters 2c4df2f070 working further on Marc2017 compatibility 2018-01-11 17:11:03 +01:00
Franz Roters 7149f9599f changes towards supporting new Marc2017 input file format
still not working
2018-01-10 17:13:25 +01:00
Zhuowen Zhao bc9c647aad assign outputID type to undifined 2018-01-09 17:25:16 -05:00
Zhuowen Zhao 8837656870 added "plastic_kinematichardeing.f90" to commercialFEM_fileList
Changed outputID type
2018-01-09 15:17:51 -05:00
ChuanlaiLiu e02a0d32be debug corrected index for stiffness tensor, phase not instance 2017-12-19 12:59:10 +08:00
Zhuowen Zhao 14c0503a7e change tag to all lowercase; make pheno hardeing only depend on its own system 2017-12-14 16:51:58 -05:00
Zhuowen Zhao dcf9e139d0 question marks on those files 2017-12-13 19:18:45 -05:00
Fengbo Han d81870dc57 output the stress-strain curve to file if yield stop criterion is used 2017-11-22 09:02:35 +01:00
Fengbo Han 2b5a536458 calculate stress and strain from the average PK stress and average deformation gradient of the whole RVE 2017-11-22 08:52:48 +01:00
Martin Diehl 8d705522ea output already disabled for worldrank !=0 2017-11-21 15:10:48 +01:00
Martin Diehl 09a66d918d (in)equality comparison for double was far too tolerant 2017-11-21 09:24:25 +01:00
Martin Diehl 37e154de65 preventing division by zero 2017-11-20 23:43:06 +01:00
Martin Diehl c81a438546 SCHMID-BOAS notation for reference 2017-11-18 22:39:13 +01:00
Zhuowen Zhao 2b8baa2f01 fixed backstress rate of change based on its own slip system evolution and its sense 2017-11-14 12:25:55 -05:00
Philip Eisenlohr e4700cda25 changed fixed_seed to random_seed for clarity 2017-11-07 14:56:28 -05:00
Zhuowen Zhao 2caf8b7ffd delete extra variable 2017-11-06 22:10:04 -05:00
Zhuowen Zhao ab11830e71 Merge branch '6-debug-replaces-openmp-preprocessor-switch' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2017-11-06 18:20:51 -05:00
Zhuowen Zhao 566cf3d855 Merge branch 'development' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2017-11-06 18:11:47 -05:00
Zhuowen Zhao 2b4a024671 renamed "tau0" to "crss0" 2017-11-06 18:11:02 -05:00
Zhuowen Zhao 23f9f03ca8 refined debug output 2017-11-06 18:09:04 -05:00
Philip Eisenlohr 426b2d0537 Merge branch 'development' into 6-debug-replaces-openmp-preprocessor-switch 2017-11-06 11:39:58 -05:00
Fengbo Han d599aec37f Merge branch 'development' into HMS 2017-10-24 11:24:54 +02:00
Fengbo Han 82758bd90f added yield criterion of facet potential 2017-10-24 11:15:34 +02:00
Martin Diehl 1c2af7bbc6 phenoplus has own brach, titanmod was not used for a long time 2017-10-08 13:41:54 +02:00
Martin Diehl eda8c24b7d Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2017-10-07 17:24:59 +02:00
Martin Diehl e32b9d9ca8 for comparison with de-facto stardard rotation definitions 2017-10-07 13:18:42 +02:00
Franz Roters 14bf4d5c7b added Marc 2017 link 2017-10-05 17:57:21 +02:00
Martin Diehl e9440ccd6f suppress warning when compiling with gfortran 2017-10-05 16:35:34 +02:00
Franz Roters 5011e201d6 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2017-10-04 17:18:57 +02:00
Franz Roters 4d4a546e2d hot fix for so far non-supported input file format of Marc 2017 2017-10-04 17:18:06 +02:00
Zhuowen Zhao 2f5d81e111 beautified debug output 2017-10-03 17:18:34 -04:00
Philip Eisenlohr 483ed745e0 replaced #ifndef _OPENMP with #ifdef DEBUG to trigger inclusion of parallelized debug statements; fixed minor bugs in debug output 2017-10-03 09:20:53 -04:00
Philip Eisenlohr 05bb3c109c Merge branch '5-flexible-dot-delta-state-relation' into 'development'
added offsetDeltaState to type-definition and respect this when applying deltaState

Closes #5

See merge request !15
2017-10-03 14:09:20 +02:00
Zhuowen Zhao 7d6d822ee0 Merge branch '5-flexible-dot-delta-state-relation' into 3-adding-plastic-constitutive-law-with-kinematic-hardening 2017-09-29 18:46:35 -04:00
Zhuowen Zhao b66a9a3527 Merge remote-tracking branch 'origin/development' into development 2017-09-29 18:40:00 -04:00
Zhuowen Zhao 81bcc72993 first shot at kinematic hardening constitutive law 2017-09-29 18:36:28 -04:00
Zhuowen Zhao d6cf3c4dd4 reduced size of aTolState from sizeState to sizeDotState
asb tolerance check is only meaningful for the state part affected by dotState
2017-09-29 18:35:36 -04:00
Zhuowen Zhao a153443239 clarified comment 2017-09-29 18:34:18 -04:00
Zhuowen Zhao fb4aadbafa moved deltaState to group of contigous memory pointers 2017-09-29 18:33:39 -04:00
Zhuowen Zhao 80bed8b8aa improved math_expand algorithm 2017-09-29 18:32:52 -04:00
Zhuowen Zhao f3292507b5 added more debugging messages and switched to DEBUG as flag 2017-09-29 18:32:07 -04:00
Philip Eisenlohr f5816a643d added offsetDeltaState to type-definition and respect this when applying deltaState 2017-09-29 17:44:10 -04:00
Martin Diehl e1bfde94f9 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2017-09-19 23:40:10 +02:00
Martin Diehl d8d42c32e7 cleaned + shortened, fixed handling of optional order parameter in math_exp33 2017-09-19 23:39:19 +02:00
Philip Eisenlohr 67ac5c7b84 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2017-09-19 08:13:08 -04:00
Philip Eisenlohr 330906d08f removed unneccesary detail from general error message, now caller gives specifics 2017-09-19 08:11:57 -04:00
Zhuowen Zhao 56f7836bf8 clarified composition of overall sizeState 2017-09-18 19:42:27 -04:00
Zhuowen Zhao ae868d3ada only one error for math_check; new "math_expand('what' by 'how')" 2017-09-14 15:25:22 -04:00
Martin Diehl f0d8d7c50a bug in eigenvalue calculation, caused wrong eigenvectors for sym33 matrices. Credits to Duancheng Ma 2017-09-06 16:20:24 +02:00
Philip Eisenlohr adc7c9d5b1 removed executable bit 2017-08-28 13:55:12 -04:00
Martin Diehl 7e0a9b75f9 simplified prime fuction 2017-08-12 06:35:44 +02:00