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_cond(imat,ntmat_,t1l,cocon,ncocon,cond) 19! 20 implicit none 21! 22! determines the thermal conductivity 23! 24 integer imat,ntmat_,id,ncocon(2,*),ncoconst,seven 25! 26 real*8 t1l,cocon(0:6,ntmat_,*),cond 27! 28 seven=7 29! 30! calculating the conductivity coefficients 31! 32 ncoconst=ncocon(1,imat) 33 if(ncoconst.eq.0) then 34 write(*,*) '*ERROR in materialdata_cond' 35 write(*,*) 36 & ' fluid conductivity is lacking' 37 call exit(201) 38 elseif(ncoconst.gt.1) then 39 write(*,*) '*ERROR in materialdata_cond' 40 write(*,*) 41 & ' conductivity for fluids must be isotropic' 42 call exit(201) 43 endif 44! 45 call ident2(cocon(0,1,imat),t1l,ncocon(2,imat),seven,id) 46 if(ncocon(2,imat).eq.0) then 47 cond=0.d0 48 continue 49 elseif(ncocon(2,imat).eq.1) then 50 cond=cocon(1,1,imat) 51 elseif(id.eq.0) then 52 cond=cocon(1,1,imat) 53 elseif(id.eq.ncocon(2,imat)) then 54 cond=cocon(1,id,imat) 55 else 56 cond=(cocon(1,id,imat)+ 57 & (cocon(1,id+1,imat)-cocon(1,id,imat))* 58 & (t1l-cocon(0,id,imat))/ 59 & (cocon(0,id+1,imat)-cocon(0,id,imat))) 60 & 61 endif 62! 63 return 64 end 65 66 67 68 69 70 71 72