Commit Graph

349 Commits

Author SHA1 Message Date
Martin Diehl 70166657a9 FILEUNIT still needed 2018-06-01 21:25:08 +02:00
Martin Diehl cd9fb1c544 adjusting variable names in isotropic to DAMASK paper 2018-06-01 21:23:25 +02:00
Martin Diehl 64270f6c95 no need to store trailing whitespace 2018-06-01 21:22:19 +02:00
Martin Diehl f843254233 Noutput not needed, using size ensures consistency 2018-06-01 20:54:40 +02:00
Martin Diehl c8cec5a121 some more functions to parse material.config from memory 2018-06-01 14:14:53 +02:00
Martin Diehl cfefcaebb8 simplifying 2018-06-01 11:33:13 +02:00
Martin Diehl 48851c00f6 sizePostResults not needed, transpose33 has no advantage over transpose: 2018-06-01 11:24:00 +02:00
Martin Diehl 5de4b47438 plastic isotropic parses material.config from main memory, not from file 2018-06-01 10:24:42 +02:00
Martin Diehl 052a0af2ac storing phase part of material.config per section in chained list 2018-06-01 10:14:46 +02:00
Martin Diehl 8867322713 data structure to store material.config
chained list data structure to store string and stringPos.
Performance-wise not 100% optimal as searching requires to parse all elements in the list.
However, secions are typically rather small (usually Order 1 to 3) and parsing from file is much slower.
Also, only the actual length of the array is stored (thats wy I preferred the list over a 2D array)
2018-06-01 10:14:12 +02:00
Martin Diehl 60e30e4a64 Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development 2018-05-28 07:03:56 +02:00
Martin Diehl e661a42a4c citation block should be separated by empty lines 2018-05-28 07:01:50 +02:00
Martin Diehl 715f2a59d0 Merge branch 'development' into release_bugfix 2018-05-27 14:10:38 +02:00
Martin Diehl 87a16b775e function to change working directory 2018-05-27 14:08:23 +02:00
Martin Diehl 39e2e8a305 len is length of string, here size is required 2018-05-27 14:07:34 +02:00
Martin Diehl e0e97abda4 citation does not match 2018-05-27 08:16:05 +02:00
Martin Diehl 60cdffacd6 integer kind need to match for array constructor
hence, a=[a,b] requires a and b of the same kind
2018-05-27 08:09:32 +02:00
Martin Diehl 466cb57d7c Merge branch 'development' into release_bugfix 2018-05-26 16:37:39 +02:00
Martin Diehl 6891a48c82 wrong label 2018-05-26 12:35:02 +02:00
Martin Diehl 896a502099 Merge branch 'development' of magit0.mpie.de:damask/DAMASK into development 2018-05-26 12:09:09 +02:00
Martin Diehl 0172987115 forgot to update default 2018-05-26 12:08:28 +02:00
Martin Diehl bcaab1c068 restored working behavior 2018-05-26 12:07:13 +02:00
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
Arka Lahiri c45d9f9ee9 Deleted performance profiling code 2018-05-17 16:33:35 +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
Arka Lahiri 4d393518ac did only compile with BUILD_TYPE=DEBUG 2018-05-09 16:54:06 +02:00
Arka Lahiri 90e9b98584 performance profiling should not be part of the source code 2018-05-09 16:35:09 +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