1 : C LINEAR SOURCE CHAT 05/01/13 01:18:32 5004 2 : c--------------------------------------------------------------------- 3 : c 4 : SUBROUTINE LINEAR (ISTEG, ITERA, NDIME, NDIM1, NSTRE, NSTR1, 5 : . ntemp, 6 : . dmatx, poiss, stran, stres, strin, 7 : . stype) 8 : c 9 : c===================================================================== 10 : c = 11 : c Evaluates total and incremental stresses assuming linear and = 12 : c elastic behaviour. = 13 : c = 14 : c Note: stran = (Sxx, Syy, Szz, Sxy, Sxz, Syz) NSTRE=6 = 15 : c stran = (Srr, Szz, S00, Srz) NSTRE=4 = 16 : c stran = (Sxx, Syy, Szz, Sxy) NSTRE=3 = 17 : c = 18 : c stres = (Sxx, Syy, Szz, Sxy, Sxz, Syz) NSTRE=6 = 19 : c stres = (Srr, Szz, S00, Srz) NSTRE=4 = 20 : c stres = (Sxx, Syy, Szz, Sxy) NSTRE=3 = 21 : c = 22 : c Author: Rui Faria Version: 93.04.20 = 23 : c -------- = 24 : c===================================================================== 25 : IMPLICIT INTEGER(I-N) 26 : integer*4 isteg, itera, ndime, ndim1, nstre, nstr1, ntemp 27 : real*8 poiss, strin, stran (nstr1), stres (nstr1), 28 : . dmatx (nstre,nstre) 29 : c 30 : character*6 stype 31 : c=====---------------------------------------------------------------- 32 : c compute the incremental stresses | 33 : c=====---------------------------------------------------------------- 33 : call opmult (ndime, ndim1, nstre, nstr1, 34 : . dmatx, stran, stres, 35 : . stype) 36 : c 37 : if (stype .eq. 'PSTRES') stres (3) = 0.0 d0 38 : if (stype .eq. 'PSTRAI') stres (3) = poiss * 39 : . (stres (1) + stres (2)) 40 : c=====---------------------------------------------------------------- 41 : c strin - stresses due to temperature | 42 : c=====---------------------------------------------------------------- 43 : c if (ntemp .eq. 1 .and. isteg .eq. 1 .and. itera .eq. 1) then 44 : c stres (1) = stres (1) + strin 45 : c stres (2) = stres (2) + strin 46 : c stres (3) = stres (3) + strin 47 : c end if 48 : c 49 : return 50 : end 51 : 52 :
© Cast3M 2003 - All rights reserved.
Disclaimer