Commit Graph

2250 Commits

Author SHA1 Message Date
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
Martin Diehl 63491bf268 for some reasons, ifort was complaining about mixed integer types.
also deleted #ifdef PETSC since we do not support the solver without PETSC
2015-03-30 20:49:17 +00:00
Pratheek Shanthraj c7418db9bd FEM : hierarchical ordering of h5 output data and more meaningful visualisation for multiple grains/crystallites/phases 2015-03-30 09:45:10 +00:00
Martin Diehl e0f1132a17 gfortran complaints: equal comparison of reals and external (MPI) functions 2015-03-29 12:54:13 +00:00
Martin Diehl 9e824d029e renamed to new names: spectralsolver and spectralfilter 2015-03-28 21:14:11 +00:00
Martin Diehl 1b2e4c0239 removed old solver 2015-03-28 17:41:40 +00:00
Martin Diehl 8a30021202 cleaned up core modul 2015-03-27 14:10:18 +00:00
Pratheek Shanthraj c8554a5641 corrected alignment of restart output byte offset. restart test working now 2015-03-27 07:17:24 +00:00
Pratheek Shanthraj 1ee81e74ea more patches to get the tests running again:
- fixed increment counting in postResults to current output scheme

- corrected node coordinates calculation

- some work on restart (restart test still not running)
2015-03-26 21:19:28 +00:00
Pratheek Shanthraj 86283359f1 PetscFinalize should be the last call before exiting 2015-03-26 13:43:18 +00:00
Pratheek Shanthraj a9cf8e8051 some patches to compile damask core module without fftw-mpi 2015-03-26 13:41:19 +00:00
Pratheek Shanthraj 2417877bdb dropped non petsc basic spectral solver 2015-03-25 16:08:41 +00:00
Pratheek Shanthraj d44fce4a76 Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me 2015-03-25 16:06:19 +00:00
Pratheek Shanthraj 37a7364a3e multi physics output now working for all solvers 2015-03-25 16:02:30 +00:00
Pratheek Shanthraj 2d4e7fb8e7 calculate dFi/dS (9x9 matrix inverse) for the analytic jacobian only when needed. should improve performance 2015-03-20 07:41:11 +00:00
Pratheek Shanthraj 2e44a846af need to keep track of temperature again since we can no longer get this from F_T 2015-03-18 18:03:18 +00:00
Pratheek Shanthraj c5a8543f17 updated petsc 2015-03-18 17:18:43 +00:00
Martin Diehl 4a0c1c1717 comments 2015-03-15 15:30:14 +00:00
Martin Diehl 12cb5e56fe fixed interpretation 'spectralOut' header. Seem to have been buggy all the time, but manifested itself only after keyword change. Also added ':' after keywords to be more robust 2015-03-15 15:29:38 +00:00
Martin Diehl ac3d82493e changed output keywords form resolution to grid and from dimension to size 2015-03-13 14:45:33 +00:00
Martin Diehl e4d70e971d cp is commercially pure ;) 2015-03-13 13:49:09 +00:00
Martin Diehl d8debc3396 updated reference geom files (used old names), updated numerics.config to reflect use of analytic jaco 2015-03-13 06:47:26 +00:00
Martin Diehl aa125bac5b aliases for grid and size not supported any more 2015-03-12 22:30:24 +00:00
Martin Diehl 00cba25a44 improved update of ip coordinates for spectral solver.
do not need to create a new fftw plan all the time, using data already defined for the convolution
2015-03-12 22:28:33 +00:00
Martin Diehl 263f997cf2 seems to be an old output 2015-03-12 08:51:56 +00:00
Martin Diehl b338f1aeb8 added some outputs 2015-03-09 23:10:36 +00:00
Su Leen Wong 6b0a2ed5a4 Changed initialization of projection matrix for transformation based on the number of slip systems 2015-03-09 12:05:38 +00:00
Martin Diehl a5ee02ff01 added some aliases for output to finally converge from 'shear_rate' to 'shearrate' and from 'accumulatedshear' to 'accumulated_shear' 2015-03-09 11:29:44 +00:00
Franz Roters 2ab90c28a5 replaced example ODF by reasonable data in fully correct format 2015-03-06 15:56:58 +00:00
Pratheek Shanthraj 35f9e91e73 some corrections to the damage models 2015-03-06 13:12:32 +00:00
Pratheek Shanthraj bbb5ff6ae9 changes related to intermediate configuration kinematics:
- switched Fi and Li from state variables to crystallite variables

- Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration

- T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi

- Updated analytic jacobian to take into account tangents wrt Fi

- Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 13:09:00 +00:00
Yannick Naunheim 4fd7e0165e more meaningful name 2015-03-05 11:07:55 +00:00
Su Leen Wong 36f14a2d1d Removed nr = ns condition 2015-02-26 14:16:29 +00:00
Su Leen Wong 2baac6fc10 Fix for 3881 because gfortran was complaining 2015-02-24 10:15:07 +00:00
Pratheek Shanthraj 3c2a95efd8 if matrix inversion error encountered during analytic tangent calculation, give a warning and return elastic tangent 2015-01-30 14:00:32 +00:00
Pratheek Shanthraj f1f4f939a3 analytic jacobian default for all solvers 2015-01-29 13:59:49 +00:00
Pratheek Shanthraj dd8458a775 updated analytic jacobian calculation to correctly take into account intermediate configuration Fi. improved convergence of Li loop in stress integration 2015-01-29 13:58:25 +00:00
Pratheek Shanthraj 0b59519a2a updated damage models:
- coupling to plasticity handled within damage module instead of plasticity module
- anisotropic models more stable
2015-01-29 13:56:09 +00:00
Martin Diehl 12c586afb9 removed unused variables 2015-01-23 13:44:16 +00:00
Martin Diehl 1f666dffc9 fixed bug: output size was not set 2015-01-23 13:08:25 +00:00
Martin Diehl e9c3137a3e disabled a faulty ntrans loop (matmul was complaining) 2015-01-23 11:23:35 +00:00
Martin Diehl 43095b1f1b added threshold_stress_slip as output (got lost once upon a time) 2015-01-22 09:34:11 +00:00
Martin Diehl 1c2dd37e28 deleted unused variables in dislo UCLA 2015-01-21 15:14:00 +00:00
David Mercier e57cb2b4f4 Wrong number of slip families in BCC material config files 2015-01-20 17:22:52 +00:00
Luv Sharma 19dc945065 Inroduced Fd in kinematic decomposition. 2015-01-19 13:39:47 +00:00
Pratheek Shanthraj 54dc4943de adaptive penalty to bound Cv in (0,1) 2015-01-16 18:31:14 +00:00
Pratheek Shanthraj f16244d240 corrected state integration 2015-01-16 17:51:43 +00:00
Pratheek Shanthraj eecb78665c removing some unused code and cleaning up 2015-01-16 17:32:47 +00:00
Pratheek Shanthraj 6411d36633 updated vacancy-damage model 2015-01-16 17:04:01 +00:00
Martin Diehl 93e50366bb fixed automatic determination of compiler 2015-01-16 11:19:30 +00:00
Luv Sharma 701f35c326 reverting the changes related to undamaged_Fi during rev. 3850/3851. 2015-01-16 00:09:22 +00:00
David Cereceda Senas 85465333ec 2015-01-15 22:42:56 +00:00
David Cereceda Senas b1ef2f9dd9 2015-01-15 22:40:17 +00:00
Martin Diehl 62ff5d1d45 gfortran/ifort is now selected depending on the compiler provided by PETSc 2015-01-15 11:35:53 +00:00
Martin Diehl ced00954fe added new model by david (LLNL/UCLA) 2015-01-15 10:56:15 +00:00
Luv Sharma 125276283b fix in analytically integrated expression. 2015-01-15 09:22:24 +00:00
Luv Sharma 19611660c2 passing undamaged Fe for calculation of stress.(Pratheek please check.) 2015-01-13 09:54:44 +00:00
Luv Sharma 9c8873f16c added functions to get degraded stiffness and undamaged intermediate deformation gradient. 2015-01-13 09:46:25 +00:00
Luv Sharma 80bb9a2e66 function to calculate mises stress. 2015-01-13 09:43:05 +00:00
Haiming Zhang c2b8627981 Currently, there is no updating for the dot state of vacancy generation, set sizeDotState = 0 2015-01-11 20:01:43 +00:00
Haiming Zhang a2552d6241 Polishing. 2015-01-11 19:10:16 +00:00
Martin Diehl e20271b89a fix for temperature in spectral solver 2015-01-08 13:26:00 +00:00
Luv Sharma 9635ae9f3d fix in integration. 2015-01-07 14:51:06 +00:00
Luv Sharma 8610531aae removed unused variable. 2015-01-04 20:13:16 +00:00
Luv Sharma 460c0d73fa Making use of extended plastic state with pointers, removed helper functions. 2015-01-04 19:26:33 +00:00
Luv Sharma e6d1b21b57 introduced extended plastic state type (including pointer/aliases) 2014-12-19 14:01:03 +00:00
Martin Diehl e4cb387571 polishing comments 2014-12-18 18:44:47 +00:00
Martin Diehl 72cb9072e1 added some debug statements 2014-12-18 18:41:02 +00:00
Martin Diehl 7fe20fe805 intent(out) specifier for getXX functions 2014-12-18 18:26:16 +00:00
Martin Diehl fdbaf00fa7 cleaned up naming and added elastic material with volume preserving stiffness 2014-12-18 18:20:11 +00:00
Martin Diehl 299518adfa cleaned tungsten parameter set, renamed to be consistent 2014-12-18 16:55:43 +00:00
Pratheek Shanthraj 3f14ebe43d helper routine to get heat generation rate needed for MARC/Abaqus 2014-12-17 13:37:13 +00:00
Franz Roters e6e72ab3a7 In mesh_marc_build_elements use mesh_NcpElems instead of mesh_Nelems (line 2906).
Fixes error when not all elements are CP elements.

