Commit Graph

2142 Commits

Author SHA1 Message Date
Su Leen Wong cf84ceddd8 Added Ndot0 for trans 2015-11-16 10:15:18 +00:00
Martin Diehl 302178ebb4 some trouble since clean was already in PETSc rules 2015-11-14 09:47:18 +00:00
Martin Diehl 9d9367ae11 still having trouble, incorporation of PETSC_ARCH not needed 2015-11-13 10:16:35 +00:00
Martin Diehl 4a52c59854 corrected/removed PETSc variables path 2015-11-12 23:13:25 +00:00
Philip Eisenlohr 900b072781 removed unused math_qDot (deprecated anyway) 2015-11-10 21:56:06 +00:00
Philip Eisenlohr 45fb209658 respect PETSC_ARCH 2015-11-10 21:18:53 +00:00
Philip Eisenlohr c5b28a22d4 small polishing 2015-11-10 21:18:13 +00:00
Su Leen Wong aefd8fa390 Added a few more trans dependent variables 2015-11-10 15:23:15 +00:00
Su Leen Wong 2778331d7b Added trans interaction matrices 2015-11-10 13:30:34 +00:00
Su Leen Wong 4e4136cd58 Consistent variable names 2015-11-10 13:01:03 +00:00
Su Leen Wong 9c91ba4608 Trans interaction matrices 2015-11-09 15:57:20 +00:00
Su Leen Wong 2b10625f0c Typo 2015-11-09 08:51:05 +00:00
Martin Diehl ff7cd17b51 not needed, PETSc is doing the job 2015-11-08 09:43:43 +00:00
Chuanlai Liu 4d9a2f8f6b introduce pointer referencing of plastic state for better readability 2015-11-06 17:00:00 +00:00
Martin Diehl 7b0c130d6f started to introduce alias name via pointer for plastic state in phenopowerlaw for more comfortable access 2015-10-30 15:48:30 +00:00
Su Leen Wong 2ae4d24fb3 Use rotated trans elasticity matrices 2015-10-27 12:32:03 +00:00
Martin Diehl 2ba15e462a corrected values and typo 2015-10-26 19:27:57 +00:00
Martin Diehl c91788d721 new intel compiler was compiling about *** in float output 2015-10-26 17:11:36 +00:00
Martin Diehl 1f356a6833 added parameters from indentation results 2015-10-22 17:45:28 +00:00
Su Leen Wong 04841cb727 List of dependent trans states 2015-10-21 13:24:26 +00:00
Su Leen Wong 20b0429517 Calculation of equilibrium separation of partials for twin and trans 2015-10-21 09:37:45 +00:00
Su Leen Wong 565d68edd7 Added some TRIP inputs 2015-10-20 15:54:01 +00:00
Su Leen Wong b9fcd9482b Renamed some input variables to denote twin 2015-10-20 13:45:06 +00:00
Franz Roters fc8a11f027 corrected typo: alias for first Piola stress should read 1stpiola NOT 1piola 2015-10-19 06:57:23 +00:00
Chen Zhang f2b29e4844 correct communication between
elements with different
phases
2015-10-14 21:04:40 +00:00
Chen Zhang 484a34b7f1 added pheno+ module
modify crystallite microstructure call
to pass orientations
2015-10-14 18:36:19 +00:00
Martin Diehl ea0b16c726 pheno plus for commercial FEM 2015-10-13 21:06:59 +00:00
Chen Zhang d6abc00218 add pheno+ module in 2015-10-13 18:52:01 +00:00
Martin Diehl 8fac635c15 made hybridIA stuff working again 2015-10-13 17:02:07 +00:00
Martin Diehl 6de8c821f8 improved consistency check on loadcase 2015-10-12 19:02:42 +00:00
Luv Sharma 9ec7c08bb6 gfortran was complaining abound mismatch 2015-10-07 13:07:21 +00:00
Martin Diehl c1a53bc8c0 introduced dummy target to check compilation 2015-10-06 18:00:17 +00:00
Martin Diehl 0f0b4eaa6b cleaned 2015-09-28 17:47:00 +00:00
Martin Diehl ad16162eb8 simplified 2015-09-24 17:50:11 +00:00
Martin Diehl c50a522dd8 some shapes not 100% correct, removed inverse laplace 2015-09-24 17:38:49 +00:00
Martin Diehl 2e70e6fded make clear that this isNaN is an intrinsic extension 2015-09-24 08:45:44 +00:00
Aritra Chakraborty b528088653 added all 13 slip families for bct(beta)-Sn
increased the allowable columns for lattice.f90 to allow for the large interaction matrix.
2015-09-22 18:42:23 +00:00
Martin Diehl 42ac7902f7 indents wrong in fortran, seeds_fromDistribution bugfix and using ASCII table facilities 2015-09-20 21:03:10 +00:00
Martin Diehl 41f09dd0f8 should make no difference, but test with ifort was failing. compierler buig? 2015-09-12 18:26:25 +00:00
Tias Maiti 03896fcd09 redefined required debug variables removed in last commit 2015-09-11 18:20:37 +00:00
Tias Maiti 840fd0e9fa report absolute value of accumulated_shear_slip in dot state 2015-09-11 17:44:26 +00:00
Martin Diehl 8716752bf8 fixed fourier convolution and div/curl calculation for even/odd grid according to Johnston 2011 (MIT, FFTW) 2015-09-11 13:05:46 +00:00
Martin Diehl 8f32d03a9e curl calculation overestimated RMS due to factor 2 instead of one for DC component and Nyquist component 2015-09-11 09:19:13 +00:00
Martin Diehl e88cedc6ae clearer naming, debug options for spectral do not work for MPI 2015-09-11 08:52:03 +00:00
Martin Diehl 751d1d7582 function statement is more clear 2015-09-10 13:01:33 +00:00
Martin Diehl 580bb76b4c pInt missing 2015-09-10 10:26:09 +00:00
Martin Diehl 427b5f4bc1 wrong results for restart 2015-09-10 10:06:14 +00:00
Martin Diehl bbe37c842e name conflict 2015-09-09 21:52:00 +00:00
Martin Diehl 34980a1d44 more logical structure when reading in 2015-09-08 19:44:32 +00:00
Martin Diehl 7d996cfa11 fixed new IO_stringPos for Abaqus 2015-09-05 16:26:55 +00:00
Martin Diehl caf0ac7e8d more verbose 2015-09-01 16:53:48 +00:00
Martin Diehl 7a3a67601f realloc lhs was not always working, fixed now and added test for new string pos function 2015-08-31 16:30:04 +00:00
Pratheek Shanthraj 197ae53553 lots of memory savings (~ stateside + 250 pReals per integration point) when using analytic jacobian 2015-08-28 10:55:38 +00:00
Martin Diehl 6848d83d13 behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name.
functions in the core module are not used, including kdTree
2015-08-28 07:38:48 +00:00
Martin Diehl 1171dc4344 MPI has a 2GB limit for writing at once, now chunking 2015-08-21 17:51:05 +00:00
Martin Diehl 8474da26d9 empty line 2015-08-18 16:31:47 +00:00
Martin Diehl bae5cfcf89 named if 2015-08-18 16:07:01 +00:00
Tias Maiti 53c574033a size of MAXCHUNKS same as maximum allocated size for “heat_time” & “heat_rate” 2015-08-18 03:20:16 +00:00
Martin Diehl f64a5e7552 output safe for 64bit integer 2015-08-14 05:03:54 +00:00
Martin Diehl 31ea4dadf1 IO_stringPos replacement not needing MAXNCHUNKS but making use of new Fortran features for dynamic allocation 2015-08-13 14:54:34 +00:00
Martin Diehl 563d9e64dd MAXNCHUNKS capitalized since it is a parameter 2015-08-06 09:24:56 +00:00
Pratheek Shanthraj 3d94ae4074 finite strain correction for respective expansions 2015-08-05 14:12:54 +00:00
Tias Maiti 09c7affbbc more efficient rate calculation in updateState when subdt is zero 2015-08-05 11:09:38 +00:00
Tias Maiti 7928ce3ac3 fixed erroneous rate calculation in updateState when time increment (subdt) was zero during initialization. 2015-08-04 21:27:30 +00:00
Tias Maiti 65d114e4f0 improved debug statements 2015-08-04 21:26:22 +00:00
Pratheek Shanthraj 905b9f8f42 petsc update (now v3.6.1) 2015-08-04 15:04:53 +00:00
Pratheek Shanthraj 5b57ea7d16 fluctuations proportional to thermal equilibrium concentration 2015-08-03 11:08:59 +00:00
Pratheek Shanthraj a01d7d8d66 corrections to crystallite_push33ToRef function 2015-08-03 11:07:19 +00:00
Pratheek Shanthraj 1ed276ce2d quadratic softening for finite dissipation energy 2015-07-30 10:29:25 +00:00
Pratheek Shanthraj c97b83c9a1 added user defined external heat source 2015-07-27 11:09:37 +00:00
Pratheek Shanthraj 4bb6664de9 tightened default thermal tolerance 2015-07-27 11:08:43 +00:00
Pratheek Shanthraj c82c4b74f1 updated config files to reflect recent changes in input parameters 2015-07-24 17:14:03 +00:00
Pratheek Shanthraj 5e09954575 remove deprecated spectral load case definition of temperature 2015-07-24 14:57:29 +00:00
Pratheek Shanthraj 7554647c8e more control over initialisation of field values. specify initial field value in the homogenisation part of the material config file using the appropriate tags 2015-07-24 14:53:50 +00:00
Pratheek Shanthraj 87d42bf447 initialise Fi correctly for initial field values away from equilibrium 2015-07-24 14:47:18 +00:00
Pratheek Shanthraj 0a5ccb3d91 possibility to set damage potential order 2015-07-24 14:43:05 +00:00
Tias Maiti 7c101cdc31 fixed issues causing compilation errors 2015-07-23 19:22:03 +00:00
Tias Maiti 7f8a8d5b0f new function to report initial thermal strain based on current temperature deviation from reference.
to be used in constitutive_initialFi, which collects from all kinematics, derives and sums Li terms to report initial Fi.
2015-07-23 15:59:25 +00:00
Philip Eisenlohr 25a03e128a added remark to include looping over sources of same type within one phase. (for instance to accommodate different activity periods.) 2015-07-16 17:20:18 +00:00
Su Leen Wong d2c748e716 TRIP input parameters for fcc to hex 2015-07-10 09:02:48 +00:00
Luv Sharma 1e9682eef9 corrected upper bound for spectral damage solver as well. 2015-07-09 15:51:34 +00:00
Shaokang Zhang d2bfd26937 1) added the error message of all MPI plans 2) deleted the suffix (MPI) for all MPI plans 2015-07-09 13:38:21 +00:00
Su Leen Wong 7a4415b28a Sanity checks for SFE 2015-07-08 15:37:16 +00:00
Su Leen Wong df5cdb6b01 Elasticity matrices for transformed phase 2015-07-08 11:58:52 +00:00
Su Leen Wong 66564c1f75 Read in elastic constants for transformed phase 2015-07-07 13:23:56 +00:00
Luv Sharma 4435244681 corrected missing initialization. newState = oldState + deltaState. But for irreversibility driving force must be compared with last converged subState. 2015-07-06 09:34:43 +00:00
Luv Sharma 3d5a85e4af corrected implementation of deltaState for isoBrittle damage 2015-07-06 08:47:40 +00:00
Martin Diehl 72b7df34a3 save some space 2015-07-02 09:45:25 +00:00
Aritra Chakraborty 74af77a81f formatted the interaction matrix for beta-Sn 2015-07-02 05:45:56 +00:00
Aritra Chakraborty 0e6ac411bb further modified to add longer lines for (bct) beta-Sn in lattice 2015-07-02 05:44:45 +00:00
Martin Diehl 6272129f01 formatted code and modified Makefile to have long lines in lattice 2015-07-02 05:09:07 +00:00
Aritra Chakraborty e61050dcde added lattice for (bct) beta- Sn. 2015-06-27 14:55:30 +00:00
Shaokang Zhang a40d247f7f lower bound for damage is residual stiffness 2015-06-22 14:27:32 +00:00
Su Leen Wong 12bc4f01ef fcc to hex orientation relationship (Shoji-Nishiyama) 2015-06-22 13:47:15 +00:00
Su Leen Wong ad17c40d92 Renamed variables to denote fcc to bcc transformation 2015-06-22 08:33:48 +00:00
Su Leen Wong 5114e0b43e Cleaning up 2015-06-22 08:10:20 +00:00
Philip Eisenlohr 3e0e187b14 clarified meaning of symEulers function in help 2015-06-19 07:29:01 +00:00
Franz Roters b1c036daab added some comments on possible code changes 2015-06-17 06:53:00 +00:00
Franz Roters 0e48fb8fc7 removed double entry in use list 2015-06-16 07:57:42 +00:00
Su Leen Wong 081bc33159 Read in c/a ratio for transformed hcp phase 2015-06-12 08:59:39 +00:00
Su Leen Wong a370372c1b transformed lattice structure should not be mandatory 2015-06-11 17:03:34 +00:00
Pratheek Shanthraj 40701cedc4 damage driving force is a history variable and damage field is strictly bounded between 0 and 1 2015-06-11 09:03:51 +00:00
Pratheek Shanthraj e8ee5d6723 moved some phase field parameters to lattice 2015-06-11 09:01:37 +00:00
Su Leen Wong b3241411f5 Added keyword to specify transformed lattice structure 2015-06-11 08:23:27 +00:00
Su Leen Wong c2781b5423 Removed unnecessary variables 2015-06-10 12:12:03 +00:00
Pratheek Shanthraj 859fb43fe5 split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options 2015-06-09 13:28:50 +00:00
Su Leen Wong b410a08995 Renamed some variables for consistency 2015-06-08 14:26:50 +00:00
Su Leen Wong 0d36a2d3e2 Removed unused variables and renamed some variables for TRIP model 2015-06-08 09:36:52 +00:00
Pratheek Shanthraj 74982294a0 added spectral thermal and damage solvers 2015-06-03 17:30:31 +00:00
Pratheek Shanthraj 3415e249ef reverting files that got mistakenly added into previous change set 2015-06-01 16:11:37 +00:00
Pratheek Shanthraj a14070bad4 changed up handling of delta states, and some bug fixes 2015-06-01 16:02:27 +00:00
Pratheek Shanthraj 8d2aa5e7be fixed typo in dislotwin label 2015-05-29 11:03:25 +00:00
Pratheek Shanthraj 8f4663985a major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs 2015-05-28 17:02:23 +00:00
Martin Diehl b58cbe424d follow up to last commit 2015-05-19 17:41:49 +00:00
Martin Diehl b07ec8232e found circular inclusion 2015-05-19 17:34:04 +00:00
Martin Diehl 71b0e283c1 clearer substructure of routines 2015-05-10 20:55:36 +00:00
Martin Diehl 4f9531e736 sliptwin is not used and will not be used as discussed in the documentation in private 2015-05-05 20:17:50 +00:00
Martin Diehl 485636b93b functions not longer needed for core module.
simplified some statements, removed double computation for direct matrix inversion
2015-05-05 06:37:59 +00:00
Martin Diehl 906c3f63a1 updated hybridIA sampling to work with new format
cropLinearODF is not working for the new format, but filterTable should be able to do the task
2015-04-27 05:30:29 +00:00
Martin Diehl 97aba96ff8 changed format of linearODF file to be more descriptive, essentially TSL OIM format with header and space separated keys 2015-04-26 11:07:39 +00:00
Martin Diehl 1f809f9d73 variable name change forgotton 2015-04-22 04:34:35 +00:00
Martin Diehl 78d913c07d clearified keywords, added raw data of rolling texture, renamed to more meaningful names 2015-04-21 19:34:50 +00:00
Pratheek Shanthraj d7df5f1934 more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature. 2015-04-21 16:11:30 +00:00
Pratheek Shanthraj 1330dae432 cleaning up, removing unused variables and variable initialisation 2015-04-21 15:16:13 +00:00
Pratheek Shanthraj 0681f67570 corrected buggy previous commit 2015-04-21 14:33:38 +00:00
Pratheek Shanthraj bec17127d0 write output information only for active homogenisations. removed unused variables 2015-04-21 14:10:34 +00:00
Pratheek Shanthraj d049eadce7 write output information only for active crystallites. removed unused variables. fixed some real == 0 comparisons 2015-04-21 12:23:00 +00:00
Pratheek Shanthraj effec313ee fixed bug in cleavage system set up for ortho lattice structures 2015-04-21 12:19:41 +00:00
Martin Diehl ea1b2191f1 deleted comment, test passed probably because fftw is system-wide available at MPIE
also, ../lib is not needed since lib is an include dir via -I
2015-04-19 09:07:16 +00:00
Philip Eisenlohr a46b9287eb restored conditional (#ifdef Spectral) inclusion of “fftw3”. Not clear how this could pass the Marc compile test..?!? (failed in my “real-world” use) 2015-04-18 15:23:37 +00:00
Philip Eisenlohr 22f9a90bcb added informative ext_msg to error 135 2015-04-18 15:22:15 +00:00
Martin Diehl eb3173c6cc ugly hack for Abaqus in prec, it seems that it links against old (10.1) ifort where IEEE_arithmetic is not avaialbe.
Took the chance and cleaned the env file and used #ifdef# statement in the wrappers because in fixed format Fortran some lines of code might be interpretated as comments.
This should not happen, but it happens for Abaqus exp
2015-04-15 18:10:56 +00:00
Martin Diehl caf6611485 unused variable 2015-04-15 18:05:27 +00:00
Martin Diehl f5762209dc changed 'range' keyword to 'limit', now using proper ASCII table (with column 'intensity') for linearODF
hybridIO_linODFsampling simplified
2015-04-15 14:37:46 +00:00
Martin Diehl e29628b459 changed keyword for texture component/ODF simulation
kewyword "origin" with possible values "center/voxelcenter" and "boundary/boundarycenter" is always needed
2015-04-14 12:45:20 +00:00
Martin Diehl 7e0b79a33c also using isNaN function 2015-04-14 11:43:31 +00:00
Martin Diehl 25665b1d80 IEEE_value (like any other function from intrinsic modules) is not allowed for parameters 2015-04-14 05:21:08 +00:00
Martin Diehl 0c7e7c66f8 integer type was use for void fraction, should be real 2015-04-13 10:36:45 +00:00
Martin Diehl 5c1804e77d improved NaN checks 2015-04-13 10:02:52 +00:00
Martin Diehl 50998bd6a4 simplified makefile and remove unused variables in RGC 2015-04-12 08:11:22 +00:00
Martin Diehl 7156e3996b gfortran noticed that these functions are external 2015-04-11 20:28:07 +00:00
Martin Diehl 0ba8f27320 took some gfortran complaints serious (unused imports, implicit castings) 2015-04-11 11:47:33 +00:00
Martin Diehl 064266c0cd updated linearODF file format (including header and keywords) 2015-04-11 10:28:10 +00:00
Martin Diehl b4456e76b3 introduced two small bugs 2015-04-11 09:25:15 +00:00
Martin Diehl 470fc2dce3 comparison for zero should be abs:
a == 0.0_pReal   becomes abs(a) <= tiny(a)
a /= 0.0_pReal   becomes abs(a) > tiny(a)
remove unused variables
2015-04-11 08:25:23 +00:00
Martin Diehl 6170209198 use statement for marc was wrong (preprocessor) 2015-04-11 05:09:15 +00:00
Martin Diehl d7b36c2c30 gfortran complaint about implicit casting,floating point comparison, and unused imported variables.
additionally changed use of math_invert33 to math_inv33 if det is not needed
2015-04-10 19:09:26 +00:00
Martin Diehl 7d0a1ebea6 math_transpose33 not used 2015-04-08 18:33:08 +00:00
Martin Diehl b6481c2513 introduced simpler multiplication and forall loops
matmul is ok for openmp, check in the web and run the state integration test.

Example program testing for new state update for rkck dot state:

program test

real, dimension(6,10) :: dotState=reshape(&
                      [1,1,1,1,1,1,1,1,1,1,&
                       2,2,2,2,2,2,2,2,2,2,&
                       3,3,3,3,3,3,3,3,3,3,&
                       4,4,4,4,4,4,4,4,4,4,&
                       5,5,5,5,5,5,5,5,5,5,&
                       6,6,6,6,6,6,6,6,6,6],[6,10])

real, dimension(10) :: residuum
real, dimension(6) :: B=2.5

integer :: i

residuum = B(1)*dotState(1,:)+& 
           B(2)*dotState(2,:)+& 
           B(3)*dotState(3,:)+& 
           B(4)*dotState(4,:)+& 
           B(5)*dotState(5,:)+& 
           B(6)*dotState(6,:)
do i =1,10
  print*,residuum(i)
enddo

residuum =  matmul(transpose(dotState),B)
do i =1,10
  print*,residuum(i)
enddo

end program test
2015-04-01 16:45:53 +00:00
Martin Diehl 693efcaa58 comment was confused by search and replace 2015-04-01 16:04:33 +00:00