1! CalculiX - A 3-dimensional finite element program 2! Copyright (C) 1998-2021 Guido Dhondt 3! 4! This program is free software; you can redistribute it and/or 5! modify it under the terms of the GNU General Public License as 6! published by the Free Software Foundation(version 2); 7! 8! 9! This program is distributed in the hope that it will be useful, 10! but WITHOUT ANY WARRANTY; without even the implied warranty of 11! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 12! GNU General Public License for more details. 13! 14! You should have received a copy of the GNU General Public License 15! along with this program; if not, write to the Free Software 16! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 17! 18 subroutine materialdata_dvi(shcon,nshcon,imat,dvi,t1l,ntmat_, 19 & ithermal) 20! 21 implicit none 22! 23! determines the dynamic viscosity 24! 25 integer imat,ntmat_,id,nshcon(*),four,ithermal 26! 27 real*8 t1l,shcon(0:3,ntmat_,*),dvi 28! 29 four=4 30! 31 if(ithermal.eq.0) then 32 dvi=shcon(2,1,imat) 33 else 34 call ident2(shcon(0,1,imat),t1l,nshcon(imat),four,id) 35 if(nshcon(imat).eq.0) then 36 continue 37 elseif(nshcon(imat).eq.1) then 38 dvi=shcon(2,1,imat) 39 elseif(id.eq.0) then 40 dvi=shcon(2,1,imat) 41 elseif(id.eq.nshcon(imat)) then 42 dvi=shcon(2,id,imat) 43 else 44 dvi=shcon(2,id,imat)+ 45 & (shcon(2,id+1,imat)-shcon(2,id,imat))* 46 & (t1l-shcon(0,id,imat))/ 47 & (shcon(0,id+1,imat)-shcon(0,id,imat)) 48 endif 49 endif 50! 51 return 52 end 53 54 55 56 57 58 59 60