Commit Graph

2187 Commits

Author SHA1 Message Date
Pratheek Shanthraj b8a45f660e reinstating changes lost in previous commit 2014-10-10 16:45:14 +00:00
Luv Sharma 9b1e55f0e4 added helper functions to communicate accumulated slip between damage and plasticity modules 2014-10-10 16:34:51 +00:00
Pratheek Shanthraj e943a3a8b0 fixed typo 2014-10-10 16:23:19 +00:00
Pratheek Shanthraj 96f036e34e updated MPI reporting in line with recent changes 2014-10-10 16:21:10 +00:00
Pratheek Shanthraj 76f5495038 fixed bug for when petsc is not defined 2014-10-10 16:02:42 +00:00
Pratheek Shanthraj eacb048c83 fixed compilation error 2014-10-10 15:58:18 +00:00
Martin Diehl d095c2484d simplified multi processor (MPI) reporting 2014-10-10 13:08:34 +00:00
Luv Sharma 1657e0f7ba started implementing a void fraction based ductile damage model. 2014-10-10 12:42:12 +00:00
Pratheek Shanthraj ba65d04bb3 removed redundant crystallite_temperature/heat and materialpoint_heat and using temperature from the thermal module 2014-10-10 12:28:57 +00:00
Pratheek Shanthraj d10f8d3af6 move heatGen to adiabatic 2014-10-10 08:58:25 +00:00
Pratheek Shanthraj 885a71c0b7 cleaned up compilation for FEM 2014-10-10 08:40:59 +00:00
Pratheek Shanthraj c8929b47d3 only output from the root processor for parallel runs 2014-10-09 20:23:06 +00:00
Pratheek Shanthraj 45bedec229 damage diffusion tensor should be in reference configuration 2014-10-09 15:56:15 +00:00
Luv Sharma 2eadb0a178 changed thermal heatGen to thermal adiabatic 2014-10-09 14:08:32 +00:00
Luv Sharma fc57c6e572 added thermal adiabatic file 2014-10-09 14:06:45 +00:00
Pratheek Shanthraj c1a5ce047f disable openMP for FEM compilation 2014-10-06 20:45:04 +00:00
Pratheek Shanthraj 2712ce0254 polishing 2014-10-06 17:01:39 +00:00
Martin Diehl 020e697b74 small polishing 2014-10-01 20:36:00 +00:00
Pratheek Shanthraj 4d7a9ce030 bug fix 2014-10-01 16:35:56 +00:00
Luv Sharma ee60ea8d5b some improvements in ductile model 2014-10-01 15:18:33 +00:00
Martin Diehl baa08d8155 made analytic tangent default for spectral solver
removed possibility to compile without PETSc (a lot of effort for little/no use)
2014-10-01 12:29:12 +00:00
Su Leen Wong 0824335b10 TRIP contribution to Lp and dotState 2014-10-01 08:11:39 +00:00
Pratheek Shanthraj 3aa5ca68c5 marc/abaqus includes for ductile damage 2014-09-30 13:25:07 +00:00
Pratheek Shanthraj 3da023a1ba thermo-mechanical example 2014-09-29 12:40:57 +00:00
Pratheek Shanthraj c7a746ddc1 constitutive dot state now takes Lp which is needed to the thermal dot state.
thermal problems should be working now
2014-09-26 20:49:25 +00:00
Pratheek Shanthraj ce3be666e8 redoing previous changes that were lost in commit 3502
please please please make sure you merge your changes to an updated version of the code before committing them
2014-09-26 18:49:40 +00:00
Pratheek Shanthraj 590eb31ed0 some more renaming or thermal related routines 2014-09-26 18:07:48 +00:00
Pratheek Shanthraj c39aac9c9d makefile to go with renaming 2014-09-26 16:10:41 +00:00
Pratheek Shanthraj 40cd991496 thermal_none to thermal_isothermal 2014-09-26 16:09:25 +00:00
Pratheek Shanthraj dabfa6d2e6 renamed get/put thermal to get/put temperature and thermal_none to thermal_isothermal 2014-09-26 16:07:26 +00:00
Pratheek Shanthraj d4076e4db9 some fixes to get ductile damage working with recent changes 2014-09-26 15:23:10 +00:00
Pratheek Shanthraj 5c3f4e3e69 some fixes to get homogenisation none working with recent changes 2014-09-26 15:16:10 +00:00
Luv Sharma d3fd4bdf90 started adding a ductile damage model 2014-09-26 10:38:13 +00:00
Pratheek Shanthraj 82932c3708 added field quantity output 2014-09-26 10:34:36 +00:00
Pratheek Shanthraj 280e563fa0 removed homogenization_sizepostresults. This information only needs to be stores per homog type in homogState and not for every element and ip 2014-09-26 10:32:56 +00:00
Pratheek Shanthraj a42fd0608a use constitutive_xxx_Noutput instead of phase_Noutput as number of outputs for more consistency with constitutive multi physics output 2014-09-26 10:25:26 +00:00
Pratheek Shanthraj 34c33cca13 changed Noutput to constitutive_nonlocal_Noutput for consistency 2014-09-26 10:21:00 +00:00
Su Leen Wong 1c38d42a0b Projection matrices for transformation 2014-09-24 12:31:16 +00:00
Luv Sharma c38b80ad49 updated module include list for abaqus interface also. 2014-09-24 11:11:17 +00:00
Luv Sharma a0b724836e typo fixed 2014-09-24 10:45:46 +00:00
Su Leen Wong 4d351ccba1 Outputs for TRIP volume fraction 2014-09-24 07:26:27 +00:00
Su Leen Wong c2ee6ff299 Renaming transformation input variables to more logical names 2014-09-23 15:38:41 +00:00
Pratheek Shanthraj 556d485ca4 typo: diffusion not difusion 2014-09-23 12:30:52 +00:00
Pratheek Shanthraj 5ce49387f6 added history variable for brittle damage to prevent healing. renamed/removed some lattice damage tensors and symmetrize them. 2014-09-23 12:22:34 +00:00
Pratheek Shanthraj d2dee58703 removed obsolete thermal and damage modules 2014-09-23 10:42:57 +00:00
Pratheek Shanthraj c2c5c0c521 subsumed constitutive_damage/thermal into constitutive 2014-09-23 10:38:20 +00:00
Pratheek Shanthraj c029554c31 better naming damage_local --> damage_brittle 2014-09-23 08:43:16 +00:00
Pratheek Shanthraj 854879d277 undid some marc/abaqus related changes 2014-09-23 08:19:15 +00:00
Pratheek Shanthraj bb7f9be6ea remove obsolete constitutive field modules 2014-09-22 20:57:12 +00:00
Pratheek Shanthraj cb4117df95 adiabatic part of possibly conductive thermal model 2014-09-22 20:55:54 +00:00
Pratheek Shanthraj 3948d99f4f local part of possibly nonlocal brittle damage model 2014-09-22 20:38:19 +00:00
Pratheek Shanthraj f6c0f7062e moved field initialization from homogenization to material since it is needed in crystallite 2014-09-22 20:36:55 +00:00
Pratheek Shanthraj 6ace12be0c removed unused variables 2014-09-22 20:34:42 +00:00
Luv Sharma e83a0fb3f7 polished field state related changes and brittle damage 2014-09-22 18:15:19 +00:00
Martin Diehl 2b35cc041e fixed wrong line continuation 2014-09-19 19:38:59 +00:00
Luv Sharma 649ec1fca8 Finished applying newstate to Homogenization. 2014-09-19 17:59:06 +00:00
Luv Sharma 4bbd6583a0 prepared for using new state structure in homogenization 2014-09-18 15:27:52 +00:00
Martin Diehl f64418f29a removed multithreaded FFTW, does not make a big time difference and will simplify linking when using PETSc 2014-09-18 15:23:11 +00:00
Chen Zhang 2d460c3b78 remove redundant debug output 2014-09-16 18:22:03 +00:00
Pratheek Shanthraj cc7064d025 changed default petsc_options for FEM and removed unused variables 2014-09-16 10:02:37 +00:00
Martin Diehl 27026c7f6a fixed missing initialization of output in RGC that got lost in Rev 3452 2014-09-16 09:50:37 +00:00
Su Leen Wong 74db2ef064 Small fix for missing constitutive dislotwin output 2014-09-16 08:02:19 +00:00
Martin Diehl 7cb7815176 polishing style 2014-09-13 10:04:44 +00:00
Luv Sharma 05d980f1e3 minor fixes 2014-09-12 15:58:03 +00:00
Su Leen Wong f8646daf85 Added projection matrix for strain-induced transformation 2014-09-12 09:10:19 +00:00
Martin Diehl 5b6f143fb0 intel fortran behaves strange for standard check:
https://software.intel.com/en-us/forums/topic/506419
worked around that
2014-09-11 16:02:05 +00:00
Pratheek Shanthraj d6080d332d updated to be compatible with latest version of petsc.
Please change your PETSC_DIR to /opt/petsc-dev and PETSC_ARCH to ifort (gfort or ifort-debug for gfortran and debugging version respectively on maws01 also available)
2014-09-11 13:28:15 +00:00
Martin Diehl 4aa2c65d1e next try to make the makefile work in all possible cases for FEM and spectral 2014-09-11 09:51:38 +00:00
Luv Sharma 7c60c45f80 corrected typos and some field related changes 2014-09-10 18:26:12 +00:00
Martin Diehl 0af667e94d set search path (VPATH) to ../private/FEM/code only in case FEM solver is the target 2014-09-10 16:27:53 +00:00
Luv Sharma 2f76ecf0e3 fixed error in use statement. 2014-09-10 15:12:14 +00:00
Luv Sharma 7b5cbf4199 now compiling thermal_adiabatic.f90 also and some field related changes 2014-09-10 15:05:28 +00:00
Luv Sharma e8f8655c98 fixed buggy initialization. Now initialising in similar way to constitutive models. 2014-09-10 14:14:03 +00:00
Su Leen Wong 8e6ea7d9c7 Added Ctrans to homogenizedC and strain-induced martensite part to basic states 2014-09-10 12:12:17 +00:00
Luv Sharma dc406a01c0 added funtions to get averged properties at integration points. 2014-09-10 08:37:12 +00:00
Martin Diehl f1e345dcdb fixed compilation for FEM 2014-09-09 15:59:41 +00:00
Franz Roters b6847ebb01 added use mesh statement needed for output of mesh_element 2014-09-09 15:18:49 +00:00
Luv Sharma 923adbc2d3 added auxillary functions: get/put 'Physics' to communicate regularised/ unregularised values between solver and constitutive physics 2014-09-05 16:31:27 +00:00
Luv Sharma a7741457b9 more work on field state. 2014-09-04 16:11:00 +00:00
Philip Eisenlohr 6612603a23 changed naming convention (keeping backward compatibility in numerics.config):
myspectralsolver --> spectral_solver (called "spectralsolver")
myfilter --> spectral_filter (called "spectralfilter")
2014-09-03 19:59:47 +00:00
Luv Sharma 79f572f869 more work on homogenisation new state and introduced field state 2014-09-03 17:12:06 +00:00
Philip Eisenlohr 2435c8d836 dropped "elasticState" (for now) 2014-09-02 20:18:37 +00:00
Philip Eisenlohr f05f47f6e3 optical polishing 2014-09-02 20:16:33 +00:00
Philip Eisenlohr e8a428613c general polishing; simplification of RKCK45 dotstate assembly for different stages. 2014-09-02 20:11:57 +00:00
Philip Eisenlohr 630d9efffd Heavy bug fixing: Butcher tableau multiplication in RK45 integration scheme, state updating of adaptive Euler, wrong assignment to evolution rate instead of state in fixed-point iterator.
Exchanged possibly non thread-safe matmul in RKCK45, simplified RK45 integration step 1--3 vs 4 handling.
(Thanks to Pratheek, Luv, and Chen for their help!)
2014-09-02 19:46:52 +00:00
Su Leen Wong 5f1e49c053 Cleaner way of calculating the Bain strain 2014-09-01 10:47:48 +00:00
Su Leen Wong 565e703f49 Minor fix for now 2014-09-01 09:30:33 +00:00
Pratheek Shanthraj fc941a6436 fixed some bugs in previous commit 2014-09-01 08:57:34 +00:00
Su Leen Wong e8a3b7ed19 Pitsch orientation relationship for phase transformation 2014-08-29 13:03:48 +00:00
Philip Eisenlohr 929ba1b04b reporting of elFE (instead of elCP) when advertised as such. 2014-08-27 15:54:50 +00:00
Philip Eisenlohr a360f14581 debug now reporting elFE in addition to elCP 2014-08-27 15:54:11 +00:00
Philip Eisenlohr 19bd2d1430 changed CPFEM_COLLECT to be 2^0=1 such that computation modes with collection step are always an odd number. (used to be 2^4 with no particular intrinsic meaning...) 2014-08-26 17:39:52 +00:00
Chen Zhang 5bc49e6441 change the debug level of cutback reporting to non-selective, of central solution from basic to extensive, of stiffness perturbation to selective.
fixed few typos and deleted excess whitespace.
2014-08-26 14:44:32 +00:00
Chen Zhang e8649ded56 change debug level for OUTDATED info from extensive to basic 2014-08-26 14:21:24 +00:00
Martin Diehl 06988353ee now compiling spectral solver even if private folder is missing 2014-08-25 06:20:25 +00:00
Martin Diehl 977392cfe7 compiler for spectral solver is now the version provided by petsc (similar to linking)
added option to compile FEM solver with same Makefile
2014-08-24 22:07:53 +00:00
Luv Sharma e555ce4827 started introducing new state structure in homogenisation 2014-08-21 17:48:20 +00:00
Luv Sharma 594c3c007b fixed wrong index 2014-08-21 12:02:09 +00:00
Pratheek Shanthraj 1e8139452c added functions for factorial, binomial and multinomial coefficients 2014-08-21 08:33:55 +00:00
Philip Eisenlohr 5c429a9939 added ..._ROOT/lib64 as possible locations for FFTW and HDF5 libraries 2014-08-20 19:13:28 +00:00
Luv Sharma ab5033adde improved strain splitting (tangent) 2014-08-19 11:27:58 +00:00
Pratheek Shanthraj 2d0c9e03c8 all constitutive models coupled with damage 2014-08-15 09:09:31 +00:00
Pratheek Shanthraj 7666013758 use subF instead of F for analytic jacobian calculation 2014-08-15 09:04:00 +00:00
Su Leen Wong 1da874a6db Modification of basic states for phase transformation 2014-08-14 13:18:33 +00:00
Su Leen Wong 591a82de26 Fixed issue with previous commit 2014-08-14 12:21:51 +00:00
Su Leen Wong 5719987f95 Calculation of eigendeformation tensor for phase transformation 2014-08-14 10:04:13 +00:00
Pratheek Shanthraj ad6467638f minor change to strain energy calculation 2014-08-14 09:59:16 +00:00
Pratheek Shanthraj 47707f86b6 fixed typo 2014-08-11 21:28:15 +00:00
Pratheek Shanthraj b36b38d7c3 introduced non-zero residual damage to prevent stiffness matrix from being ill-conditioned 2014-08-11 17:13:06 +00:00
Luv Sharma 33b4c0d908 minor typo 2014-08-11 13:33:17 +00:00
Pratheek Shanthraj 584beee94e dTdFe bug fix 2014-08-11 08:34:44 +00:00
Pratheek Shanthraj 8f444b0695 minor changes to multiphysicsOut 2014-08-10 11:14:43 +00:00
Pratheek Shanthraj 0a01d5f3a1 hardening does not depend on damage 2014-08-10 10:47:12 +00:00
Pratheek Shanthraj b6ccbc0fe9 reworked energy splitting for damage 2014-08-10 10:45:07 +00:00
Pratheek Shanthraj ec71d77038 corrections to analytic dPdF calculation 2014-08-10 10:27:35 +00:00
Pratheek Shanthraj cfab313e3f corrected incorrectly set dLp/dS for plasticity none 2014-08-10 10:06:03 +00:00
Pratheek Shanthraj 3700e132b0 cleaned up strain energy splitting. coupled damage to plasticity for J2 and phenopowerlaw 2014-08-08 22:29:38 +00:00
Pratheek Shanthraj 9c40f187a0 corrections to plastic work calculation and volumetric strain split 2014-08-08 17:01:20 +00:00
Luv Sharma c4b9ddf45d added strain decomposition to uncouple damage from negative hydrostatic part 2014-08-08 15:37:32 +00:00
Pratheek Shanthraj 10d2cfebd3 changed naming scheme for helper functions. damage and thermal dependent stress calculation more modularised 2014-08-08 14:02:20 +00:00
Pratheek Shanthraj da5f3d8f85 removed some unused variables 2014-08-08 13:54:08 +00:00
Martin Diehl f55dc4b639 added kmc example material.config
added elastic test for mkc, titanmod, and J2
2014-08-08 13:51:11 +00:00
Martin Diehl 0ecdd2692b started to introduce state for phase transformation (adding one set of additional basic states) 2014-08-08 12:20:43 +00:00
Martin Diehl dafccad7f5 bufgix for last commit, Franz please rather check this 2014-08-08 12:16:51 +00:00
Martin Diehl 1a0e0f4fa2 fixed wrong? homogenizedC calculation 2014-08-08 11:13:04 +00:00
Martin Diehl ea66511160 added dislokmc (dislotwin modification) 2014-08-08 11:04:40 +00:00
Pratheek Shanthraj 9d615b3d3b added getDamage and getTemperature helper functions 2014-08-08 07:39:17 +00:00
Pratheek Shanthraj 2f5d9376dc call IO_error if inversion error in analytic jacobian calculation 2014-08-07 21:16:37 +00:00
Pratheek Shanthraj 98a297cd9e removed old inaccurate analytic jacobian with exact jacobian. seems to work as well as perturbed jacobian and better than the old analytic jacobian for some simple tests. 2014-08-07 21:08:34 +00:00
Martin Diehl 0fc5e3717c finished updating to new state of dislokmc 2014-08-07 15:50:52 +00:00
Martin Diehl f0ade8d0a3 started to add David Cerecedas Model 2014-08-07 14:59:28 +00:00
Martin Diehl fcfe26f7d3 added lineaODF example file 2014-08-07 09:29:50 +00:00
Pratheek Shanthraj 8a683e2371 reactivated plastic work contribution to damage. default state integration tolerance is 0.001. use aTol_damage tag in material.config to specify your own tolerance 2014-08-06 11:37:47 +00:00
Martin Diehl fa45ef456e fixed commercial FEM interfaces to work with damage_local 2014-08-04 18:18:28 +00:00
Martin Diehl a787d66763 updated restart test, deviations seems to be related to tolerance in stress BC 2014-08-04 17:50:01 +00:00
Pratheek Shanthraj 6371d0d2bc slight modifications to gradient damage 2014-08-01 15:15:02 +00:00
Pratheek Shanthraj 3d83ac64fd added code for local damage 2014-08-01 14:54:57 +00:00
Martin Diehl 2e377056ed missing adjustment of interface for call to CPFEM_general related to last commit 2014-07-24 12:26:01 +00:00
Martin Diehl e7962e6dd1 disabled calculation and storage of Cauchy stress and derived quantities for spectral solver and own FEM:
command: /usr/bin/time -v
OLD DAMASK terminated on:
Date:               24/07/2014
Time:               12:03:58
STOP 0
	Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom"
	User time (seconds): 12230.57
	System time (seconds): 45.98
	Percent of CPU this job got: 353%
	Elapsed (wall clock) time (h:mm:ss or m:ss): 57:49.24
	Average shared text size (kbytes): 0
	Average unshared data size (kbytes): 0
	Average stack size (kbytes): 0
	Average total size (kbytes): 0
	Maximum resident set size (kbytes): 243392
	Average resident set size (kbytes): 0
	Major (requiring I/O) page faults: 0
	Minor (reclaiming a frame) page faults: 13839
	Voluntary context switches: 452541
	Involuntary context switches: 2233168
	Swaps: 0
	File system inputs: 0
	File system outputs: 556112
	Socket messages sent: 0
	Socket messages received: 0
	Signals delivered: 0
	Page size (bytes): 4096
	Exit status: 0

