Martin Diehl
7ebd5e29e3
line too long
2016-01-17 19:42:24 +00:00
Martin Diehl
12653bf1a4
more rename g -> c/ipc, wrong ipc iteration for crystallite_integrateStateRK4 and single run fixed
2016-01-17 17:56:24 +00:00
Martin Diehl
b10774aae9
unused variable
2016-01-17 15:06:35 +00:00
Martin Diehl
cfd6579024
split CPFEM module into one for spectral solver (also suitable for PETSc FEM) and one for commercial FEM codes as all the ping pong handling, cut back detection etc. is not needed for spectral.
...
Needs certainly a better name, or might be superfluous at all if crystallite, constitutive, and homogenization take care of their data for restart write/read and forwarding
2016-01-17 15:03:54 +00:00
Martin Diehl
0799570a03
renamed g -> c/ipc
2016-01-17 14:50:33 +00:00
Martin Diehl
1caec2e608
calc done is used only internally
2016-01-17 14:15:38 +00:00
Martin Diehl
3d2187cc07
moved some variables from FEsolving to more appropriate CPFEM module
2016-01-17 13:29:42 +00:00
Martin Diehl
ff840ae8ad
remove superfluous code and commented
2016-01-17 11:14:06 +00:00
Martin Diehl
11df75dfb2
added documentation and changed some names
2016-01-16 17:27:19 +00:00
Martin Diehl
83a6d00c60
remove superfluous variables
2016-01-16 07:06:34 +00:00
Martin Diehl
208196d213
introduced pointers as names. thx to Chuanlai for doing the work
2016-01-15 15:26:24 +00:00
Martin Diehl
117c8edfbd
corrected comment
2016-01-15 14:03:24 +00:00
Philip Eisenlohr
684992bf50
exchanged mappingConstitutive array for more descriptive "phaseAt and phasememberAt"
2016-01-15 00:19:44 +00:00
Martin Diehl
4b10e4792e
adjusted function for polar decomposition to actual need (no return of U)
2016-01-12 21:39:31 +00:00
Martin Diehl
f090a1b216
generalized dimension of calls to lapack, for performance reason special 3x3 variants will follow
2016-01-12 11:00:23 +00:00
Martin Diehl
2eafefe652
removed norm functions from math in favor of intrinsic and simplified (mostly by using existing functions, merge intrinsic and array constructors)
2016-01-10 13:34:26 +00:00
Martin Diehl
519cd29c6f
using norm2 function as L2 norm (vector, math_norm3) and Frobenius norm (tensor, math_norm33)
2016-01-09 16:01:30 +00:00
Tias Maiti
53d07215a9
reverting back to converged Lp as the initial guess for next iteration
2016-01-09 15:03:18 +00:00
Martin Diehl
885b9095c4
option not supported by slightly older gfortran versions. thx to Luv for pointing out
2016-01-09 12:33:33 +00:00
Martin Diehl
5378dbdcc9
cleaned
2016-01-09 12:12:31 +00:00
Martin Diehl
b5ab94a0f1
added new plasticity module
2016-01-09 12:12:05 +00:00
Tias Maiti
7ba031adb3
fixed wrong Li and dLi_dTstar_3333 initialization when Tstar is zero.
2016-01-09 03:41:56 +00:00
Tias Maiti
a30b596e7f
combined J2 and I2 isotropic mechanical response. Default is identical to old plastic_j2. /dilatation/ flag adds (spherical) Li calculation based on hydrostatic pressure.
2016-01-08 19:45:20 +00:00
Tias Maiti
6cf92913a3
use current Lp for next iteration instead of converged Lp inherited from end of former increment as starting guess
2016-01-08 19:01:53 +00:00
Tias Maiti
1d1672ebfc
added math_spherical33 to calculate hydrostatic part of a 2nd order tensor
2016-01-08 18:57:37 +00:00
Martin Diehl
cbb6b7a689
adding fpe report for gfortran, some compiler opions are now used for linking (Debug + Optimization)
2016-01-07 12:58:19 +00:00
Martin Diehl
1857e47f75
div by zero in output
2016-01-07 11:48:30 +00:00
Martin Diehl
86885cb2ba
wrong initialization order, thanks to N Grilli for reporting the bug
2016-01-06 16:46:37 +00:00
Martin Diehl
268c7f359c
wrong index in output
2015-12-23 21:15:49 +00:00
Martin Diehl
74086fd0fb
correct initialization
2015-12-23 16:13:08 +00:00
Martin Diehl
7127403bd9
usused variables
2015-12-23 11:25:25 +00:00
Martin Diehl
dfe09f11f5
cleaner handling of NaN
2015-12-22 10:03:15 +00:00
Martin Diehl
ee47c9302a
small readability improvement
2015-12-22 09:22:57 +00:00
Martin Diehl
e127bd8db4
should reflect name of parameter
2015-12-21 07:27:21 +00:00
Chen Zhang
2a2f558416
update physics in phenoplus module
2015-12-17 16:11:13 +00:00
Martin Diehl
5f37ade722
followed some suggestions from gfortran
2015-12-16 22:42:08 +00:00
Martin Diehl
5901af1212
removing code fragments after introducing pointer state
2015-12-16 06:15:21 +00:00
Martin Diehl
9baec1a0ff
get rid of overly long file names 1/2
2015-12-15 20:45:54 +00:00
Martin Diehl
c710709b1a
example Makefile typically does not work, deleted
2015-12-15 14:34:03 +00:00
Martin Diehl
37aa559bd4
fixed missing structure type in material.config, small changes in the other code
2015-12-15 14:30:17 +00:00
Martin Diehl
88999c847a
reverted handling of NaN in converged reason, seems to be a compiler bug
2015-12-14 20:04:59 +00:00
Martin Diehl
dc990a88d6
unused variables
2015-12-14 19:48:11 +00:00
Pratheek Shanthraj
3a7f4bf43a
added new set of accelerated spectral methods based on derivative approximations. use spectral_derivative to set the kind of derivative you like. valid options are:
...
‘continuous’: computed continuous derivatives (standard spectral method)
‘central_difference’: central difference approximation of derivatives
‘fwbw_difference’: forward difference for gradient and backward difference divergence
2015-12-14 18:12:09 +00:00
Martin Diehl
da913ec3c7
bug in dot state, thx to Haiming for pointing out
2015-12-11 18:36:58 +00:00
Martin Diehl
76fe37e496
commented and respected compiler warning
2015-12-08 18:10:06 +00:00
Martin Diehl
4174e6f9e2
removed dislokmc, use disloucla for tungsten simulations
2015-12-05 22:07:03 +00:00
Martin Diehl
6ecbd3ae56
makefile had problems with static flag, mentat spectral Box was still using table.croak
2015-12-02 20:02:19 +00:00
Martin Diehl
0e76ef85c1
checking return code in util.exececute, cleaning makefile, deleted dummy docs
2015-12-02 13:21:17 +00:00
Martin Diehl
22a9d65825
spectral solver accepting NaN in SNES converged reason as this might happen for homogeneous microstructures
2015-12-01 22:36:19 +00:00
Martin Diehl
3365edc6cb
simplified and updated abaqus
2015-12-01 15:13:33 +00:00
Martin Diehl
4eb43908c3
some simplifications as discussed with IBF colleagues
2015-11-25 20:55:17 +00:00
Su Leen Wong
a6cc799375
Removed unnecessary variables
2015-11-17 16:55:23 +00:00
Su Leen Wong
51059abaf0
Martensite volume fraction evolution
2015-11-17 16:30:06 +00:00
Su Leen Wong
a8b157a87c
Did not initialize a variable
2015-11-17 15:20:36 +00:00
Su Leen Wong
4a8477c2c3
Characteristic transformation shear
2015-11-17 14:50:59 +00:00
Su Leen Wong
cf84ceddd8
Added Ndot0 for trans
2015-11-16 10:15:18 +00:00
Martin Diehl
302178ebb4
some trouble since clean was already in PETSc rules
2015-11-14 09:47:18 +00:00
Martin Diehl
9d9367ae11
still having trouble, incorporation of PETSC_ARCH not needed
2015-11-13 10:16:35 +00:00
Martin Diehl
4a52c59854
corrected/removed PETSc variables path
2015-11-12 23:13:25 +00:00
Philip Eisenlohr
900b072781
removed unused math_qDot (deprecated anyway)
2015-11-10 21:56:06 +00:00
Philip Eisenlohr
45fb209658
respect PETSC_ARCH
2015-11-10 21:18:53 +00:00
Philip Eisenlohr
c5b28a22d4
small polishing
2015-11-10 21:18:13 +00:00
Su Leen Wong
aefd8fa390
Added a few more trans dependent variables
2015-11-10 15:23:15 +00:00
Su Leen Wong
2778331d7b
Added trans interaction matrices
2015-11-10 13:30:34 +00:00
Su Leen Wong
4e4136cd58
Consistent variable names
2015-11-10 13:01:03 +00:00
Su Leen Wong
9c91ba4608
Trans interaction matrices
2015-11-09 15:57:20 +00:00
Su Leen Wong
2b10625f0c
Typo
2015-11-09 08:51:05 +00:00
Martin Diehl
ff7cd17b51
not needed, PETSc is doing the job
2015-11-08 09:43:43 +00:00
Chuanlai Liu
4d9a2f8f6b
introduce pointer referencing of plastic state for better readability
2015-11-06 17:00:00 +00:00
Martin Diehl
7b0c130d6f
started to introduce alias name via pointer for plastic state in phenopowerlaw for more comfortable access
2015-10-30 15:48:30 +00:00
Su Leen Wong
2ae4d24fb3
Use rotated trans elasticity matrices
2015-10-27 12:32:03 +00:00
Martin Diehl
2ba15e462a
corrected values and typo
2015-10-26 19:27:57 +00:00
Martin Diehl
c91788d721
new intel compiler was compiling about *** in float output
2015-10-26 17:11:36 +00:00
Martin Diehl
1f356a6833
added parameters from indentation results
2015-10-22 17:45:28 +00:00
Su Leen Wong
04841cb727
List of dependent trans states
2015-10-21 13:24:26 +00:00
Su Leen Wong
20b0429517
Calculation of equilibrium separation of partials for twin and trans
2015-10-21 09:37:45 +00:00
Su Leen Wong
565d68edd7
Added some TRIP inputs
2015-10-20 15:54:01 +00:00
Su Leen Wong
b9fcd9482b
Renamed some input variables to denote twin
2015-10-20 13:45:06 +00:00
Franz Roters
fc8a11f027
corrected typo: alias for first Piola stress should read 1stpiola NOT 1piola
2015-10-19 06:57:23 +00:00
Chen Zhang
f2b29e4844
correct communication between
...
elements with different
phases
2015-10-14 21:04:40 +00:00
Chen Zhang
484a34b7f1
added pheno+ module
...
modify crystallite microstructure call
to pass orientations
2015-10-14 18:36:19 +00:00
Martin Diehl
ea0b16c726
pheno plus for commercial FEM
2015-10-13 21:06:59 +00:00
Chen Zhang
d6abc00218
add pheno+ module in
2015-10-13 18:52:01 +00:00
Martin Diehl
8fac635c15
made hybridIA stuff working again
2015-10-13 17:02:07 +00:00
Martin Diehl
6de8c821f8
improved consistency check on loadcase
2015-10-12 19:02:42 +00:00
Luv Sharma
9ec7c08bb6
gfortran was complaining abound mismatch
2015-10-07 13:07:21 +00:00
Martin Diehl
c1a53bc8c0
introduced dummy target to check compilation
2015-10-06 18:00:17 +00:00
Martin Diehl
0f0b4eaa6b
cleaned
2015-09-28 17:47:00 +00:00
Martin Diehl
ad16162eb8
simplified
2015-09-24 17:50:11 +00:00
Martin Diehl
c50a522dd8
some shapes not 100% correct, removed inverse laplace
2015-09-24 17:38:49 +00:00
Martin Diehl
2e70e6fded
make clear that this isNaN is an intrinsic extension
2015-09-24 08:45:44 +00:00
Aritra Chakraborty
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
Martin Diehl
42ac7902f7
indents wrong in fortran, seeds_fromDistribution bugfix and using ASCII table facilities
2015-09-20 21:03:10 +00:00
Martin Diehl
41f09dd0f8
should make no difference, but test with ifort was failing. compierler buig?
2015-09-12 18:26:25 +00:00
Tias Maiti
03896fcd09
redefined required debug variables removed in last commit
2015-09-11 18:20:37 +00:00
Tias Maiti
840fd0e9fa
report absolute value of accumulated_shear_slip in dot state
2015-09-11 17:44:26 +00:00
Martin Diehl
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
Martin Diehl
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
Martin Diehl
e88cedc6ae
clearer naming, debug options for spectral do not work for MPI
2015-09-11 08:52:03 +00:00
Martin Diehl
751d1d7582
function statement is more clear
2015-09-10 13:01:33 +00:00
Martin Diehl
580bb76b4c
pInt missing
2015-09-10 10:26:09 +00:00
Martin Diehl
427b5f4bc1
wrong results for restart
2015-09-10 10:06:14 +00:00
Martin Diehl
bbe37c842e
name conflict
2015-09-09 21:52:00 +00:00
Martin Diehl
34980a1d44
more logical structure when reading in
2015-09-08 19:44:32 +00:00
Martin Diehl
7d996cfa11
fixed new IO_stringPos for Abaqus
2015-09-05 16:26:55 +00:00
Martin Diehl
caf0ac7e8d
more verbose
2015-09-01 16:53:48 +00:00
Martin Diehl
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
Pratheek Shanthraj
197ae53553
lots of memory savings (~ stateside + 250 pReals per integration point) when using analytic jacobian
2015-08-28 10:55:38 +00:00
Martin Diehl
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
Martin Diehl
1171dc4344
MPI has a 2GB limit for writing at once, now chunking
2015-08-21 17:51:05 +00:00
Martin Diehl
8474da26d9
empty line
2015-08-18 16:31:47 +00:00
Martin Diehl
bae5cfcf89
named if
2015-08-18 16:07:01 +00:00
Tias Maiti
53c574033a
size of MAXCHUNKS same as maximum allocated size for “heat_time” & “heat_rate”
2015-08-18 03:20:16 +00:00
Martin Diehl
f64a5e7552
output safe for 64bit integer
2015-08-14 05:03:54 +00:00
Martin Diehl
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
Martin Diehl
563d9e64dd
MAXNCHUNKS capitalized since it is a parameter
2015-08-06 09:24:56 +00:00
Pratheek Shanthraj
3d94ae4074
finite strain correction for respective expansions
2015-08-05 14:12:54 +00:00
Tias Maiti
09c7affbbc
more efficient rate calculation in updateState when subdt is zero
2015-08-05 11:09:38 +00:00
Tias Maiti
7928ce3ac3
fixed erroneous rate calculation in updateState when time increment (subdt) was zero during initialization.
2015-08-04 21:27:30 +00:00
Tias Maiti
65d114e4f0
improved debug statements
2015-08-04 21:26:22 +00:00
Pratheek Shanthraj
905b9f8f42
petsc update (now v3.6.1)
2015-08-04 15:04:53 +00:00
Pratheek Shanthraj
5b57ea7d16
fluctuations proportional to thermal equilibrium concentration
2015-08-03 11:08:59 +00:00
Pratheek Shanthraj
a01d7d8d66
corrections to crystallite_push33ToRef function
2015-08-03 11:07:19 +00:00
Pratheek Shanthraj
1ed276ce2d
quadratic softening for finite dissipation energy
2015-07-30 10:29:25 +00:00
Pratheek Shanthraj
c97b83c9a1
added user defined external heat source
2015-07-27 11:09:37 +00:00
Pratheek Shanthraj
4bb6664de9
tightened default thermal tolerance
2015-07-27 11:08:43 +00:00
Pratheek Shanthraj
c82c4b74f1
updated config files to reflect recent changes in input parameters
2015-07-24 17:14:03 +00:00
Pratheek Shanthraj
5e09954575
remove deprecated spectral load case definition of temperature
2015-07-24 14:57:29 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
87d42bf447
initialise Fi correctly for initial field values away from equilibrium
2015-07-24 14:47:18 +00:00
Pratheek Shanthraj
0a5ccb3d91
possibility to set damage potential order
2015-07-24 14:43:05 +00:00
Tias Maiti
7c101cdc31
fixed issues causing compilation errors
2015-07-23 19:22:03 +00:00
Tias Maiti
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
Philip Eisenlohr
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
Su Leen Wong
d2c748e716
TRIP input parameters for fcc to hex
2015-07-10 09:02:48 +00:00
Luv Sharma
1e9682eef9
corrected upper bound for spectral damage solver as well.
2015-07-09 15:51:34 +00:00
Shaokang Zhang
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
Su Leen Wong
7a4415b28a
Sanity checks for SFE
2015-07-08 15:37:16 +00:00
Su Leen Wong
df5cdb6b01
Elasticity matrices for transformed phase
2015-07-08 11:58:52 +00:00
Su Leen Wong
66564c1f75
Read in elastic constants for transformed phase
2015-07-07 13:23:56 +00:00
Luv Sharma
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
Luv Sharma
3d5a85e4af
corrected implementation of deltaState for isoBrittle damage
2015-07-06 08:47:40 +00:00
Martin Diehl
72b7df34a3
save some space
2015-07-02 09:45:25 +00:00
Aritra Chakraborty
74af77a81f
formatted the interaction matrix for beta-Sn
2015-07-02 05:45:56 +00:00
Aritra Chakraborty
0e6ac411bb
further modified to add longer lines for (bct) beta-Sn in lattice
2015-07-02 05:44:45 +00:00
Martin Diehl
6272129f01
formatted code and modified Makefile to have long lines in lattice
2015-07-02 05:09:07 +00:00
Aritra Chakraborty
e61050dcde
added lattice for (bct) beta- Sn.
2015-06-27 14:55:30 +00:00
Shaokang Zhang
a40d247f7f
lower bound for damage is residual stiffness
2015-06-22 14:27:32 +00:00
Su Leen Wong
12bc4f01ef
fcc to hex orientation relationship (Shoji-Nishiyama)
2015-06-22 13:47:15 +00:00
Su Leen Wong
ad17c40d92
Renamed variables to denote fcc to bcc transformation
2015-06-22 08:33:48 +00:00
Su Leen Wong
5114e0b43e
Cleaning up
2015-06-22 08:10:20 +00:00
Philip Eisenlohr
3e0e187b14
clarified meaning of symEulers function in help
2015-06-19 07:29:01 +00:00
Franz Roters
b1c036daab
added some comments on possible code changes
2015-06-17 06:53:00 +00:00
Franz Roters
0e48fb8fc7
removed double entry in use list
2015-06-16 07:57:42 +00:00
Su Leen Wong
081bc33159
Read in c/a ratio for transformed hcp phase
2015-06-12 08:59:39 +00:00
Su Leen Wong
a370372c1b
transformed lattice structure should not be mandatory
2015-06-11 17:03:34 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
e8ee5d6723
moved some phase field parameters to lattice
2015-06-11 09:01:37 +00:00
Su Leen Wong
b3241411f5
Added keyword to specify transformed lattice structure
2015-06-11 08:23:27 +00:00
Su Leen Wong
c2781b5423
Removed unnecessary variables
2015-06-10 12:12:03 +00:00
Pratheek Shanthraj
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
Su Leen Wong
b410a08995
Renamed some variables for consistency
2015-06-08 14:26:50 +00:00
Su Leen Wong
0d36a2d3e2
Removed unused variables and renamed some variables for TRIP model
2015-06-08 09:36:52 +00:00
Pratheek Shanthraj
74982294a0
added spectral thermal and damage solvers
2015-06-03 17:30:31 +00:00
Pratheek Shanthraj
3415e249ef
reverting files that got mistakenly added into previous change set
2015-06-01 16:11:37 +00:00
Pratheek Shanthraj
a14070bad4
changed up handling of delta states, and some bug fixes
2015-06-01 16:02:27 +00:00
Pratheek Shanthraj
8d2aa5e7be
fixed typo in dislotwin label
2015-05-29 11:03:25 +00:00
Pratheek Shanthraj
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
Martin Diehl
b58cbe424d
follow up to last commit
2015-05-19 17:41:49 +00:00
Martin Diehl
b07ec8232e
found circular inclusion
2015-05-19 17:34:04 +00:00
Martin Diehl
71b0e283c1
clearer substructure of routines
2015-05-10 20:55:36 +00:00
Martin Diehl
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
Martin Diehl
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
Martin Diehl
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
Martin Diehl
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
Martin Diehl
1f809f9d73
variable name change forgotton
2015-04-22 04:34:35 +00:00
Martin Diehl
78d913c07d
clearified keywords, added raw data of rolling texture, renamed to more meaningful names
2015-04-21 19:34:50 +00:00
Pratheek Shanthraj
d7df5f1934
more cleaning up. added pure attribute to some functions. differentiating between ‘local’ and ‘nonlocal’ temperature.
2015-04-21 16:11:30 +00:00
Pratheek Shanthraj
1330dae432
cleaning up, removing unused variables and variable initialisation
2015-04-21 15:16:13 +00:00
Pratheek Shanthraj
0681f67570
corrected buggy previous commit
2015-04-21 14:33:38 +00:00
Pratheek Shanthraj
bec17127d0
write output information only for active homogenisations. removed unused variables
2015-04-21 14:10:34 +00:00
Pratheek Shanthraj
d049eadce7
write output information only for active crystallites. removed unused variables. fixed some real == 0 comparisons
2015-04-21 12:23:00 +00:00
Pratheek Shanthraj
effec313ee
fixed bug in cleavage system set up for ortho lattice structures
2015-04-21 12:19:41 +00:00
Martin Diehl
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
Philip Eisenlohr
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
Philip Eisenlohr
22f9a90bcb
added informative ext_msg to error 135
2015-04-18 15:22:15 +00:00
Martin Diehl
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
Martin Diehl
caf6611485
unused variable
2015-04-15 18:05:27 +00:00
Martin Diehl
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
Martin Diehl
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
Martin Diehl
7e0b79a33c
also using isNaN function
2015-04-14 11:43:31 +00:00
Martin Diehl
25665b1d80
IEEE_value (like any other function from intrinsic modules) is not allowed for parameters
2015-04-14 05:21:08 +00:00
Martin Diehl
0c7e7c66f8
integer type was use for void fraction, should be real
2015-04-13 10:36:45 +00:00
Martin Diehl
5c1804e77d
improved NaN checks
2015-04-13 10:02:52 +00:00
Martin Diehl
50998bd6a4
simplified makefile and remove unused variables in RGC
2015-04-12 08:11:22 +00:00
Martin Diehl
7156e3996b
gfortran noticed that these functions are external
2015-04-11 20:28:07 +00:00
Martin Diehl
0ba8f27320
took some gfortran complaints serious (unused imports, implicit castings)
2015-04-11 11:47:33 +00:00
Martin Diehl
064266c0cd
updated linearODF file format (including header and keywords)
2015-04-11 10:28:10 +00:00
Martin Diehl
b4456e76b3
introduced two small bugs
2015-04-11 09:25:15 +00:00
Martin Diehl
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
Martin Diehl
6170209198
use statement for marc was wrong (preprocessor)
2015-04-11 05:09:15 +00:00
Martin Diehl
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
Martin Diehl
7d0a1ebea6
math_transpose33 not used
2015-04-08 18:33:08 +00:00
Martin Diehl
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
Martin Diehl
693efcaa58
comment was confused by search and replace
2015-04-01 16:04:33 +00:00
Martin Diehl
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
Pratheek Shanthraj
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
Martin Diehl
e0f1132a17
gfortran complaints: equal comparison of reals and external (MPI) functions
2015-03-29 12:54:13 +00:00
Martin Diehl
9e824d029e
renamed to new names: spectralsolver and spectralfilter
2015-03-28 21:14:11 +00:00
Martin Diehl
1b2e4c0239
removed old solver
2015-03-28 17:41:40 +00:00
Martin Diehl
8a30021202
cleaned up core modul
2015-03-27 14:10:18 +00:00
Pratheek Shanthraj
c8554a5641
corrected alignment of restart output byte offset. restart test working now
2015-03-27 07:17:24 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
86283359f1
PetscFinalize should be the last call before exiting
2015-03-26 13:43:18 +00:00
Pratheek Shanthraj
a9cf8e8051
some patches to compile damask core module without fftw-mpi
2015-03-26 13:41:19 +00:00
Pratheek Shanthraj
2417877bdb
dropped non petsc basic spectral solver
2015-03-25 16:08:41 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
37a7364a3e
multi physics output now working for all solvers
2015-03-25 16:02:30 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
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
Pratheek Shanthraj
c5a8543f17
updated petsc
2015-03-18 17:18:43 +00:00
Martin Diehl
4a0c1c1717
comments
2015-03-15 15:30:14 +00:00
Martin Diehl
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
Martin Diehl
ac3d82493e
changed output keywords form resolution to grid and from dimension to size
2015-03-13 14:45:33 +00:00
Martin Diehl
e4d70e971d
cp is commercially pure ;)
2015-03-13 13:49:09 +00:00
Martin Diehl
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
Martin Diehl
aa125bac5b
aliases for grid and size not supported any more
2015-03-12 22:30:24 +00:00
Martin Diehl
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
Martin Diehl
263f997cf2
seems to be an old output
2015-03-12 08:51:56 +00:00
Martin Diehl
b338f1aeb8
added some outputs
2015-03-09 23:10:36 +00:00
Su Leen Wong
6b0a2ed5a4
Changed initialization of projection matrix for transformation based on the number of slip systems
2015-03-09 12:05:38 +00:00
Martin Diehl
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
Franz Roters
2ab90c28a5
replaced example ODF by reasonable data in fully correct format
2015-03-06 15:56:58 +00:00
Pratheek Shanthraj
35f9e91e73
some corrections to the damage models
2015-03-06 13:12:32 +00:00
Pratheek Shanthraj
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
Yannick Naunheim
4fd7e0165e
more meaningful name
2015-03-05 11:07:55 +00:00
Su Leen Wong
36f14a2d1d
Removed nr = ns condition
2015-02-26 14:16:29 +00:00
Su Leen Wong
2baac6fc10
Fix for 3881 because gfortran was complaining
2015-02-24 10:15:07 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
f1f4f939a3
analytic jacobian default for all solvers
2015-01-29 13:59:49 +00:00
Pratheek Shanthraj
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
Pratheek Shanthraj
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
Martin Diehl
12c586afb9
removed unused variables
2015-01-23 13:44:16 +00:00
Martin Diehl
1f666dffc9
fixed bug: output size was not set
2015-01-23 13:08:25 +00:00
Martin Diehl
e9c3137a3e
disabled a faulty ntrans loop (matmul was complaining)
2015-01-23 11:23:35 +00:00
Martin Diehl
43095b1f1b
added threshold_stress_slip as output (got lost once upon a time)
2015-01-22 09:34:11 +00:00
Martin Diehl
1c2dd37e28
deleted unused variables in dislo UCLA
2015-01-21 15:14:00 +00:00
David Mercier
e57cb2b4f4
Wrong number of slip families in BCC material config files
2015-01-20 17:22:52 +00:00
Luv Sharma
19dc945065
Inroduced Fd in kinematic decomposition.
2015-01-19 13:39:47 +00:00
Pratheek Shanthraj
54dc4943de
adaptive penalty to bound Cv in (0,1)
2015-01-16 18:31:14 +00:00