Commit Graph

1264 Commits

Author SHA1 Message Date
Martin Diehl 81531097f1 changes on convergence tolerances of AL and Polarisation, switched back to immediate correction of stress bc but only when last two average stresses are close to each other (cosine decay) 2013-07-30 15:32:55 +00:00
Pratheek Shanthraj 03d8f14a98 introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms.
dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!).

not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 16:25:37 +00:00
Martin Diehl 216063fe36 added Polarisation scheme as solver 2013-07-24 13:06:16 +00:00
Franz Roters 4f7cd76b29 renamed texture_rotation into texture_transformation
changed keyword for material.config to axes in line with geom_fromAng.py
former keyword rotation is still recognized for compatibility
2013-07-24 11:09:39 +00:00
Franz Roters f896b4a01e updated Marc/Mentat integration once more to minimze interference with original instal
added support for Marc2013, which is now default in apply_DAMASK_modifications
2013-07-24 08:41:58 +00:00
Martin Diehl ebe8361af0 decoupled compatibility/equilibrium calculation from stress BC correction. needs more iteration for mildly contrasted materials, but there the basic scheme is better suited anyway. but now converges better for highly contrasted VEs 2013-07-23 17:42:15 +00:00
Christoph Kords 3749cd5564 changed "math_AxisAngleToR" to "math_axisAngleToR" to be consistent with spelling in math.f90 2013-07-22 09:43:35 +00:00
Martin Diehl bbb0803ce4 itmin was always +1 2013-07-17 19:58:48 +00:00
Martin Diehl 74a7840a59 wrong intent(in) statement fixed 2013-07-15 10:13:28 +00:00
Martin Diehl 8210b14e90 comments for doxygen 2013-07-12 06:57:15 +00:00
Martin Diehl 4407edd802 checked convergence reporting for AL and BasicPETSc, seems to be ok now. 2013-07-08 15:48:13 +00:00
Martin Diehl 27c64ba4ad corrected line breaking for /echo/ tag 2013-07-08 09:46:38 +00:00
Franz Roters 1db24d2a72 removed duploicate writes
changed last occurences of old mode syntax
one more change in the direction of a plain mode, i.e. when there is no collect cycle
2013-07-04 11:59:45 +00:00
Pratheek Shanthraj 8f146ad385 fixed bug in cosine filter 2013-07-02 04:38:18 +00:00
Philip Eisenlohr 8ccdfb27f3 reworked the hexagonal slip system order and adjusted interaction matrices to a logic that facilitates later system additions without altering existing structure. 2013-07-01 14:01:42 +00:00
Franz Roters 4a291dc372 implementation of twin nucleation criteria according to Davids PhD thesis
applies for fcc crystal structure only
2013-07-01 13:06:01 +00:00
Martin Diehl 6ec57db0c6 fixed writing out of initial configuration for spectral solver in case of ngrains != 1, small formating improvements in other files 2013-07-01 06:40:09 +00:00
Martin Diehl 40ace5c666 more documentation and parameters capitalization unified and using ipc in all modules (sometimes called gr)
nonlocal: only missing line continuation in string fixed
2013-07-01 06:10:42 +00:00
Philip Eisenlohr f408ea3381 corrected error msg typos 2013-06-30 00:17:58 +00:00
Martin Diehl 5f20c86905 small documentation and output fixes 2013-06-28 18:59:21 +00:00
Martin Diehl ef2c2af471 doxygen documentation for J2 2013-06-28 18:58:10 +00:00
Philip Eisenlohr fd8d85896a error msg in case of unknown keyword was not properly trimmed --> 64kB of spaces... 2013-06-27 16:41:00 +00:00
Philip Eisenlohr f9f76fa158 missed these three files in last commit!
[added possibility for multi-level inclusion of files in *.config and loadcase files.
include subfiles by stating

{path/to/include}
]
2013-06-26 19:19:41 +00:00
Philip Eisenlohr 98528f9a89 added possibility for multi-level inclusion of files in *.config and loadcase files.
include subfiles by stating

{path/to/include}
2013-06-26 19:19:00 +00:00
Martin Diehl 0f1bb4f555 corrected doxygen documentation tags 2013-06-25 06:08:33 +00:00
Harsha Phukan aea8ab144a changed all file reads consistently to 65536 bytes length 2013-06-24 13:33:30 +00:00
Franz Roters 9a1ad7e089 added third value for slip twin interaction in fcc dislotwin example 2013-06-24 12:53:56 +00:00
Pratheek Shanthraj 7b2e48df68 new 'gradient' filter based on greens function from gradient elasticity theory. work in progress 2013-06-21 17:26:05 +00:00
Philip Eisenlohr 81cdd2b1cc added second type prismatic slip system for hexagonal lattices as 6th slip family.
consequently, 4 new interaction types (21,22,23,24) available for slip/twin and twin/slip.