NEW DAMASK terminated on:
Date:               24/07/2014
Time:               14:12:16
STOP 0
	Command being timed: "DAMASK_spectral -l tensionX.load -g 20grains16x16x16.geom"
	User time (seconds): 11948.80
	System time (seconds): 37.28
	Percent of CPU this job got: 365%
	Elapsed (wall clock) time (h:mm:ss or m:ss): 54:42.97
	Average shared text size (kbytes): 0
	Average unshared data size (kbytes): 0
	Average stack size (kbytes): 0
	Average total size (kbytes): 0
	Maximum resident set size (kbytes): 235280
	Average resident set size (kbytes): 0
	Major (requiring I/O) page faults: 0
	Minor (reclaiming a frame) page faults: 14200
	Voluntary context switches: 362509
	Involuntary context switches: 1808383
	Swaps: 0
	File system inputs: 0
	File system outputs: 556192
	Socket messages sent: 0
	Socket messages received: 0
	Signals delivered: 0
	Page size (bytes): 4096
	Exit status: 0
2014-07-24 12:19:15 +00:00
Martin Diehl 7afc21f48f fixed reading in of state for restarting 2014-07-24 08:38:52 +00:00
Martin Diehl a8c0ca66fc corrected linking: use PETSC MPI (when available) rather then system wide version.
little polishing (e.g. debug output for new state)
2014-07-23 13:26:05 +00:00
Martin Diehl ab9f36cfe7 fixed bug in J2 related to new state.
make clean now deletes old marc compilations
2014-07-23 08:23:56 +00:00
Martin Diehl 4e3a202f02 set local option as default for plastic state 2014-07-22 17:38:49 +00:00
Su Leen Wong 1707f7d367 updated to read in TRIP parameters 2014-07-22 07:43:03 +00:00
Martin Diehl 0252fea3d7 fixed small bug in newstate 2014-07-10 08:47:00 +00:00
Martin Diehl 146d151357 disabled multiphysics output because homogenization output was not accessible anymore 2014-07-09 13:34:33 +00:00
Luv Sharma e16eb379ad some more consistency 2014-07-09 12:47:42 +00:00
Luv Sharma a1162e6fb4 removed some redundant parts 2014-07-08 14:58:23 +00:00
Pratheek Shanthraj 4af3459424 fixed bug in symmetrizing 2nd order tensor 2014-07-08 01:02:29 +00:00
Pratheek Shanthraj a68b216773 reworked damage. cleaned up homogenisation. removed VI stuff and using projection instead 2014-07-07 17:50:40 +00:00
Luv Sharma 7dfda9d8af remved bug in phenopower law state initialization and a minor change in Makefile 2014-07-07 14:21:58 +00:00
Pratheek Shanthraj cd1d286611 updated makefile for FEM solver and minor changes to spectral makefile 2014-07-04 10:46:23 +00:00
Pratheek Shanthraj e71241a627 reverted change in commit 3263. correct index was being used in use of mapping. changed to material_phase for clarity 2014-07-04 09:46:43 +00:00
Luv Sharma a1468b2693 for consistency, pushed assignment of sizePostResults from aggregators to individual thermal/damage label. Also fixed wrong index while use of mapping 2014-07-04 09:22:57 +00:00
Pratheek Shanthraj 9cc20ad751 fixed bug in counting outputs and changed damage formulation 2014-07-04 08:26:59 +00:00
Pratheek Shanthraj 546984b075 fixed bug in phase instance initialisation 2014-07-04 08:26:25 +00:00
Luv Sharma a2cd7815e9 corrected none variable in pheno case 2014-07-03 14:02:20 +00:00
Pratheek Shanthraj 3ee743ce66 corrected dependencies for thermal and damage modules in makefile 2014-07-03 13:32:21 +00:00
Martin Diehl 27d861decc fixed calculation of post results size introduced with new state 2014-07-03 13:17:29 +00:00
Luv Sharma d5952138e3 corrected wrong dimensions of aTolstate and inclusion order of files in abaqus and marc interfaces 2014-07-03 12:52:33 +00:00
Luv Sharma 499eb7cdbe included damage and thermal files in Abaqus and Marc intefaces 2014-07-03 10:59:27 +00:00
Luv Sharma 4812939b21 fixed error concerning state indices check 2014-07-02 14:54:43 +00:00
Martin Diehl 8fa2dcffbd changed to new state, please report bugs to Luv or Martin 2014-07-02 12:27:39 +00:00
Pratheek Shanthraj ad452508b0 cleaned up some errors with post results size 2014-07-01 23:02:24 +00:00
Pratheek Shanthraj ac36190f10 no need to allocate a large array for constitutive(thermal,damage)_sizePostResults (changes apply only to new state) 2014-06-30 14:47:30 +00:00
Martin Diehl 9be7c0d5f3 fixed floating point exception in nonlocal (calculation of dUpper at several places) 2014-06-26 20:01:38 +00:00
Luv Sharma 4c7d6c4464 fixed error with Noutput counter and few other newstate related changes (nonlocal) 2014-06-26 13:53:12 +00:00
Pratheek Shanthraj 57a865c886 changed pointer attribute to allocatable in p_vec and p_intvec 2014-06-25 16:08:49 +00:00
Pratheek Shanthraj 03d1bb6b8b changed thermal and damage defaults from undefined to none 2014-06-25 16:06:21 +00:00
Pratheek Shanthraj baaf104238 changed pointer attribute to allocatable in tState. since the pointers are only used for allocation, allocatables should do the same job more efficiently 2014-06-25 14:22:08 +00:00
Pratheek Shanthraj 8b11b1b6b0 added thermal and damage effects to elastic hooke's law 2014-06-24 23:22:52 +00:00
Pratheek Shanthraj 676cba13a8 modified state integrators to perform staggered integration of plasticity damage and thermal internal variables 2014-06-24 23:21:25 +00:00
Pratheek Shanthraj 6dc7f1bac9 thermal and damage modules now synched with DAMASK 2014-06-24 23:19:21 +00:00
Pratheek Shanthraj ab9a9f4f58 thermal and damage modules now linked to DAMASK 2014-06-24 23:18:39 +00:00
Pratheek Shanthraj d594996e7c added mesh dependent characteristic length as a numerical parameter needed for gradient and phase field problems on unstructured meshes 2014-06-24 23:18:07 +00:00
Pratheek Shanthraj 5d88a78206 added thermal and damage modules as examples of multi physics modules. only works with new state layout and still under testing.
damage_none: does nothing
damage_gradient: interacts with solver to solve gradient damage problems

