Commit Graph

1300 Commits

Author SHA1 Message Date
Martin Diehl c399a06c97 tried to implement an inverse reconstruction, meaning to calculate an average deformation gradient from the 8 cornes of a node in reference and current configuration. Not working yet.
polishing, mainly in fft reconstruction.
added fftw library file and fortran file for type specification
added folder references with literature for visualization/geometry reconstruction
2011-02-21 16:30:18 +00:00
Philip Eisenlohr 34336bc659 fixed bug when no argument was supplied 2011-02-18 09:02:29 +00:00
Martin Diehl 438bf95105 added ft-based reconstruction of deformed configuration to postprocessingMath.f90 and postprocessingMath.pyf
also added function to calculate average of tensor
some polishing in mpie_spectral.f90, added sanity check to see im resolution is power of 2
2011-02-14 17:21:31 +00:00
Philip Eisenlohr 874d5067c7 switches drawing of servo links temporarily off to speed up processing. 2011-02-11 11:29:58 +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
Philip Eisenlohr 65ae7de42b fixed bug for missing file name 2011-02-09 09:18:21 +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
Philip Eisenlohr 6511b4b5a2 added possibility to select predefined colorscheme from config file 2011-02-04 09:36:26 +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 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
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 f84ac44f9a now compiles Fortran sources on the fly 2011-01-26 15:17:03 +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
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 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 a662a28a97 added support for spectral method result files 2011-01-12 16:55:56 +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 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
Martin Diehl 2f47d4b09b added conversion of 'resolution' in result file to integer array 2010-12-02 14:53:17 +00:00
Martin Diehl ab82dbc99c continued working on reading header of output file 2010-11-30 14:26:30 +00:00
Martin Diehl 9d32bc7ae6 added output of complete mesh to voronoi fast.f90
first version of post processing for spectral method: spectral_post.py
2010-11-29 10:08:53 +00:00
Philip Eisenlohr b2cc34349f corrected mistake in nodeCoordinate assignment, bug-fix for option " --filter 'x,y,z' " 2010-11-12 08:56:04 +00:00
Philip Eisenlohr be265aef37 slight improvement of remaining time estimate... 2010-11-02 16:55:56 +00:00
Philip Eisenlohr 796bffee2e 1) speed-up of elementTensor output
2) added estimate of remaining time
2010-11-02 15:45:23 +00:00
Martin Diehl 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.
2010-10-31 16:06:10 +00:00
Philip Eisenlohr 74743c0856 little typo fixed... 2010-10-26 18:26:04 +00:00
Philip Eisenlohr 487466dd88 scripts are now aware of symbolic links to themselves... 2010-10-26 16:26:55 +00:00
Philip Eisenlohr 18344b4e65 forgot linebreaks and executable flag... 2010-10-26 15:32:03 +00:00
Philip Eisenlohr 5f64cacc74 generate alias in ~/bin for processing scripts 2010-10-26 15:02:20 +00:00
Philip Eisenlohr 5be6604fd6 included definitions of "material" and "geometry" properties after mesh generation... 2010-10-26 13:45:23 +00:00
Philip Eisenlohr 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.
2010-10-25 17:57:57 +00:00
Philip Eisenlohr 3a4743f1d3 fixed bug in grain naming scheme. occurred for gaps in consecutive numberings... 2010-09-28 16:56:49 +00:00
Philip Eisenlohr 23746e7e78 chopped line beginnings
introduced check for PIL, if this fails, no image will be output
2010-09-23 07:59:52 +00:00
Philip Eisenlohr d8a4cd1739 added optional argument --baseindex to specify what colorset should be influenced. identify colors, for instance, start at idx 214... 2010-09-13 07:31:27 +00:00
Alankar Alankar 006dce7d8a printing euler angles (and defgrads and orientations) for individual ips now works 2010-09-08 11:55:57 +00:00
Philip Eisenlohr 4a31b175da 1) introduced default colors
2) complaints are more verbose
2010-09-06 16:51:11 +00:00
Philip Eisenlohr 3b671faae7 added some flexibility (and bug-fixing) to the pre/post processing. the file 'MSCpath' contains the LOCAL path to the MSC installation (/msc in our case) 2010-08-17 18:21:22 +00:00
Philip Eisenlohr 94efd57663 started to put Marc/Mentat related scripts to corresponding post/pre processing folders. We should move development from my Code folder to here (sorry, lost history then...) 2010-08-16 20:47:27 +00:00