223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity
4,336 Commits 5 Branches 20 Tags 292 MiB
Commit Graph

4336 Commits

Author SHA1 Message Date
Philip Eisenlohr 4bcc80911e changed grid to float from (wrong) integer type. 
row and column probability density normalization now takes place BEFORE limiting data range.
 2015-04-13 19:55:28 +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 d008b566b4 ensure to bring all seeds within target range 2015-04-10 20:47:20 +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 5526d2a721 updated test according to improvments 2015-04-10 18:46:26 +00:00
Martin Diehl 9d3cd8f997 added options do define number of maximum seeds to move simulateneously 2015-04-10 18:08:17 +00:00
Martin Diehl 4aa63e1ab9 handling of seeds files without microstructure (but eulerangles) 2015-04-10 16:59:08 +00:00
Martin Diehl 09bfd1e18e more verbose error handling for euler angles limits 
ensuring that ROUNDED euler angles are not above the limits
 2015-04-10 14:50:23 +00:00
Martin Diehl e01e270300 improved search for double microstructures (was N^2 bogosort) in geom_fromAng and introduced it in geom_fromVPSC 2015-04-10 12:57:25 +00:00
Martin Diehl b08141d2b3 less strict tolerances for basic PETSC compared to basic 2015-04-09 10:07:44 +00:00
Tias Maiti 5257a2161f added options for x and y normalization 2015-04-09 06:45:21 +00:00
Martin Diehl d8260b12df updated ignore list 2015-04-08 18:49:47 +00:00
Martin Diehl 57a2c90d30 corrected according to vtk specification and remove an unused option 2015-04-08 18:48:35 +00:00
Martin Diehl 7d0a1ebea6 math_transpose33 not used 2015-04-08 18:33:08 +00:00
Martin Diehl db657a7a9e hdf5 not needed (either from PETSc or h5py for python) 2015-04-08 18:04:01 +00:00
Haiming Zhang 976647b9e4 add the convergence criterion, the method is : 
1. calculate the L2 norm of the residual of all the stress points
2. store the relative errors of the L2 norm
3. if the standard deviation of the relative errors of last five fittings is less than 0.05, that it is considered that the relative errors is stabilized, so the fitting is finished.
 2015-04-08 17:48:26 +00:00
Haiming Zhang b77768fd4d all the criteria support plane stress 2015-04-08 16:24:20 +00:00
Martin Diehl 0544706c7b corrected help for legacy format 2015-04-08 12:51:40 +00:00
Martin Diehl 6e650a42b4 some comments for discussion 2015-04-07 07:42:46 +00:00
Haiming Zhang ab89c84ec8 Correct the message of last change 
1. rename BBC2003 as BBC2000;
2. add the real BBC2003 yield criterion, BBC2003 works; 
3. now the BBC family yield criteria include: BBC2000, BBC2003, BBC2005
 2015-04-03 21:01:37 +00:00
Haiming Zhang 412028f33f 1. rename BBC2003 as BBC2000; 
2. add the real BBC2003 yield criterion, BBC2003 works;
3. now the BBC family yield criteria include: BBC2000, BBC2003, BBC2008
 2015-04-03 20:59:23 +00:00
Haiming Zhang cd2a744db9 add Yld2000 yield criterion, Yld2000 works 2015-04-03 16:20:03 +00:00
Haiming Zhang 2b9964bba5 add Barlat-Lian 1989 yield criterion, Barlat1989 works. 2015-04-03 12:31:13 +00:00
Haiming Zhang 0097d39873 fix some type errors; add more comments 2015-04-02 19:32:48 +00:00
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