223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Commit Graph

  • 6731659aa4 added hint that extreme values of stress and tangent belong to previous cycle Christoph Kords 2011-02-08 12:43:30 +0000
  • 0a30fb739a now reporting extreme values of stress and tangent, that were returned to marc, in output file Christoph Kords 2011-02-08 10:25:51 +0000
  • 7a7ca1aab7 mpie_spectral and mpie_interface: changed workingdir from pathToLoadFile to pathToGeomFile mpie_spectral and numerics: added switch to prevent pre calculation of gamma_hat. slower, but saves memory 3Dvisualize: started to add support for gmsh (not fully working yet) reconstruct: new version of f2py/Fortran subroutines for output of results from spectral method Martin Diehl 2011-02-07 14:35:42 +0000
  • faba13f7fd * now able to introduce some scatter in the initial single dislocation density (only for nonlocal model!) ; setting the parameter "rhoSglScatter" to a positive value generates a gauss distribution for the dislocation density with standard deviation equal to "rhoSglScatter" * dislocation stress calculation is only done for nonlocal constitution ("/nonlocal/" keyword is present in material.config) Christoph Kords 2011-02-04 15:41:32 +0000
  • 6511b4b5a2 added possibility to select predefined colorscheme from config file Philip Eisenlohr 2011-02-04 09:36:26 +0000
  • b362755457 grain orientation shuffling was skipped for /elemhomo/. no clue why, now dropped... Alankar Alankar 2011-02-03 17:25:58 +0000
  • 05350d573e fixed a bug in new --sort --separation option now keeps data point location as intact as meaningful after applying averaging Philip Eisenlohr 2011-02-01 18:32:20 +0000
  • 669893b397 binary \n will not longer cause problems Martin Diehl 2011-02-01 18:25:40 +0000
  • 9d2b1aa942 specified line endings and activated $Id$ Philip Eisenlohr 2011-02-01 16:43:58 +0000
  • 14d3d0ba99 less output but can now handle 'n/a' in data files without crashing Philip Eisenlohr 2011-02-01 16:43:00 +0000
  • ebf201e89c sorry---wrong file version uploaded before. Philip Eisenlohr 2011-02-01 13:24:19 +0000
  • 958c51fc26 added "--sort" option extending the existing sorting according due to "--separation". This is, for instance, useful for separating /elementhomogeneous/ results by element and later listing the results along x,y,z... Philip Eisenlohr 2011-02-01 13:13:05 +0000
  • d6edb64929 produce VTK visualizations from postResults files Philip Eisenlohr 2011-02-01 10:48:44 +0000
  • fec2c14a4e removed hard-coded tolerances from mpie_spectral and put them to numerics/IO removed storage of full cauchy stres field from mpie_spectral.f90, only average is stored now added cauchy stress and von mises equivalent calculation to spectral post. Martin Diehl 2011-01-31 17:07:42 +0000
  • b72d75ed05 corrected indices in flux density output Christoph Kords 2011-01-28 12:49:58 +0000
  • 0cf1404db0 changed from element 117 (hexahedral with reduced integration) to element 7 (hexahedral with full integration). one-IP elements give trouble in the ping-pong scheme by being frequently/always outdated. this is likely connected to their finite element technology such that even high stiffness values do not prevent the solver from altering displacements in subsequent cycles. remedy is to go for full integration combined with /elementhomogeneous/ phase... Philip Eisenlohr 2011-01-27 14:32:29 +0000
  • c00a42a9a4 fixed last remaining tensor outputs to their transposed versions Philip Eisenlohr 2011-01-27 14:29:06 +0000
  • f84ac44f9a now compiles Fortran sources on the fly Philip Eisenlohr 2011-01-26 15:17:03 +0000
  • aa7d5afd94 changed ifort options back to previous ones as -fast is not working Martin Diehl 2011-01-26 13:33:44 +0000
  • 17b8205e3f reconstruction of geometry is now working. It is directly implemented in python (file spectral_post.py). reconstruction in fortran is not working (file reconstruct.f90) due to some problems with f2py Martin Diehl 2011-01-26 12:56:52 +0000
  • 05d4d5fef2 * kinetics now according to Mohles; enables solid solution hardening * changed example in material.config accordingly Christoph Kords 2011-01-26 10:17:42 +0000
  • 6ac06af517 sorry, forgot to alter the most important part: post_results..! eventually, the t16 has defgrads and friends in correct format: du_i/dx_j is listed as linear array (1,1),(1,2),(1,3),(2,1),... Philip Eisenlohr 2011-01-24 16:23:37 +0000
  • 38409ebc71 included both new voronoi_xxx scripts Philip Eisenlohr 2011-01-24 16:21:49 +0000
  • ee7a8ad52a general speed up spectral format records physical coordinates tessellation now based on phys coords instead of discretization. Philip Eisenlohr 2011-01-24 16:21:17 +0000
  • 4ee40df5ba transposed writing into t16 file (and screen/debug) did not work with just stating array(1:3,:)... fixed by using math_transpose3x3 (and similar) instead. Philip Eisenlohr 2011-01-24 15:02:19 +0000
  • 794aeb253b crystallite and CPFEM now print/output tensors in ij notation (used to be, particularly in the t16, ji, i.e. transposed...) Philip Eisenlohr 2011-01-20 19:25:45 +0000
  • e721f2b654 sight restructuring of algorithm -- same functionality Philip Eisenlohr 2011-01-20 19:22:22 +0000
  • 52f06009be rewrote voronoi tessellation. the VE generation is now divided into two parts: -generation of random positions (in a 1x1x1 parameter space) and random angels -actual voronoi tessellation for the generation of the positions, a resolution can be specified to ensure that the tessellation works at that resolution. The tessellation will always work on larger resolutions but not necessarily at smaller ones. Martin Diehl 2011-01-20 17:21:38 +0000
  • 5990d94ca3 added description and $ID$ Philip Eisenlohr 2011-01-19 15:14:05 +0000
  • 581955f68d parsePostFile screwed up the position of outputs for a grain number larger than 1 Christoph Kords 2011-01-14 13:15:19 +0000
  • b546d3bb91 some corrections regarding Multicore-support on mpie_spectral.f90 and makefile slight changes on mpie_spectral.f90 regarding file output Martin Diehl 2011-01-12 17:02:42 +0000
  • a662a28a97 added support for spectral method result files Martin Diehl 2011-01-12 16:55:56 +0000
  • 062c6b3d54 * allow zero cutoff radius * added missing error codes Christoph Kords 2011-01-12 12:36:48 +0000
  • d92f992eb7 if nonlocal ip has local neighbor assume perfect transparency of the (perhaps virtual) phase boundary. this is only a quick fix! need to rethink the way constitutive and real boundaries are treated! Christoph Kords 2011-01-11 18:00:29 +0000
  • 6a3dd78259 fluxes are disabled when keyword "/nonlocal/" is omitted in the material.config (this is a way to use the constitutive_nonlocal model as a local dislocation based model) Christoph Kords 2011-01-11 14:55:36 +0000
  • 322cbd2597 changed tools for voronoi tessellation to match new specification (file extension and resolution). Also renamed the files, because the space characters cause trouble Martin Diehl 2011-01-07 14:37:05 +0000
  • ccc6aac10b did not delete mpie_spectral2.f90 properly Martin Diehl 2011-01-07 13:11:32 +0000
  • 71fb5eedf2 changed makefile: material subroutines are now compiled using openmp, enabled linux threads for fftw3 (also did some slight changes to enable that in the fortran source code) renamed mpie_spectral2.f90 to mpie_spectral2d.f90 (testing file, not properly working at the moment) changed file extension and variable names in mpie_spectral.f90 and mpie_spectral_interface.f90 from "mesh" to "geom". Removed direct output from mpie_spectral.f90, all output is now base on materialpoint_results(:,1,:) Martin Diehl 2011-01-07 12:56:45 +0000
  • 9447330ab0 ATTENTION repository has been moved to new server! Now at musvn1:/home/svn/repos/cpfem Franz Roters 2011-01-05 16:11:51 +0000
  • 20c00288b8 ---had some problems with svn, hope everything is ok now--- added new tools to generate colormaps for paraview and gmsh, written in python removed old fortran colormap generator. Martin Diehl 2011-01-05 14:53:31 +0000
  • 3b0eeb9956 added files for geometry reconstruction from given deformation gradient. Martin Diehl 2010-12-22 10:59:54 +0000
  • 8d6ddcf28b restructured (or structured at all) reading of header, first try of reading in deformation gradient Martin Diehl 2010-12-09 15:53:52 +0000
  • ac676623ca reading of header is now working Martin Diehl 2010-12-03 14:57:28 +0000
  • 2f47d4b09b added conversion of 'resolution' in result file to integer array Martin Diehl 2010-12-02 14:53:17 +0000
  • 096a4803b0 some polishing concerning openMP initialization number of threads is now written to *.out file Franz Roters 2010-12-02 11:04:29 +0000
  • ab82dbc99c continued working on reading header of output file Martin Diehl 2010-11-30 14:26:30 +0000
  • 9d32bc7ae6 added output of complete mesh to voronoi fast.f90 first version of post processing for spectral method: spectral_post.py Martin Diehl 2010-11-29 10:08:53 +0000
  • 36f27edd75 1. Adding more comments 2. Changing 'subroutine' into 'function' for optimizing the computation time Denny Tjahjanto 2010-11-26 11:50:20 +0000
  • 6b199e7830 allow 2digit values to be read in for MPIE_NUM_THREADS Christoph Kords 2010-11-25 10:34:55 +0000
  • 4216a3e066 corrected name of include file Franz Roters 2010-11-24 07:55:08 +0000
  • c07a91d19e removed -save option from the compile statement using OpenMP Franz Roters 2010-11-23 16:25:58 +0000
  • ea6a6c9503 * environment variable MPIE_NUM_THREADS now determines number of threads used in parallel regions * default value of the OMP_NUM_THREADS variable has to be restored at the end of mpie subroutine, since marc also seems to use and change(!) this * usage: "export MPIE_NUM_THREADS=<number of threads>" to set variable in shell, then restart mentat and compile with option 3 (at the moment this does only work on ws 6, since all other workstations use compiler option "-save"; this puts all local variables by default in static memory, which is a killer for parallelization!) Christoph Kords 2010-11-19 17:45:27 +0000
  • b9b226548c * enclose wall time measurement for LpAndItsTangent and count of leapfrogbreaks by CRITICAL construct * better use SINGLE (having an implicit barrier at the end) instead of MASTER construct * deleted all explicit BARRIERs after do loops since parallel loop construct implies barrier at the end Christoph Kords 2010-11-19 17:29:29 +0000
  • 15728d3305 enclose wall time measurement for dotState and dotTemperature function by CRITICAL construct Christoph Kords 2010-11-19 17:19:03 +0000
  • b2cc34349f corrected mistake in nodeCoordinate assignment, bug-fix for option " --filter 'x,y,z' " Philip Eisenlohr 2010-11-12 08:56:04 +0000
  • 4f267f4375 parallelization now works properly (not yet tested extensively though): * had to add some BARRIER constructs * only the master thread is allowed to increase the state counter Christoph Kords 2010-11-11 16:16:05 +0000
  • 14594dc9e2 * OpenMP does not like unbounded array subscripts in function calls * some write statements were still not enclosed by CRITICAL construct, some needed a MASTER construct Christoph Kords 2010-11-11 13:14:53 +0000
  • f3a7e323ea pretty printing of hex slip/twin systems Philip Eisenlohr 2010-11-09 13:20:33 +0000
  • bca951325d * cpfem will now produce a warning if all components of the stiffness tensor are close to zero Christoph Kords 2010-11-04 18:18:01 +0000
  • 54748bb6a5 cycleCounter reset was 0 now -1 ==> occurence of jacobian calculation is now correct, i.e. always calc at start of new increment Philip Eisenlohr 2010-11-04 18:15:50 +0000
  • e561358a90 presentation to illustrate Marc/Mentat interface for restarting capability (thx to Denny!) Philip Eisenlohr 2010-11-04 15:31:29 +0000
  • 405d5529e7 * checked and corrected parallelization of code. compiles now successfully, but simply aborts computation with first parallel directive without any comment :-((( * also put a call to constitutive_microstructure at the start of each crystallite_integration subroutine like it was before. need that for nonlocal model in case of crystallite cutback Christoph Kords 2010-11-03 17:22:48 +0000
  • 0dd99cb965 lattice: (re)introduced _symmetryType function to replace unsafe lookup array Philip Eisenlohr 2010-11-03 14:58:11 +0000
  • 763c20b302 Introducing the capability to restart jobs that crashed in the middle of sims. At the moment, this feature is exclusive for Marc. Denny Tjahjanto 2010-11-03 14:39:18 +0000
  • be265aef37 slight improvement of remaining time estimate... Philip Eisenlohr 2010-11-02 16:55:56 +0000
  • 796bffee2e 1) speed-up of elementTensor output 2) added estimate of remaining time Philip Eisenlohr 2010-11-02 15:45:23 +0000
  • e80e055c75 added three small programs useful for spectral method (fortran sourececode) voronoi fast.f90 and voronoi small memory.f90 are two variants to do a voronoi tessellation and write the result to a mesh file that can be interpreted by mpie_spectral.f90. Difference is the memory management resulting in one faster and one memory saving version. voronoi fast.f90 has also the ability to write out a file for Ricardo Lebensohns spectral code. colormap.f90 is a simple code that can be used to generate colormaps for gmsh. Martin Diehl 2010-10-31 16:06:10 +0000
  • 4c20daa00d state loop count was wrong for adaptive euler integration Christoph Kords 2010-10-28 12:22:17 +0000
  • 17812c1f9e added raw output, gmsh output is still included Martin Diehl 2010-10-27 17:15:49 +0000
  • 9320d48305 in RKCK45 state integration: dotState was not saved in last Runge-Kutta step Christoph Kords 2010-10-27 08:48:04 +0000
  • 74743c0856 little typo fixed... Philip Eisenlohr 2010-10-26 18:26:04 +0000
  • 487466dd88 scripts are now aware of symbolic links to themselves... Philip Eisenlohr 2010-10-26 16:26:55 +0000
  • 18344b4e65 forgot linebreaks and executable flag... Philip Eisenlohr 2010-10-26 15:32:03 +0000
  • 5f64cacc74 generate alias in ~/bin for processing scripts Philip Eisenlohr 2010-10-26 15:02:20 +0000
  • f5f3ef5ecc * have to reset dotState to zero before call to collectDotState (only important for nonlocal calculation) * in Fixed Point Iteration: update dependent states after state preguess was missing; on the other hand, the first call to constitutive_microstructure was obsolete Christoph Kords 2010-10-26 14:04:33 +0000
  • 4439282239 no change to previous version, just forget to mention one important change of the last revision: * fluxes are now again calculated and distributed only! by the originating material point. this means that the central MP might change the dotState of its neighbor. have to see whether locks slow down parallel computation Christoph Kords 2010-10-26 13:49:58 +0000
  • 5be6604fd6 included definitions of "material" and "geometry" properties after mesh generation... Philip Eisenlohr 2010-10-26 13:45:23 +0000
  • 62d06001ea * corrected compatibility for screws (always positive) * detection of grain boundary in constitutive_nonlocal_microstructure with the help of transmissivity * enforce positive densities in constitutive_nonlocal_microstructure (needed because dotState does not create cutbacks for negative densities anymore) * reset single mobile densities below certain threshold to zero (also done in constitutive_nonlocal_microstructure) * constitutive_nonlocal_kinetics only gets local state variable as input, no need for the entire array here * dv_dtau is always positive * multiplication is only active when there is already some initial density of the respective type Christoph Kords 2010-10-26 13:42:18 +0000
  • 366d52bd71 * former "relevantRho" and "relevantResistance" is renamed to "atol_rho" and "atol_resistance" and now used as an absolute tolerance for the state residuum. before it was used rather as significant state, so whenever the state dropped below that value it was considered converged. (In dislotwin and titanmod constitutive law there is only one value atol_rho which is used for all the states, though the state array consists not only of densities. We need further parameters here!) * somewhat simplified convergence check for adaptive euler and runge-kutta Christoph Kords 2010-10-26 13:16:37 +0000
  • d965f14f90 generate hexahedral mesh in Mentat based on a data file of Ricardo Lebensohn's "txfft" format (phi1, Phi, phi2, x, y, z, grainId, phaseId) Only grainId is used, physical dimension and subdivisions need to be provided on command line. improvement possible by parsing x, y, z columns to get those directly. Philip Eisenlohr 2010-10-25 17:57:57 +0000
  • 26cb618b6d changed calculation to small strain/cauchy stress Martin Diehl 2010-10-20 11:27:10 +0000
  • da35eaa8a1 changed spectral algorithm to version proposed by Ricardo 2001 (using tau_field instead of stress_field). Still using largestrain-formulation Martin Diehl 2010-10-20 09:27:47 +0000
  • 3837dad51e fixed bug in mpie_spectral.f90 concerning spectral method (now deltaF hat get the right value) added mpie_spectral2.f90, a version that should get the new algorithm proposed in 2010. until now, it is the same as mpie_spectral.f90 (large strain formulation by suquet et al) but with c2c, c2c FFT Martin Diehl 2010-10-20 08:59:00 +0000
  • e49de75fe3 stiffness calculation now only needs a single function call to the respective integrator method instead of one call per grain (which seems to heavily slow down the computation). also no special treatment for non-local material points anymore. Christoph Kords 2010-10-15 14:57:13 +0000
  • fffe731447 * had to correct nonlocal slip system compatibility and flux calculation. last update gave wrong results. Christoph Kords 2010-10-15 13:19:26 +0000
  • 61f8a5fcbe corrected UpperCase "SolidSolutionStrength" in config file parsing. Philip Eisenlohr 2010-10-15 12:59:00 +0000
  • 3c502561ee added matrix multiplication 3333x33 to math.f90 added some parameters for spectral method to numerics.f90 (tolerance) changed error message concerning spectral method in IO.f90 corrected calculation of stress BC in mpie_spectral.f90 Martin Diehl 2010-10-13 16:04:44 +0000
  • 724960686f * in nonlocal model: dislocation flux now with valid description also for large angle grain boundaries; transfer of dislocations from one slip system to the other according to new variable "constitutive_nonlocal_compatibility" which depends on the angle between slip plane normals and slip directions and is updated once per cycle in function "crystallite_orientations" * reactivated debugging functionality for "non-standard" integration methods Christoph Kords 2010-10-12 13:08:54 +0000
  • efd92d9b51 allocation of special state variables for "non-standard integration" was erroneous. used to produce segmentation fault for integration method 4 and 5. Christoph Kords 2010-10-12 12:11:24 +0000
  • fce7590c17 * now able to choose method for state integration (integrator and integratorStiffness in numerics.config); standard is "Fixed Point Iteration", which is basically the same as before; others available are "Explicit Euler", "AdaptiveEuler", "Classical Runge-Kutta" and "Runge-Kutta Cash-Karp" Christoph Kords 2010-10-01 12:18:49 +0000
  • 9927cd7adb changed format of resolution. instead of reading exponent of two, now the number given equals directly to the number of FPs. Allows us to use arbitrary number of FPs. example: OLD: reading 5 results in 2^5=32 FPs NEW: reading 33 results in 33 FPs Martin Diehl 2010-10-01 10:42:15 +0000
  • 2cc7bc38d7 Fixed a missing ','. Assigned plane and direction names. Alankar Alankar 2010-10-01 06:32:44 +0000
  • 1ae252c1ce altered hexagonal slip systems (needs comparison to former set, which is potentially wrong) Philip Eisenlohr 2010-09-30 17:31:19 +0000
  • b85f9b0770 material.config: changed texture of microstructure 1 from 1 to 3, because 1 needs a ODF file which is not in the repository and makes the one_element_model.proc complain. Claudio Zambaldi 2010-09-30 12:55:38 +0000
  • a8dfdc6487 added _mul33xx33 double contraction: A_ij * B_ij = scalar Philip Eisenlohr 2010-09-30 09:32:49 +0000
  • cb4b4afc88 corrected header for math_qRot Christoph Kords 2010-09-30 08:46:58 +0000
  • b259fbd9c6 1) added distribution of leapfrog breaks 2) lattice_symmetryType is now a function (former lookup array was buggy) 3) stricter check of state var values (>0!) and memory deallocation done Philip Eisenlohr 2010-09-30 07:31:53 +0000
  • e562df35a9 Added more state dependent variables in post results. Alankar Alankar 2010-09-29 06:35:08 +0000
  • 3a4743f1d3 fixed bug in grain naming scheme. occurred for gaps in consecutive numberings... Philip Eisenlohr 2010-09-28 16:56:49 +0000
  • 5a51082977 polishing, changed calculation of stress BC. prevent updating of gamma hat in case of perfect plasticity Martin Diehl 2010-09-24 13:27:53 +0000