Commit Graph

1318 Commits

Author SHA1 Message Date
Martin Diehl 82faf74363 added doxygen documentation and unified variable names and some common parts of the code 2013-10-08 16:27:26 +00:00
Philip Eisenlohr 0a7052a7da fixed two errors in hexagonal slipslip and twintwin interaction matrix 2013-10-08 13:52:19 +00:00
Philip Eisenlohr 42cd4eb38f fixed some e12.6 string formats to avoid compile time warnings 2013-10-08 13:51:36 +00:00
Philip Eisenlohr a2169e338f fixed type-casting problems (mainly with Lib_VTK_IO) 2013-09-26 17:21:46 +00:00
Martin Diehl ca2a8d0a03 fixed typos in Abaqus related functions 2013-09-23 09:01:46 +00:00
Martin Diehl 595c8860a2 renamed myInstance -> matID, myStructure -> structID to be consistent with other constitutive models 2013-09-20 17:03:11 +00:00
Martin Diehl ad32aa3ced fixed calculation of CurlRMS when grid(1)==1 2013-09-20 16:17:25 +00:00
Martin Diehl a0f304021f fixed restarting for Polarisation and AL, updated corresponding test and polished output in python (some statements are screwed up when making them Python3 compatible) 2013-09-20 14:22:37 +00:00
Martin Diehl 40c399bf2e fixed bug in last commit and remove one more goto marker 2013-09-19 15:28:55 +00:00
David Mercier ce926aaf1f Minor contribution :
- add SVN keyword to the material.config file for cp-Ti
- add comments in lattice file
2013-09-19 09:17:09 +00:00
Martin Diehl ca0b3b7b81 changed name of some variables to be in accordance with constitutive_none,j2,phenopowerlaw 2013-09-19 07:46:01 +00:00
Martin Diehl 70c888f6ea fixed a typo in documentation 2013-09-18 19:54:39 +00:00
Martin Diehl 56d21286a6 fixed comment on values changed automatically 2013-09-18 14:33:40 +00:00
Martin Diehl 3eec0ecdcd introduced #EOF# as module wide parameter in IO.f90
renamed read/write binary file to read/write real/int
removed suffix job from readFile functions as the name of the model is needed as an argument
2013-09-18 14:07:55 +00:00
Martin Diehl 19a353edb4 made error in case of invalid precision more verbose 2013-09-18 13:59:42 +00:00
David Mercier a8eb6a985e - Addition of an example of a material.config file for hcp material (cp-Ti).
- New Error Message #214 when elastic constants are not correctly defined.
- Error message when the c/a ratio for hcp is higher than 2 (no physical sense).
2013-09-17 12:20:50 +00:00
Martin Diehl cb16b30409 set output order of spectral parameters to more reasonable one 2013-09-17 09:19:16 +00:00
Martin Diehl 3bf0ed84c8 some comments and minor improvements.
removed tol_gravityNodePos from prec because it's not used
tol_math_check is now used by spectral driver to check rotation
2013-09-14 10:59:35 +00:00
Christoph Kords 5f973a21c0 require all hardening coefficients/ non schmid coefficients to be specified in material.config, do not silently ignore any
warning message now more meaningful (?) when using less than maximum number of slip families
2013-09-12 14:47:09 +00:00
Martin Diehl ddfc40afe6 improve restart test, deletion of results was messed up. fixed restart for BasicPETSc, AL and Polarization currently don't run when restarting 2013-09-05 12:07:58 +00:00
Christoph Kords 4ebf46d2ac nonSchmidProjection was reset to zero each time during instance loop in init, so entries from previous instances got lost; fixed that 2013-08-23 13:53:03 +00:00
Christoph Kords 1f77594ca0 dislocation velocity was not correctly stored to state in last version 2013-08-21 13:54:36 +00:00
Christoph Kords 4f8664baa3 in kinetics: non schmid stresses only influence peierls mechanism, but not solid solution hardening; as a result the derivative of the velocity with respect to the resolved stress has to be split into a Schmid and a non-Schmid part 2013-08-21 12:21:52 +00:00
Christoph Kords ec377a6e8e removed obsolete "dead zone" scaling 2013-08-21 08:39:43 +00:00
Christoph Kords f706ba3ff9 rearranged arguments of "constitutive_nonlocal_kinetics", got "tauTreshold" as argument, not complete state
don't call "constitutive_nonlocal_kinetics" twice for edges in case of nonSchmid behavior, but just call once and copy results from positive to negative edges
2013-08-21 07:55:34 +00:00
Christoph Kords 277972edbd "dv_dtau" argument in "constitutive_nonlocal_kinetics" not optional anymore
renamed "screwStressProjection" to "nonSchmidProjection"
2013-08-21 07:23:02 +00:00
Christoph Kords 39c05e0d9e for bcc: exchanged 2nd and third hardening coefficient to be consistent with fcc; now 2nd is coplanar and 3rd is collinear; added newer reference concerning the bcc hardening coefficients 2013-08-14 11:47:03 +00:00
Christoph Kords e909b825cc correction of hardening coefficients also enabled for bcc
for fcc: now applied to all hardening coefficients, not only hirth, lomer, and glissile
2013-08-14 11:44:30 +00:00
Martin Diehl b8d5f3124f added citation on phenopowerlaw aluminum parameters 2013-08-09 21:24:42 +00:00
Pratheek Shanthraj e62b760a6e added F_avg = F_aim in boundary condition convergence check 2013-08-09 16:25:13 +00:00
Martin Diehl 822fa3d97c forgotten line continuation 2013-08-09 04:16:16 +00:00
Martin Diehl c93373a4ac removed long lines 2013-08-08 19:33:21 +00:00
Martin Diehl 3bb5a5e7b5 line too long in DAMASK_spectral_solverPolarisation.f90, reordered curl calculation for clearer look, fixed test for Marc 2013 2013-08-08 18:10:58 +00:00
Franz Roters 591a3023e8 made plain mode work with Abaqus
ATTENTION plain mode (Marc and Avaqus) does not support writung out the deformed geometries
2013-08-08 15:00:22 +00:00
Christoph Kords 3fcc00910c fixed bug that occured under marc (and abaqus?) with openmp:
"lastLovl" (for marc) and "lastStep" (for abaqus) have to be global module variables in FEsolving.f90 instead of local variables in the respective DAMASK_* files; otherwise they are initialized with zero each time hypela2/umat is called.
2013-08-08 12:53:03 +00:00
Martin Diehl 7f3c24e1f1 fixed some buggy last minute changes in the last commit 2013-08-08 09:13:29 +00:00
Pratheek Shanthraj 23e8fcbc51 updated config file with new convergence parameters 2013-08-08 09:12:12 +00:00
Martin Diehl ac92b90e0b fixed al and pol solver, now checking for div(p) = curl(f) = 0 2013-08-07 17:20:05 +00:00
Christoph Kords f2d6703148 multiplication now sensitive of crystal structure, yet multiplication in bcc identical to fcc except for the contribution of edges which is completely dropped for bcc 2013-08-05 09:26:37 +00:00
Christoph Kords 813ed130b0 nonSchmid behavior now ready to use for bcc with phenopowerlaw and nonlocal constitutive model.
nonSchmid tensors according to Koester,Ma,Hartmaier,2012.
extended lattice_Sslip with non-Schmid tensors (the full non-symmetric tensors are required for the tangent and cannot be retrieved from the symmetrized Mandel notation)
2013-08-05 09:23:21 +00:00
Christoph Kords 9f0fe873ba plain mode (no ping pong) now seems to work for marc, further testing needed
regular sequence of computation modes is
inc 0: 1,1,1,…
inc 1: 0,0,0,…
          7,1,1,…
          0,0,0,…
          1,1,1,…
