Commit Graph

7436 Commits

Author SHA1 Message Date
Luc Hantcherli b476aa4028 CPFEM.f90 contains the following modifications:
-Temperature is now defined as a CPFEM_Temperature, and filled for each new call of CPFEM_general according the value output by MARC
-constitutive_Microstructure(state_new,CPFEM_Temperature,grain,CPFEM_in,cp_en) is used in the state loop and before entering in the stress loop. This subroutine computes all parameters,arrays or matrices required from the state, i.e. forest or parallel or mobile dislocation densities from the statistically stored dislocation densities. Its definition is made in module constitutive.
2007-10-16 08:33:59 +00:00
Luc Hantcherli ab0fddf518 Module crystal contains a so-called interaction type matrix for each crystal structure. This matrix helps for computing the hardening matrix in case of a phenomenological constitutive or the dislocation interaction intensity matrix in case of a dislocation-based constitutive.
Definition can be discuss.
Corresponding changes should be carried on constitutive.f90
2007-10-16 08:23:38 +00:00
Luc Hantcherli ccafbf9a0f Temperature T(1) is now transmitted to the CPFEM_general subroutine. Corresponding changes should be carried on CPFEM.f90 2007-10-16 08:18:57 +00:00
Philip Eisenlohr e3dcaa6915 outer crystal system loop now runs up to _MaxCrystalStructure (used to be hardcoded 3) 2007-10-15 14:36:25 +00:00
Philip Eisenlohr 5a82f22681 resubmit of IO_countContinousIntValues (did not check fortran syntax -- should be fine now) 2007-10-15 13:57:50 +00:00
Philip Eisenlohr 08a9985179 added function IO_countContinousIntValues
to read number of items in set definitions:
123 to 423
or something like
1 3 5 6 7 9 10 c
12 14
2007-10-15 13:55:52 +00:00
Luc Hantcherli 80e76fcd3f Added in an extra folder: modified subroutines with dislocation as state variables (version 0.1) 2007-10-15 12:05:47 +00:00
Luc Hantcherli 73694790f8 2007-10-15 12:03:52 +00:00
Franz Roters 532d3c2a16 added use crystal statement before call of crystal_init 2007-10-15 10:16:51 +00:00
Luc Hantcherli 64f3f6ccee crystal_init() included in CPFEM.f90 2007-10-15 09:30:57 +00:00
Philip Eisenlohr b0ef9e4a94 added element types
11 (4node 2D linear full) and
27 (8node 2D quadratic full)
2007-10-12 13:48:29 +00:00
Franz Roters 4ce8270415 corrected some variable names 2007-10-11 11:42:37 +00:00
Luc Hantcherli a6fef91c6b Extraction crystal structure information from constitutive.f90. Creation of crystal.f90 2007-10-11 11:36:09 +00:00
William Counts 7bc27e02ca keyword "table" appeared more than once -- restricted interpretation to header 2007-09-28 14:56:26 +00:00
Philip Eisenlohr 892f33591c adapted to new default table style in 2007R1
hypoelastic part of input file needs skipping an additional line if new style tables are active
2007-09-28 14:16:03 +00:00
Franz Roters 28fc9cc221 put recalculation of elastic stress back in (was accidentally deleted)
changed check of residuum improvement for small residuum
2007-08-13 07:12:24 +00:00
Franz Roters 85df3b95ee changed restriction of dTstar, its now based on the residuum, if it does not improve 0.5*dTstar is used instead 2007-08-07 08:02:31 +00:00
Franz Roters a5f49a0578 changed first elastic guess for small Tstar
changed residuum to only consider deviatoric part of Tstar
changed scaling when dTstar becomes too large to proportionally scaling all components of dTstar
2007-07-20 13:32:44 +00:00
Franz Roters a1c3df84ab changed quite a few things, mainly using last result as initial guess now, thinks work a lot better now but still not perfect 2007-06-06 14:38:06 +00:00
Franz Roters 0efe7756bc included fix in constitutive_post_results to handle case of constitutive_Nresults==0 2007-05-22 06:36:52 +00:00
Franz Roters 6f3f5ee23c corrected error in deMandeling of jacoby matrix 2007-05-22 06:34:58 +00:00
Franz Roters 09c47f15e3 changed constitutive to open file using IO_open_file and use fileunit 200 instead of 1
changed constitutive to skip empty lines when parsing files
activated call quit in IO_error
2007-05-16 14:36:03 +00:00
Franz Roters e1ae35c5e4 corrected typo in constitutive_Sslip_v
corrected demandle of jacoby in hypela2
2007-05-15 08:47:44 +00:00
Franz Roters 49886a072c corrected CPFEM_stressIP to use FFN1 at the end of the time increment
calculate for ncycle==0 and not only inc==0 in hypela2
2007-05-14 12:13:36 +00:00
Franz Roters 6aece52d1c changed Calculation of the tangent of Lp to always use state(i) and not constitutive_state_new 2007-05-11 13:22:16 +00:00
Franz Roters e5749cbff0 changed scatter check for gauss components to if(noise==0) 2007-04-26 14:46:54 +00:00
Philip Eisenlohr 09aad5ccb4 changed _Xvalue() return codes on error to huge() 2007-04-26 12:40:06 +00:00
Franz Roters 11cc6d1ee1 math_sampleGaussOri did not treat the special case of zero noise (scatter) 2007-04-26 11:46:58 +00:00
Franz Roters 6fd473eaa7 scatter was not transformed to rad when read from material file 2007-04-26 11:42:16 +00:00
Franz Roters 1fd722e4cf changed format for reading float from('F') to * in IO_floatValue 2007-04-25 14:38:22 +00:00
Franz Roters 033a6a03b3 determine number of grains as texture_Ngrains(mesh_element(4,cp_en)) instead of constitutive_Ngrains (which is not used anymore) in CPFEM_stressIP 2007-04-25 13:58:10 +00:00
Franz Roters 1747c0cc98 corrected calculation of constitutive_maxNstatevars 2007-04-25 13:50:28 +00:00
Franz Roters 03f244995e corrected mesh_build_elements (unit) so that state variables are read correctly from input file 2007-04-25 07:33:24 +00:00
Franz Roters 3a4ec4c2d1 removed several unused variables 2007-04-24 06:49:13 +00:00
Franz Roters c0ebb4d846 reordered subroutine calls in mesh_init to allocate nodemap and elementmap BEFORE they are used 2007-04-24 06:29:22 +00:00
Philip Eisenlohr 83d2dea8c6 CPFEM_general call now uses CPFEM_en instead of cp_en -- the required mapping is not in place on first call! 2007-04-23 13:25:09 +00:00
Philip Eisenlohr 471d53e3ba _general needs CPFEM_en NOT cp_en..!
stress loop counter was not been reset
matmul(A,B) now replaced by stored AB
2007-04-23 13:23:03 +00:00
Philip Eisenlohr 9e0b9a7096 removed syntax error in _LpAndItsTangent (enforcing symmetry was flawed...) 2007-04-17 08:02:27 +00:00
Luc Hantcherli 0f6ff82f30 Checked indices for Jacobi calculation 2007-04-17 07:58:53 +00:00
Luc Hantcherli c76ab95b21 Added in LpandIstTangent 0.5*(transpose+untranspose) 2007-04-17 07:53:06 +00:00
Philip Eisenlohr b26704c23f added missing use statement in _init 2007-04-13 15:37:10 +00:00
Philip Eisenlohr 4718abfc70 agreed on wrong calc of scaled state residual 2007-04-13 14:22:55 +00:00
Philip Eisenlohr 09495e852d added necessary use statements 2007-04-13 14:20:59 +00:00
Philip Eisenlohr bbfb86fb8f added necessary use statements 2007-04-13 14:20:33 +00:00
Luc Hantcherli d90162dc38 Corrected Rstress(i)/state_new(i) into Rstate(i)/state_new(i) 2007-04-13 11:30:49 +00:00
Luc Hantcherli eb292917b9 untied interlinked constitutive and CPFEM modules
now constitutive stores EulerAngles and CPFEM generates Fp_old from these
2007-04-11 15:28:46 +00:00
Luc Hantcherli cb4a85319d corrected small syntax errors 2007-04-11 14:55:06 +00:00
Luc Hantcherli 3c5d490148 added _plot_results function to tell about constitutive results (state + current shear per slip system) 2007-04-11 14:52:47 +00:00
Luc Hantcherli f9f3e2bd9b result processing now in stressIP
i.e. no extra polar decompositions performed any more
2007-04-11 14:51:49 +00:00
Philip Eisenlohr 9704a4e83c syntax polishing 2007-04-11 13:21:22 +00:00