Commit Graph

3807 Commits

Author SHA1 Message Date
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
Martin Diehl 5f986f3d54 updated to contain latest python improvements 2014-09-10 16:14:37 +00:00
Luv Sharma e61e6c4ca7 Now checking field type to set number of active fields. 2014-09-10 15:14:31 +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
Pratheek Shanthraj b2787f8e66 added support for symmetric quadrature rules unto order 5 2014-09-10 09:25:45 +00:00
Luv Sharma dc406a01c0 added funtions to get averged properties at integration points. 2014-09-10 08:37:12 +00:00
Pratheek Shanthraj 9b4c1fc43d added DAMASK_FEM.exe to bin_link list 2014-09-10 08:09:35 +00:00
Martin Diehl f1e345dcdb fixed compilation for FEM 2014-09-09 15:59:41 +00:00
Pratheek Shanthraj d2d52ee666 updated FEM examples 2014-09-09 15:52:54 +00:00
Pratheek Shanthraj 0c5d3a88ca updated to native Petsc FEM and quadrature definitions 2014-09-09 15:39:39 +00:00
Franz Roters b6847ebb01 added use mesh statement needed for output of mesh_element 2014-09-09 15:18:49 +00:00
Philip Eisenlohr 5a159dcdfa script to add scalar data and/or RGB color tuples to rectilinear VTK file (*.vtr) 2014-09-05 20:28:39 +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
Tias Maiti e286cf4c74 fixed bug "theTable" instead of "table" 2014-09-03 14:53:12 +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
Tias Maiti b566eed5a8 fixed error in last commit 2014-08-28 18:50:48 +00:00
Tias Maiti 5d96f62b44 fixed bug "theTable" instead of "table" 2014-08-28 17:37:27 +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
Chen Zhang 5059b5ce2d polishing 2014-08-27 15:19:04 +00:00
Test User 3697bdb06e Added report of SUCCESSFUL test for Revision 3423 2014-08-27 02:44:13 +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
Luv Sharma ff390f96f4 added quadrature info for 2d elements also 2014-08-26 13:14:46 +00:00
Martin Diehl 8d53ec09eb updated test to last changes in cli interface 2014-08-26 06:14:15 +00:00
Martin Diehl 6cdcdecdfa post processing scripts are working again.
@Philip: definining one offset for shiftData and one factor for scaleData for tensor/vector/special>1 seems to be more reasonable than per component
2014-08-25 21:13:39 +00:00
Luv Sharma fa45d376a1 added quadrature points and weight for dim=3 order=4,5 elements 2014-08-25 16:08:57 +00:00
Martin Diehl d93c40a3dd fixed small bug ('theTable' instead of 'table') and simplified 2014-08-25 12:53:11 +00:00
Martin Diehl 1ba7cbb046 fixed small bug ('theTable' instead of 'table') and simplified 2014-08-25 12:39:12 +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
Martin Diehl c7837e1f66 changed manual input/output close to class method 2014-08-22 16:58:53 +00:00
Philip Eisenlohr 28a7c4c727 corrected (probable) bug in disorientation calculation. 2014-08-22 15:45:03 +00:00
Philip Eisenlohr d5e2e42f21 file handle closing can now be done from within class.
methods:
.input_close()
.output_close()
2014-08-22 15:37:46 +00:00
Philip Eisenlohr 99e3c9a397 fixed (potentially) buggy rotation of tensors. convention is:
with A v = b
v' = R v
A' = R A R^T
such that
A' v' = R A R^T  R v = R A v (= R b = b', i.e. old result transformed into other coordinate basis)
2014-08-22 15:32:48 +00:00
Philip Eisenlohr df27115472 extended "scalar" capabilities to "special" having arbitrary dimension (defaults to 1, i.e. former scalar case). Breaks backward compatibility since "--scalar" is now "--special". However, "-s" still works as before... 2014-08-22 15:25:51 +00:00