d94236b2d6Declaration of filename etc. was missing. Fixed now.
Alankar Alankar
2010-09-23 09:13:46 +0000
83a921eec2crystallite: polishing rest: added "$ID" and line ending settings...
Philip Eisenlohr
2010-09-23 08:05:50 +0000
a467fb4ae3forgot "$ID" and setting line endings...
Philip Eisenlohr
2010-09-23 08:03:00 +0000
89136afadedebug now has config file, too
Philip Eisenlohr
2010-09-23 08:01:41 +0000
23746e7e78chopped line beginnings introduced check for PIL, if this fails, no image will be output
Philip Eisenlohr
2010-09-23 07:59:52 +0000
6ea8623f65added functions math_mul33x3_complex math_Plain99to3333 to math.f90. mpie_spectral.f90 uses both functions math_Plain99to3333 is used for inversion of c0 math_mul33x3_complex is used for equilibrium check in fourier space also did some cleaning up on mpie_spectral.f90
Martin Diehl
2010-09-22 12:04:43 +0000
80618c9814Cleaned-up, rearranged state variables
Alankar Alankar
2010-09-22 09:00:40 +0000
7548858ffaSome useful changes in damper for quicker convergence.
Alankar Alankar
2010-09-22 08:54:36 +0000
604992a9e1fixes bugs in screen output, changed gmsh output, now the deformed mesh is written out
Martin Diehl
2010-09-22 08:51:34 +0000
6920d8d939cleaning up, added new output to gmsh and to screen, prepared calculation of deformed mesh new convergence criteria: divergence of stressfield (in fourier space)
Martin Diehl
2010-09-21 15:38:25 +0000
8488cafe3anow consistent with numerics.f90 values
Alankar Alankar
2010-09-21 14:37:52 +0000
9a20f742ealeaner memory requirement to store states and their rates when calculating the crystallite stiffness
Philip Eisenlohr
2010-09-13 09:13:25 +0000
d8a4cd1739added optional argument --baseindex to specify what colorset should be influenced. identify colors, for instance, start at idx 214...
Philip Eisenlohr
2010-09-13 07:31:27 +0000
006dce7d8aprinting euler angles (and defgrads and orientations) for individual ips now works
Alankar Alankar
2010-09-08 11:55:57 +0000
1693bfca47changed calculation of defgrad from cauchy-stress to pk-stress, now working in large-strain-formulation. output of msh-file changed to deformed configuration, removed output of defgrad as a field
Martin Diehl
2010-09-07 16:37:55 +0000
3c5f38643dBrushed up accountability of twinning to Lp
Alankar Alankar
2010-09-07 14:44:37 +0000
7d4c7f7fa7corrected bug in debug: stressloop info was erroneous mesh init now showing much more ip statistics in verbose mode
Christoph Kords
2010-09-07 09:06:02 +0000
4a31b175da1) introduced default colors 2) complaints are more verbose
Philip Eisenlohr
2010-09-06 16:51:11 +0000
a3c92061f91) terminallyIll was reset before FE did cutback --> useless extra calculations of same problem over and over... 2) local stiffness calculation is now standard for non-local grains 3) stressLoopDistribution discriminates between (a) central solution and (b) stiffness perturbation 4) debugger is switched on as standard... (but verboseDebugger not!)
Philip Eisenlohr
2010-09-06 16:06:41 +0000
00922705ebnow using current average stiffness for reference material stiffness
Martin Diehl
2010-09-06 10:00:59 +0000
c19524f2641) fixed terminallyIll bug occurring with Ngrains=1 homogenization 2) brushed up the output to be more easily readable/understandable
Philip Eisenlohr
2010-09-01 21:04:02 +0000
4f29e8c2feIntroduced twinning. Twinning based on accumulated shear on twin systems. No slip-twin interactions yet.
Alankar Alankar
2010-09-01 10:07:52 +0000
a5c228fd02changed back transform to complex-to-real, removed redundant variables, reduced size of arrays where possible
Martin Diehl
2010-09-01 08:05:11 +0000
85febf0803added fftw3 as fft(library will not versioned, should be in a linkable folder) , did some corrections on the code, splitted main file up (allows use of makefile), added makefile
Martin Diehl
2010-08-27 16:39:38 +0000
05568ef73bcalculation of state damper was erroneous for use with nonlocal and other constitution at the same time
Christoph Kords
2010-08-26 12:54:31 +0000
197d7ea8d9rearranged output statements in crystallite init for better readability in output file
Christoph Kords
2010-08-26 09:52:40 +0000
ec4269bdd8Added 'dislocation-type' dependent interactions. Ideal for hexagonal and BCC structures.
Alankar Alankar
2010-08-23 11:36:51 +0000
da53656c9aadded C3D20R (type 9) to list of possible non-local elements (Christoph, should we invent a reasonable error number for this? "-1" may not be optimal...)
Philip Eisenlohr
2010-08-19 23:00:26 +0000
c352146fe6exchanged forall (warning causers) to nested do loops
Philip Eisenlohr
2010-08-19 21:35:38 +0000
af2437c0afremoved hardwired installation directory from submitX scripts and added corresponding flexibility into make script
Philip Eisenlohr
2010-08-17 21:22:24 +0000
3b671faae7added some flexibility (and bug-fixing) to the pre/post processing. the file 'MSCpath' contains the LOCAL path to the MSC installation (/msc in our case)
Philip Eisenlohr
2010-08-17 18:21:22 +0000
f78b07448fsome modifications in line with Davids analytical model input some additional variables via material.config
Franz Roters
2010-08-17 14:23:55 +0000
bb9899e7de1) added element type 57 (C3D20R) reduced quadratic hexahedral 2) need discussion about the present omission of the FE_nodesAtIP for type 8 --> added comment into source code
Philip Eisenlohr
2010-08-16 22:53:24 +0000
94efd57663started to put Marc/Mentat related scripts to corresponding post/pre processing folders. We should move development from my Code folder to here (sorry, lost history then...)
Philip Eisenlohr
2010-08-16 20:47:27 +0000
d1ed15fbdfpath to MSC is now command line option(al)
Philip Eisenlohr
2010-08-16 20:36:02 +0000
9f071c3e01Added kink-pair pre-factor for screw dislocation mobilities on non-basal slip systems in hex
Alankar Alankar
2010-08-06 15:53:45 +0000
298cecbfeccorrected syntax errors (long lines, line continuation by \) and logical mistake in mpie_cpfem_marc / abq_std which went unnoticed in ifort. Using SunStudio f90 surfaced those...
Philip Eisenlohr
2010-08-03 23:47:00 +0000
2c4f5c37dcsource code check to find lines exceeding 132 char limit of Fortran
Philip Eisenlohr
2010-08-03 23:45:10 +0000
c7e5cb9806reworked verbosity and structure of installation script
Philip Eisenlohr
2010-07-30 19:39:54 +0000
60decea760Calibration (Prismatic slip) of the new dislocation density based CPFE model for single crystal alpha-Ti (Titanmod). Formulated edge and screw dislocations mobilities in presence of kink pairs. Updated comments. The model shows the plateau in stress evolution in prismatic orientation which is observed in experiments.
Alankar Alankar
2010-07-15 07:16:15 +0000
49926d5d66F77 --> F90 polishing/condensation as far as possible. next step is change of Re/Im FFT to Re only...
Martin Diehl
2010-07-13 15:29:26 +0000
4d110126daadopted improvements done by Arun Prakash.
Philip Eisenlohr
2010-07-13 10:26:07 +0000
eb7830dc8fNew Titanium CPFE model based on dislocation density. Running but calibration and changes needed.
Alankar Alankar
2010-07-13 08:19:25 +0000
397db06fb3cpfe_general now returns also Kirchhoff-stress P and dPdF
Martin Diehl
2010-07-07 09:58:18 +0000
80016f8429added P and dPdF for the call to cpfem_general, made some comments in the code, aligned some loops
Martin Diehl
2010-07-07 09:10:54 +0000
f9834bc612changed types of integers and reals in fourier transform, changed some variables and cleaned up code to make it easier to understand
Martin Diehl
2010-07-05 16:01:36 +0000
c3e222dbbderrs added etc.
Ricardo Lebensohn
2010-07-05 11:33:48 +0000
4f76eada31started suggestions for F77 --> F90 style
Philip Eisenlohr
2010-07-02 17:15:53 +0000
d6ba9d54b6changed declaration of two variables from real to int, spectral method is now working. started to commend Ricardos code, layoutet loops and removed some redundant variables. until now, no error calculation is done. at the moment calculations are in an infinite loop
Martin Diehl
2010-07-02 14:10:36 +0000
661bb97800corrected some problems in building mesh
Martin Diehl
2010-07-01 15:20:39 +0000
2ec39c936cadded first version of spectral method, changed corrected functions to get correct path to files
Martin Diehl
2010-07-01 15:20:06 +0000
c3e51e5330added functions to read out resolution, dimension and homogenization from mesh file
Martin Diehl
2010-06-25 11:31:05 +0000
eb0cc24a8eno need for material parameter "a" anymore, dislocation stress calculation now works without "core spreading"
Christoph Kords
2010-06-21 16:01:05 +0000
7d6e52067bdislocation stress based on dislocation density gradients
Christoph Kords
2010-06-21 15:58:56 +0000
6d874e2c1fnonlocal stiffness calculation: rather perturb all components at once (and optionally decrease the frequency of the Jacobian update with the iJaco parameter) than perturbing only a single component per cycle
Christoph Kords
2010-06-17 06:32:56 +0000
1c72439350first shot at init with CPFEM_general better would be to manually read the mesh resolution, set up Spectral solution strategy and call CPFEM_general just during BVP solution.
Philip Eisenlohr
2010-06-10 15:32:06 +0000
8ea52ac495mpie_spectral now can read loadcase data adapted error codes
Philip Eisenlohr
2010-06-10 14:51:10 +0000
386ca7ebc3added some basic I/O functions to mpie_spectral.f90
Martin Diehl
2010-06-10 08:50:04 +0000
740db98090- use upwind differences scheme for fluxes instead of central differences - introduced possibility to "reflect" fluxes at free surfaces in order to have quasi-periodic fluxes
Christoph Kords
2010-06-09 08:56:00 +0000
11b98fbfa6added ID functionality and correct line endings
Philip Eisenlohr
2010-06-08 10:10:57 +0000
82687d31c8first shot at mpie_spectral wrapper
Philip Eisenlohr
2010-06-08 10:08:15 +0000
7e78932976renamed "mpie_cpfem_init" subroutine within mpie_interface module to more general name "mpie_interface_init"
Philip Eisenlohr
2010-06-08 09:35:13 +0000
71c6179ee5dislocation stress had wrong sign in some components
Christoph Kords
2010-06-08 09:34:23 +0000
3ce1882d29when calculating dislocation stress at grain boundaries, densities are now extrapolated similarly to like it was already done at free surfaces
Christoph Kords
2010-06-07 16:01:37 +0000
3989cc2688dislocation stress now considers dislocation segments of length V^(1/3) and ensures non-singular solution for points that are collinear to the dislocation line by means of a core spreading parameter "a". stress is scaled by a reference cutoff radius "r" to avoid mesh-dependent stress values.
Christoph Kords
2010-06-07 14:32:23 +0000
5283ba5546if it is pReal it should be 1.0 not just 1
Franz Roters
2010-05-31 07:12:45 +0000
c1fed9d197polishing of verbosity
Philip Eisenlohr
2010-05-28 16:05:16 +0000
76b392e12fdislocation stress calculation in constitutive_nonlocal_microstructure: corrected mapping from lattice to slip system frame and the neighboring dislocation density formulation
Christoph Kords
2010-05-28 13:12:36 +0000
0d69623b74grpmf...
Philip Eisenlohr
2010-05-27 14:32:16 +0000
89f09ed784use testMe in this subdir to check for compilation errors.
Philip Eisenlohr
2010-05-27 14:31:04 +0000
920bccc556single IP, reduced integration element for 3D
Philip Eisenlohr
2010-05-26 17:30:16 +0000
bb73a74325there was a line break missing so that it did not compile
Franz Roters
2010-05-26 17:10:26 +0000
9932c90c73restoring the crystallite.f90, IO.f90 and material.f90 to the version 575. the version 576 of these files contain some uncleared debugging statements.
Denny Tjahjanto
2010-05-26 16:34:44 +0000
d114a600c3fixed some errors in QuaternionToEuler, RToQuaternion to get the correct texture representation in the output.
Denny Tjahjanto
2010-05-26 15:52:54 +0000
e5f0af638efixed calculation of dislocation stress for elements with a single ip and no neighbors; used to get segmentation fault
Christoph Kords
2010-05-26 09:17:10 +0000
0c5bc83469in nonlocal_microstructure: at the model surface the excess density is now extrapolated from opposite neighbor instead of assuming zero density; this results in better modeling of the dislocation stress at the surface restructured nonlocal_dotstate, to be able to easily switch on and off particular effects in the microstructure evolution nonlocal_dotstate now enforces a cutback when single density runs the risk of becoming negative; in the case of a state already below the relevantState dotState is set to zero introduced two new output variables: rho_dot_edge and rho_dot_screw
Christoph Kords
2010-05-21 08:51:15 +0000
dba4ae7ef1new numerics parameter "defgradTolerance" replaces "relevantStrain" in cpfem, since "relevantStrain" is also used in crystallite but in a different context. "defgradTolerance" gets the same standard value as "relevantStrain" so calculations should be consistent with older ones
Christoph Kords
2010-05-20 14:55:11 +0000
3a23dc4bbdfirst call to hypela2 in inc 0 is done with computation mode 2 (=plain calc), but since rev. 539 the second call was done with computation mode 6 (=only return values) instead of 2. This doesnt work, because computation mode 2 might calculate stress and tangent for all materialpoints, but the mapping from lattice to current configuration is only done for the element/IP that triggers the call.
Christoph Kords
2010-05-20 09:14:59 +0000
51da3b6242sorry, in rev 568 some variables were not declared in crystallite_postResults
Christoph Kords
2010-05-18 12:36:09 +0000
472684c30emake element type 8 (= hex8 reduced integration) available for nonlocal constitutive law
Christoph Kords
2010-05-18 08:13:23 +0000
19b6fae88badded a missing exclamation mark in front of an $OMP directive
Christoph Kords
2010-05-18 08:07:55 +0000
d96780967badded some output variables in crystallite: F, Fe, Fp, Ee, P, S
Christoph Kords
2010-05-18 07:57:13 +0000
15328e51cfLike this Euler angles AND grain rotation should be OK.
Claudio Zambaldi
2010-05-11 15:06:21 +0000
483584ceb3improved version of general I/O functions renamed module cpfem_interface to mpie_interface
Franz Roters
2010-05-11 06:57:15 +0000
e810e5cfa1changed fileopen routines in IO.f90 to use solver dependant functions defined in the interface routines DID NOT test ABAQUS versions yet, not shure whether (V)GETOUTDIR returns path with or without terminating slash
Franz Roters
2010-05-10 15:02:59 +0000
0416c9a616reworked build_ipNeighborhood added C3D8R reduced integration hexahedral element
Philip Eisenlohr
2010-05-10 14:54:59 +0000
2405a51042crystallite: - grainrotation calculation now is done with symmetryID 0, i.e. without symmetry reduction since we want the absolute misorientation. - While math has everything in radians, post results eulerangles and axisangle are given in degrees. And: grain rotation seems OK after the previous changes in math module.
Claudio Zambaldi
2010-05-07 12:01:46 +0000
97f206d0bcnew errors regarding Spectral method input file parsing mesh setup from Spectral input file $Id$ in FEsolving activated
Philip Eisenlohr
2010-05-06 16:40:47 +0000
5a0963fb81added variable working dir (if apply_script is called from somewhere else) and made the script executable
Philip Eisenlohr
2010-05-06 14:17:53 +0000
e4067f2d26debugged zoo of rotation operations and transformations
Philip Eisenlohr
2010-05-06 14:07:21 +0000
4954228b3aforgot to change this file accordingly (main.msb revmoved)
Franz Roters
2010-05-06 13:19:45 +0000
f1950cacc9corrected apply script removed binary menu file, it is now created by the apply script
Franz Roters
2010-05-06 13:17:28 +0000
f23b4d16d1updated section with dislotwin paramters
Franz Roters
2010-05-06 09:01:54 +0000
383a36b226added comment on reformating to common block files added modifications necessary for running the subroutine with Marc/Mentat2010
Franz Roters
2010-05-06 08:50:54 +0000
393a67e587now stating instance when sanity checks run into trouble
Alankar Alankar
2010-05-05 13:36:59 +0000
1ee5df41fdSorry, that was the wrong file, now the correct one...
Claudio Zambaldi
2010-05-04 17:34:51 +0000
48c554811bA simple one-element model that you can use to try out changed code before commits. Hopefully, the start of a testing suite.
Claudio Zambaldi
2010-05-04 17:32:17 +0000
0a36527f61fixed bugs in math_quaternionDisorientation though running, there still seem to be some problems with the quaternion functions, since crystallite_grainRotation gives some dubious results!
Christoph Kords
2010-05-04 16:02:05 +0000
Alankar Alankar