Martin Diehl
ddfc40afe6
improve restart test, deletion of results was messed up. fixed restart for BasicPETSc, AL and Polarization currently don't run when restarting
2013-09-05 12:07:58 +00:00
Martin Diehl
7f3c24e1f1
fixed some buggy last minute changes in the last commit
2013-08-08 09:13:29 +00:00
Martin Diehl
ac92b90e0b
fixed al and pol solver, now checking for div(p) = curl(f) = 0
2013-08-07 17:20:05 +00:00
Martin Diehl
81531097f1
changes on convergence tolerances of AL and Polarisation, switched back to immediate correction of stress bc but only when last two average stresses are close to each other (cosine decay)
2013-07-30 15:32:55 +00:00
Pratheek Shanthraj
03d8f14a98
introduced dynamics for high strain rate simulations. so far only works for materials with homogeneous density. the density is read in as an optional argument in the load file using the keyword 'den' or 'density'. if density > 0 then inertial terms are added to the stress equilibrium. an implicit euler time discretisation is used to calculate the inertial terms.
...
dynamic problems are length scale dependent so make sure the geometry size is meaningful (1 m^3 or N m^3 for example is too big for a RVE!).
not specifying a density (density = 0) will perform a quasi static simulation as before.
2013-07-26 16:25:37 +00:00
Martin Diehl
4407edd802
checked convergence reporting for AL and BasicPETSc, seems to be ok now.
2013-07-08 15:48:13 +00:00
Franz Roters
a78f91311e
added svn properties: line ending: LF and keywords: ID
2013-06-14 09:49:33 +00:00
Ding Wang
ce9f92a884
polishing output
2013-06-11 10:27:38 +00:00
Martin Diehl
872f6a9d90
introduced possibility so specify deformation gradient aim at end of load case, rate will be calculated using difference between start of load case and aim. Needed for cyclic loading.
...
Use keyword "f" for this behavior, don't use it as short name for freq any more!
2013-05-13 09:44:23 +00:00
Martin Diehl
85d4a37d95
moved public data res,size and homog from mesh to DAMASK_spectral_utilities (as grid and geomSize)
2013-05-08 15:52:29 +00:00
Martin Diehl
39a70e8a19
fixed bug in coordinate calculation when no average F was given (scaled wrong by ncp_elems**2)
...
other changes: just polishing + some more comments
2013-03-27 12:28:55 +00:00
Franz Roters
02f092492e
updated copyright header
2013-03-22 17:35:05 +00:00
Pratheek Shanthraj
be655ae536
Changed AL solver errors to absolute. Stress boundary error now max of relative and absolute error
2013-03-22 14:46:55 +00:00
Martin Diehl
1779510dda
remove some unused variables, set default value for divergence calculation to 1.e-5 and made calculation independent of dim and res
2013-03-12 18:48:28 +00:00
Pratheek Shanthraj
1f4d7c2ca4
changed reference stiffness to min max avg from volume avg. using terminally ill information in convergence check
2013-03-06 14:31:13 +00:00
Pratheek Shanthraj
07ff1f7c93
reporting current time for most module inits
2013-02-25 16:34:59 +00:00
Martin Diehl
cc9eb685fe
fixed wrong setting of coordinates in spectral solvers and removed old python scripts for setup of code and processing
2013-02-05 12:31:44 +00:00
Martin Diehl
96577b18fb
introduced error code for run_test.py
2013-01-23 18:33:46 +00:00
Martin Diehl
9ee8108b6b
added doxygen documentation to material.f90 and marked read-only quantities as protected where possible, removed substituted "call flush" by "flush"
2013-01-18 11:30:52 +00:00
Martin Diehl
c6a79d2b3d
divergence_correction for basic solver variants has now 3 possibilities:
...
0: uncorrected, slope per sidelength (physical dimension) e = res/dim
1: corrected by sidelength, slope per unitlength e = res/1
2: corrected such that distance between FPs e = 1
alway regarding the medium length of x,y,z direction
2013-01-16 10:40:53 +00:00
Martin Diehl
645b2605e1
corrected reporting to statistics file (*.sta)
2013-01-11 10:40:16 +00:00
Martin Diehl
6e3e06bed7
corrected PETSc macro expansion
2013-01-10 15:36:55 +00:00
Martin Diehl
563b1f5e4b
added some warning and explicit size of arrays
2013-01-10 13:33:43 +00:00
Martin Diehl
0d5e91ac87
corrected reporting of basic PETSc variant and unified reporting of all spectral solvers
...
improved warning and error in IO, now able to report correctly elements up to 9 digits
2013-01-09 22:19:32 +00:00
Martin Diehl
55b88e47b7
fixed wrong temperature when using spectral solver
2013-01-09 18:08:08 +00:00
Martin Diehl
ce7a0571fd
fixed bug in forwarding fields for AL solver
2013-01-08 10:12:03 +00:00
Martin Diehl
7a43d1b6ad
added funtionality to specify working directory to spectral solver
2013-01-02 17:02:12 +00:00
Martin Diehl
1baf8dea5d
load case rotation no working for Basic PETSc solver
2012-12-15 23:52:06 +00:00
Martin Diehl
b1f28ba5e7
made restart working for all solvers and added corresponding test
...
corrected error handling for FFT based function in case of odd numbers etc.
2012-12-14 17:30:22 +00:00
Martin Diehl
7b87987751
introduced option for regridding to numerics.f90, working as follows:
...
first cut back is tried if material point model oder BVP solver does not converge.
If no regridding is enabled after max cut back, in case of non-converged material point the simulation stops and in case of non-converged BVP solver it continues.
set regridMode to 2 to enable regridding if BVP solver OR materialPoint model do not converge,
set regridMode to 1 to enable regridding if materialPoint model do not converge, non-converged BVP solver will be ignored as in the standard case.
For regridding, the load case need to have a restart freq set.
enabled restarting for Basic PETSc variant
2012-12-14 15:18:04 +00:00
Martin Diehl
d559653722
small improvements on tests
2012-12-13 20:20:04 +00:00
Martin Diehl
5ea0139678
simplified the rate calculation interface in DAMASK_spectral_utilities.f90, corrected bug in rotation parsing and added debug option "rotation" to spectral solver options. This will show the current average stress and deformation additionally in lab coordinate system
2012-11-28 15:04:05 +00:00
Martin Diehl
7e683ca7c3
removed coordinates from call to CPFEM_general and introduced direct storage of coordinates to mesh_ipCoordinates
2012-11-14 14:38:10 +00:00
Martin Diehl
70c4e11742
added comments and structured the code. Temperature is not longer stored for each point but only to simplify restart behavior as long as it is not fully supported by the constitutive models
2012-11-12 14:14:39 +00:00
Martin Diehl
b1c3c57412
some more finetuning, still not running on gfortran 4.5.0
2012-11-08 21:33:58 +00:00
Martin Diehl
ee5903c767
replaced real 'guessmode' with logical 'guess', caused trouble on gfortran 6.3 with ultra optimization
2012-11-08 19:32:00 +00:00
Martin Diehl
bbcffa668b
some small changes in including PETSc *.h90 files to prevent warnings.
...
It seems that a bug in old gfortran 4.5.0 causes the problems, at least with 4.7.3 it seems to work after fixing the bugs in IO.f90
2012-11-07 13:11:41 +00:00
Martin Diehl
3ada4897fb
reordered PETSc initialization (now first part done in interface, setting of parameters done in numerics), removed unnecessary includes for PETSc
...
still not running with gfortran, use at own risk!
2012-11-06 16:00:51 +00:00
Martin Diehl
13b55275b1
documented utilities and structured, worked on the restart capabilities of the new basic solver
2012-10-24 11:31:40 +00:00
Martin Diehl
b9f97ba5da
unified naming scheme, fixed a bug in the new basic scheme, and added a statistic file similar to abaqus to store the information on needed cut backs and iterations for each inc
2012-10-19 08:44:21 +00:00