thermal_none: does nothing
thermai_adiabatic: local heating only
thermal_conduction: interacts with conduction solver to solve coupled heat transfer problems
2014-06-24 22:59:16 +00:00
Pratheek Shanthraj 0c66204904 reference temperature for thermal problems parsed and stored in lattice 2014-06-24 22:54:19 +00:00
Pratheek Shanthraj 81729dab2a changed output parsing slightly to not give an error for undefined (for plasticity) parameters or outputs. get Noutput for plasticity from constitutive_xxx_noutput and not from phase_Noutput which is the total number of outputs (plasticity+thermal....) 2014-06-24 22:53:25 +00:00
Pratheek Shanthraj 18b1245874 added damage and thermal state vectors and phase instances. only works with new state layout. 2014-06-24 13:51:17 +00:00
Pratheek Shanthraj 9cdee47500 now parsing thermal conductivity, thermal expansion coeff and damage tensor from material config if present and symmetrizing according to lattice structure. setting up for multi physics 2014-06-24 13:48:29 +00:00
Pratheek Shanthraj c1fab338c3 new state with constitutive none was giving segfaults. fixed. 2014-06-24 13:44:50 +00:00
Luv Sharma 6ad6c7ad82 removed few bugs from nonlocal (newstate) 2014-06-24 09:24:59 +00:00
Luv Sharma 9ef7d2261e removed few bugs from nonlocal (newstate) 2014-06-24 09:24:19 +00:00
Pratheek Shanthraj 6347adcb65 fixed error in previous commit 2014-06-23 22:09:15 +00:00
Pratheek Shanthraj d3356bde4e moved mappings for new state layout to material 2014-06-22 18:58:29 +00:00
Martin Diehl 257bf2f9e4 renamed to fit naming scheme 2014-06-18 14:36:29 +00:00
Pratheek Shanthraj bf24bd83f2 fixed error check for CoverA ratio 2014-06-18 13:26:16 +00:00
Martin Diehl eaa3e32f21 improved warning and error messages 2014-06-18 09:10:16 +00:00
Martin Diehl 2160b37b1e fixed warning message in constitutive 2014-06-17 15:24:44 +00:00
Martin Diehl 3533138936 removed delta state for J2 and phenopowerlaw when using new state 2014-06-17 06:54:49 +00:00
Martin Diehl 5bd80123b4 fixed two severe bugs in new state, now a little bit faster (399 vs. 378 sec) and more memory efficicient (251312 vs. 236832 kbytes) for Phenopowerlaw 16x16x16 example (short load only) 2014-06-16 19:19:38 +00:00
Luv Sharma 6da07dd437 new state nonlocal output fixed 2014-06-16 14:28:11 +00:00
Martin Diehl 488927a231 worked on phenopowerlaw test, made material.config writer case insensitive, polished lattice, removed deprecated part of documentation 2014-06-16 13:11:26 +00:00
Luv Sharma e8be3b871b commented out the logical 'nonlocal' from newState structure(might be used later) 2014-06-14 11:50:46 +00:00
Luv Sharma 5755ed547b applied new State for nonlocal model and a trivial changes in crystallite 2014-06-13 20:53:17 +00:00
Martin Diehl da8f10fe6f consistency checks for lattice type improved 2014-06-13 08:57:00 +00:00
Martin Diehl bc5456b6c7 added dummy homogenization 2014-06-12 13:58:54 +00:00
Luv Sharma 0ca7f01186 fixed previous incomplete commit and few other trivial changes in newstate 2014-06-11 16:52:18 +00:00
Luv Sharma 8130f4b21f changes to make newstate aware of damage variables 2014-06-11 16:32:09 +00:00
Luv Sharma 401b31c951 new state related changes 2014-06-11 12:27:41 +00:00
Luv Sharma a54bb6ab24 prepared new State for dislotwin 2014-06-11 12:11:14 +00:00
Martin Diehl f5ca6b5b36 made (homogeneous) temperature working for spectral solver 2014-06-11 08:19:07 +00:00
Martin Diehl 6327a0c471 removed debug statement 2014-06-10 18:00:13 +00:00
Pratheek Shanthraj 01a0a0d0bb reverted some changes from previous commit 2014-06-06 08:57:24 +00:00
Pratheek Shanthraj 639ca89133 DAMASK (except nonlocal) now sees and interacts with FEM solver 2014-06-06 00:38:29 +00:00
Martin Diehl 4d3b09030c fixed array out of bounds error due to misplaced #ifdef 2014-06-05 15:22:35 +00:00
Luv Sharma bc1750f6f3 prepared new state for dislotwin 2014-06-03 13:46:42 +00:00
Martin Diehl 164252213b more work on the new state 2014-05-27 14:46:03 +00:00
Martin Diehl f7e574d7ab fixed missing sizeState initialization causing problem with output (only for new state) 2014-05-27 12:10:16 +00:00
Martin Diehl 8f6200d01a small bug fix for last commit 2014-05-23 23:17:22 +00:00
Martin Diehl fcb89f7e75 polishing (partly redoing changes that got lost when going to rev 3000) 2014-05-23 17:13:08 +00:00
Martin Diehl a284e7e6c0 more changes related to new state 2014-05-22 15:24:12 +00:00
Martin Diehl adfa634d05 improved on new state, fixed wrong indices in output 2014-05-22 15:16:05 +00:00
Martin Diehl 052aee6d04 new state copying/restoring 2014-05-22 12:07:50 +00:00
Martin Diehl 8480fc706d some changes on the new state, J2 seems to work already 2014-05-22 10:37:57 +00:00
Martin Diehl d9e8e8fc10 added consitency check: At least one phase must be present 2014-05-21 13:03:59 +00:00
Martin Diehl 8044c5b737 improved on IO_read for recursion case 2014-05-21 10:03:57 +00:00
Philip Eisenlohr 5d21ef5a84 fixed selective debugging logic flaw in crystallite_stressAndItsTangent 2014-05-19 16:43:32 +00:00
Martin Diehl 880f0a0c48 remove OMP Critical for output in regions not parallel 2014-05-16 14:01:27 +00:00
Martin Diehl 8be013b6ea removed unneeded omp pragmas 2014-05-16 10:23:10 +00:00
Philip Eisenlohr ae7adcaa2b error 100 now reports iostat return value (as "element") 2014-05-15 13:08:02 +00:00
Philip Eisenlohr 88f94cbbec increased debug threshold for some outputs (decreased verbosity level at "basic") 2014-05-15 13:01:54 +00:00
Martin Diehl e62c5fdc05 added consistency check in material.f90: Microstructure index in geometry must not exceed number of sections in material.config
better error messages in case of recursive file input in IO.f90. also supports absolute path now
2014-05-15 09:40:43 +00:00
Martin Diehl 715e7fd918 some small iprovements: reading in geometry in mesh, readability of math, output formatting in fesolving and debug, hdf5 in constitutive 2014-05-15 08:52:16 +00:00
Martin Diehl 4da866b29b exchanged race condition prevention from waiting random amount of time to omp critical 2014-05-14 13:57:25 +00:00
Martin Diehl 4bfced1a70 fixed small allocation flaws for new state 2014-05-12 13:00:37 +00:00
Martin Diehl 838a8be1b9 removed forgotten debug statements from last commit
fixed material.config for abaqus example
2014-05-12 09:28:32 +00:00
Luv Sharma d820a5aaa6 fixed some errors in changes related to newstate 2014-05-12 00:44:44 +00:00
Martin Diehl ee31bb1cae some changes related to new state 2014-05-09 09:04:09 +00:00
Martin Diehl 1298f6ea5e updated material.config of tests, had too many slip systems defined.
simplified none homogenization
2014-05-08 17:44:28 +00:00
Martin Diehl 986926aaf2 some changes related to new state
fixed assigning of wrong output in J2 model
2014-05-08 14:55:19 +00:00
Martin Diehl 72c7cfc7af still some linking issues 2014-04-30 21:13:06 +00:00
Martin Diehl e8df11b62f splitted link and compile options 2014-04-30 14:13:21 +00:00
Martin Diehl 8cc3bf4049 added forgotten linker options 2014-04-30 13:22:08 +00:00
Christoph Kords c027e8e1a8 added comment and citation for non-Schmid matrices for bcc 2014-04-30 09:56:30 +00:00
Martin Diehl 5b5a6ceaac makefile modified not to link against wrong libraries when using PETSc, reading in of values in phenopowerlaw improved, HDF5 in mesh 2014-04-29 17:50:59 +00:00
Martin Diehl 6bce04df48 changed linker command slightly to ensure dynamic linking and link now agains non-parallel (but threadsafe) versions of LAPACK
(compile_CoreModule.py, Makefile, abaqus_v6.env, mod_MarcMentat)

