2ba15e462a
corrected values and typo
2015-10-26 19:27:57 +00:00
c91788d721
new intel compiler was compiling about *** in float output
2015-10-26 17:11:36 +00:00
1f356a6833
added parameters from indentation results
2015-10-22 17:45:28 +00:00
04841cb727
List of dependent trans states
2015-10-21 13:24:26 +00:00
20b0429517
Calculation of equilibrium separation of partials for twin and trans
2015-10-21 09:37:45 +00:00
565d68edd7
Added some TRIP inputs
2015-10-20 15:54:01 +00:00
b9fcd9482b
Renamed some input variables to denote twin
2015-10-20 13:45:06 +00:00
fc8a11f027
corrected typo: alias for first Piola stress should read 1stpiola NOT 1piola
2015-10-19 06:57:23 +00:00
f2b29e4844
correct communication between
...
elements with different
phases
2015-10-14 21:04:40 +00:00
484a34b7f1
added pheno+ module
...
modify crystallite microstructure call
to pass orientations
2015-10-14 18:36:19 +00:00
ea0b16c726
pheno plus for commercial FEM
2015-10-13 21:06:59 +00:00
d6abc00218
add pheno+ module in
2015-10-13 18:52:01 +00:00
8fac635c15
made hybridIA stuff working again
2015-10-13 17:02:07 +00:00
6de8c821f8
improved consistency check on loadcase
2015-10-12 19:02:42 +00:00
9ec7c08bb6
gfortran was complaining abound mismatch
2015-10-07 13:07:21 +00:00
c1a53bc8c0
introduced dummy target to check compilation
2015-10-06 18:00:17 +00:00
0f0b4eaa6b
cleaned
2015-09-28 17:47:00 +00:00
ad16162eb8
simplified
2015-09-24 17:50:11 +00:00
c50a522dd8
some shapes not 100% correct, removed inverse laplace
2015-09-24 17:38:49 +00:00
2e70e6fded
make clear that this isNaN is an intrinsic extension
2015-09-24 08:45:44 +00:00
b528088653
added all 13 slip families for bct(beta)-Sn
...
increased the allowable columns for lattice.f90 to allow for the large interaction matrix.
2015-09-22 18:42:23 +00:00
42ac7902f7
indents wrong in fortran, seeds_fromDistribution bugfix and using ASCII table facilities
2015-09-20 21:03:10 +00:00
41f09dd0f8
should make no difference, but test with ifort was failing. compierler buig?
2015-09-12 18:26:25 +00:00
03896fcd09
redefined required debug variables removed in last commit
2015-09-11 18:20:37 +00:00
840fd0e9fa
report absolute value of accumulated_shear_slip in dot state
2015-09-11 17:44:26 +00:00
8716752bf8
fixed fourier convolution and div/curl calculation for even/odd grid according to Johnston 2011 (MIT, FFTW)
2015-09-11 13:05:46 +00:00
8f32d03a9e
curl calculation overestimated RMS due to factor 2 instead of one for DC component and Nyquist component
2015-09-11 09:19:13 +00:00
e88cedc6ae
clearer naming, debug options for spectral do not work for MPI
2015-09-11 08:52:03 +00:00
751d1d7582
function statement is more clear
2015-09-10 13:01:33 +00:00
580bb76b4c
pInt missing
2015-09-10 10:26:09 +00:00
427b5f4bc1
wrong results for restart
2015-09-10 10:06:14 +00:00
bbe37c842e
name conflict
2015-09-09 21:52:00 +00:00
34980a1d44
more logical structure when reading in
2015-09-08 19:44:32 +00:00
7d996cfa11
fixed new IO_stringPos for Abaqus
2015-09-05 16:26:55 +00:00
caf0ac7e8d
more verbose
2015-09-01 16:53:48 +00:00
7a3a67601f
realloc lhs was not always working, fixed now and added test for new string pos function
2015-08-31 16:30:04 +00:00
197ae53553
lots of memory savings (~ stateside + 250 pReals per integration point) when using analytic jacobian
2015-08-28 10:55:38 +00:00
6848d83d13
behavior of IO_stringPos was unsafe in some situations, also took the chance to rename the variables used in its context to a more meaningful name.
...
functions in the core module are not used, including kdTree
2015-08-28 07:38:48 +00:00
1171dc4344
MPI has a 2GB limit for writing at once, now chunking
2015-08-21 17:51:05 +00:00
8474da26d9
empty line
2015-08-18 16:31:47 +00:00
bae5cfcf89
named if
2015-08-18 16:07:01 +00:00
53c574033a
size of MAXCHUNKS same as maximum allocated size for “heat_time” & “heat_rate”
2015-08-18 03:20:16 +00:00
f64a5e7552
output safe for 64bit integer
2015-08-14 05:03:54 +00:00
31ea4dadf1
IO_stringPos replacement not needing MAXNCHUNKS but making use of new Fortran features for dynamic allocation
2015-08-13 14:54:34 +00:00
563d9e64dd
MAXNCHUNKS capitalized since it is a parameter
2015-08-06 09:24:56 +00:00
3d94ae4074
finite strain correction for respective expansions
2015-08-05 14:12:54 +00:00
09c7affbbc
more efficient rate calculation in updateState when subdt is zero
2015-08-05 11:09:38 +00:00
7928ce3ac3
fixed erroneous rate calculation in updateState when time increment (subdt) was zero during initialization.
2015-08-04 21:27:30 +00:00
65d114e4f0
improved debug statements
2015-08-04 21:26:22 +00:00
905b9f8f42
petsc update (now v3.6.1)
2015-08-04 15:04:53 +00:00
5b57ea7d16
fluctuations proportional to thermal equilibrium concentration
2015-08-03 11:08:59 +00:00
a01d7d8d66
corrections to crystallite_push33ToRef function
2015-08-03 11:07:19 +00:00
1ed276ce2d
quadratic softening for finite dissipation energy
2015-07-30 10:29:25 +00:00
c97b83c9a1
added user defined external heat source
2015-07-27 11:09:37 +00:00
4bb6664de9
tightened default thermal tolerance
2015-07-27 11:08:43 +00:00
c82c4b74f1
updated config files to reflect recent changes in input parameters
2015-07-24 17:14:03 +00:00
5e09954575
remove deprecated spectral load case definition of temperature
2015-07-24 14:57:29 +00:00
7554647c8e
more control over initialisation of field values. specify initial field value in the homogenisation part of the material config file using the appropriate tags
2015-07-24 14:53:50 +00:00
87d42bf447
initialise Fi correctly for initial field values away from equilibrium
2015-07-24 14:47:18 +00:00
0a5ccb3d91
possibility to set damage potential order
2015-07-24 14:43:05 +00:00
7c101cdc31
fixed issues causing compilation errors
2015-07-23 19:22:03 +00:00
7f8a8d5b0f
new function to report initial thermal strain based on current temperature deviation from reference.
...
to be used in constitutive_initialFi, which collects from all kinematics, derives and sums Li terms to report initial Fi.
2015-07-23 15:59:25 +00:00
25a03e128a
added remark to include looping over sources of same type within one phase. (for instance to accommodate different activity periods.)
2015-07-16 17:20:18 +00:00
d2c748e716
TRIP input parameters for fcc to hex
2015-07-10 09:02:48 +00:00
1e9682eef9
corrected upper bound for spectral damage solver as well.
2015-07-09 15:51:34 +00:00
d2bfd26937
1) added the error message of all MPI plans 2) deleted the suffix (MPI) for all MPI plans
2015-07-09 13:38:21 +00:00
7a4415b28a
Sanity checks for SFE
2015-07-08 15:37:16 +00:00
df5cdb6b01
Elasticity matrices for transformed phase
2015-07-08 11:58:52 +00:00
66564c1f75
Read in elastic constants for transformed phase
2015-07-07 13:23:56 +00:00
4435244681
corrected missing initialization. newState = oldState + deltaState. But for irreversibility driving force must be compared with last converged subState.
2015-07-06 09:34:43 +00:00
3d5a85e4af
corrected implementation of deltaState for isoBrittle damage
2015-07-06 08:47:40 +00:00
72b7df34a3
save some space
2015-07-02 09:45:25 +00:00
74af77a81f
formatted the interaction matrix for beta-Sn
2015-07-02 05:45:56 +00:00
0e6ac411bb
further modified to add longer lines for (bct) beta-Sn in lattice
2015-07-02 05:44:45 +00:00
6272129f01
formatted code and modified Makefile to have long lines in lattice
2015-07-02 05:09:07 +00:00
e61050dcde
added lattice for (bct) beta- Sn.
2015-06-27 14:55:30 +00:00
a40d247f7f
lower bound for damage is residual stiffness
2015-06-22 14:27:32 +00:00
12bc4f01ef
fcc to hex orientation relationship (Shoji-Nishiyama)
2015-06-22 13:47:15 +00:00
ad17c40d92
Renamed variables to denote fcc to bcc transformation
2015-06-22 08:33:48 +00:00
5114e0b43e
Cleaning up
2015-06-22 08:10:20 +00:00
3e0e187b14
clarified meaning of symEulers function in help
2015-06-19 07:29:01 +00:00
b1c036daab
added some comments on possible code changes
2015-06-17 06:53:00 +00:00
0e48fb8fc7
removed double entry in use list
2015-06-16 07:57:42 +00:00
081bc33159
Read in c/a ratio for transformed hcp phase
2015-06-12 08:59:39 +00:00
a370372c1b
transformed lattice structure should not be mandatory
2015-06-11 17:03:34 +00:00
40701cedc4
damage driving force is a history variable and damage field is strictly bounded between 0 and 1
2015-06-11 09:03:51 +00:00
e8ee5d6723
moved some phase field parameters to lattice
2015-06-11 09:01:37 +00:00
b3241411f5
Added keyword to specify transformed lattice structure
2015-06-11 08:23:27 +00:00
c2781b5423
Removed unnecessary variables
2015-06-10 12:12:03 +00:00
859fb43fe5
split petsc_options into petsc_defaultOptions and petsc_options, so you only need to specify the petsc options you intend to overwrite or append to the default options
2015-06-09 13:28:50 +00:00
b410a08995
Renamed some variables for consistency
2015-06-08 14:26:50 +00:00
0d36a2d3e2
Removed unused variables and renamed some variables for TRIP model
2015-06-08 09:36:52 +00:00
74982294a0
added spectral thermal and damage solvers
2015-06-03 17:30:31 +00:00
3415e249ef
reverting files that got mistakenly added into previous change set
2015-06-01 16:11:37 +00:00
a14070bad4
changed up handling of delta states, and some bug fixes
2015-06-01 16:02:27 +00:00
8d2aa5e7be
fixed typo in dislotwin label
2015-05-29 11:03:25 +00:00
8f4663985a
major restructuring of multi field handling in DAMASK and added some example config files for multi field simulations. please report bugs
2015-05-28 17:02:23 +00:00
b58cbe424d
follow up to last commit
2015-05-19 17:41:49 +00:00
b07ec8232e
found circular inclusion
2015-05-19 17:34:04 +00:00
71b0e283c1
clearer substructure of routines
2015-05-10 20:55:36 +00:00
4f9531e736
sliptwin is not used and will not be used as discussed in the documentation in private
2015-05-05 20:17:50 +00:00
485636b93b
functions not longer needed for core module.
...
simplified some statements, removed double computation for direct matrix inversion
2015-05-05 06:37:59 +00:00
906c3f63a1
updated hybridIA sampling to work with new format
...
cropLinearODF is not working for the new format, but filterTable should be able to do the task
2015-04-27 05:30:29 +00:00
97aba96ff8
changed format of linearODF file to be more descriptive, essentially TSL OIM format with header and space separated keys
2015-04-26 11:07:39 +00:00
1f809f9d73
variable name change forgotton
2015-04-22 04:34:35 +00:00
78d913c07d
clearified keywords, added raw data of rolling texture, renamed to more meaningful names
2015-04-21 19:34:50 +00:00
d7df5f1934
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 16:11:30 +00:00
1330dae432
cleaning up, removing unused variables and variable initialisation
2015-04-21 15:16:13 +00:00
0681f67570
corrected buggy previous commit
2015-04-21 14:33:38 +00:00
bec17127d0
write output information only for active homogenisations. removed unused variables
2015-04-21 14:10:34 +00:00
d049eadce7
write output information only for active crystallites. removed unused variables. fixed some real == 0 comparisons
2015-04-21 12:23:00 +00:00
effec313ee
fixed bug in cleavage system set up for ortho lattice structures
2015-04-21 12:19:41 +00:00
ea1b2191f1
deleted comment, test passed probably because fftw is system-wide available at MPIE
...
also, ../lib is not needed since lib is an include dir via -I
2015-04-19 09:07:16 +00:00
a46b9287eb
restored conditional (#ifdef Spectral) inclusion of “fftw3”. Not clear how this could pass the Marc compile test..?!? (failed in my “real-world” use)
2015-04-18 15:23:37 +00:00
22f9a90bcb
added informative ext_msg to error 135
2015-04-18 15:22:15 +00:00
eb3173c6cc
ugly hack for Abaqus in prec, it seems that it links against old (10.1) ifort where IEEE_arithmetic is not avaialbe.
...
Took the chance and cleaned the env file and used #ifdef# statement in the wrappers because in fixed format Fortran some lines of code might be interpretated as comments.
This should not happen, but it happens for Abaqus exp
2015-04-15 18:10:56 +00:00
caf6611485
unused variable
2015-04-15 18:05:27 +00:00
f5762209dc
changed 'range' keyword to 'limit', now using proper ASCII table (with column 'intensity') for linearODF
...
hybridIO_linODFsampling simplified
2015-04-15 14:37:46 +00:00
e29628b459
changed keyword for texture component/ODF simulation
...
kewyword "origin" with possible values "center/voxelcenter" and "boundary/boundarycenter" is always needed
2015-04-14 12:45:20 +00:00
7e0b79a33c
also using isNaN function
2015-04-14 11:43:31 +00:00
25665b1d80
IEEE_value (like any other function from intrinsic modules) is not allowed for parameters
2015-04-14 05:21:08 +00:00
0c7e7c66f8
integer type was use for void fraction, should be real
2015-04-13 10:36:45 +00:00
5c1804e77d
improved NaN checks
2015-04-13 10:02:52 +00:00
50998bd6a4
simplified makefile and remove unused variables in RGC
2015-04-12 08:11:22 +00:00
7156e3996b
gfortran noticed that these functions are external
2015-04-11 20:28:07 +00:00
0ba8f27320
took some gfortran complaints serious (unused imports, implicit castings)
2015-04-11 11:47:33 +00:00
064266c0cd
updated linearODF file format (including header and keywords)
2015-04-11 10:28:10 +00:00
b4456e76b3
introduced two small bugs
2015-04-11 09:25:15 +00:00
470fc2dce3
comparison for zero should be abs:
...
a == 0.0_pReal becomes abs(a) <= tiny(a)
a /= 0.0_pReal becomes abs(a) > tiny(a)
remove unused variables
2015-04-11 08:25:23 +00:00
6170209198
use statement for marc was wrong (preprocessor)
2015-04-11 05:09:15 +00:00
d7b36c2c30
gfortran complaint about implicit casting,floating point comparison, and unused imported variables.
...
additionally changed use of math_invert33 to math_inv33 if det is not needed
2015-04-10 19:09:26 +00:00
7d0a1ebea6
math_transpose33 not used
2015-04-08 18:33:08 +00:00
b6481c2513
introduced simpler multiplication and forall loops
...
matmul is ok for openmp, check in the web and run the state integration test.
Example program testing for new state update for rkck dot state:
program test
real, dimension(6,10) :: dotState=reshape(&
[1,1,1,1,1,1,1,1,1,1,&
2,2,2,2,2,2,2,2,2,2,&
3,3,3,3,3,3,3,3,3,3,&
4,4,4,4,4,4,4,4,4,4,&
5,5,5,5,5,5,5,5,5,5,&
6,6,6,6,6,6,6,6,6,6],[6,10])
real, dimension(10) :: residuum
real, dimension(6) :: B=2.5
integer :: i
residuum = B(1)*dotState(1,:)+&
B(2)*dotState(2,:)+&
B(3)*dotState(3,:)+&
B(4)*dotState(4,:)+&
B(5)*dotState(5,:)+&
B(6)*dotState(6,:)
do i =1,10
print*,residuum(i)
enddo
residuum = matmul(transpose(dotState),B)
do i =1,10
print*,residuum(i)
enddo
end program test
2015-04-01 16:45:53 +00:00
693efcaa58
comment was confused by search and replace
2015-04-01 16:04:33 +00:00
63491bf268
for some reasons, ifort was complaining about mixed integer types.
...
also deleted #ifdef PETSC since we do not support the solver without PETSC
2015-03-30 20:49:17 +00:00
c7418db9bd
FEM : hierarchical ordering of h5 output data and more meaningful visualisation for multiple grains/crystallites/phases
2015-03-30 09:45:10 +00:00
e0f1132a17
gfortran complaints: equal comparison of reals and external (MPI) functions
2015-03-29 12:54:13 +00:00
9e824d029e
renamed to new names: spectralsolver and spectralfilter
2015-03-28 21:14:11 +00:00
1b2e4c0239
removed old solver
2015-03-28 17:41:40 +00:00
8a30021202
cleaned up core modul
2015-03-27 14:10:18 +00:00
c8554a5641
corrected alignment of restart output byte offset. restart test working now
2015-03-27 07:17:24 +00:00
1ee81e74ea
more patches to get the tests running again:
...
- fixed increment counting in postResults to current output scheme
- corrected node coordinates calculation
- some work on restart (restart test still not running)
2015-03-26 21:19:28 +00:00
86283359f1
PetscFinalize should be the last call before exiting
2015-03-26 13:43:18 +00:00
a9cf8e8051
some patches to compile damask core module without fftw-mpi
2015-03-26 13:41:19 +00:00
2417877bdb
dropped non petsc basic spectral solver
2015-03-25 16:08:41 +00:00
d44fce4a76
Spectral solver now fully parallel (parallel IO, domain decomposition, FFTs and restart). Working but not extensively tested so please report bugs to me
2015-03-25 16:06:19 +00:00
37a7364a3e
multi physics output now working for all solvers
2015-03-25 16:02:30 +00:00
2d4e7fb8e7
calculate dFi/dS (9x9 matrix inverse) for the analytic jacobian only when needed. should improve performance
2015-03-20 07:41:11 +00:00
2e44a846af
need to keep track of temperature again since we can no longer get this from F_T
2015-03-18 18:03:18 +00:00
c5a8543f17
updated petsc
2015-03-18 17:18:43 +00:00
4a0c1c1717
comments
2015-03-15 15:30:14 +00:00
12cb5e56fe
fixed interpretation 'spectralOut' header. Seem to have been buggy all the time, but manifested itself only after keyword change. Also added ':' after keywords to be more robust
2015-03-15 15:29:38 +00:00
ac3d82493e
changed output keywords form resolution to grid and from dimension to size
2015-03-13 14:45:33 +00:00
e4d70e971d
cp is commercially pure ;)
2015-03-13 13:49:09 +00:00
d8debc3396
updated reference geom files (used old names), updated numerics.config to reflect use of analytic jaco
2015-03-13 06:47:26 +00:00
aa125bac5b
aliases for grid and size not supported any more
2015-03-12 22:30:24 +00:00
00cba25a44
improved update of ip coordinates for spectral solver.
...
do not need to create a new fftw plan all the time, using data already defined for the convolution
2015-03-12 22:28:33 +00:00
263f997cf2
seems to be an old output
2015-03-12 08:51:56 +00:00
b338f1aeb8
added some outputs
2015-03-09 23:10:36 +00:00
6b0a2ed5a4
Changed initialization of projection matrix for transformation based on the number of slip systems
2015-03-09 12:05:38 +00:00
a5ee02ff01
added some aliases for output to finally converge from 'shear_rate' to 'shearrate' and from 'accumulatedshear' to 'accumulated_shear'
2015-03-09 11:29:44 +00:00
2ab90c28a5
replaced example ODF by reasonable data in fully correct format
2015-03-06 15:56:58 +00:00
35f9e91e73
some corrections to the damage models
2015-03-06 13:12:32 +00:00
bbb5ff6ae9
changes related to intermediate configuration kinematics:
...
- switched Fi and Li from state variables to crystallite variables
- Lp and Li are now work conjugate with the corresponding mandel stresses defined in their respective configuration
- T, Lp and Li need to return tangents wrt Fi arising from the convection of the material frame due to Fi
- Updated analytic jacobian to take into account tangents wrt Fi
- Updated Lp and Li residual jacobians to take into account tangents wrt Fi
2015-03-06 13:09:00 +00:00
4fd7e0165e
more meaningful name
2015-03-05 11:07:55 +00:00
36f14a2d1d
Removed nr = ns condition
2015-02-26 14:16:29 +00:00
2baac6fc10
Fix for 3881 because gfortran was complaining
2015-02-24 10:15:07 +00:00
3c2a95efd8
if matrix inversion error encountered during analytic tangent calculation, give a warning and return elastic tangent
2015-01-30 14:00:32 +00:00
f1f4f939a3
analytic jacobian default for all solvers
2015-01-29 13:59:49 +00:00
dd8458a775
updated analytic jacobian calculation to correctly take into account intermediate configuration Fi. improved convergence of Li loop in stress integration
2015-01-29 13:58:25 +00:00
0b59519a2a
updated damage models:
...
- coupling to plasticity handled within damage module instead of plasticity module
- anisotropic models more stable
2015-01-29 13:56:09 +00:00
12c586afb9
removed unused variables
2015-01-23 13:44:16 +00:00
1f666dffc9
fixed bug: output size was not set
2015-01-23 13:08:25 +00:00
e9c3137a3e
disabled a faulty ntrans loop (matmul was complaining)
2015-01-23 11:23:35 +00:00
43095b1f1b
added threshold_stress_slip as output (got lost once upon a time)
2015-01-22 09:34:11 +00:00
1c2dd37e28
deleted unused variables in dislo UCLA
2015-01-21 15:14:00 +00:00
e57cb2b4f4
Wrong number of slip families in BCC material config files
2015-01-20 17:22:52 +00:00
19dc945065
Inroduced Fd in kinematic decomposition.
2015-01-19 13:39:47 +00:00
54dc4943de
adaptive penalty to bound Cv in (0,1)
2015-01-16 18:31:14 +00:00
f16244d240
corrected state integration
2015-01-16 17:51:43 +00:00
eecb78665c
removing some unused code and cleaning up
2015-01-16 17:32:47 +00:00
6411d36633
updated vacancy-damage model
2015-01-16 17:04:01 +00:00
93e50366bb
fixed automatic determination of compiler
2015-01-16 11:19:30 +00:00
701f35c326
reverting the changes related to undamaged_Fi during rev. 3850/3851.
2015-01-16 00:09:22 +00:00
85465333ec
2015-01-15 22:42:56 +00:00
b1ef2f9dd9
2015-01-15 22:40:17 +00:00
62ff5d1d45
gfortran/ifort is now selected depending on the compiler provided by PETSc
2015-01-15 11:35:53 +00:00
ced00954fe
added new model by david (LLNL/UCLA)
2015-01-15 10:56:15 +00:00
125276283b
fix in analytically integrated expression.
2015-01-15 09:22:24 +00:00
19611660c2
passing undamaged Fe for calculation of stress.(Pratheek please check.)
2015-01-13 09:54:44 +00:00
9c8873f16c
added functions to get degraded stiffness and undamaged intermediate deformation gradient.
2015-01-13 09:46:25 +00:00
80bb9a2e66
function to calculate mises stress.
2015-01-13 09:43:05 +00:00
c2b8627981
Currently, there is no updating for the dot state of vacancy generation, set sizeDotState = 0
2015-01-11 20:01:43 +00:00
a2552d6241
Polishing.
2015-01-11 19:10:16 +00:00
e20271b89a
fix for temperature in spectral solver
2015-01-08 13:26:00 +00:00
9635ae9f3d
fix in integration.
2015-01-07 14:51:06 +00:00
8610531aae
removed unused variable.
2015-01-04 20:13:16 +00:00
460c0d73fa
Making use of extended plastic state with pointers, removed helper functions.
2015-01-04 19:26:33 +00:00
e6d1b21b57
introduced extended plastic state type (including pointer/aliases)
2014-12-19 14:01:03 +00:00
e4cb387571
polishing comments
2014-12-18 18:44:47 +00:00
72cb9072e1
added some debug statements
2014-12-18 18:41:02 +00:00
7fe20fe805
intent(out) specifier for getXX functions
2014-12-18 18:26:16 +00:00
fdbaf00fa7
cleaned up naming and added elastic material with volume preserving stiffness
2014-12-18 18:20:11 +00:00
299518adfa
cleaned tungsten parameter set, renamed to be consistent
2014-12-18 16:55:43 +00:00
3f14ebe43d
helper routine to get heat generation rate needed for MARC/Abaqus
2014-12-17 13:37:13 +00:00
e6e72ab3a7
In mesh_marc_build_elements use mesh_NcpElems instead of mesh_Nelems (line 2906).
...
Fixes error when not all elements are CP elements.
Some formating and comment improvements.
2014-12-15 11:51:32 +00:00
7c8f2ede31
fixed missing initialization causing a NaN in gdot and crashes
2014-12-11 15:04:16 +00:00
5bc93efd63
Fix bug about vacancy generation in subroutine constitutive_microstructure
2014-12-11 13:07:15 +00:00
bdb6bf8090
recasting vacancy dot state as microstructure state
2014-12-10 13:40:30 +00:00
39a3a5b5a0
cleaned up FEM examples
2014-12-09 18:59:37 +00:00
3064c14ff6
fig a bug of energy barrier of vacancy generation.
2014-12-09 18:56:19 +00:00
3f8678c7c4
plastic dot states evaluated at effective stress not damaged stress
2014-12-09 18:23:48 +00:00
55e2de6ffd
made dotstate and state pointers (targets are not allowed in derived types) and introduced plastic state having aliases to parts of the state.
2014-12-09 12:12:53 +00:00
d1e1191df4
corrected bug: renamed file list for commercial FEM solvers to start with small letter
2014-12-08 16:36:22 +00:00
2458aa9c3e
BBar stabilisation now optional (default is off)
2014-12-08 16:27:23 +00:00
c34252c71c
changed constitutive_* prefix for plasticity laws to plastic_* to be consistent with novel thermal_, damage_, vacancy_, etc. laws.
2014-12-08 15:55:30 +00:00
0bd4ab48f8
numerics with respect to vacancy diffusion.
2014-12-08 15:37:44 +00:00
db0e02b94d
changed to lower canse convention
2014-12-08 15:37:02 +00:00
67c3eb3a2e
bounding the analytically integrated damage variable
2014-12-08 14:05:38 +00:00
347dac74c6
added vacancy state integrators
2014-12-08 10:45:12 +00:00
9301b56aa3
initialize restoration points of vacancy concentration at material points, postResults of VC at material points.
2014-12-08 09:27:32 +00:00
62f254c9de
age results and write restart data for vacancy concentration
2014-12-08 09:09:54 +00:00
55c2b729a0
while node coordinates are always 3D displacements are passed into hypela2 as 2D array (dispt) for 2D simulations
2014-12-05 07:35:28 +00:00
334ade5260
Every line in a recipe must begin with a tab character.
...
http://stackoverflow.com/questions/14576237/errors-in-make-file-missing-separator-stop
2014-12-03 13:37:04 +00:00
9d12e026e8
renumbered some errors, fixed non-existing errors
2014-12-03 00:42:35 +00:00
01a8c500ca
detabbing…
2014-12-03 00:41:12 +00:00
fa6f784102
deleted unused stuff
2014-12-02 21:49:33 +00:00
82942d1618
bug fix for previous commit
2014-12-02 17:23:47 +00:00
b796e46f98
residualStiffness is not just an FEM parameter
2014-12-02 17:17:35 +00:00
2fc264d90a
forgot file in previous commit
2014-11-28 20:07:22 +00:00
37dd574790
recast power law and viscous dot state variables as derived quantities using analytical solutions of the evolution equation
2014-11-28 20:06:24 +00:00
8d64e513da
improved commenting and variable names.
...
values for unused slip/twin families are ignored in material.config (do not need to be numeric anymore…)
2014-11-26 21:23:02 +00:00
f8b29bdaf0
fixed bug in type definition for ‘drivingForce’
2014-11-26 18:40:31 +00:00
54de033c45
using analytic solution for local/nonlocal coupled damage (as opposed to having crystallite iterate to find a consistent solution)
2014-11-26 17:29:45 +00:00
61bb9a373e
all damage models consistently semi-explicit
2014-11-26 11:32:50 +00:00
6e90260817
updated local damage post results for newer models
2014-11-26 11:32:19 +00:00
a91ebafa0e
semi-explicit damage evolution
2014-11-26 10:36:14 +00:00
646f1c48dc
typo, IO_warning and not IO_error
2014-11-26 10:16:30 +00:00
9714a172d1
updated to get local damage from microstructure state
2014-11-26 00:42:20 +00:00
c89ef9de61
relaxed some sanity checks
2014-11-26 00:29:54 +00:00
26c01f64e5
updated file list
2014-11-26 00:06:16 +00:00
32ee3a4a66
corrected constitutive_getVacancyPotentialDrivingForce function
2014-11-26 00:04:06 +00:00
36dc59b09f
reverted back to rate independent damage evolution. updated example config files to recent changes
2014-11-25 23:56:52 +00:00
a4584fecd2
potential driving force for conservative cahn hilliard vacancy diffusion
2014-11-25 20:13:33 +00:00
20437b2ae0
added brittle(elastic energy release rate)/ductile(vacancy condensation) phase field damage model coupled to vacancy concentration
...
simplified vacancy_generation
2014-11-25 17:23:37 +00:00
65b22f9a40
calculate analytic jacobian from xxx0 and not xxxsub0
2014-11-25 17:20:52 +00:00
e6edc3d8cd
Some sanity checks
2014-11-25 13:20:12 +00:00
296146490f
changed some default parameters for FEM
2014-11-21 14:20:54 +00:00
3d4b2d2f52
tolerance for damage can be any positive value. Set to 0.001 if found negative.
2014-11-21 13:12:54 +00:00
ee1d409bad
critical shear can be set to zero for family not in use.
2014-11-21 09:03:54 +00:00
9e6ab9866f
Read in TRIP inputs, changed 'totalvolfrac' to 'totalvolfrac_twin'
2014-11-21 08:54:20 +00:00
36c2cc4347
added some sanity checks in damage modules.
2014-11-20 13:48:07 +00:00
f538469892
reverted name change: In current PETSc release, its DM_ instead of DMDA_!
2014-11-19 05:35:10 +00:00
73a13ddffc
added more comments
2014-11-18 15:36:17 +00:00
b3a96d40c5
Started adding TRIP model to phenopowerlaw
2014-11-18 10:46:02 +00:00
23d68114e3
cleaned up stress integrations and some documentation
2014-11-13 12:53:20 +00:00
b407e6ac80
complement of pore nucleation.
2014-11-12 17:55:02 +00:00
fc9f4835c3
switched order of Lp-Li nesting in stress integration loop for better convergence. temperature is integrated during stress integration so it need not be a dot state
2014-11-12 16:40:50 +00:00
e89199d119
correction to thermal velocity gradient
2014-11-11 23:14:45 +00:00
7567582781
updated constitutive in line with changes to damage
2014-11-11 22:46:41 +00:00
5ff3c882b5
thermodynamically consistent treatment of nonlocality over all damage models. switched from power law to exponential with viscous drag type rate formulation for damage evolution
2014-11-11 22:44:11 +00:00
f0f04a25bf
Pore nucleation, and the decrease of vacancy concentration due the formation of pores.
2014-11-11 16:03:06 +00:00
e0e0f1503b
svn properties for new files
2014-11-07 20:55:16 +00:00
9544fe55a6
FEM solver for vacancy diffusion
2014-11-07 20:40:54 +00:00
eb1028eb6f
dislokmc coupled to ductile damage
2014-11-07 12:28:18 +00:00
646f632977
no need to pass fileunit during init of none types
2014-11-07 12:15:28 +00:00
c4e54793fc
missing factor of 2 in damage driving force
2014-11-07 12:10:10 +00:00
dc0745800f
reworked anisotropic brittle damage
2014-11-07 12:08:56 +00:00
69b150319f
fully introduced non-schmid behavior
2014-11-07 11:23:34 +00:00
eee863956f
added flush(6) for debugging
2014-11-06 21:05:35 +00:00
4fcd1821cf
reworked anisotropic ductile model
2014-11-06 17:52:43 +00:00
9c474796ac
cleaning up and meaningful renaming
2014-11-06 17:41:34 +00:00
a72a9ed8ec
Add the generation rate of vacancy induced by dislocation movement.
2014-11-06 16:53:34 +00:00
29951a9e1a
named 99 representation according to other constitutive laws
2014-11-06 11:49:37 +00:00
b443ac6321
nicer preprocessor statements
2014-11-06 11:47:27 +00:00
95ec8b7b7b
removed imported module variables that are not longer needed, commented loops and simplified calculation for rates, i.e. do not store per system anymore but only once
2014-11-06 08:27:48 +00:00
1242aab2bf
unified calculation for slip rate for output to avoid doubling of code
...
introduced calculation of negative and positive stress (three times in the code) for non-schmid behavior.
Still, non-schmid behavior is not active at all
2014-11-05 20:53:25 +00:00
6fc3908c71
prepared for merging postResults rate calculation in one loop
2014-11-05 19:39:23 +00:00
f3b7b5bb96
corrected substraction by volume fraction of twin for LpAndItsTangent and removed the substraction in postResults (like for dislowtin)
...
cleaned up and unified notation in calculation of slip rates
2014-11-05 19:11:09 +00:00
accb571c53
corrected substraction by sum of twinned and transformed value fraction of Lp and derivative.
...
removed unused variables
2014-11-05 17:52:49 +00:00
a61c3059ef
missing file in previous commit
2014-11-05 17:47:46 +00:00
f2211260e3
introduced possibly nonlocal anisotropic ductile damage model
2014-11-05 17:41:08 +00:00
94ba7e6246
added accumulated shear as a dot state for cleaner integration with ductile damage models
2014-11-05 17:07:31 +00:00
ade918246c
reduced number of state variables. some cleaning up
2014-11-05 17:05:40 +00:00
a21f4d80db
removed some leftovers of the shearbanding part
2014-11-05 09:53:08 +00:00
c2799585d9
introduced non schmid tensor calculation in dot state, also not active yet
2014-11-05 09:35:49 +00:00
c66fc95dda
made a bunch of arrays simple scalars to improve readability, commented loops
2014-11-05 09:24:08 +00:00
816169fbec
continued work on integration of non schmid behavior. not active yet
2014-11-05 09:03:04 +00:00
18b49fbac8
line in OMP statement too long
2014-11-04 21:32:16 +00:00
e36880e3dc
fixed typo 'orthorhombic'
2014-11-03 15:41:05 +00:00
76442942b1
damage def grad updated in stress integration loop and not as a dot state
2014-11-03 10:43:36 +00:00
0ad917dceb
nicer FEM parameters
2014-11-03 10:36:07 +00:00
261bffa745
switched to smoother traction separation law for better convergence. anisotropic brittle now working.
...
lowered default max iterations for FEM
2014-11-01 17:51:35 +00:00
9c94eb4be6
OMP bug fix from previous commit
2014-11-01 11:43:43 +00:00
f9cec3079e
overhaul of kinematics from FeFp to FeFiFp decomposition
2014-10-31 19:03:08 +00:00
26a609dbee
nicer permutation of cleavage planes
2014-10-31 18:37:37 +00:00
629ff8bdef
reverted back to r 3644 to get rid of buggy behavior
2014-10-31 10:39:02 +00:00
ad4bd86fd2
Polishing
...
Modify
texture: do i = 1_pInt,myNconstituents
grain = sum(NgrainsOfConstituent(1_pInt:i-1_pInt))
as:
grain = 0_pInt
texture: do i = 1_pInt,myNconstituents
grain = grain + NgrainsOfConstituent(i-1_pInt)
First try of developing DAMASK, :)
2014-10-30 09:51:06 +00:00
b0469854c8
brittle anisotropic damage on cleavage planes
2014-10-28 20:57:12 +00:00
a6f88c0e37
introduced crystallographic cleavage systems for each lattice type
...
FCC: (100) and (111) planes
BCC: (100) and (110) planes
HCP: (0001) planes
iso/orthotropic: (100) planes
2014-10-28 18:05:51 +00:00
Martin Diehl