Some formating and comment improvements.
2014-12-15 11:51:32 +00:00
Martin Diehl 7c8f2ede31 fixed missing initialization causing a NaN in gdot and crashes 2014-12-11 15:04:16 +00:00
Haiming Zhang 5bc93efd63 Fix bug about vacancy generation in subroutine constitutive_microstructure 2014-12-11 13:07:15 +00:00
Pratheek Shanthraj bdb6bf8090 recasting vacancy dot state as microstructure state 2014-12-10 13:40:30 +00:00
Pratheek Shanthraj 39a3a5b5a0 cleaned up FEM examples 2014-12-09 18:59:37 +00:00
Haiming Zhang 3064c14ff6 fig a bug of energy barrier of vacancy generation. 2014-12-09 18:56:19 +00:00
Pratheek Shanthraj 3f8678c7c4 plastic dot states evaluated at effective stress not damaged stress 2014-12-09 18:23:48 +00:00
Martin Diehl 55e2de6ffd made dotstate and state pointers (targets are not allowed in derived types) and introduced plastic state having aliases to parts of the state. 2014-12-09 12:12:53 +00:00
David Mercier d1e1191df4 corrected bug: renamed file list for commercial FEM solvers to start with small letter 2014-12-08 16:36:22 +00:00
Pratheek Shanthraj 2458aa9c3e BBar stabilisation now optional (default is off) 2014-12-08 16:27:23 +00:00
Philip Eisenlohr c34252c71c changed constitutive_* prefix for plasticity laws to plastic_* to be consistent with novel thermal_, damage_, vacancy_, etc. laws. 2014-12-08 15:55:30 +00:00
Haiming Zhang 0bd4ab48f8 numerics with respect to vacancy diffusion. 2014-12-08 15:37:44 +00:00
Martin Diehl db0e02b94d changed to lower canse convention 2014-12-08 15:37:02 +00:00
Pratheek Shanthraj 67c3eb3a2e bounding the analytically integrated damage variable 2014-12-08 14:05:38 +00:00
Pratheek Shanthraj 347dac74c6 added vacancy state integrators 2014-12-08 10:45:12 +00:00
Haiming Zhang 9301b56aa3 initialize restoration points of vacancy concentration at material points, postResults of VC at material points. 2014-12-08 09:27:32 +00:00
Haiming Zhang 62f254c9de age results and write restart data for vacancy concentration 2014-12-08 09:09:54 +00:00
Franz Roters 55c2b729a0 while node coordinates are always 3D displacements are passed into hypela2 as 2D array (dispt) for 2D simulations 2014-12-05 07:35:28 +00:00
Luv Sharma 334ade5260 Every line in a recipe must begin with a tab character.
http://stackoverflow.com/questions/14576237/errors-in-make-file-missing-separator-stop
2014-12-03 13:37:04 +00:00
Philip Eisenlohr 9d12e026e8 renumbered some errors, fixed non-existing errors 2014-12-03 00:42:35 +00:00
Philip Eisenlohr 01a8c500ca detabbing… 2014-12-03 00:41:12 +00:00
Martin Diehl fa6f784102 deleted unused stuff 2014-12-02 21:49:33 +00:00
Pratheek Shanthraj 82942d1618 bug fix for previous commit 2014-12-02 17:23:47 +00:00
Pratheek Shanthraj b796e46f98 residualStiffness is not just an FEM parameter 2014-12-02 17:17:35 +00:00
Pratheek Shanthraj 2fc264d90a forgot file in previous commit 2014-11-28 20:07:22 +00:00
Pratheek Shanthraj 37dd574790 recast power law and viscous dot state variables as derived quantities using analytical solutions of the evolution equation 2014-11-28 20:06:24 +00:00
Philip Eisenlohr 8d64e513da improved commenting and variable names.
values for unused slip/twin families are ignored in material.config (do not need to be numeric anymore…)
2014-11-26 21:23:02 +00:00
Tias Maiti f8b29bdaf0 fixed bug in type definition for ‘drivingForce’ 2014-11-26 18:40:31 +00:00
Pratheek Shanthraj 54de033c45 using analytic solution for local/nonlocal coupled damage (as opposed to having crystallite iterate to find a consistent solution) 2014-11-26 17:29:45 +00:00
Pratheek Shanthraj 61bb9a373e all damage models consistently semi-explicit 2014-11-26 11:32:50 +00:00
Pratheek Shanthraj 6e90260817 updated local damage post results for newer models 2014-11-26 11:32:19 +00:00
Pratheek Shanthraj a91ebafa0e semi-explicit damage evolution 2014-11-26 10:36:14 +00:00
Luv Sharma 646f1c48dc typo, IO_warning and not IO_error 2014-11-26 10:16:30 +00:00
Pratheek Shanthraj 9714a172d1 updated to get local damage from microstructure state 2014-11-26 00:42:20 +00:00
Pratheek Shanthraj c89ef9de61 relaxed some sanity checks 2014-11-26 00:29:54 +00:00
Pratheek Shanthraj 26c01f64e5 updated file list 2014-11-26 00:06:16 +00:00
Pratheek Shanthraj 32ee3a4a66 corrected constitutive_getVacancyPotentialDrivingForce function 2014-11-26 00:04:06 +00:00
Pratheek Shanthraj 36dc59b09f reverted back to rate independent damage evolution. updated example config files to recent changes 2014-11-25 23:56:52 +00:00
Pratheek Shanthraj a4584fecd2 potential driving force for conservative cahn hilliard vacancy diffusion 2014-11-25 20:13:33 +00:00
Pratheek Shanthraj 20437b2ae0 added brittle(elastic energy release rate)/ductile(vacancy condensation) phase field damage model coupled to vacancy concentration
simplified vacancy_generation
2014-11-25 17:23:37 +00:00
Pratheek Shanthraj 65b22f9a40 calculate analytic jacobian from xxx0 and not xxxsub0 2014-11-25 17:20:52 +00:00
Su Leen Wong e6edc3d8cd Some sanity checks 2014-11-25 13:20:12 +00:00
Pratheek Shanthraj 296146490f changed some default parameters for FEM 2014-11-21 14:20:54 +00:00
Luv Sharma 3d4b2d2f52 tolerance for damage can be any positive value. Set to 0.001 if found negative. 2014-11-21 13:12:54 +00:00
Luv Sharma ee1d409bad critical shear can be set to zero for family not in use. 2014-11-21 09:03:54 +00:00
Su Leen Wong 9e6ab9866f Read in TRIP inputs, changed 'totalvolfrac' to 'totalvolfrac_twin' 2014-11-21 08:54:20 +00:00
Luv Sharma 36c2cc4347 added some sanity checks in damage modules. 2014-11-20 13:48:07 +00:00
Martin Diehl f538469892 reverted name change: In current PETSc release, its DM_ instead of DMDA_! 2014-11-19 05:35:10 +00:00
Philip Eisenlohr 73a13ddffc added more comments 2014-11-18 15:36:17 +00:00
Su Leen Wong b3a96d40c5 Started adding TRIP model to phenopowerlaw 2014-11-18 10:46:02 +00:00
Pratheek Shanthraj 23d68114e3 cleaned up stress integrations and some documentation 2014-11-13 12:53:20 +00:00
Haiming Zhang b407e6ac80 complement of pore nucleation. 2014-11-12 17:55:02 +00:00
Pratheek Shanthraj fc9f4835c3 switched order of Lp-Li nesting in stress integration loop for better convergence. temperature is integrated during stress integration so it need not be a dot state 2014-11-12 16:40:50 +00:00
Pratheek Shanthraj e89199d119 correction to thermal velocity gradient 2014-11-11 23:14:45 +00:00
Pratheek Shanthraj 7567582781 updated constitutive in line with changes to damage 2014-11-11 22:46:41 +00:00
Pratheek Shanthraj 5ff3c882b5 thermodynamically consistent treatment of nonlocality over all damage models. switched from power law to exponential with viscous drag type rate formulation for damage evolution 2014-11-11 22:44:11 +00:00
Haiming Zhang f0f04a25bf Pore nucleation, and the decrease of vacancy concentration due the formation of pores. 2014-11-11 16:03:06 +00:00
Pratheek Shanthraj e0e0f1503b svn properties for new files 2014-11-07 20:55:16 +00:00
Pratheek Shanthraj 9544fe55a6 FEM solver for vacancy diffusion 2014-11-07 20:40:54 +00:00
Pratheek Shanthraj eb1028eb6f dislokmc coupled to ductile damage 2014-11-07 12:28:18 +00:00
Pratheek Shanthraj 646f632977 no need to pass fileunit during init of none types 2014-11-07 12:15:28 +00:00
Pratheek Shanthraj c4e54793fc missing factor of 2 in damage driving force 2014-11-07 12:10:10 +00:00
Pratheek Shanthraj dc0745800f reworked anisotropic brittle damage 2014-11-07 12:08:56 +00:00
Martin Diehl 69b150319f fully introduced non-schmid behavior 2014-11-07 11:23:34 +00:00
Martin Diehl eee863956f added flush(6) for debugging 2014-11-06 21:05:35 +00:00
Pratheek Shanthraj 4fcd1821cf reworked anisotropic ductile model 2014-11-06 17:52:43 +00:00
Pratheek Shanthraj 9c474796ac cleaning up and meaningful renaming 2014-11-06 17:41:34 +00:00
Haiming Zhang a72a9ed8ec Add the generation rate of vacancy induced by dislocation movement. 2014-11-06 16:53:34 +00:00
Martin Diehl 29951a9e1a named 99 representation according to other constitutive laws 2014-11-06 11:49:37 +00:00
Martin Diehl b443ac6321 nicer preprocessor statements 2014-11-06 11:47:27 +00:00
Martin Diehl 95ec8b7b7b removed imported module variables that are not longer needed, commented loops and simplified calculation for rates, i.e. do not store per system anymore but only once 2014-11-06 08:27:48 +00:00
Martin Diehl 1242aab2bf unified calculation for slip rate for output to avoid doubling of code
introduced calculation of negative and positive stress (three times in the code) for non-schmid behavior. 
Still, non-schmid behavior is not active at all
2014-11-05 20:53:25 +00:00
Martin Diehl 6fc3908c71 prepared for merging postResults rate calculation in one loop 2014-11-05 19:39:23 +00:00
Martin Diehl f3b7b5bb96 corrected substraction by volume fraction of twin for LpAndItsTangent and removed the substraction in postResults (like for dislowtin)
cleaned up and unified notation in calculation of slip rates
2014-11-05 19:11:09 +00:00
Martin Diehl accb571c53 corrected substraction by sum of twinned and transformed value fraction of Lp and derivative.
removed unused variables
2014-11-05 17:52:49 +00:00
Pratheek Shanthraj a61c3059ef missing file in previous commit 2014-11-05 17:47:46 +00:00
Pratheek Shanthraj f2211260e3 introduced possibly nonlocal anisotropic ductile damage model 2014-11-05 17:41:08 +00:00
Pratheek Shanthraj 94ba7e6246 added accumulated shear as a dot state for cleaner integration with ductile damage models 2014-11-05 17:07:31 +00:00
Pratheek Shanthraj ade918246c reduced number of state variables. some cleaning up 2014-11-05 17:05:40 +00:00
Martin Diehl a21f4d80db removed some leftovers of the shearbanding part 2014-11-05 09:53:08 +00:00
Martin Diehl c2799585d9 introduced non schmid tensor calculation in dot state, also not active yet 2014-11-05 09:35:49 +00:00
Martin Diehl c66fc95dda made a bunch of arrays simple scalars to improve readability, commented loops 2014-11-05 09:24:08 +00:00
Martin Diehl 816169fbec continued work on integration of non schmid behavior. not active yet 2014-11-05 09:03:04 +00:00
Martin Diehl 18b49fbac8 line in OMP statement too long 2014-11-04 21:32:16 +00:00
Luv Sharma e36880e3dc fixed typo 'orthorhombic' 2014-11-03 15:41:05 +00:00
Pratheek Shanthraj 76442942b1 damage def grad updated in stress integration loop and not as a dot state 2014-11-03 10:43:36 +00:00
Pratheek Shanthraj 0ad917dceb nicer FEM parameters 2014-11-03 10:36:07 +00:00
Pratheek Shanthraj 261bffa745 switched to smoother traction separation law for better convergence. anisotropic brittle now working.
lowered default max iterations for FEM
2014-11-01 17:51:35 +00:00
Pratheek Shanthraj 9c94eb4be6 OMP bug fix from previous commit 2014-11-01 11:43:43 +00:00
Pratheek Shanthraj f9cec3079e overhaul of kinematics from FeFp to FeFiFp decomposition 2014-10-31 19:03:08 +00:00
Pratheek Shanthraj 26a609dbee nicer permutation of cleavage planes 2014-10-31 18:37:37 +00:00
Martin Diehl 629ff8bdef reverted back to r 3644 to get rid of buggy behavior 2014-10-31 10:39:02 +00:00
Haiming Zhang ad4bd86fd2 Polishing
Modify
       texture: do i = 1_pInt,myNconstituents
         grain            = sum(NgrainsOfConstituent(1_pInt:i-1_pInt))
