7d84a0911eremoved unused variables
Franz Roters
2011-04-13 14:16:22 +0000
b4678112cdsorry, last version did not compile. now with correct encoding.
Christoph Kords
2011-04-13 11:51:46 +0000
dd2057bac5changed header to unicode utf-8 encoding
Christoph Kords
2011-04-13 11:24:36 +0000
a9774df43badded comment on need to install AMD Core Math Library
Franz Roters
2011-04-13 09:06:58 +0000
a007c7c4b3removed instructions "how to apply to your program" at end of document
Philip Eisenlohr
2011-04-12 20:13:46 +0000
482be626e0major recoding. now with tiny memory footprint and better guessing of remaining time.
Philip Eisenlohr
2011-04-12 17:46:35 +0000
ef4fc9d0eeadded <crystallite> output "texture" to list of examples
Philip Eisenlohr
2011-04-12 15:33:29 +0000
c1b8391110changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
Franz Roters
2011-04-07 07:20:28 +0000
e4cbdc5743now also works for 2D quadratic triangular element type 125
Nan Jia
2011-04-06 13:11:25 +0000
4d0fefd231security factor of 1.2 for CFL condition in constitutive_nonlocal_dotState
Christoph Kords
2011-04-06 11:47:24 +0000
a989aa979ecorrected increment offsetting for spectral file type. should now be consistent with new (0+N) mpie_spectral output scheme
Philip Eisenlohr
2011-04-06 10:00:19 +0000
51763ed93ecorrected number of incs reported in spectralOut-file
Philip Eisenlohr
2011-04-06 09:58:17 +0000
7c7c929455added comments on origin of files
Franz Roters
2011-04-06 09:40:39 +0000
4ea1fe436badded check for CFL condition in constitutive_nonlocal_dotState
Christoph Kords
2011-04-06 09:07:36 +0000
e00d073ee3added new 2D triangle elements added some sourcecode commenting on internal database formats corrected database to allow for ipVol and Area calc in 2D element cases
Philip Eisenlohr
2011-04-06 08:35:37 +0000
a200bae6cetwo-dimensional plane strain elements with 3 IPs (internal element type 10)
Philip Eisenlohr
2011-04-04 16:10:34 +0000
fffff202b02D mesh option now also alters corresponding job option
Philip Eisenlohr
2011-04-04 15:47:09 +0000
fcdb805225added copyright text to all f90 (free) format files
Franz Roters
2011-04-04 14:09:54 +0000
ed47c09424added file containing copyright text to be added to each file ATTENTION uses f90 comment marker (!)
Franz Roters
2011-04-04 14:07:42 +0000
73accf79b6added option for 2D mesh (3D extrusion step is skipped)
Philip Eisenlohr
2011-04-04 12:12:01 +0000
6ca8942f62choose correct output format for debug_verbosity: integer instead of logical
Christoph Kords
2011-04-04 08:38:55 +0000
790dbed1e4* removed last remnants of old debugger * replaced "dble" intrinsic function by "real" with pReal kind in constitutive_nonlocal.f90 * removed useless line breaks in output of state in CPFEM.f90
Christoph Kords
2011-04-04 08:34:52 +0000
06fc83ac14Final version of titanmod
Alankar Alankar
2011-03-31 09:21:43 +0000
a546f8dbd2changed wording of input queries to (hopefully) clearer statements
Philip Eisenlohr
2011-03-30 15:48:18 +0000
fc471d73f3added Arun to list of authors
Franz Roters
2011-03-30 14:39:28 +0000
dc8961f728added first files for release under GPL 3.0
Franz Roters
2011-03-30 14:38:33 +0000
c36968fc68distance for periodic images in internal stress calculation was incorrect
Christoph Kords
2011-03-30 11:57:41 +0000
08da2b5e4dadded minimum abaqus_v6.env file containing all changes necessary for the use of the MPIE subroutine changed comments in interface routines to refelect this
Franz Roters
2011-03-29 13:49:50 +0000
e9e6abbdafintroduced variables for compiler options added -fpp switch to invoke precompiler
Franz Roters
2011-03-29 07:55:38 +0000
314ca3fe7fInternal stress calculation in nonlocal model: instead of integration of excess gradients (->Bayley) we now sum up contributions from adjacent superdislocations within a certain radius R. When periodicity is used, also periodic images are considered in stress evaluation.
Christoph Kords
2011-03-29 07:34:33 +0000
3d51dd36fa* Introduced preprocessor directives in order to suppress compilation of most write statements when using openmp. This tremendously improves efficiency of parallelization.
Christoph Kords
2011-03-29 07:27:19 +0000
1b3dba4113added -fpp switch to all compiler calls to explicitly use the precompiler removed -save from OpenMP compiler calls to use OpenMP auto behaviour
Franz Roters
2011-03-28 12:18:15 +0000
5c05c5d068corrected mixed up array size for hardeningMatrices (twin/slip)
Philip Eisenlohr
2011-03-24 17:20:35 +0000
fa323f3f0bremoved division by zero from debug_reset(), uses huge() now
Franz Roters
2011-03-24 07:38:56 +0000
78d4e7d1dcmpie_cpfem_abaqus_std.f: first call is with kinc==1 mesh: allow multiple element sets with identical material
Franz Roters
2011-03-23 16:20:12 +0000
8467c91dabmoved setup scripts back to subdirs where they are needed
Franz Roters
2011-03-22 15:42:53 +0000
31a570768bnew error 237: singularity in internal stress calculation
Christoph Kords
2011-03-22 13:58:42 +0000
7be2edb10eadded "texture " as possible crystallite output => reports ID for later discrimination...
Nan Jia
2011-03-22 13:40:27 +0000
b8b75edd7dno changes to previous version, just wanted to add some comment to the previous commit: subroutine "constitutive_nonlocal_dotState" only affects the current dotState, so no need to pass the whole array as argument
Christoph Kords
2011-03-21 15:23:36 +0000
9b1e1e35eadotState
Christoph Kords
2011-03-21 15:18:09 +0000
462eda6736split parallel region in integrateStateFPI into two
Christoph Kords
2011-03-21 14:30:10 +0000
5ee73dee72age dotState before resetting it to zero, otherwise the statedamping does not work.
Christoph Kords
2011-03-21 13:06:11 +0000
8bef254be7changed name to include *.py ending to make Windows file association work again
Franz Roters
2011-03-21 11:58:24 +0000
fc6f2ae68fNo need to call constitutive_microstructure at the beginning of each crystallite step, since it's already been done at the end of the previous step. Just do it once after initialization to start with correct values for the dependent state variables.
Christoph Kords
2011-03-21 10:35:42 +0000
11138e3ee2debugging output is now controlled by the "verbosity" parameter in the debug.config ranging from 0 (=almost no output) to 8 (=very detailed output) 0 : only version infos and all from "hypela2"/"umat" 1 : basic outputs from "CPFEM.f90", basic output from initialization routines, debug_info 2 : extensive outputs from "CPFEM.f90", extensive output from initialization routines 3 : basic outputs from "homogenization.f90" 4 : extensive outputs from "homogenization.f90" 5 : basic outputs from "crystallite.f90" 6 : extensive outputs from "crystallite.f90" 7 : basic outputs from the constitutive files 8 : extensive outputs from the constitutive files
Christoph Kords
2011-03-21 10:31:17 +0000
9d7ede7e03extreme values of stress and jacobian now recorded in CPFEM_general. variable declaration and generation of output moved to debug module.
Christoph Kords
2011-03-17 13:13:13 +0000
235266b169openmp parallelization working again (at least for j2 and nonlocal constitutive model). In order to keep it like that, please follow these simple rules:
Christoph Kords
2011-03-17 10:46:17 +0000
6ac2b4cf88deleted useless storage of data resulting in partial speed-up.
Philip Eisenlohr
2011-03-11 16:32:01 +0000
97a0146672deactivated one debug statement at the end of crystallite_updateState(g,i,e) as it leads to division by zero if constitutive_state(g,i,e)%p(1:mySize) contains zero values statement should be changed to check for zero values
Franz Roters
2011-03-11 10:05:30 +0000
1f0569230bnew installation/setup dir to hold all top-level make_*-scripts
Philip Eisenlohr
2011-03-11 09:20:37 +0000
41c0bf098cstarting to substructure manual
Philip Eisenlohr
2011-03-11 09:19:31 +0000
bdc17f7d72extended comment on modification of abaqus_v6.env
Franz Roters
2011-03-11 09:05:09 +0000
d860015d57first shot at reviving the manual...
Philip Eisenlohr
2011-03-11 08:43:29 +0000
621c1198a4Marc/Mentat related things are only executed for filetype marc. this has been too lax before, causing trouble.
Philip Eisenlohr
2011-03-10 09:45:57 +0000
6475235fb4included "-heap-arrays 500000000" in remark on how to change compile command
Franz Roters
2011-03-08 08:18:04 +0000
18e668340fAbaqus interfaces dapted to newest version of the rest of the code does not yet check for restart, have to check out first how restarts work in Abaqus explicit code not checked as we do not have a license
Franz Roters
2011-03-07 12:49:27 +0000
654b5ea987set type for acos_arg
Franz Roters
2011-03-04 14:57:22 +0000
2975c99144now the maybe not perfect quartenions caused trouble in math_QuaternionToEuler
Franz Roters
2011-03-03 14:23:39 +0000
e7c7ccdcddimproper rotation matrix from pdecomposition could cause negative arguments in squareroot
Franz Roters
2011-03-03 10:47:07 +0000
6667e94238abaqus routines VGETOUTDIR and GETOUTDIR do not return the trailing slash, which we assume to be there, so it has to be appended
Franz Roters
2011-03-03 10:29:59 +0000
6b63d488cbneed py_post, not py_mentat
Christoph Kords
2011-02-28 10:00:08 +0000
9312aed545set executable bit
Philip Eisenlohr
2011-02-27 21:51:56 +0000
138e9c43ebsome comments (talking to Ricardo)
Martin Diehl
2011-02-25 16:24:23 +0000
a6d1e5c911now deleted old _single files
Martin Diehl
2011-02-25 12:49:18 +0000
8cae4d609asome polishing for single precision version. Now only prec and mpie_spectral are needed in different versions
Martin Diehl
2011-02-25 12:41:46 +0000
ad4706673b* removed calculations for dipole formation/dissociation by stress change, since it is not used anyways; also removed associated constitutive outputs from material.config * removed input variables in constitutive_collectDotState and constitutive_postResults that are not needed anymore (because of recent changes in constitutive_nonlocal)
Christoph Kords
2011-02-25 09:53:20 +0000
cd5407b08bremoved all math functions only for double precision by the more flexible counterpart, e.g. "dsqrt --> sqrt", "dsin --> sin". Should not cause any harm, as long as "implicit none" is used. Now it is possible to compile a single precision spectral solver/crystal plasticity by replacing mesh.f90 and prec.f90 with mesh_single.f90 and prec_single.f90. For the spectral method, just call "make precision=single" instead of "make". Use "make clean" evertime you switch precision
Martin Diehl
2011-02-25 09:25:53 +0000
2f503f5cdbadded "Lp" as crystallite output
Christoph Kords
2011-02-25 08:15:26 +0000
e022810e66fixed problem in internal stress calculation for periodic neighborhood
Christoph Kords
2011-02-25 08:10:11 +0000
f525c02dedadded single precision libraries for FFTW First try of implement single precision crystal plasticity, not working yet. polishing text about geometry construction. polishing postResults, still having problems concerning machines without MSC installation
Martin Diehl
2011-02-24 18:39:57 +0000
e49e5e13af* in flux calculation: take care of special case for single element model * do not overwrite "constitutive_nonlocal_rhoDotFlux" when doing a state perturbation, so we can use it for output
Christoph Kords
2011-02-24 10:01:41 +0000
8aef9b1c13check for twin nodes only in direction of the surface normal. this solves a problem with models of single element thickness and periodic fluxes.
Christoph Kords
2011-02-24 09:26:30 +0000
fed53a7e60in subroutine "mesh_marc_get_mpieOptions": need at least one additional string after keywords "$mpie" and "periodic", not two as it used to be
Christoph Kords
2011-02-23 12:33:51 +0000
430966eb20sorry, forgot to commit the numerics.f90, too
Christoph Kords
2011-02-23 12:30:52 +0000
96d3682d5e* new global integer variable "numerics_integrationMode" (1 indicating integration of central solution, 2 indicating integration of perturbed state) * combined "integrator" and "integratorStiffness" in new global variable "numerics_integrator"
Christoph Kords
2011-02-23 08:29:51 +0000
d835380bc0* "constitutive_nonlocal_flux" is set to zero in "constitutive_nonlocal_dotState", not in "constitutive_nonlocal_microstructure" * dislocation flux and internal stress calculation now consistent with new definition of slip system lattice according to paper (polarity of screws inverted)
Christoph Kords
2011-02-23 08:08:06 +0000
3c944e79fadid some things Philip told me to do (postResults for spectral files is now working without MSC installation)
Martin Diehl
2011-02-22 16:18:47 +0000
003b0def46fixed a stupido error
Philip Eisenlohr
2011-02-22 16:06:24 +0000
333445be8bchanged meshing routine from python to fortran (f2py)
Martin Diehl
2011-02-22 16:00:22 +0000
2dc3b35725no need anymore for msc python interface when dealing with spectralOut files
Philip Eisenlohr
2011-02-22 15:57:27 +0000
8a5e28d5a6now using fft reconstruction in 3Dvisualize, linear python code for reconstruction is removed
Martin Diehl
2011-02-21 21:03:21 +0000
c399a06c97tried to implement an inverse reconstruction, meaning to calculate an average deformation gradient from the 8 cornes of a node in reference and current configuration. Not working yet. polishing, mainly in fft reconstruction. added fftw library file and fortran file for type specification added folder references with literature for visualization/geometry reconstruction
Martin Diehl
2011-02-21 16:30:18 +0000
8dd1a694a3extended IO to cope with different name for solverJob and Model polishing, added error codes added FFTW library files
Martin Diehl
2011-02-21 14:37:38 +0000
34336bc659fixed bug when no argument was supplied
Philip Eisenlohr
2011-02-18 09:02:29 +0000
b21c153761restructured constitutive_nonlocal_compatibility and incorporated the "surfaceTransmissivity" in the compatibility calculation
Christoph Kords
2011-02-16 16:38:18 +0000
24d33bf2ff* added a new material parameter "surfaceTransmissivity" (default value 1.0) which allows to change the transmissivity of the material surface between 0 and 1 * now complaining when encountering an unknown nonlocal parameter in material.config * use same error ID for all material parameters out of bounds * symmetric flux calculation in side dotState can now be omitted (because of new treatment of periodicity) * switching back to "local flux balance" (add leaving and entering fluxes at central MP, don't touch neighbor) instead of "flux distribution" (subtract leaving fluxes from central MP and add them at neighboring MP). This has the advantage that there is almost no need for CRITICAL statements in parallelization, so hopefully this results in some speed up.
Christoph Kords
2011-02-16 16:35:38 +0000
8f626c8989Now able to have real periodicity for fluxes (fluxes leaving the model on one side enter on the other side). To enable this feature one has to add the following somewhere in the marc input file: $mpie periodic x y z for having periodicity in all directions $mpie periodic z x for having periodicity in x and z direction etc. Note that this only works for regular meshes!!!
Christoph Kords
2011-02-16 16:23:08 +0000
c39e04644bSorting of hex slip-systems done in update 659 was missing. Included again.
Alankar Alankar
2011-02-15 12:21:54 +0000
438bf95105added ft-based reconstruction of deformed configuration to postprocessingMath.f90 and postprocessingMath.pyf also added function to calculate average of tensor some polishing in mpie_spectral.f90, added sanity check to see im resolution is power of 2
Martin Diehl
2011-02-14 17:21:31 +0000
cdeb786721added default path of ACML BLAS now also works on msuws2
Franz Roters
2011-02-14 15:03:45 +0000
bf2082a15bmarc run scripts now link against ACML (AMD Core Math Library) to enable use of BLAS routines ACML has certainly to be installed first (currently on msuws1 only) and is included for this purpose
Franz Roters
2011-02-14 14:25:16 +0000
874d5067c7switches drawing of servo links temporarily off to speed up processing.
Philip Eisenlohr
2011-02-11 11:29:58 +0000
171f5375f1shearrate output now with sign (no |abs| as before)
Philip Eisenlohr
2011-02-11 10:19:41 +0000
3d7fad6ba9mpie_spectral.f90: changed fourier transormation, now using the advanced interface to do the tranformation both ways with just one call. now also using the same variable for pk-stress in both domains and the change of deformation gradient in both domains.
Martin Diehl
2011-02-09 17:47:28 +0000
803e1a8c05In order to help preventing further memory leaks all array sections now have an explicit instead of assumed shape, e.g. Fe(1:3,1:3) instead of Fe(:,:).
Christoph Kords
2011-02-09 13:12:46 +0000
65ae7de42bfixed bug for missing file name
Philip Eisenlohr
2011-02-09 09:18:21 +0000
466ced911e* check convergence for nonlocals also in stiffness calculation mode
Christoph Kords
2011-02-09 08:39:07 +0000
921967c2ba* found memory leak in nonlocal model! "constitutive_nonlocal_rhoDotFlux" has to have upper index bound of 10 for second dimension, not 8. * need to recalculate dislocation velocity in postResults, otherwise we take values of last perturbed state! So the following outputs were up to now showing the perturbed state: shearrate, dislocation velocity, all density rates!
Christoph Kords
2011-02-09 08:28:47 +0000
17a6c0e59benhanced comment on dislocation kinetics
Christoph Kords
2011-02-08 12:48:50 +0000