openmp is not working for Abaqus exp, probably because the ifort 14 is not supported.
Test is now using a version without openmp
2014-04-29 16:43:59 +00:00
Martin Diehl 01416b7ab6 set num threads back and forth between abaqus and DAMASK (like marc) 2014-04-29 16:11:13 +00:00
Martin Diehl 3fdf602acc modified reading in of values (realizes, if phase is not used in the model at all), started to include HDF5 output and new state 2014-04-15 10:49:24 +00:00
Martin Diehl 9b220c5c9a modified reading in of values (realizes, if phase is not used in the model at all), started to include HDF5 output and new state 2014-04-15 10:09:20 +00:00
Martin Diehl 6314271136 polishing, HDF5 related changes 2014-04-15 09:49:50 +00:00
Martin Diehl 6f2f04fa49 added and modified a bunch of HDF5 related functions/subroutines/variables 2014-04-15 09:43:35 +00:00
Martin Diehl 302423d266 introduced new state structure, compile with STATE=NEW 2014-04-15 09:20:38 +00:00
Martin Diehl eca211588e new compiler is more sensitive, needed to disable complaints in case of DEBUG=ON 2014-04-11 08:00:20 +00:00
Martin Diehl ac3bd2d725 space in frot caused compilation problem, fixed now 2014-04-07 14:35:40 +00:00
Martin Diehl e62bf8b1b7 some minor polishing on the way 2014-04-04 14:40:30 +00:00
Martin Diehl c712f30635 some hdf5 related functions added 2014-04-04 07:33:13 +00:00
Martin Diehl c365e1995c remove zeros for non existing slip/twin families 2014-04-02 10:59:14 +00:00
Martin Diehl 4aacca9f63 changed encodig of abaqus files, wrong directory in deploy me 2014-04-02 06:32:43 +00:00
Martin Diehl 342ca0be47 tested new scripts to update shebang, all files got same shebang (and for python files encoding) 2014-04-01 18:41:14 +00:00
Martin Diehl 436fa63912 corrected svn properties (executable, Id etc.) 2014-04-01 16:34:04 +00:00
David Mercier 11e138150a Update material config file for cp-Ti (alpha) 2014-03-31 16:49:05 +00:00
Martin Diehl 84ce6e429a changed handling of non-converged BVP solution: By default, exit if no solution is found to prevent "user errors"
Removed ambiguous "regridMode" keyword which was used to trigger this behavior in favor of  "continueCalculation", Set this to 1 to get back old behavior, e.g. report non-converged increments.
2014-03-31 10:04:11 +00:00
Franz Roters 7b27606000 modified constitutive description in line with other dislocation density based models 2014-03-30 15:04:06 +00:00
Martin Diehl 4bb5439994 removed copyright from *.f90 and *.f files, will be added using the deployMe script for release/nightly build. Now based on presence of svn:keywords property set (usually "Id"). Similar thing will be done for .py, .sh, and .config files but more involved because they might be executable 2014-03-29 08:20:36 +00:00
Martin Diehl aae3f95f76 added script to add Copyright information (complete header) to f90 files having the svn property "MPIE" set.
Development version don't need an header anymore
2014-03-26 14:57:33 +00:00
Martin Diehl 62497b648a removed forgotten debug statement 2014-03-26 09:09:28 +00:00
Martin Diehl 58ad750b2a small changes regarding restart and HDF5 2014-03-26 08:41:45 +00:00
Martin Diehl 99fe2dfccd small changes to HDF5, put link to homepage into README not to have double effort 2014-03-25 17:21:47 +00:00
Martin Diehl 2659ee51d4 restart is working for all solvers, but slight deviations still occur. Reason is most probably, that for CPFEM init, a 0s timestep is required which seems to alter the state a tiny bit, leading to small deviations. For AL and Polarization, this is even more severe since the regular call to CPFEM has F-residual as the current deformation gradient, however the init step uses only F. 2014-03-25 15:44:16 +00:00
Martin Diehl c7994db5ef writing out f_aim and f_aim_last inc for restart, was calculated as average before 2014-03-24 09:59:30 +00:00
Franz Roters 04b2011195 reintroduced dislocation hardening for bcc 2014-03-17 15:42:46 +00:00
Martin Diehl bac0e24d21 added parameter to dislotwin 'dipoleformationfactor'. set it to 0.0 to turn hardening due to dipole formation off (like in the updated tungsten example). standard behavior (e.g. a value of 1) is recovered if no value is given 2014-03-13 23:50:55 +00:00
Martin Diehl 5d4017bbfc introduced homogenization_none to substitute isostrain with ngrains 1.
cleaned up homogenization isostrain (has no state)
2014-03-13 23:20:50 +00:00
Duancheng Ma 670c0caabc 2014-03-13 12:12:19 +00:00
Martin Diehl 3f7a389ff7 changed state parsing for local models (and for delta_state) such that only the needed part of the state array (for the given material point) is used 2014-03-13 06:43:49 +00:00
Martin Diehl 7c79b31f6c remove not needed use statement 2014-03-13 05:50:56 +00:00
Martin Diehl e1a31457ba fixed potential array mismatch in assignment 2014-03-13 05:49:07 +00:00
Martin Diehl fced0168f0 added a little bit more HDF5 functionality
made abbreviations in  configure options consistently capitals
2014-03-12 16:51:01 +00:00
Martin Diehl 9afc1e3a15 fixed wrong omp statement 2014-03-12 15:29:14 +00:00
Martin Diehl 3aea8b39e9 added some HDF5 functionality (needs to be activated with preprocessor makro) 2014-03-12 07:33:51 +00:00
Martin Diehl ff1b1c1a50 fixed bug introduced with lattice_structure change nonlocal, but with DEBUG=ON OPTIMIZATION=OFF there is an FPE. Division by zero? Marked in the code
forgot to commit dislotwin last time, now seems to work
2014-03-11 23:55:40 +00:00
Martin Diehl ef8fbf4dda introduced case in dislotwin for bcc (peierls stress as critical stress)
renamed some parameters, now exponents (p,q) for slip are per family, shear band got own ones. exponent for twin (r) now per family