as:
       grain = 0_pInt
       texture: do i = 1_pInt,myNconstituents
         grain            = grain + NgrainsOfConstituent(i-1_pInt)

First try of developing DAMASK, :)
2014-10-30 09:51:06 +00:00
Pratheek Shanthraj b0469854c8 brittle anisotropic damage on cleavage planes 2014-10-28 20:57:12 +00:00
Pratheek Shanthraj a6f88c0e37 introduced crystallographic cleavage systems for each lattice type
FCC: (100) and (111) planes
BCC: (100) and (110) planes
HCP: (0001) planes
iso/orthotropic: (100) planes
2014-10-28 18:05:51 +00:00
Pratheek Shanthraj ceb294b77a labels should always be in lower case 2014-10-28 18:02:18 +00:00
Pratheek Shanthraj 16fc93d7db deleting and renaming files consistently with svn 2014-10-28 16:41:29 +00:00
Luv Sharma df58327dff renamed damage modules/functions etc.
brittle -> isoBrittle ; ductile -> isoDuctile; anisotropic-> anisoBrittle
2014-10-28 10:49:12 +00:00
Pratheek Shanthraj 14d71eb35b cleaner interface for ductile damage models 2014-10-28 02:42:25 +00:00
Pratheek Shanthraj ebd285f565 reworking damage routines 2014-10-28 01:18:10 +00:00
Luv Sharma 88159d1458 added function to return degraded stiffness where Fe..Fd..Fp decomposition not implemented 2014-10-27 15:33:35 +00:00
Su Leen Wong 53ae245e4f Added martensite lamellar term to mean free path of slip 2014-10-27 15:14:45 +00:00
Luv Sharma b65ccb0181 changed isotropic ductile damage to be inline with Fe..Fd..Fp decomposition 2014-10-27 14:15:25 +00:00
Martin Diehl bec0af7b06 started to introduce non-schmid behavior and disabled shear banding in dislokmc 2014-10-27 09:19:36 +00:00
Noriki Fujita aa608573b9 fixed bug about internal energy calculation 2014-10-24 07:45:00 +00:00
Su Leen Wong dfd5be70b5 Tolerance for TRIP volume fraction 2014-10-23 13:51:18 +00:00
Noriki Fujita 0fb61d8f4a Fixed SDV output, added internal energy output for ABAQUS Explicit 2014-10-23 09:28:21 +00:00
Noriki Fujita 32c4eeb0ff Added thermo-mechanical element type for ABAQUS 2014-10-23 08:17:37 +00:00
Pratheek Shanthraj cb81624b74 history variable is initial state not current state 2014-10-22 10:29:40 +00:00
Pratheek Shanthraj b54782046d increased size of petsc_options string 2014-10-21 16:25:53 +00:00
Martin Diehl 78eb90a754 added missing ticks 2014-10-21 15:35:14 +00:00
Martin Diehl 5c4e99f965 created new file containing all files to include for commercial FEM to simplify adding/renaming/removing modules 2014-10-21 15:31:28 +00:00
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
Martin Diehl f9a1e71207 added source 2014-10-20 08:17:09 +00:00
Pratheek Shanthraj fec71e2986 missing variable definitions 2014-10-16 08:15:57 +00:00
Pratheek Shanthraj 898aa0e9bf use slip system families for anisotropic ductile damage 2014-10-16 02:11:14 +00:00
Pratheek Shanthraj 4a445a9f54 fracture strain energy and plastic strain indexed by instance and not phase 2014-10-16 01:50:47 +00:00
Luv Sharma 92d5c89d4b started implementing anisotropic damage 2014-10-15 15:02:29 +00:00
Pratheek Shanthraj 32d558d1d6 set svn Id property for recently added files 2014-10-15 13:48:19 +00:00
Martin Diehl f9f258d639 fixed bug in initializtion of phase instance, added svn properties for FEM code 2014-10-15 12:22:09 +00:00
Pratheek Shanthraj 1f795d674c phase_thermalInstance was not being counted 2014-10-14 07:54:25 +00:00
Luv Sharma 674185a8de big fixes 2014-10-14 03:30:59 +00:00
Luv Sharma 4e3f8245a7 some improvements for void nucleation 2014-10-14 01:10:42 +00:00
Luv Sharma fc6c61471f pvec still used in FEZoo 2014-10-14 00:33:38 +00:00
Pratheek Shanthraj 3bd152db2d temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use. 2014-10-13 18:08:54 +00:00
Pratheek Shanthraj a057c540e0 temperature initialised to value from load case file 2014-10-13 17:54:27 +00:00
Pratheek Shanthraj 45dffc704c removed unused function arguments 2014-10-13 17:52:33 +00:00
Martin Diehl 75d8cacfbc removed old state pointer (p_vec) 2014-10-13 12:31:04 +00:00
Martin Diehl eeac38dd1a polishing 2014-10-13 08:44:49 +00:00
Martin Diehl 8f1a2cd2ae core module does not compile because it does not include PETSc 2014-10-11 16:54:59 +00:00
Pratheek Shanthraj bc9cf40371 vacancy diffusion proportional to accumulated slip (pipe diffusion along dislocations) 2014-10-11 10:39:44 +00:00
Pratheek Shanthraj c0584b47e1 modified damage diffusion tensor for brittle damage so this is now doing griffith's fracture criterion 2014-10-11 10:09:36 +00:00
Pratheek Shanthraj 616a44f898 - added get slip rate helper function
- added vacancy generation rate
- cleaned up some bugs
2014-10-11 09:45:30 +00:00
Pratheek Shanthraj dd16851ab7 implemented possibly diffusive vacancy physics. to be coupled with micro void nucleation and ductile damage 2014-10-10 20:55:09 +00:00
Pratheek Shanthraj b8a45f660e reinstating changes lost in previous commit 2014-10-10 16:45:14 +00:00
Luv Sharma 9b1e55f0e4 added helper functions to communicate accumulated slip between damage and plasticity modules 2014-10-10 16:34:51 +00:00
Pratheek Shanthraj e943a3a8b0 fixed typo 2014-10-10 16:23:19 +00:00
Pratheek Shanthraj 96f036e34e updated MPI reporting in line with recent changes 2014-10-10 16:21:10 +00:00
Pratheek Shanthraj 76f5495038 fixed bug for when petsc is not defined 2014-10-10 16:02:42 +00:00
Pratheek Shanthraj eacb048c83 fixed compilation error 2014-10-10 15:58:18 +00:00
Martin Diehl d095c2484d simplified multi processor (MPI) reporting 2014-10-10 13:08:34 +00:00
Luv Sharma 1657e0f7ba started implementing a void fraction based ductile damage model. 2014-10-10 12:42:12 +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 d10f8d3af6 move heatGen to adiabatic 2014-10-10 08:58:25 +00:00
Pratheek Shanthraj 885a71c0b7 cleaned up compilation for FEM 2014-10-10 08:40:59 +00:00
Pratheek Shanthraj c8929b47d3 only output from the root processor for parallel runs 2014-10-09 20:23:06 +00:00
Pratheek Shanthraj 45bedec229 damage diffusion tensor should be in reference configuration 2014-10-09 15:56:15 +00:00
Luv Sharma 2eadb0a178 changed thermal heatGen to thermal adiabatic 2014-10-09 14:08:32 +00:00
Luv Sharma fc57c6e572 added thermal adiabatic file 2014-10-09 14:06:45 +00:00
Pratheek Shanthraj c1a5ce047f disable openMP for FEM compilation 2014-10-06 20:45:04 +00:00
Pratheek Shanthraj 2712ce0254 polishing 2014-10-06 17:01:39 +00:00
Martin Diehl 020e697b74 small polishing 2014-10-01 20:36:00 +00:00
Pratheek Shanthraj 4d7a9ce030 bug fix 2014-10-01 16:35:56 +00:00
Luv Sharma ee60ea8d5b some improvements in ductile model 2014-10-01 15:18:33 +00:00
Martin Diehl baa08d8155 made analytic tangent default for spectral solver
removed possibility to compile without PETSc (a lot of effort for little/no use)
2014-10-01 12:29:12 +00:00
Su Leen Wong 0824335b10 TRIP contribution to Lp and dotState 2014-10-01 08:11:39 +00:00
Pratheek Shanthraj 3aa5ca68c5 marc/abaqus includes for ductile damage 2014-09-30 13:25:07 +00:00
Pratheek Shanthraj 3da023a1ba thermo-mechanical example 2014-09-29 12:40:57 +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 ce3be666e8 redoing previous changes that were lost in commit 3502
please please please make sure you merge your changes to an updated version of the code before committing them
2014-09-26 18:49:40 +00:00
Pratheek Shanthraj 590eb31ed0 some more renaming or thermal related routines 2014-09-26 18:07:48 +00:00
Pratheek Shanthraj c39aac9c9d makefile to go with renaming 2014-09-26 16:10:41 +00:00
Pratheek Shanthraj 40cd991496 thermal_none to thermal_isothermal 2014-09-26 16:09:25 +00:00
Pratheek Shanthraj dabfa6d2e6 renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal 2014-09-26 16:07:26 +00:00
Pratheek Shanthraj d4076e4db9 some fixes to get ductile damage working with recent changes 2014-09-26 15:23:10 +00:00
Pratheek Shanthraj 5c3f4e3e69 some fixes to get homogenisation none working with recent changes 2014-09-26 15:16:10 +00:00
Luv Sharma d3fd4bdf90 started adding a ductile damage model 2014-09-26 10:38:13 +00:00
Pratheek Shanthraj 82932c3708 added field quantity output 2014-09-26 10:34:36 +00:00
Pratheek Shanthraj 280e563fa0 removed homogenization_sizepostresults. This information only needs to be stores per homog type in homogState and not for every element and ip 2014-09-26 10:32:56 +00:00
Pratheek Shanthraj a42fd0608a use constitutive_xxx_Noutput instead of phase_Noutput as number of outputs for more consistency with constitutive multi physics output 2014-09-26 10:25:26 +00:00
Pratheek Shanthraj 34c33cca13 changed Noutput to constitutive_nonlocal_Noutput for consistency 2014-09-26 10:21:00 +00:00
Su Leen Wong 1c38d42a0b Projection matrices for transformation 2014-09-24 12:31:16 +00:00
Luv Sharma c38b80ad49 updated module include list for abaqus interface also. 2014-09-24 11:11:17 +00:00
Luv Sharma a0b724836e typo fixed 2014-09-24 10:45:46 +00:00
Su Leen Wong 4d351ccba1 Outputs for TRIP volume fraction 2014-09-24 07:26:27 +00:00
Su Leen Wong c2ee6ff299 Renaming transformation input variables to more logical names 2014-09-23 15:38:41 +00:00
Pratheek Shanthraj 556d485ca4 typo: diffusion not difusion 2014-09-23 12:30:52 +00:00
Pratheek Shanthraj 5ce49387f6 added history variable for brittle damage to prevent healing. renamed/removed some lattice damage tensors and symmetrize them. 2014-09-23 12:22:34 +00:00
Pratheek Shanthraj d2dee58703 removed obsolete thermal and damage modules 2014-09-23 10:42:57 +00:00
Pratheek Shanthraj c2c5c0c521 subsumed constitutive_damage/thermal into constitutive 2014-09-23 10:38:20 +00:00
Pratheek Shanthraj c029554c31 better naming damage_local --> damage_brittle 2014-09-23 08:43:16 +00:00
Pratheek Shanthraj 854879d277 undid some marc/abaqus related changes 2014-09-23 08:19:15 +00:00
Pratheek Shanthraj bb7f9be6ea remove obsolete constitutive field modules 2014-09-22 20:57:12 +00:00
Pratheek Shanthraj cb4117df95 adiabatic part of possibly conductive thermal model 2014-09-22 20:55:54 +00:00
Pratheek Shanthraj 3948d99f4f local part of possibly nonlocal brittle damage model 2014-09-22 20:38:19 +00:00
Pratheek Shanthraj f6c0f7062e moved field initialization from homogenization to material since it is needed in crystallite 2014-09-22 20:36:55 +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 2b35cc041e fixed wrong line continuation 2014-09-19 19:38:59 +00:00
Luv Sharma 649ec1fca8 Finished applying newstate to Homogenization. 2014-09-19 17:59:06 +00:00
Luv Sharma 4bbd6583a0 prepared for using new state structure in homogenization 2014-09-18 15:27:52 +00:00
Martin Diehl f64418f29a removed multithreaded FFTW, does not make a big time difference and will simplify linking when using PETSc 2014-09-18 15:23:11 +00:00
Chen Zhang 2d460c3b78 remove redundant debug output 2014-09-16 18:22:03 +00:00
Pratheek Shanthraj cc7064d025 changed default petsc_options for FEM and removed unused variables 2014-09-16 10:02:37 +00:00
Martin Diehl 27026c7f6a fixed missing initialization of output in RGC that got lost in Rev 3452 2014-09-16 09:50:37 +00:00
Su Leen Wong 74db2ef064 Small fix for missing constitutive dislotwin output 2014-09-16 08:02:19 +00:00
Martin Diehl 7cb7815176 polishing style 2014-09-13 10:04:44 +00:00
Luv Sharma 05d980f1e3 minor fixes 2014-09-12 15:58:03 +00:00
Su Leen Wong f8646daf85 Added projection matrix for strain-induced transformation 2014-09-12 09:10:19 +00:00
Martin Diehl 5b6f143fb0 intel fortran behaves strange for standard check:
https://software.intel.com/en-us/forums/topic/506419
worked around that
2014-09-11 16:02:05 +00:00
Pratheek Shanthraj d6080d332d updated to be compatible with latest version of petsc.
Please change your PETSC_DIR to /opt/petsc-dev and PETSC_ARCH to ifort (gfort or ifort-debug for gfortran and debugging version respectively on maws01 also available)
2014-09-11 13:28:15 +00:00
Martin Diehl 4aa2c65d1e next try to make the makefile work in all possible cases for FEM and spectral 2014-09-11 09:51:38 +00:00
Luv Sharma 7c60c45f80 corrected typos and some field related changes 2014-09-10 18:26:12 +00:00
Martin Diehl 0af667e94d set search path (VPATH) to ../private/FEM/code only in case FEM solver is the target 2014-09-10 16:27:53 +00:00
Luv Sharma 2f76ecf0e3 fixed error in use statement. 2014-09-10 15:12:14 +00:00
Luv Sharma 7b5cbf4199 now compiling thermal_adiabatic.f90 also and some field related changes 2014-09-10 15:05:28 +00:00
Luv Sharma e8f8655c98 fixed buggy initialization. Now initialising in similar way to constitutive models. 2014-09-10 14:14:03 +00:00
Su Leen Wong 8e6ea7d9c7 Added Ctrans to homogenizedC and strain-induced martensite part to basic states 2014-09-10 12:12:17 +00:00
Luv Sharma dc406a01c0 added funtions to get averged properties at integration points. 2014-09-10 08:37:12 +00:00
Martin Diehl f1e345dcdb fixed compilation for FEM 2014-09-09 15:59:41 +00:00
Franz Roters b6847ebb01 added use mesh statement needed for output of mesh_element 2014-09-09 15:18:49 +00:00
Luv Sharma 923adbc2d3 added auxillary functions: get/put 'Physics' to communicate regularised/ unregularised values between solver and constitutive physics 2014-09-05 16:31:27 +00:00
Luv Sharma a7741457b9 more work on field state. 2014-09-04 16:11:00 +00:00
Philip Eisenlohr 6612603a23 changed naming convention (keeping backward compatibility in numerics.config):
myspectralsolver --> spectral_solver (called "spectralsolver")
myfilter --> spectral_filter (called "spectralfilter")
2014-09-03 19:59:47 +00:00
Luv Sharma 79f572f869 more work on homogenisation new state and introduced field state 2014-09-03 17:12:06 +00:00
Philip Eisenlohr 2435c8d836 dropped "elasticState" (for now) 2014-09-02 20:18:37 +00:00
Philip Eisenlohr f05f47f6e3 optical polishing 2014-09-02 20:16:33 +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
Su Leen Wong 5f1e49c053 Cleaner way of calculating the Bain strain 2014-09-01 10:47:48 +00:00
Su Leen Wong 565e703f49 Minor fix for now 2014-09-01 09:30:33 +00:00