Commit Graph

1131 Commits

Author SHA1 Message Date
Christoph Kords 3b788720e5 added crystallite output "neighboringelement" and "neighboringip"
some polishing in crystallite_postResults
2013-05-08 12:02:30 +00:00
Martin Diehl b2a3c9235b fixed some warnings issued by gfortran (type conversion, not needed use statements) 2013-05-08 09:23:47 +00:00
Christoph Kords 6e64b5c476 better not allow for deads to contribute to flux, since not clear at which interface they got stuck 2013-05-07 18:47:17 +00:00
Christoph Kords 1c2c9b8436 ip neighborhood was screwed up for spectral solver as of rev 2332 due to new "mesh_spectral_ipNeighborhood": used to have fastest index in z direction although it should be in x direction!
also removed output of node twins in "mesh_tell_statistics" when using spectral solver, since not allocated
2013-05-07 17:50:58 +00:00
Christoph Kords fb5f94812f corrected format strings in mesh_read_meshfile and mesh_write_meshfile 2013-05-07 13:27:51 +00:00
Christoph Kords c7f987a3c1 "unitlength" parameter, which determines the physical size of the mesh, now available as a global mesh variable "mesh_unitlength" and written to the .mesh file during init.
Hence, it is available to the marc_deformedGeometry.py script via "mesh_init_postprocessing()" and "mesh_get_unitlength()", so no need for setting the scaling of the displacement vectors explicitly through an option; now displacements and nodal positions are always consistent
2013-05-07 13:06:29 +00:00
Franz Roters d2416abdb4 introduced new keyword 'rotation' to specify texture rotation
ATTENTION: so far it is read but NOT yet applied!!!
waiting for updated code for orientation assignement
2013-05-02 08:35:37 +00:00
Martin Diehl b98ce3d22d fixed out of range bug introduced with singleRun changes in rev. 2356 2013-04-30 12:14:07 +00:00
Martin Diehl 2a5dfbbfa0 removed #define Marc and substituted #ifdef Marc by #ifdef Marc4DAMASK which is defined by either the submit scripts or DAMASK_marcXXXX.f90 2013-04-30 09:49:30 +00:00
Philip Eisenlohr e912c74da5 added possibility to select "mapping" of multiple grains.
'average', 'avg', 'mean' --> result averaged over Ngrains.
'parallel', 'sum' --> result added up from Ngrains.
2013-04-30 09:40:06 +00:00
Franz Roters 0d95d33646 use disk shape for twin also in initialization 2013-04-30 07:34:05 +00:00
Martin Diehl 533fcec96f named some loops and made BUTCHER TABLEAUs parameters in crystallite_integrateStateRKCK45 2013-04-29 11:17:30 +00:00
Christoph Kords cc72c0b17b added two new scripts:
"marc_extractData" is a simple variant of the "postResult" script that extracts all data from a macro *.t16 file (quadratic elements not yet supported) and writes it into an asciitable
"addDataToGeometry" searches for corresponding *.txt and *.vtk file in a given directory and adds the data from the *.txt file as SCALARS to the *.vtk file

first script needed an additional function in mesh that returns the corresponding node of a given ip in a specific element
2013-04-28 20:56:17 +00:00
Martin Diehl 9fb6d8167b simplified logic for single run, now not checking for extra arguments but simply checking if start and end elem/IP are the same 2013-04-26 13:23:36 +00:00
Martin Diehl 60739abb0d corrected statistics reporting in case of cutback 2013-04-26 09:13:39 +00:00
Martin Diehl f18c6c5903 fixed bug causing no fluctuation in case of using with 3Dvisualize 2013-04-25 14:43:20 +00:00
Martin Diehl d05ebff400 introduce commented new output with stream (no markers at begin/end) and with grid/size instead of res/dim 2013-04-25 14:03:54 +00:00
Martin Diehl 6e1cbb427d parallel execution currently doesn't work, now disabled for explicit 2013-04-24 12:32:25 +00:00
Christoph Kords 650b7ef4ac added python script that produces vtk files with the (node based and ip based) deformed mesh from marc output file 2013-04-23 18:30:56 +00:00
Christoph Kords 51552b027d fixed bug: initialize computation mode to zero, used to be unset for lovl=4 and not(timinc < theDelta .and. theInc == inc), which unpredictably caused trouble with marc 2013-04-23 14:42:22 +00:00
Martin Diehl 5837b25e27 accidently commited constitutive.f90 wasn't working, fixed now 2013-04-23 07:37:14 +00:00
Martin Diehl d5e4a354e6 misplaced space corrected 2013-04-22 15:35:15 +00:00
Franz Roters b5b6bf56d4 reintroduced Calculation of forest projections for edge dislocations, got lost after rev.1908 2013-04-22 15:17:24 +00:00
Martin Diehl 8da7544978 marc doesn't need error prone setup_code any more, moved Marc includes from code/include to lib/includeMarc
DAMASK_marcXXXX.f90 now is not a copy any more but contains one definition and an include statement