fixed bugs in lattice ("empty" interactions should be 1 not zero)

see both Phase examples to check what is needed
2014-03-11 17:40:59 +00:00
Franz Roters 6781d7d8a5 last line got lost
as it is fcc it should actually be only one slip system family
2014-03-11 10:49:48 +00:00
Martin Diehl c2e81414e3 updated TWIP steel data set to have only 2 slip system families, added data set for tungsten 2014-03-11 10:34:36 +00:00
Martin Diehl 721183e702 removed some forgotten debug statements 2014-03-08 22:51:32 +00:00
Martin Diehl 2b589c3d71 moved reading in of lattice type and elastic constants to lattice module
removed structure type for hex, fcc, bcc, now defining slip/twin systems for each phase found in material.config
constitutive modules will only be initialized if needed
homogenizedC function is only needed for models incorporating twinning in a physical way (titanmod and dislotwin)
2014-03-08 20:50:31 +00:00
Martin Diehl 32493675d6 fixed bug in constitutive introduced in rev 2988 causing wrong names in *.outputConstitutive 2014-03-05 08:06:21 +00:00
Christoph Kords 0265978941 renamed "maxNmatIDs" accordingly to "maxNinstances" 2014-03-04 13:47:04 +00:00
Martin Diehl a0d75ee05e moved quaternion disorientation to lattice because it requires the knowledge of the lattice structure. 2014-02-28 13:28:27 +00:00
Martin Diehl 92bf1565fc formatting 2014-02-28 13:05:11 +00:00
Martin Diehl b0f191c88c also rename matID/i to instance like in the other constitutive models 2014-02-28 13:03:21 +00:00
Martin Diehl b9722b866c improved error handling when getting number of threads 2014-02-28 10:30:07 +00:00
Martin Diehl 650b71ffa9 renamed instance consistently to "instance" as a preparation for new structure / elastic matrix handling 2014-02-28 10:18:40 +00:00
Martin Diehl fe13bb16fe added some consistency check and remove unnamed constants 2014-02-27 18:57:46 +00:00
Martin Diehl ad203d32c1 circular dependencie (copy and paste error) removed 2014-02-26 12:04:52 +00:00
Martin Diehl 35be468b24 fixed detection of working directory for ifort 2014-02-25 21:10:55 +00:00
Martin Diehl a554bc5267 now ignoring c/a ratio keyword 2014-02-19 07:40:21 +00:00
Martin Diehl 69992638a8 added missing covera_ration in hex test. was not necessary before for constitutive_none, now it is 2014-02-17 09:45:54 +00:00
Martin Diehl 8bbba55506 function version seems to be quasei-standard, not subroutine for getcwd 2014-02-14 19:59:37 +00:00
Martin Diehl 4efb1e5617 fixed bug: wrong initialization of lattice structure, had effect on j2 and phenopowerlaw 2014-02-14 18:05:59 +00:00
Franz Roters d3d4931061 corrected checking for ~/.damask/damask.conf 2014-02-14 15:58:01 +00:00
Franz Roters 4399fb9535 improved comment on Euler angles, now includes unit 2014-02-10 15:48:00 +00:00
Martin Diehl d8311d46b3 .. last commit message was missing because I wanted to abort, so:
move reading in of lattice structure to lattice (elastic constants, c/a ratio, lattice structure). Still needs a little bit of documentation, the lattice structure names are not easy to understand.
tested it for j2 and phenopowerlaw and seems ok, but before continuation with the other constitutives models might simplifiy a little bit.
2014-02-10 14:38:05 +00:00
Martin Diehl d45aea4467 moved reading in of lattice stru 2014-02-10 14:31:19 +00:00
Martin Diehl 20bea80f0c changed order to have PETSc variables before others, fixed problem when MPI and FFTW were installed in standard location but PETSc build own MPI 2014-02-09 18:30:31 +00:00
Martin Diehl 33ab0ff0fe introduce enums, remove gdot_twin and tau_twin from possible outputs because they weren't output anyway 2014-02-09 17:08:13 +00:00
Martin Diehl a07c738778 polished a little bit 2014-02-08 10:48:09 +00:00
Martin Diehl 8741f2368d first (not too complicated) workaround for the shortcomings of the new intel compiler. still crashes, but later ;) 2014-02-06 17:48:01 +00:00
Martin Diehl bc4cc20c55 slightly polished makefile and numerics, recompile fftw with icc; still not running with ifort 14.01 2014-02-06 10:41:34 +00:00
Martin Diehl 1ea48bb8ff added accumulatedshear as output for twin and slip 2014-02-05 15:26:22 +00:00
Martin Diehl 83e4bc3247 updated test for restart, now also checking if number of increments is the same
fixed output of test class
fixed small bug? in asciitable.py, @philip: please check
2014-02-03 19:30:28 +00:00
Martin Diehl 1444739d60 prevent ifort from complaining too loud about old syntax in openmpi files 2014-02-03 16:47:16 +00:00
Christoph Kords caca65148d fixed bug in restart: stiffness values that were read from file used to be overwritten by "utilities_constitutiveResponse"
Still not fixed in polarization and AL solver !!!
2014-02-03 15:57:04 +00:00
Christoph Kords b9b87a785c added crystallite_F to debugging output 2014-02-03 13:19:49 +00:00
Franz Roters c5b31b373c no longer needed 2014-02-03 10:38:59 +00:00
Christoph Kords aa4eee2b50 fixed bug in math_rotate_forward3333 2014-01-31 14:28:20 +00:00
Martin Diehl c7e59821d1 fixed compilation bug when using IMKL single core 2014-01-30 15:38:10 +00:00
Martin Diehl 986591682c remove equals sign in rpath definition to work with mac, changed linker options in compile_CoreModule to work with linux again, please let me know if mac doesn't work anymore (it should) 2014-01-30 13:07:44 +00:00
Martin Diehl 9e5a2d8e10 fixed runpath 2014-01-29 19:52:53 +00:00
Martin Diehl f82e5ec26e some small changes on the installation scripts, IMKL now works with gfortran, more verbosity and more sanity checks 2014-01-29 12:20:56 +00:00
Martin Diehl ee99216edc polishing compile_CoreModule, bug fix for Makefile 2014-01-28 07:34:13 +00:00
Martin Diehl 02441bef0f introduced rpath to tell executable where libraries are located 2014-01-27 14:10:49 +00:00
Pratheek Shanthraj b9e55b4495 moved setting of CPFEM forwarding flag to the forwarding subroutine 2014-01-27 11:05:36 +00:00
Martin Diehl 656e3cc6ca some more small changes, including bug fix in configure + relative linking in Makefile for spectral solver 2014-01-27 10:12:29 +00:00
Martin Diehl 4850e5987d remove LD_LIBRARY_PATH (http://linuxmafia.com/faq/Admin/ld-lib-path.html) and fixed wrong named xxx_ROOT 2014-01-25 00:15:04 +00:00
Martin Diehl d66b41c686 fixed installation, will test now 2014-01-24 12:30:27 +00:00
Martin Diehl c99ca8b0ee fixed magnesium parameters 2014-01-23 13:33:09 +00:00
Martin Diehl e6d88e4458 prevent twin volume fraction from going above 1.0 2014-01-22 15:47:49 +00:00
Christoph Kords 391d83fa2b moved rate sensitivity block to where it belongs to, namely the analytic jacobian calculation 2014-01-22 15:34:10 +00:00
Pratheek Shanthraj 4b7890f133 fixed wrong indices in analytic jacobian calculation following christoph's changes. convergence is restored at least for the spectral solver 2014-01-22 10:16:55 +00:00
Christoph Kords 3dfdbaff5b Fixed wrong indices in tangents dT_dFe and dFe_dLp, which however luckily did not have any effect in the perturbed stiffness since they were transposed such that the double contraction of both remained unchanged.
In contrast, the analytical jacobian will probably be affected by this change!
@Pratheek: Can you check with me how this can be fixed?
2014-01-22 08:38:13 +00:00
Christoph Kords c10d0d15cb crystallite_integrateStateFPI: convergence check always ran about all e,i,g; now checking only those specified in FEsolving_execElem, FEsolving_execIP 2014-01-21 18:45:41 +00:00
Christoph Kords ff648492c1 perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode
(missed some part in last commit)
2014-01-21 18:28:21 +00:00
Christoph Kords 029826d84f perturbation calculation in crystallite_stressAndItsTangent: only backup and restore those chunks of crystallite_Fe, _Fp etc that are actually used in order to avoid pointless data copying in non ping pong mode 2014-01-21 16:35:12 +00:00
Christoph Kords 351aa6ba29 setting todo and converged flag only for those e,i,g that are to be executed (FE_elemExec etc)
used to spoil convergence checking in non ping pong mode and lead to slow performance
2014-01-21 16:09:00 +00:00
Philip Eisenlohr 74cbfb6d48 added few comments 2014-01-17 01:38:35 +00:00
Christoph Kords 77cba87bf3 in crystallite_integrateStress: pass full 2nd PK stress to LpAndItsTangent instead of only the deviatoric part; otherwise constitutive laws that depend on the hydrostatic pressure do not work; constitutive laws that require a deviatoric stress need to calculate it internally (as is e.g. already done in the J2 model) 2014-01-16 10:36:40 +00:00
Christoph Kords a011b0a4a2 corrected tangent of plastic velocity gradient dLp_dT for slip (twin volume fraction was missing) 2014-01-16 09:53:45 +00:00
Christoph Kords 9223cf52b4 corrected non-Schmid effect: used to calculate resolvedstress based on Mandel notation which can be incorrect for nonsymmetric non-Schmid tensors 2014-01-16 09:50:40 +00:00
Christoph Kords b24dd62431 corrected tangent of plastic velocity gradient dLp_dT for slip (twin volume fraction was missing) 2014-01-16 09:29:21 +00:00
Christoph Kords 563c4ac94b corrected tangent of plastic velocity gradient dLp_dT, which improves convergence a lot 2014-01-16 09:18:26 +00:00
Martin Diehl 278f854e41 introduced functionality to get output results by point from the HDF5 file 2014-01-15 22:17:32 +00:00
Martin Diehl 3fffb2ae3e added parameter set for Magensium 2014-01-15 12:29:44 +00:00
Martin Diehl 9f3a17b1b6 fixed bug in hexagonal slip system definition 2014-01-15 12:08:46 +00:00
Franz Roters 98f677916c fixed crash in RtoEuler for PHI=180° 2014-01-15 08:32:43 +00:00
Martin Diehl 449c66c70f did some testing... 2014-01-14 19:32:55 +00:00
Martin Diehl 01c8c434f7 finished work on the new installation scheme (if no bugs are found ;)) 2014-01-14 19:03:41 +00:00
Martin Diehl aab0208c65 more modifications, only configure needs modifications 2014-01-12 17:29:00 +00:00
Martin Diehl 00eb9b2cf6 forgot Makefile 2014-01-11 17:18:59 +00:00
Christoph Kords a1927b60a7 allow flux between neighboring material points if one of both got a velocity of zero 2014-01-08 16:35:10 +00:00
Christoph Kords b41593f279 blocked density does no longer contribute to dislocation flux; used to generate source term for mobile density if burgers vector had any non-zero component perpendicular to blocking boundary. 2014-01-08 16:31:50 +00:00
Christoph Kords 629295242f added two new outputs: "rho_dot_flux_mobile" and "rho_dot_sgl_mobile"
corrected outputs of dislocation density rates: take into account that blocked dislocation density is a signed quantity
2014-01-07 15:55:28 +00:00
Martin Diehl 79307de6a7 added some code parts for HDF5 in Fortran 2013-12-30 10:06:01 +00:00
Martin Diehl 5b4b132176 introduced preprocessor statements in mesh to hide unneeded variables as a preparation for HDF5 2013-12-27 20:03:28 +00:00
Martin Diehl 831e905b16 added missing lp output to crystallite 2013-12-20 11:13:12 +00:00
Pratheek Shanthraj de3fc70a02 moved restart writing to the forwarding routine 2013-12-20 10:49:14 +00:00
Martin Diehl 98d5a99a43 reverted unwanted changes to Makefile (configure) 2013-12-20 10:31:33 +00:00
Martin Diehl d595f59f0e improved sanity checks, mainly for RGC 2013-12-20 08:36:15 +00:00
Martin Diehl af075aaa9d found bug (array index out of bounds), Philip/Franz please check if correct (Gfortran was complaining) 2013-12-20 08:33:46 +00:00
Franz Roters 543639bb8f added support for Marc/Mentat2013.1 2013-12-19 13:37:35 +00:00
Martin Diehl 61c6839723 fixed over-sensitive error in homogeniztion and ifort option in makefile 2013-12-19 08:49:47 +00:00
Pratheek Shanthraj 6f6bd86078 fixed bug in heat generation and analytic jacobian calculations 2013-12-18 10:19:48 +00:00
Pratheek Shanthraj f3ea92f50c fixed error in previous commit 2013-12-18 09:38:44 +00:00
Pratheek Shanthraj 61981617d7 separated forwarding and solution subroutines for better control at the load step looping level 2013-12-18 09:35:05 +00:00
Pratheek Shanthraj ff6211b78c rolled back phase field changes 2013-12-18 09:09:32 +00:00
Martin Diehl 6fa9ed8f48 homogenization: added enums and sourced allocation for RGC, some higher verbosity for isostrain 2013-12-18 07:28:01 +00:00
Martin Diehl c4a592867f updated to PETSc 3.4.3, will not compile on all workstations except for maws0X 2013-12-17 15:37:14 +00:00
Martin Diehl 676e621af2 new version of Lib_VTK_IO 2013-12-17 13:24:34 +00:00
Martin Diehl 7b9d76cdcc added missing allocation of outputID caused crashing 2013-12-17 08:02:21 +00:00
Martin Diehl ba89df7784 continued with sourced allocation and EOL improvements 2013-12-16 11:58:03 +00:00
Martin Diehl 118dff9f23 fixed wrong initialized symmetry value 2013-12-16 10:56:56 +00:00
Martin Diehl d5bdcf5a6d some small improvements, still something wrong, probably in crystallite.f90 2013-12-13 14:14:17 +00:00
Martin Diehl 84d51d743c fixed small issues, ATTENTION: Still a bug introduced with rev. 2783 not fixed 2013-12-13 13:19:17 +00:00
Martin Diehl 7ba3cc0ab7 removed symlink to marc 2010 2013-12-13 10:04:25 +00:00
Martin Diehl acad3d4a33 using sourced allocation 2013-12-12 23:03:37 +00:00
Martin Diehl da337962af fixed bug for Abaqus input parsing 2013-12-12 22:29:40 +00:00
Martin Diehl 65ae979920 indroduced sourced allocation, enums where applicable (some parts still missing). fixed bug when having recursive file input 2013-12-12 17:09:59 +00:00
Christoph Kords 9b9f4dd624 adapted changes from rev 2776 for constitutive_nonlocal: improved usage of "enum", fixed bug in file reading, usage of "source" specifier for "allocate" 2013-12-12 10:16:50 +00:00
Christoph Kords dd77690a68 fixed bug introduced in rev2777: renamed all remaining occurrences of variable "ncomponents" to "nconstituents" 2013-12-12 10:06:39 +00:00
Martin Diehl 102712d91f added enum for dislotwin output, fixed bug when using recursive file input function 2013-12-11 23:42:33 +00:00
Martin Diehl 95d6430b09 improved on enums (introduced them for output in j2), fixed wrong recursion depth bug in IO, fixed a bug in reading in for none, j2, phenopowerlaw 2013-12-11 22:03:09 +00:00
Martin Diehl 7885ebaf8f added 'reset' flag for recursive function IO_read, need to reset stack when stopping exectution (will be used by constitutive and homogenization) 2013-12-11 16:49:20 +00:00
Franz Roters efa7a6dae2 added example material.config again 2013-12-05 15:11:30 +00:00
Franz Roters bf6e749abd adopted deployMe.sh to new svn server
editorial changes in LICENSE and CPFEM.f90
2013-12-05 14:31:28 +00:00
Christoph Kords 61ebfb4c85 removed obsolete output 'boundarylayer' from material.config 2013-11-29 12:30:06 +00:00
Christoph Kords b0695c0904 enumerators for case switch of outputs as already used in phenopowerlaw
removed obsolete output 'boundarylayer'
2013-11-29 12:29:47 +00:00
Christoph Kords f56cdafccb polishing indentation levels 2013-11-29 12:25:36 +00:00
Martin Diehl 2d624d022f fixed small bug in detecting the lattice structure. only first 3 letters count (e.g. 'ort' and 'orthography' will give the orthorombic lattice type) 2013-11-28 08:56:02 +00:00
Martin Diehl 345e5cd47c split up material.config and added a small helper script to create wiki content from all material.configS
(Work done by Luv, just his TortoiseSVN is not yet installed)
2013-11-27 17:49:34 +00:00
Martin Diehl 0bc99a9622 change latticeName to latticeID 2013-11-27 16:20:27 +00:00
Martin Diehl 6c2ab1f1f8 added ID (integer, enum) case switching for phenopowerlaw and isostrain 2013-11-27 11:39:28 +00:00
Martin Diehl 080edffd87 file was missing in last commit 2013-11-27 08:05:23 +00:00
Martin Diehl 2f7efa2055 introduced named constants for homogenization_type, phase_plasticitiy, and phase_elasticity to replace expensive string comparison. Labels for different schemes need to be stored in material.f90 because they are needed when parsing the config file 2013-11-27 08:04:05 +00:00
Martin Diehl c9ce8d2c52 some improvements on HDF5 test, polishing math 2013-11-22 10:31:52 +00:00
Christoph Kords e8de8f69a8 renaming of "maxNinstance" to be consistent with naming scheme introduced in rev2662 2013-11-21 13:40:14 +00:00
Christoph Kords 3cc6781eaf renaming of "instance" and thelike to be consistent with naming scheme introduced in rev2662 2013-11-21 13:35:43 +00:00
Christoph Kords 883669bd77 in preguess of FPintegrator: resetting of previous dotstates can be done in same loop as collectdotstate
polishing: indentation level and capitalization
2013-11-21 10:58:41 +00:00
Martin Diehl 9083aa53fb fixed long lines due to makro expansion 2013-11-19 16:06:53 +00:00
Pratheek Shanthraj 9c34c1b871 error reporting for phase fields only if there are active phase fields 2013-11-14 10:36:49 +00:00
Pratheek Shanthraj f51c600189 convergence check for phase field only if there are active phase fields 2013-11-14 10:32:41 +00:00
Pratheek Shanthraj 3b113f6a62 more errors :) 2013-11-13 19:38:25 +00:00
Pratheek Shanthraj 8444b4116d corrected error in previous commit 2013-11-13 19:26:36 +00:00
Pratheek Shanthraj 50db944c0c reworked phase field interface to damask spectral solvers. now specify 'thermal a b c' or 'fracture a b c' to activate either phase field where a b c are the initial value, diffusion coefficient and mobility respectively.
Right now only thermal and fracture phase fields implemented and only in the basic petsc solver
2013-11-13 19:21:35 +00:00
Pratheek Shanthraj 724380c624 Coupled phase field implementation for the basic petsc solver. needs clearing up… work in progress 2013-11-12 00:44:23 +00:00
Martin Diehl 9ac8846e41 compiling with hdf5 is now working for DAMASK (if HDF5=ON) and for playground example in private 2013-11-07 12:52:36 +00:00
Philip Eisenlohr c857381881 added "-d", "--directory" as possible cmd switches for working directory.
fixed (sometimes?) non-working check for valid working directory path.
2013-10-24 18:56:56 +00:00
Martin Diehl ea980f428e fixed long line coming from PETSc macro expansion 2013-10-23 17:52:58 +00:00
Pratheek Shanthraj 57cf472982 Improved convergence checking when using newton solver (does not oversolve the problem like before) 2013-10-23 15:22:12 +00:00
Philip Eisenlohr eb0d2b7e24 fixed bug that (sometimes, always?) prevented reading in of debug.config, when variable 'line' was still '#EOF#' from former IO_read... 2013-10-23 11:23:41 +00:00
Philip Eisenlohr bdafc6f1b8 condensed error msg numbering,
introduced new error for "microstructure without any constituent"
2013-10-23 11:21:48 +00:00
Martin Diehl 205eb66ba8 bugfix for last commit, too many variables declared as private 2013-10-23 10:29:38 +00:00
Martin Diehl 96b03841e8 redefined macro for intel compiler, should work with windows (but Nicolo had some problems)
made configure verbose
2013-10-23 08:38:14 +00:00
Martin Diehl 7b28dedfe3 removed temperature calculation and debugging statements 2013-10-22 15:52:58 +00:00
Martin Diehl ea664688f8 introduced dummy temperature calculation. 2013-10-18 20:56:10 +00:00
Martin Diehl 3ecc8103f0 introduced dummy heat calculation, available as an crystallite output
temperature, avgP and avgF are available as homogenization output. move crystallite output of ipcoordinates to homogenization
2013-10-18 18:57:28 +00:00
Martin Diehl 1c7ed337d7 corrected linebreak in comment 2013-10-17 12:52:46 +00:00
Martin Diehl 03e26af37b small bufix, didn't compile 2013-10-16 20:44:04 +00:00
Martin Diehl dc95c82d4a removed temperature integration and corresponding data structures and debugging options
temperature is stored in crystallite, but homogeneous on one IP (not an component (grain) quantity and an input value parsed in by the BVP solver.
introduced heat, a component (grain) quantity which is homogenized before returned to the heat transfer solver.
went ahead with removal of dummy functions in homogenization and constitutive, this time mainly reduced function signatures to reflect actually needed quantities.
2013-10-16 13:04:59 +00:00
Martin Diehl 6a1c40d540 introduced error when element/IP/component(grain) selected for debugging does not exist 2013-10-16 12:38:00 +00:00
Philip Eisenlohr 0a292d8d8b fixed final bug in hex twintwin matrix 2013-10-15 22:02:39 +00:00
Franz Roters 4d184b1ee4 improved comments 2013-10-14 14:35:41 +00:00
Martin Diehl 295d0cd28e removed a bunch of dummy functions (delta state for all constitutive_XXX except nonlocal, microstructure for none, j2 and phenopowerla), additional simplifications for none, averageBurgers is only relevant for RGC 2013-10-14 10:54:45 +00:00
Martin Diehl 029502c12e fixed typo causing compilation error 2013-10-14 06:11:50 +00:00
Martin Diehl b651f334fe set useful default values to remove dummy functions in individual homogenization schemes (stateInit, averageTemperature, and updateState (isostrain only)) 2013-10-11 16:01:53 +00:00
Franz Roters 80cb78c1fd corrected activation of slipbands; both velocity and resistance have to be non zero 2013-10-11 13:52:00 +00:00
Martin Diehl 3b5a5ab812 changed error in case of number of found interaction types is less than expected to warning 2013-10-11 09:17:03 +00:00
Martin Diehl bb0b35e466 made compilation_info.f90 windows compatible and added abaqus_v6_windows.env as a template for using DAMASK with windows.
Deleted old fftw and ACML installation instruction (can be found on damask.mpie.de)
2013-10-10 09:38:31 +00:00
Martin Diehl ed47d25adc remove dotTemperature (returned always 0 anyway) to make it easier to replace it with a heating rate funtction
remove potentially dangerous short circuiting in homogenization files
added doxygen comment for nonlocal (giving the functions a name)
fixed wrong definition of maxnchungs introcuded into dislotwin with last commit, causes trouble with hex
2013-10-09 06:12:16 +00:00
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
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