Commit Graph

4354 Commits

Author SHA1 Message Date
Chen Zhang 43b665ba48 add getAverageOrientation() to Orientation class as class method. 2015-04-02 19:15:09 +00:00
Martin Diehl 86f39de462 remove non-utf symbol 2015-04-02 16:13:00 +00:00
Haiming Zhang d8a99b23bf 1.compact all the criteria into a single Class;
2.the first thorough check of the script;
3.add the option of exponents for all non-quadratic yield criteria, now the user 1)can specify the exponent, for example, m=6 for Barlat 1991, or 2)see the exponent as an undetermined parameters;
4.add the pre-specified bounds for all criteria;
5.add the user defined equivalent stress for some anisotropic yield criteria
2015-04-02 12:08:55 +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 32f9e7f9d8 now giving just a 'Remark' when h5py is not found 2015-04-01 14:43:34 +00:00
Martin Diehl 92a67e1a9d geom_fromAng.py crashed when found empty line (eg. at the end of the file), now fixed 2015-04-01 14:42:25 +00:00
Martin Diehl a1956436be fixed bug introduced last time when searching for keywords in spectralOut header
made reading in of header more robust by explicitly using the information given by fortran specifies the data length
introduced a legacy mode again, now it is used to read old spectralOut files prior to rev 4017
2015-04-01 10:02:21 +00:00
Martin Diehl 3a42a18b8d Abaqus example only runs with perturbed jacobian 2015-03-31 17:58:12 +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 8b1b397be6 cleaning up 2015-03-30 13:32:25 +00:00
Pratheek Shanthraj 41b2d68856 some polishing. usage now similar to postResults script 2015-03-30 12:56:39 +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
Tias Maiti f5801601d8 fixed UnboundLocalError (calling variable before assignment) in “inversePole” function 2015-03-29 19:16:45 +00:00
Philip Eisenlohr 38943c1e94 improved help 2015-03-29 15:37:40 +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 913c5347a5 removed handling for old legacy format (we might introduce a new one) 2015-03-29 09:08:42 +00:00
Martin Diehl 9e824d029e renamed to new names: spectralsolver and spectralfilter 2015-03-28 21:14:11 +00:00
Martin Diehl d733db7638 removed old solver 2015-03-28 17:47:59 +00:00
Martin Diehl 1b2e4c0239 removed old solver 2015-03-28 17:41:40 +00:00
Martin Diehl 66d4ff097f removed obsolete comment 2015-03-28 11:44:17 +00:00
Martin Diehl f54ff567a0 improved style 2015-03-28 11:36:52 +00:00
Tias Maiti 0e7325b27d added function — “inversePole” axis rotated according to orientation (using crystal symmetry to ensure location falls into SST) 2015-03-28 07:43:49 +00:00
Haiming Zhang 72c2ead277 polishing, replace variables calculation with arrays calculation. 2015-03-27 17:58:49 +00:00
Martin Diehl 151d4ff75b fixed bug (crash) in case (meaningless) strings found after part or section key.
DAMASK itself is fine with it, and it can be used for commenting.
thanks to B. Bode for figuring out
2015-03-27 16:23:21 +00:00
Martin Diehl 8a30021202 cleaned up core modul 2015-03-27 14:10:18 +00:00
Martin Diehl 004375cc79 removed table rewind causing problem with STDIN/STDOUT 2015-03-27 13:57:34 +00:00
Pratheek Shanthraj 5e61ed0382 removed some no longer used tests for the basic spectral solver 2015-03-27 07:46:51 +00:00
Pratheek Shanthraj e9a1ffd82e switched from binary to h5 for fem output. improved post processing script to identify results by the variable name instead of offset number 2015-03-27 07:31:58 +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
Martin Diehl 69b8e02a3a removed unwanted ASCII table data rewind causing problem for STDIN 2015-03-26 21:15:59 +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
Martin Diehl 643f7fe934 some comments 2015-03-26 05:04:12 +00:00
Pratheek Shanthraj bf108ee09f updated to current petsc directories in the FEM makefiles 2015-03-25 16:23:37 +00:00
Pratheek Shanthraj 6e6f96f926 updated spectral example to show current way of starting a parallel simulation 2015-03-25 16:16:01 +00:00
Pratheek Shanthraj fe4373a9f9 updated post processing script for MPI IO spectralOut file 2015-03-25 16:09:56 +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
Test User 8b09d78ed1 Added report of SUCCESSFUL test for Revision 4012 2015-03-25 02:19:36 +00:00
Martin Diehl 1a3bfc7d1f updated according to new convention for ip and nodes 2015-03-24 13:52:40 +00:00
Test User df566f6466 Added report of SUCCESSFUL test for Revision 4010 2015-03-24 02:27:00 +00:00
Martin Diehl d03aef882b adjusted to new location 2015-03-20 15:55:08 +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
Test User 8b04854c23 Added report of SUCCESSFUL test for Revision 4005 2015-03-17 04:32:27 +00:00
Haiming Zhang 226381586c compact the code of Yld2004-18p and Karafillis-Boyce, use array instead of components to calculate the partial derivative. 2015-03-16 16:42:26 +00:00