and so on

after a cutback the computation modes follow as
inc 5: 8,0,0,…
          1,1,1,…
          0,0,0,…
          1,1,1,…
2013-08-02 16:19:45 +00:00
Christoph Kords fa3be9a967 corrected lastLovl introduced in rev2593 2013-08-02 15:44:28 +00:00
Martin Diehl 5b80ef3a4c fixed wrong error calculation in Polarisation and AL solvers, and strictened tolerances for tests 2013-08-02 13:55:44 +00:00
Christoph Kords 8d6b840802 removed variable lastMode from FE_solving, used to be used for detection of first call after ping pong; now this is done by checking for a change in the lovl (macro) or a change in the calculation step (abaqus) 2013-08-02 13:28:50 +00:00
Martin Diehl 74791a6686 corrected invalid error for mismatch_p 2013-08-02 11:48:24 +00:00
Christoph Kords 3d687be8b0 Changed order of blocks associated with different tasks:
backup jacobian
restore jacobian
age results
collect 
calc
This becomes important if multiple tasks are requested by one call to CPFEM_general. In plain mode this could, e.g., be 
1) backup jaco, 2) age results, 3) calc
2013-08-02 11:36:51 +00:00
Christoph Kords 7d2206356e parallelExecution flag now passed as input variable by CPFEM_general; flag is set on the solver level (DAMASK_abaqus_exp.f, DAMASK_marc.f90, etc.) 2013-08-02 11:20:11 +00:00
Franz Roters 4f9dbfa193 last commit also solves problem with infinite loop in case of a cutback in homogenization
corrected indentation of a comment
2013-08-02 07:59:55 +00:00
Franz Roters a08d263865 remnoved some unused fuctions from use statements 2013-08-02 06:18:41 +00:00
Christoph Kords a7f0f4ed79 moved some small part of the code that remembers maximum and minimum stress/stiffness values to different location 2013-08-01 16:27:37 +00:00
Christoph Kords 728facd451 corrected indentation level of parts of the code
moved debugging output of stress and stiffness to different position in code
2013-08-01 16:10:56 +00:00
Martin Diehl ac2ca43cfc added svn properties 2013-08-01 09:13:46 +00:00
Martin Diehl 2dc985f558 set tolerances to useful values and renamed them.
allowed "keyword = value" in numerics.config (simply removing "=" from line)
2013-07-31 15:03:38 +00:00
Martin Diehl b5a11f9d31 fixed bug in AL, weakened correction slightly 2013-07-30 21:04:41 +00:00
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