1! 2! CalculiX - A 3-dimensional finite element program 3! Copyright (C) 1998-2021 Guido Dhondt 4! 5! This program is free software; you can redistribute it and/or 6! modify it under the terms of the GNU General Public License as 7! published by the Free Software Foundation(version 2); 8! 9! 10! This program is distributed in the hope that it will be useful, 11! but WITHOUT ANY WARRANTY; without even the implied warranty of 12! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 13! GNU General Public License for more details. 14! 15! You should have received a copy of the GNU General Public License 16! along with this program; if not, write to the Free Software 17! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 18! 19 subroutine materialdata_cfd1(nef,vel,shcon,nshcon,ielmatf, 20 & ntmat_,mi,cvel,physcon,ithermal,umel,rhcon,nrhcon,nefa,nefb) 21! 22! calculation of material properties at elements centers 23! (incompressible fluids) 24! 25 implicit none 26! 27 integer nef,i,imat,ntmat_,mi(*),ielmatf(mi(3),*),ithermal(*), 28 & nshcon(2,*),nrhcon(*),nefa,nefb 29! 30 real*8 t1l,vel(nef,0:7),shcon(0:3,ntmat_,*),cvel(*), 31 & physcon(*),umel(*),rhcon(0:1,ntmat_,*) 32! 33! 34! 35 do i=nefa,nefb 36 t1l=vel(i,0) 37 imat=ielmatf(1,i) 38! 39! density 40! 41 call materialdata_rho(rhcon,nrhcon,imat,vel(i,5),t1l,ntmat_, 42 & ithermal) 43! 44! heat capacity at constant volume 45! (for liquids: =heat capacity at constant pressure) 46! 47 call materialdata_cp_sec(imat,ntmat_,t1l,shcon,nshcon,cvel(i), 48 & physcon) 49! 50! dynamic viscosity 51! 52 call materialdata_dvi(shcon,nshcon,imat,umel(i),t1l,ntmat_, 53 & ithermal) 54 enddo 55! 56 return 57 end 58