Commit Graph

858 Commits

Author SHA1 Message Date
Philip Eisenlohr 171f5375f1 shearrate output now with sign (no |abs| as before) 2011-02-11 10:19:41 +00:00
Martin Diehl 3d7fad6ba9 mpie_spectral.f90: changed fourier transormation, now using the advanced interface to do the tranformation both ways with just one call. now also using the same variable for pk-stress in both domains and the change of deformation gradient in both domains.
postprocessing: renamed name of python/f2py modul from "reconstruct" to "postprocessingMath", added some numerical operations to use for postprocessing.
2011-02-09 17:47:28 +00:00
Christoph Kords 803e1a8c05 In order to help preventing further memory leaks all array sections now have an explicit instead of assumed shape, e.g. Fe(1:3,1:3) instead of Fe(:,:). 2011-02-09 13:12:46 +00:00
Philip Eisenlohr 65ae7de42b fixed bug for missing file name 2011-02-09 09:18:21 +00:00
Christoph Kords 466ced911e * check convergence for nonlocals also in stiffness calculation mode 2011-02-09 08:39:07 +00:00
Christoph Kords 921967c2ba * found memory leak in nonlocal model! "constitutive_nonlocal_rhoDotFlux" has to have upper index bound of 10 for second dimension, not 8.
* need to recalculate dislocation velocity in postResults, otherwise we take values of last perturbed state! So the following outputs were up to now showing the perturbed state: shearrate, dislocation velocity, all density rates!
2011-02-09 08:28:47 +00:00
Christoph Kords 17a6c0e59b enhanced comment on dislocation kinetics 2011-02-08 12:48:50 +00:00
Christoph Kords 6731659aa4 added hint that extreme values of stress and tangent belong to previous cycle 2011-02-08 12:43:30 +00:00
Christoph Kords 0a30fb739a now reporting extreme values of stress and tangent, that were returned to marc, in output file 2011-02-08 10:25:51 +00:00
Martin Diehl 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
2011-02-07 14:35:42 +00:00
Christoph Kords 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)
2011-02-04 15:41:32 +00:00
Philip Eisenlohr 6511b4b5a2 added possibility to select predefined colorscheme from config file 2011-02-04 09:36:26 +00:00
Alankar Alankar b362755457 grain orientation shuffling was skipped for /elemhomo/. no clue why, now dropped... 2011-02-03 17:25:58 +00:00
Philip Eisenlohr 05350d573e fixed a bug in new --sort
--separation option now keeps data point location as intact as meaningful after applying averaging
2011-02-01 18:32:20 +00:00
Martin Diehl 669893b397 binary \n will not longer cause problems 2011-02-01 18:25:40 +00:00
Philip Eisenlohr 9d2b1aa942 specified line endings and activated $Id$ 2011-02-01 16:43:58 +00:00
Philip Eisenlohr 14d3d0ba99 less output but can now handle 'n/a' in data files without crashing 2011-02-01 16:43:00 +00:00
Philip Eisenlohr ebf201e89c sorry---wrong file version uploaded before.
now correct: 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...

slight improvement in file-error reporting
2011-02-01 13:24:19 +00:00
Philip Eisenlohr 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... 2011-02-01 13:13:05 +00:00
Philip Eisenlohr d6edb64929 produce VTK visualizations from postResults files 2011-02-01 10:48:44 +00:00
Martin Diehl 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.
2011-01-31 17:07:42 +00:00
Christoph Kords b72d75ed05 corrected indices in flux density output 2011-01-28 12:49:58 +00:00
Philip Eisenlohr 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... 2011-01-27 14:32:29 +00:00
Philip Eisenlohr c00a42a9a4 fixed last remaining tensor outputs to their transposed versions 2011-01-27 14:29:06 +00:00
Philip Eisenlohr f84ac44f9a now compiles Fortran sources on the fly 2011-01-26 15:17:03 +00:00
Martin Diehl aa7d5afd94 changed ifort options back to previous ones as -fast is not working 2011-01-26 13:33:44 +00:00
Martin Diehl 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
2011-01-26 12:56:52 +00:00
Christoph Kords 05d4d5fef2 * kinetics now according to Mohles; enables solid solution hardening
* changed example in material.config accordingly
2011-01-26 10:17:42 +00:00
Philip Eisenlohr 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),... 2011-01-24 16:23:37 +00:00
Philip Eisenlohr 38409ebc71 included both new voronoi_xxx scripts 2011-01-24 16:21:49 +00:00
Philip Eisenlohr ee7a8ad52a general speed up
spectral format records physical coordinates
tessellation now based on phys coords instead of discretization.
2011-01-24 16:21:17 +00:00
Philip Eisenlohr 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.
2011-01-24 15:02:19 +00:00
Philip Eisenlohr 794aeb253b crystallite and CPFEM now print/output tensors in ij notation (used to be, particularly in the t16, ji, i.e. transposed...)
math reports random seed used for rnd() generation
2011-01-20 19:25:45 +00:00
Philip Eisenlohr e721f2b654 sight restructuring of algorithm -- same functionality 2011-01-20 19:22:22 +00:00
Martin Diehl 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.
2011-01-20 17:21:38 +00:00
Philip Eisenlohr 5990d94ca3 added description and $ID$ 2011-01-19 15:14:05 +00:00
Christoph Kords 581955f68d parsePostFile screwed up the position of outputs for a grain number larger than 1 2011-01-14 13:15:19 +00:00
Martin Diehl b546d3bb91 some corrections regarding Multicore-support on mpie_spectral.f90 and makefile
slight changes on mpie_spectral.f90 regarding file output
2011-01-12 17:02:42 +00:00
Martin Diehl a662a28a97 added support for spectral method result files 2011-01-12 16:55:56 +00:00
Christoph Kords 062c6b3d54 * allow zero cutoff radius
* added missing error codes
2011-01-12 12:36:48 +00:00
Christoph Kords 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! 2011-01-11 18:00:29 +00:00
Christoph Kords 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) 2011-01-11 14:55:36 +00:00
Martin Diehl 322cbd2597 changed tools for voronoi tessellation to match new specification (file extension and resolution). Also renamed the files, because the space characters cause trouble
improved reconstruct.f90 and spectral_post.py, both files work now under linux

added make_reconstruct.py, small shell script for using f2py
2011-01-07 14:37:05 +00:00
Martin Diehl ccc6aac10b did not delete mpie_spectral2.f90 properly 2011-01-07 13:11:32 +00:00
Martin Diehl 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,:)
2011-01-07 12:56:45 +00:00
Franz Roters 9447330ab0 ATTENTION repository has been moved to new server!
Now at musvn1:/home/svn/repos/cpfem
2011-01-05 16:11:51 +00:00
Martin Diehl 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.

removed test.py (not longer needed) and the python module reconstruct.pyd (not running under linux)
2011-01-05 14:53:31 +00:00
Martin Diehl 3b0eeb9956 added files for geometry reconstruction from given deformation gradient.
reconstruct.f90 is a fortran source file with comments to use with f2py
reconstruct.pyd is the compiled python module

test.py is a test file to check if the reconstruction works (unfortunately i does not at the moment, probably because there are type conflicts between fortran and python)
2010-12-22 10:59:54 +00:00
Martin Diehl 8d6ddcf28b restructured (or structured at all) reading of header, first try of reading in deformation gradient 2010-12-09 15:53:52 +00:00
Martin Diehl ac676623ca reading of header is now working 2010-12-03 14:57:28 +00:00