Philip Eisenlohr
ef4fc9d0ee
added <crystallite> output "texture" to list of examples
2011-04-12 15:33:29 +00:00
Franz Roters
c1b8391110
changed enconding of all source files to UTF-8 without BOM (signature) Codepage 65001
2011-04-07 07:20:28 +00:00
Nan Jia
e4cbdc5743
now also works for 2D quadratic triangular element type 125
2011-04-06 13:11:25 +00:00
Christoph Kords
4d0fefd231
security factor of 1.2 for CFL condition in constitutive_nonlocal_dotState
2011-04-06 11:47:24 +00:00
Philip Eisenlohr
a989aa979e
corrected increment offsetting for spectral file type. should now be consistent with new (0+N) mpie_spectral output scheme
2011-04-06 10:00:19 +00:00
Philip Eisenlohr
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
Franz Roters
7c7c929455
added comments on origin of files
2011-04-06 09:40:39 +00:00
Christoph Kords
4ea1fe436b
added check for CFL condition in constitutive_nonlocal_dotState
2011-04-06 09:07:36 +00:00
Philip Eisenlohr
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
Philip Eisenlohr
a200bae6ce
two-dimensional plane strain elements with 3 IPs (internal element type 10)
2011-04-04 16:10:34 +00:00
Philip Eisenlohr
fffff202b0
2D mesh option now also alters corresponding job option
2011-04-04 15:47:09 +00:00
Franz Roters
fcdb805225
added copyright text to all f90 (free) format files
2011-04-04 14:09:54 +00:00
Franz Roters
ed47c09424
added file containing copyright text to be added to each file
...
ATTENTION uses f90 comment marker (!)
2011-04-04 14:07:42 +00:00
Philip Eisenlohr
73accf79b6
added option for 2D mesh (3D extrusion step is skipped)
2011-04-04 12:12:01 +00:00
Christoph Kords
6ca8942f62
choose correct output format for debug_verbosity: integer instead of logical
2011-04-04 08:38:55 +00:00
Christoph Kords
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
Alankar Alankar
06fc83ac14
Final version of titanmod
2011-03-31 09:21:43 +00:00
Philip Eisenlohr
a546f8dbd2
changed wording of input queries to (hopefully) clearer statements
2011-03-30 15:48:18 +00:00
Franz Roters
fc471d73f3
added Arun to list of authors
2011-03-30 14:39:28 +00:00
Franz Roters
dc8961f728
added first files for release under GPL 3.0
2011-03-30 14:38:33 +00:00
Christoph Kords
c36968fc68
distance for periodic images in internal stress calculation was incorrect
2011-03-30 11:57:41 +00:00
Franz Roters
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
Franz Roters
e9e6abbdaf
introduced variables for compiler options
...
added -fpp switch to invoke precompiler
2011-03-29 07:55:38 +00:00
Christoph Kords
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
Christoph Kords
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
Franz Roters
1b3dba4113
added -fpp switch to all compiler calls to explicitly use the precompiler
...
removed -save from OpenMP compiler calls to use OpenMP auto behaviour
2011-03-28 12:18:15 +00:00
Philip Eisenlohr
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
Franz Roters
fa323f3f0b
removed division by zero from debug_reset(), uses huge() now
2011-03-24 07:38:56 +00:00
Franz Roters
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
Franz Roters
8467c91dab
moved setup scripts back to subdirs where they are needed
2011-03-22 15:42:53 +00:00
Christoph Kords
31a570768b
new error 237: singularity in internal stress calculation
2011-03-22 13:58:42 +00:00
Nan Jia
7be2edb10e
added "texture
...
" as possible crystallite output => reports ID for later discrimination...
2011-03-22 13:40:27 +00:00
Christoph Kords
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
Christoph Kords
9b1e1e35ea
dotState
2011-03-21 15:18:09 +00:00
Christoph Kords
462eda6736
split parallel region in integrateStateFPI into two
2011-03-21 14:30:10 +00:00
Christoph Kords
5ee73dee72
age dotState before resetting it to zero, otherwise the statedamping does not work.
2011-03-21 13:06:11 +00:00
Franz Roters
8bef254be7
changed name to include *.py ending to make Windows file association work again
2011-03-21 11:58:24 +00:00
Christoph Kords
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
Christoph Kords
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
Christoph Kords
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
Christoph Kords
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
Philip Eisenlohr
6ac2b4cf88
deleted useless storage of data resulting in partial speed-up.
...
large datasets (>1 mio elements) cannot be processed at the moment---investigation required, why each element eats about 20kB of memory?
2011-03-11 16:32:01 +00:00
Franz Roters
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
Philip Eisenlohr
1f0569230b
new installation/setup dir to hold all top-level make_*-scripts
2011-03-11 09:20:37 +00:00
Philip Eisenlohr
41c0bf098c
starting to substructure manual
2011-03-11 09:19:31 +00:00
Franz Roters
bdc17f7d72
extended comment on modification of abaqus_v6.env
2011-03-11 09:05:09 +00:00
Philip Eisenlohr
d860015d57
first shot at reviving the manual...
2011-03-11 08:43:29 +00:00
Philip Eisenlohr
621c1198a4
Marc/Mentat related things are only executed for filetype marc. this has been too lax before, causing trouble.
...
filetype is now lowercased right from the start, keeping effort low downstream...
2011-03-10 09:45:57 +00:00
Franz Roters
6475235fb4
included "-heap-arrays 500000000" in remark on how to change compile command
2011-03-08 08:18:04 +00:00
Franz Roters
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