660c9212f8
introduced error 667 for out of bounds user result requesting
2011-05-04 16:02:18 +00:00
affd383ef8
simplified disorientation calculation. performed FEM check that "grainrotation" reflects the *active* rotation from the starting orientation to the current one with respect to the lab coord system (as expected).
2011-04-14 19:09:44 +00:00
5814262f55
moved some write statements within stateIntegrationEuler
2011-04-14 09:35:56 +00:00
349f022100
reinstantiated calculation of grainrotation (still needs check for correctness, i.e. active or passive rotation output)
2011-04-13 16:08:52 +00:00
7d84a0911e
removed unused variables
2011-04-13 14:16:22 +00:00
b4678112cd
sorry, last version did not compile. now with correct encoding.
2011-04-13 11:51:46 +00:00
dd2057bac5
changed header to unicode utf-8 encoding
2011-04-13 11:24:36 +00:00
ef4fc9d0ee
added <crystallite> output "texture" to list of examples
2011-04-12 15:33:29 +00:00
c1b8391110
changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
2011-04-07 07:20:28 +00:00
4d0fefd231
security factor of 1.2 for CFL condition in constitutive_nonlocal_dotState
2011-04-06 11:47:24 +00:00
51763ed93e
corrected number of incs reported in spectralOut-file
...
inc 0 contains undeformed results
plus lots of typographic polishing
2011-04-06 09:58:17 +00:00
7c7c929455
added comments on origin of files
2011-04-06 09:40:39 +00:00
4ea1fe436b
added check for CFL condition in constitutive_nonlocal_dotState
2011-04-06 09:07:36 +00:00
e00d073ee3
added 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
2011-04-06 08:35:37 +00:00
fcdb805225
added copyright text to all f90 (free) format files
2011-04-04 14:09:54 +00:00
6ca8942f62
choose correct output format for debug_verbosity: integer instead of logical
2011-04-04 08:38:55 +00:00
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
2011-04-04 08:34:52 +00:00
06fc83ac14
Final version of titanmod
2011-03-31 09:21:43 +00:00
fc471d73f3
added Arun to list of authors
2011-03-30 14:39:28 +00:00
c36968fc68
distance for periodic images in internal stress calculation was incorrect
2011-03-30 11:57:41 +00:00
08da2b5e4d
added minimum abaqus_v6.env file containing all changes necessary for the use of the MPIE subroutine
...
changed comments in interface routines to refelect this
2011-03-29 13:49:50 +00:00
e9e6abbdaf
introduced variables for compiler options
...
added -fpp switch to invoke precompiler
2011-03-29 07:55:38 +00:00
314ca3fe7f
Internal 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.
2011-03-29 07:34:33 +00:00
3d51dd36fa
* Introduced preprocessor directives in order to suppress compilation of most write statements when using openmp. This tremendously improves efficiency of parallelization.
...
* Also added some more openmp directives to increase percentage of parallelized code.
* "implicit none" was missing in two subroutines of homogenization and constitutive.
2011-03-29 07:27:19 +00:00
5c05c5d068
corrected mixed up array size for hardeningMatrices (twin/slip)
...
added explicit array range 1:lattice_maxNslip/twinFamily
2011-03-24 17:20:35 +00:00
fa323f3f0b
removed division by zero from debug_reset(), uses huge() now
2011-03-24 07:38:56 +00:00
78d4e7d1dc
mpie_cpfem_abaqus_std.f: first call is with kinc==1
...
mesh: allow multiple element sets with identical material
2011-03-23 16:20:12 +00:00
8467c91dab
moved setup scripts back to subdirs where they are needed
2011-03-22 15:42:53 +00:00
31a570768b
new error 237: singularity in internal stress calculation
2011-03-22 13:58:42 +00:00
7be2edb10e
added "texture
...
" as possible crystallite output => reports ID for later discrimination...
2011-03-22 13:40:27 +00:00
b8b75edd7d
no 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
2011-03-21 15:23:36 +00:00
9b1e1e35ea
dotState
2011-03-21 15:18:09 +00:00
462eda6736
split parallel region in integrateStateFPI into two
2011-03-21 14:30:10 +00:00
5ee73dee72
age dotState before resetting it to zero, otherwise the statedamping does not work.
2011-03-21 13:06:11 +00:00
fc6f2ae68f
No 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.
2011-03-21 10:35:42 +00:00
11138e3ee2
debugging 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
If verbosity is equal to zero, all counters in debug are not set during calculation (e.g. debug_StressLoopDistribution or debug_cumDotStateTicks). This might speed up parallel calculation, because all these need critical statements which extremely slow down parallel computation.
2011-03-21 10:31:17 +00:00
9d7ede7e03
extreme values of stress and jacobian now recorded in CPFEM_general. variable declaration and generation of output moved to debug module.
2011-03-17 13:13:13 +00:00
235266b169
openmp parallelization working again (at least for j2 and nonlocal constitutive model).
...
In order to keep it like that, please follow these simple rules:
DON'T use implicit array subscripts:
example: real, dimension(3,3) :: A,B
A(:,2) = B(:,1) <--- DON'T USE
A(1:3,2) = B(1:3,1) <--- BETTER USE
In many cases the use of explicit array subscripts is inevitable for parallelization. Additionally, it is an easy means to prevent memory leaks.
Enclose all write statements with the following:
!$OMP CRITICAL (write2out)
<your write statement>
!$OMP END CRITICAL (write2out)
Whenever you change something in the code and are not sure if it affects parallelization and leads to nonconforming behavior, please ask me and/or Franz to check this.
2011-03-17 10:46:17 +00:00
97a0146672
deactivated 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
2011-03-11 10:05:30 +00:00
1f0569230b
new installation/setup dir to hold all top-level make_*-scripts
2011-03-11 09:20:37 +00:00
bdc17f7d72
extended comment on modification of abaqus_v6.env
2011-03-11 09:05:09 +00:00
6475235fb4
included "-heap-arrays 500000000" in remark on how to change compile command
2011-03-08 08:18:04 +00:00
18e668340f
Abaqus 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
2011-03-07 12:49:27 +00:00
654b5ea987
set type for acos_arg
2011-03-04 14:57:22 +00:00
2975c99144
now the maybe not perfect quartenions caused trouble in math_QuaternionToEuler
2011-03-03 14:23:39 +00:00
e7c7ccdcdd
improper rotation matrix from pdecomposition could cause negative arguments in squareroot
2011-03-03 10:47:07 +00:00
6667e94238
abaqus routines VGETOUTDIR and GETOUTDIR do not return the trailing slash, which we assume to be there, so it has to be appended
2011-03-03 10:29:59 +00:00
138e9c43eb
some comments (talking to Ricardo)
2011-02-25 16:24:23 +00:00
a6d1e5c911
now deleted old _single files
2011-02-25 12:49:18 +00:00
8cae4d609a
some polishing for single precision version.
...
Now only prec and mpie_spectral are needed in different versions
2011-02-25 12:41:46 +00:00
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)
2011-02-25 09:53:20 +00:00
cd5407b08b
removed 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
2011-02-25 09:25:53 +00:00
2f503f5cdb
added "Lp" as crystallite output
2011-02-25 08:15:26 +00:00
e022810e66
fixed problem in internal stress calculation for periodic neighborhood
2011-02-25 08:10:11 +00:00
f525c02ded
added 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
2011-02-24 18:39:57 +00:00
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
2011-02-24 10:01:41 +00:00
8aef9b1c13
check for twin nodes only in direction of the surface normal. this solves a problem with models of single element thickness and periodic fluxes.
2011-02-24 09:26:30 +00:00
fed53a7e60
in subroutine "mesh_marc_get_mpieOptions": need at least one additional string after keywords "$mpie" and "periodic", not two as it used to be
2011-02-23 12:33:51 +00:00
430966eb20
sorry, forgot to commit the numerics.f90, too
2011-02-23 12:30:52 +00:00
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"
2011-02-23 08:29:51 +00:00
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)
2011-02-23 08:08:06 +00:00
8dd1a694a3
extended IO to cope with different name for solverJob and Model
...
polishing, added error codes
added FFTW library files
2011-02-21 14:37:38 +00:00
b21c153761
restructured constitutive_nonlocal_compatibility and incorporated the "surfaceTransmissivity" in the compatibility calculation
2011-02-16 16:38:18 +00:00
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.
2011-02-16 16:35:38 +00:00
8f626c8989
Now 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!!!
2011-02-16 16:23:08 +00:00
89ed6f5f66
Did not compile. Fixed missing ,
2011-02-15 12:29:01 +00:00
c39e04644b
Sorting of hex slip-systems done in update 659 was missing. Included again.
2011-02-15 12:21:54 +00:00
438bf95105
added 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
2011-02-14 17:21:31 +00:00
171f5375f1
shearrate output now with sign (no |abs| as before)
2011-02-11 10:19:41 +00:00
3d7fad6ba9
mpie_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.
...
postprocessing: renamed name of python/f2py modul from "reconstruct" to "postprocessingMath", added some numerical operations to use for postprocessing.
2011-02-09 17:47:28 +00:00
803e1a8c05
In 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(:,:).
2011-02-09 13:12:46 +00:00
466ced911e
* check convergence for nonlocals also in stiffness calculation mode
2011-02-09 08:39:07 +00:00
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!
2011-02-09 08:28:47 +00:00
17a6c0e59b
enhanced comment on dislocation kinetics
2011-02-08 12:48:50 +00:00
6731659aa4
added hint that extreme values of stress and tangent belong to previous cycle
2011-02-08 12:43:30 +00:00
0a30fb739a
now reporting extreme values of stress and tangent, that were returned to marc, in output file
2011-02-08 10:25:51 +00:00
7a7ca1aab7
mpie_spectral and mpie_interface: changed workingdir from pathToLoadFile to pathToGeomFile
...
mpie_spectral and numerics: added switch to prevent pre calculation of gamma_hat. slower, but saves memory
3Dvisualize: started to add support for gmsh (not fully working yet)
reconstruct: new version of f2py/Fortran subroutines for output of results from spectral method
2011-02-07 14:35:42 +00:00
faba13f7fd
* now able to introduce some scatter in the initial single dislocation density (only for nonlocal model!) ; setting the parameter "rhoSglScatter" to a positive value generates a gauss distribution for the dislocation density with standard deviation equal to "rhoSglScatter"
...
* dislocation stress calculation is only done for nonlocal constitution ("/nonlocal/" keyword is present in material.config)
2011-02-04 15:41:32 +00:00
b362755457
grain orientation shuffling was skipped for /elemhomo/. no clue why, now dropped...
2011-02-03 17:25:58 +00:00
fec2c14a4e
removed hard-coded tolerances from mpie_spectral and put them to numerics/IO
...
removed storage of full cauchy stres field from mpie_spectral.f90, only average is stored now
added cauchy stress and von mises equivalent calculation to spectral post.
2011-01-31 17:07:42 +00:00
b72d75ed05
corrected indices in flux density output
2011-01-28 12:49:58 +00:00
c00a42a9a4
fixed last remaining tensor outputs to their transposed versions
2011-01-27 14:29:06 +00:00
aa7d5afd94
changed ifort options back to previous ones as -fast is not working
2011-01-26 13:33:44 +00:00
05d4d5fef2
* kinetics now according to Mohles; enables solid solution hardening
...
* changed example in material.config accordingly
2011-01-26 10:17:42 +00:00
6ac06af517
sorry, forgot to alter the most important part: post_results..! eventually, the t16 has defgrads and friends in correct format: du_i/dx_j is listed as linear array (1,1),(1,2),(1,3),(2,1),...
2011-01-24 16:23:37 +00:00
4ee40df5ba
transposed writing into t16 file (and screen/debug) did not work with just stating array(1:3,:)...
...
fixed by using math_transpose3x3 (and similar) instead.
2011-01-24 15:02:19 +00:00
794aeb253b
crystallite and CPFEM now print/output tensors in ij notation (used to be, particularly in the t16, ji, i.e. transposed...)
...
math reports random seed used for rnd() generation
2011-01-20 19:25:45 +00:00
b546d3bb91
some corrections regarding Multicore-support on mpie_spectral.f90 and makefile
...
slight changes on mpie_spectral.f90 regarding file output
2011-01-12 17:02:42 +00:00
062c6b3d54
* allow zero cutoff radius
...
* added missing error codes
2011-01-12 12:36:48 +00:00
d92f992eb7
if nonlocal ip has local neighbor assume perfect transparency of the (perhaps virtual) phase boundary. this is only a quick fix! need to rethink the way constitutive and real boundaries are treated!
2011-01-11 18:00:29 +00:00
6a3dd78259
fluxes are disabled when keyword "/nonlocal/" is omitted in the material.config (this is a way to use the constitutive_nonlocal model as a local dislocation based model)
2011-01-11 14:55:36 +00:00
ccc6aac10b
did not delete mpie_spectral2.f90 properly
2011-01-07 13:11:32 +00:00
71fb5eedf2
changed makefile: material subroutines are now compiled using openmp, enabled linux threads for fftw3 (also did some slight changes to enable that in the fortran source code)
...
renamed mpie_spectral2.f90 to mpie_spectral2d.f90 (testing file, not properly working at the moment)
changed file extension and variable names in mpie_spectral.f90 and mpie_spectral_interface.f90 from "mesh" to "geom". Removed direct output from mpie_spectral.f90, all output is now base on materialpoint_results(:,1,:)
2011-01-07 12:56:45 +00:00
9447330ab0
ATTENTION repository has been moved to new server!
...
Now at musvn1:/home/svn/repos/cpfem
2011-01-05 16:11:51 +00:00
096a4803b0
some polishing concerning openMP initialization
...
number of threads is now written to *.out file
2010-12-02 11:04:29 +00:00
36f27edd75
1. Adding more comments
...
2. Changing 'subroutine' into 'function' for optimizing the computation time
2010-11-26 11:50:20 +00:00
6b199e7830
allow 2digit values to be read in for MPIE_NUM_THREADS
2010-11-25 10:34:55 +00:00
ea6a6c9503
* environment variable MPIE_NUM_THREADS now determines number of threads used in parallel regions
...
* default value of the OMP_NUM_THREADS variable has to be restored at the end of mpie subroutine, since marc also seems to use and change(!) this
* usage: "export MPIE_NUM_THREADS=<number of threads>" to set variable in shell, then restart mentat and compile with option 3 (at the moment this does only work on ws 6, since all other workstations use compiler option "-save"; this puts all local variables by default in static memory, which is a killer for parallelization!)
2010-11-19 17:45:27 +00:00
b9b226548c
* enclose wall time measurement for LpAndItsTangent and count of leapfrogbreaks by CRITICAL construct
...
* better use SINGLE (having an implicit barrier at the end) instead of MASTER construct
* deleted all explicit BARRIERs after do loops since parallel loop construct implies barrier at the end
2010-11-19 17:29:29 +00:00
15728d3305
enclose wall time measurement for dotState and dotTemperature function by CRITICAL construct
2010-11-19 17:19:03 +00:00
Philip Eisenlohr