setup_code of spectral solver is now replaced by the standard 1) configure 2) make 3) make install
2013-04-22 14:41:33 +00:00
Christoph Kords a7706090d2 neighborhood relations were screwed up for some elementtypes as of version 2318, might have caused segmentation faults with nonlocal model, now fixed 2013-04-22 13:35:35 +00:00
Philip Eisenlohr 49ef6d7193 spectral_build_neighborhood now "hidden" in spectral compile part. (should not interfere with Marc/Abaqus compile anymore) 2013-04-22 09:01:58 +00:00
Christoph Kords ca99274a5d fluxes now always periodic for spectral solver, no need to explicitly state this in the geom file by use of "periodic" keyword 2013-04-22 08:32:51 +00:00
Christoph Kords 60fa91be9a calculation of cellnode positions now available in python via damask.core
first call damask.core.mesh.mesh_init_postprocessing(meshfilename) to initialize all necessary mesh variables
then damask.core.mesh.mesh_build_cellnodes(nodes) calculates the cellnode positions for a given list of node positions
the meshfile that is needed for the init is created automatically by mesh_init in DAMASK

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
 ./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-21 18:48:59 +00:00
Philip Eisenlohr d203e8dece removed extra comma (might have caused issues?) 2013-04-19 16:05:13 +00:00
Philip Eisenlohr e8f238a413 added specialized ipNeighborhood calculation for spectral solver case 2013-04-19 16:04:20 +00:00
Christoph Kords 12760026a4 last version still did not compile: IO_write_jobFile was not used by mesh_init 2013-04-19 13:02:37 +00:00
Christoph Kords 55bc741d1d last version did not compile: fixed small typo in mesh_init 2013-04-19 12:59:45 +00:00
Martin Diehl 4bc5e6717b fixed small bug when calculating mesh_Ncells 2013-04-19 12:41:06 +00:00
Christoph Kords 103c770ee6 fixed bug in cell generation: "FE_Ncellnodes" was not correct for elem type 17 and 136
too much memory allocated for mesh_cell
2013-04-18 17:56:14 +00:00
Martin Diehl 6b46a9c338 introduced output of initial geometry (vtk file) to mesh 2013-04-18 16:40:49 +00:00
Martin Diehl 8b2d9d8155 further improvements on implementing the plain mode 2013-04-16 17:07:27 +00:00
Martin Diehl 2ba986be2a introduced new dummy module "libs" for inclusion of further library sources. contains kdtree2 and Lib_VTK_IO at and is enabled for the spectral solver only at the moment 2013-04-16 12:45:12 +00:00
Christoph Kords 2a6ac44496 Major restructuring of mesh.f90 due to introduction of "cells". A cell is the volume (area in 2D elements) that surrounds an ip. It is described by up to 8 cell nodes and is used to calculate ip volume, area, and face normals. It can also be used to generate a mesh for e.g. paraview that contains ip output as cell centered data.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
 ./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-15 08:13:20 +00:00
Christoph Kords eaccf88e12 fixed bug that was introduced in the element type mapping in version 2305 2013-04-10 13:22:02 +00:00
Christoph Kords 9b1fd579d4 corrected phase numbering in example material.config 2013-04-10 12:09:09 +00:00
Martin Diehl ddcc795461 allowed arbitrary grid (formerly know as resolution), i.e. any integer>1 in each dimension is possible.
removed square root correction for divergence (doesn't make sense)
fixed FFT debug of spectral solver
2013-04-10 10:19:16 +00:00
Christoph Kords 4338cd13bc - added element type 54 for macro (2d 8node 4ip reduced integration)
- automatically raising error for unknown element type inside mapping function
2013-04-10 09:38:40 +00:00
Christoph Kords 338e338c9a - fixed bug in volume calculation of 2d elements
- subroutine "mesh_build_ipVolumes" uses openmp parallelization
2013-04-10 08:24:53 +00:00
Christoph Kords d80d416c32 2d elements now have a volume corresponding to a thickness of 1; used to have zero volume
now also 2d elements can be used with nonlocal model
2013-04-09 18:07:30 +00:00
Christoph Kords 5b508b5ee4 no need to exclude any geometry type from usage of nonlocal model 2013-04-09 13:10:31 +00:00
Martin Diehl 924d943edc simplified interface of CPFEM_general (removed P and F, made cs and dcs/dE optional)
commented and cleaned up the marc interface.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
For marc simulations, run
 ./code/setup/setup_code.sh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2013-04-09 10:08:00 +00:00
Martin Diehl 7b5a202e8c enabled alternative (and soon standard) keywords grid (resolution) and size (dimension).
removed errors for odd resolution
2013-04-08 14:22:32 +00:00
Christoph Kords ad879ebcf9 forgot to delete all occurrences of obsolete "constitutive_nonlocal_ accumulatedShear" in last commit 2013-04-08 08:51:12 +00:00
Christoph Kords 10e50bf41d improved accuracy of "accumulated shear": now added as a state variable, which facilitates integration (dotstate=shearrate); needs absolute tolerance value for state integration in material.config 2013-04-04 13:37:14 +00:00
Christoph Kords 3c25723139 - make sure that formerly blocked density is not able to move to neighbor until it is remobilized
- added some explicit indices instead of ":"
2013-04-03 16:29:48 +00:00