interaction of slip/slip got rearranged and now has 42 types...
list of former --> new:
6	7
7	8
8	9
9	10
10	11
11	13
12	14
13	15
14	16
15	18
16	19
17	20
18	22
19	23
20	25
21	28
22	29
23	30
24	31
25	33
26	34
27	35
28	37
29	38
30	40
2013-06-14 10:20:31 +00:00
Franz Roters a78f91311e added svn properties: line ending: LF and keywords: ID 2013-06-14 09:49:33 +00:00
Pratheek Shanthraj 1c25403b5b changed default values of AL spectral solver tolerances. should've been done along with commit 2254. 2013-06-12 10:06:04 +00:00
Christoph Kords daf89806b1 fixed bug in deltaState; used to lead to unpredictable values in accumulated shear, since deltaState was not defined at this state position 2013-06-11 21:03:23 +00:00
Martin Diehl 296a5040b9 fixed bug in initialization when reading in material.config.
Fortran does not short-circuit, e.g. (if a>0 .and. b(a) ==c) might cause an out-of-bounds error when a=-1 it'll try to access b(-1)!
2013-06-11 20:16:40 +00:00
Martin Diehl 6f7740a243 made doxygen working for MSC.Marc again, small changes on the todo-statements for doxygen, they don't work on single lines of codes but only on module/variables 2013-06-11 16:35:04 +00:00
Ding Wang ce9f92a884 polishing output 2013-06-11 10:27:38 +00:00
Martin Diehl 1b53bad9a4 removed doubled equal sign in compiler name definition 2013-06-11 09:06:25 +00:00
Martin Diehl a7b6093edc fixed bug in include_linux64 for Marc 2012, tests should run again after updating the copy in /msc/....
polished two other versions of this file
unified output format for /echo/ in material.f90
2013-06-11 07:28:08 +00:00
Christoph Kords 55735b9580 grainrotation output now uses corrected rotation math 2013-06-07 09:35:00 +00:00
Martin Diehl 5d7f7f1bc5 reverted accidently commited changes to Makefile.
one public statement in math.f90 still prevented compilation of core module, fixed that
2013-06-07 07:27:12 +00:00
Martin Diehl e6ca36c7d7 removed old DAMASK_marcXXXX.f90 source files.
If you just use one Marc version, simply select DAMASK_marc.f90, this matches for all releases. 
DAMAKS_marc.marc compile with Marc 2010 does not work with Marc 2012. In case you need more than one Marc version in parallel, run symLink_Code.py in /code/setup to generate symlinks for Marc version 2010,2011,2012
2013-06-06 16:04:40 +00:00
Martin Diehl 06b977ff61 added public statements to fftw functions and constants after whole math module is private now 2013-06-06 09:25:27 +00:00
Christoph Kords 622e2dcf15 Fixed bug in math rotations: passive and active rotations were not clearly distinguished and partly mixed up; yet, luckily, only resulted in wrong output of euler angles.
Now fixed with following convention:
Rotation conversions do not switch implicitly from active to passive or vice versa EXCEPT when converting to or from any "Euler" type (Euler angles or Euler axis/angles). Those functions expect a passive rotation as input if converting to Euler type and return a passive rotation if converting from Euler type.
2013-06-05 19:10:37 +00:00
Philip Eisenlohr 84c43741a6 reworked distribution of microstructure constituents. now each material point represents "as good as possible" the fractional content of constituents.
removed error for volume fractions not equalling 1.

implemented capability to rotate the texture given in material.config.
2013-05-29 17:23:49 +00:00
Martin Diehl a21dd816c7 init reporting of constitutive_*.f90 had less/to much spaces, renamed label to LABEL because it is a parameter.
removed debug output of geom_fromEuclideanDistance.py
2013-05-28 17:31:55 +00:00
Christoph Kords a6361738fb (again!) fixed bug that was introduced in rev 2085: messed up indices of stiffness matrix in calculation of Poisson's ratio 2013-05-26 14:31:31 +00:00
Christoph Kords ab997f19ac fixed bug in nonlocal state init that was introduced in rev 2438 and used to reset states from other constitutions 2013-05-24 13:43:44 +00:00
Christoph Kords c0539d2383 replaced all remaining occurrences of state indices 2013-05-24 11:48:34 +00:00
Christoph Kords e2d970ce57 combined some state indices to an array with a more generic name 2013-05-24 09:02:30 +00:00
Christoph Kords ad4a66cb4e use named state variable indices in "deltaState" 2013-05-23 21:46:15 +00:00
Christoph Kords 332b207270 use named state variable indices in "kinetics" and "microstructure" 2013-05-23 21:40:00 +00:00
Christoph Kords 5900076d72 rename state variable index for critical resolved shear stress from "iTau" to "iTauF" 2013-05-23 21:32:36 +00:00
Christoph Kords ad2a245cdf use named state variable indices in "microstructure" 2013-05-23 21:28:45 +00:00
Christoph Kords c02645c996 use named state variable indices in "aTolState" 2013-05-23 21:13:56 +00:00
Christoph Kords f04a531e9b use named state variable indices in "stateInit" 2013-05-23 21:10:31 +00:00
Christoph Kords 65524b7f84 corrected named indices for state variables; however, not yet used 2013-05-23 20:30:06 +00:00
Christoph Kords e249366ded introduced named indices for state variables; however, not yet used 2013-05-23 20:15:23 +00:00
Christoph Kords 72129cae00 removed prefix "constitutive_nonlocal_" from all private variables 2013-05-23 19:56:36 +00:00
Christoph Kords b559666436 replaced FE_maxNipNeighbors with mesh_maxNipNeighbors, which is sufficient for memory allocation, since this gives the maximum of neighbors in all cp elements 2013-05-23 17:46:21 +00:00
Christoph Kords 77e9b80a16 bug fix for calculation of fluxdensity and slipnormal/direction in constitutive_nonlocal_postResults: indices for Fe were missing 2013-05-23 14:52:57 +00:00
Christoph Kords 41faa60323 fixed bug in backstress calculation for elements with periodic neighbors 2013-05-23 12:36:48 +00:00
Christoph Kords 3aaf60cffe polishing 2013-05-23 12:25:56 +00:00
Christoph Kords ec2503253e added output of slip direction and slip normal
changed name of fluxdensity outputs
2013-05-23 08:19:36 +00:00
Pratheek Shanthraj 399a0218c7 bug fix on commit 2424 that caused compilation errors 2013-05-22 09:40:23 +00:00
Christoph Kords e758ee32be improved accuracy of accumulated shear output: take shear rate from LpAndItsTangent for time integration 2013-05-21 10:04:52 +00:00
Nader Zaafarani 882c3e155f Adding two titles for accumulated shear due to "slip" and "twin" 2013-05-19 20:34:40 +00:00
Nader Zaafarani a104fda3b9 Adding accumulated shear due to "slip" and "twin" as two outputs 2013-05-19 20:32:06 +00:00
Christoph Kords 3759024b42 forgot to set neighboring_rho (also reverting changes made in rev 2384) 2013-05-19 19:23:31 +00:00
Nader Zaafarani 871df5427d Correction for shear_rate_twin in postresults: characteristic twin shear was missing 2013-05-19 15:45:34 +00:00
Martin Diehl c7ba8a2a9b introduced plain mode (no ping pong) again and added test for MSC.MArc 2012 as prove that its working 2013-05-17 17:52:46 +00:00
Martin Diehl 4416efe89b doxygen documentation for constitutive_none 2013-05-17 17:09:42 +00:00
Martin Diehl 7833f79f75 replaced numerics_unitlength with mesh_unitlength 2013-05-17 16:52:19 +00:00
Christoph Kords 0504fb3b21 reverted changes made in rev2384 2013-05-17 12:54:47 +00:00
Ding Wang 57b591ff39 corrected PETSc debug and Capitalization 2013-05-16 08:55:10 +00:00
David Mercier 190a59b707 set Lib_VTK_IO-master.zip and related files back to initial version, new version is not working with MSC.Marc 2010 2013-05-15 16:43:18 +00:00
Christoph Kords 7ae786e802 fixed bug in mesh_spectral_build_ipneighborhood that was introduced in rev2390: loop counter were mixed up 2013-05-14 09:45:22 +00:00
Martin Diehl 364267afa2 fixed bug when having . at end of directory name, added check to see if the specified working directory exists 2013-05-13 14:10:48 +00:00
Martin Diehl 872f6a9d90 introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading.
Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 09:44:23 +00:00
Christoph Kords 6fac7c347c fixed bug in constitutive_nonlocal_kinetics that was introduced in rev2088 along with the non-Schmid behavior: nonSchmid matrix of last slip system was used for all slip systems; led to extremely bad convergence due to flawed dLp_dT 2013-05-10 22:29:12 +00:00
Martin Diehl b83b557af2 updated Lib_VTK_IO again. Contains small bug fix ans support for Base64 encoding 2013-05-10 12:57:00 +00:00
Christoph Kords d9cc42daef slight change in flux formulation: no flux to neighbor, not only if there is a sign change in the velocity, but also if the neighboring velocity is exactly zero 2013-05-10 08:42:42 +00:00
Martin Diehl cfac73043d some lines in crystallite.f90 exceeded limit of 132 characters
fixed bug in Lib_VTK_IO.f90 (informed author as well)
added zip file containing source of Lib_VTK_IO.f90
2013-05-08 17:48:00 +00:00
Martin Diehl 755e0e2440 vtk files are now correctly written to current working directory.
updated Lib_VTK_IO.f90
2013-05-08 16:10:21 +00:00
Martin Diehl 85d4a37d95 moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize) 2013-05-08 15:52:29 +00:00
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