/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | networkstiff.f | 19 subroutine networkstiff(voldl,s,imat,konl,mi,ntmat_,shcon, argument 29 real*8 voldl(0:mi(2),9),s(60,60),gastemp,shcon(0:3,ntmat_,*), local 34 gastemp=(voldl(0,1)+voldl(0,3))/2.d0 39 if(voldl(1,2).gt.0.d0) then 40 s(3,1)=-cp*voldl(1,2) 43 s(1,1)=-cp*voldl(1,2)
|
H A D | calcstressheatfluxfem.f | 38 & shcon(0:3,ntmat_,*),voldl(0:mi(2),20),vold(0:mi(2),*),temp, 78 voldl(i2,i1)=vold(i2,konl(i1)) 88 voldl(i2,i1)=vold(i2,konl(i1)) 123 vkl(j1,k1)=vkl(j1,k1)+shp(k1,i1)*voldl(j1,i1) 136 dtem(j1)=dtem(j1)+shp(j1,i1)*voldl(0,i1) 146 temp=temp+shp(4,i1)*voldl(0,i1) 183 xkin=xkin+shp(4,i1)*voldl(5,i1) 184 xtuf=xtuf+shp(4,i1)*voldl(6,i1)
|
H A D | springstiff_n2f_th.f | 19 subroutine springstiff_n2f_th(xl,voldl,s,imat,elcon,nelcon, argument 35 & al(3),s(60,60),voldl(0:mi(2),10),pl(3,10),xn(3),dm, 57 pl(j,i)=xl(j,i)+voldl(j,i) 151 t1l=t1l+ratio(j)*voldl(0,j) 153 dtemp=t1l-voldl(0,nope) 154 tmean=(voldl(0,nope)+t1l)/2.d0 171 temp(1)=voldl(0,nope)
|
H A D | springdamp_f2f.f | 19 subroutine springdamp_f2f(xl,elas,voldl,s,imat,elcon, argument 33 & voldl(0:mi(2),19),pl(3,19),xn(3),c3,elcon(0:ncmat_,ntmat_,*), local 56 pl(j,i)=xl(j,i)+voldl(j,i) 64 pl(j,i)=xl(j,i)+voldl(j,i)+clearini(j,i-nopem,jfaces) 95 voldl(k,nopep)=0.d0 98 voldl(k,nopep)=voldl(k,nopep)+shp2s(4,j)*voldl(k,nopem+j)
|
H A D | springstiff_f2f.f | 19 subroutine springstiff_f2f(xl,elas,voldl,s,imat,elcon,nelcon, argument 37 & al(3),s(60,60),voldl(0:mi(2),19),pl(3,19),xn(3), local 66 pl(j,i)=xl(j,i)+voldl(j,i) 74 pl(j,i)=xl(j,i)+voldl(j,i)+clearini(j,i-nopem,jfaces) 102 voldl(k,nopep)=0.d0 105 voldl(k,nopep)=voldl(k,nopep)+shp2s(4,j)*voldl(k,nopem+j) 250 alnew(i)=voldl(i,nopep) 252 alnew(i)=alnew(i)-shp2m(4,j)*voldl(i,j)
|
H A D | e_c3d_v1rhs.f | 152 voldl(i2,i1)=vold(i2,konl(i1)) 203 temp=temp+shp(4,i1)*voldl(0,i1) 205 vel(j1)=vel(j1)+shp(4,i1)*voldl(j1,i1) 207 dtem(j1)=dtem(j1)+shp(j1,i1)*voldl(0,i1) 211 dpress(j1)=dpress(j1)+shp(j1,i1)*voldl(4,i1) 289 xkin=xkin+shp(4,i1)*voldl(5,i1) 290 xtuf=xtuf+shp(4,i1)*voldl(6,i1) 788 temp=temp+shp(4,i1)*voldl(0,i1) 790 vel(j1)=vel(j1)+shp(4,i1)*voldl(j1,i1) 861 xkin=xkin+shp(4,i1)*voldl(5,i1) [all …]
|
H A D | springstiff_f2f_th.f | 19 subroutine springstiff_f2f_th(xl,voldl,s,imat,elcon,nelcon, argument 39 & al(3),s(60,60),voldl(0:mi(2),19),pl(3,19),xn(3), local 66 pl(j,i)=xl(j,i)+voldl(j,i) 72 pl(j,i)=xl(j,i)+voldl(j,i)+clearini(j,i-nopem,jfaces) 110 t1ls=t1ls+shp2s(4,j)*voldl(0,nopem+j) 140 t1lm=t1lm+shp2m(4,j)*voldl(0,j)
|
H A D | dashdamp.f | 19 subroutine dashdamp(xl,elas,konl,voldl,s,imat,elcon,nelcon, argument 32 real*8 xl(3,9),elas(21),s(60,60),voldl(3,9),xn(3),dd, local 48 pl(j,i)=xl(j,i)+voldl(j,i)
|
H A D | e_damp.f | 49 & t0(*),t1(*),voldl(3,20),vold(0:mi(2),*),xload(2,*),omx, local 74 voldl(i2,i1)=vold(i2,konl(i1)) 101 call dashdamp(xl,elas,konl,voldl,s,imat,elcon,nelcon,
|
H A D | advecforc.f | 19 subroutine advecforc(nope,voldl,ithermal,xl,nelemload,nelemadvec, argument 41 real*8 tl2(9),voldl(0:mi(2),20),xl(3,9),sinktemp,xi,et,weight, local 64 tl2(i)=voldl(0,i) 73 xl2(j,i)=xl(j,i)+voldl(j,i)
|
H A D | printoutfacefem.f | 43 & voldl(0:mi(2),8),xl2(3,8),xsj2(3), 156 voldl(i2,i1)=vold(i2,konl(i1)) 268 temp=temp+shp(4,i1)*voldl(0,i1) 269 pres=pres+shp(4,i1)*voldl(4,i1) 272 vkl(j1,k1)=vkl(j1,k1)+shp(k1,i1)*voldl(j1,i1)
|
H A D | advecstiff.f | 19 subroutine advecstiff(nope,voldl,ithermal,xl,nelemload,nelemadvec, argument 41 real*8 tl2(9),voldl(0:mi(2),9),xl(3,9),sinktemp,xi,et,weight, local 64 tl2(i)=voldl(0,i) 73 xl2(j,i)=xl(j,i)+voldl(j,i)
|
H A D | springstiff_n2f.f | 19 subroutine springstiff_n2f(xl,elas,konl,voldl,s,imat,elcon,nelcon, argument 39 & al(3),s(60,60),voldl(0:mi(2),10),pl(3,10),xn(3),dm, 62 pl(j,i)=xl(j,i)+voldl(j,i) 140 overclosure=-dd-xn(1)*(voldl(1,2)-voldl(1,1)) 141 & -xn(2)*(voldl(2,2)-voldl(2,1)) 142 & -xn(3)*(voldl(3,2)-voldl(3,1)) 575 alnew(i)=voldl(i,nope) 577 alnew(i)=alnew(i)-ratio(j)*voldl(i,j)
|
H A D | printoutface.f | 43 & voldl(0:mi(2),20),cocon(0:6,ntmat_,*),xl2(3,8),xsj2(3), 234 voldl(i2,i1)=vold(i2,konl(i1)) 244 xl(j,i)=xl(j,i)+voldl(j,i) 440 temp=temp+shp(4,i1)*voldl(0,i1) 441 pres=pres+shp(4,i1)*voldl(4,i1) 445 & shp(k1,i1)*voldl(j1,i1) 489 temp=temp+shp(4,i1)*voldl(0,i1) 491 vkl(0,k1)=vkl(0,k1)+shp(k1,i1)*voldl(0,i1)
|
H A D | springdamp_n2f.f | 19 subroutine springdamp_n2f(xl,elas,voldl,s,imat,elcon, argument 33 & al(3),s(60,60),voldl(0:mi(2),10),pl(3,10),xn(3),dm, local 48 pl(j,i)=xl(j,i)+voldl(j,i) 279 alnew(i)=voldl(i,nope) 281 alnew(i)=alnew(i)-ratio(j)*voldl(i,j)
|
H A D | mafilldm.f | 54 & p2(3),ad(*),au(*),bodyf(3),voldl(0:mi(2),26),clearini(3,9,*), local 128 voldl(j,k)=vold(j,konl(k)) 147 call springdamp_n2f(xl,elas,voldl,s,imat,elcon, 154 call springdamp_f2f(xl,elas,voldl,s,imat,elcon,
|
H A D | e_c3d_th.f | 72 & voldl(0:mi(2),26),springarea(2,*),plkcon(0:2*npmat_,ntmat_,*), local 272 voldl(i2,i1)=vold(i2,konl(i1)) 284 call springstiff_n2f_th(xl,voldl,s,imat,elcon, 292 call springstiff_f2f_th(xl,voldl,s,imat,elcon, 305 call advecstiff(nope,voldl,ithermal,xl,nelemload, 319 call networkstiff(voldl,s,imat,konl,mi,ntmat_,shcon,
|
H A D | mafilldmss.f | 64 & ff(60),xl(3,26),voldl(0:mi(2),26), local 158 voldl(j,k)=vold(j,konl(k)) 177 call springdamp_n2f(xl,elas,voldl,s,imat,elcon, 184 call springdamp_f2f(xl,elas,voldl,s,imat,elcon,
|
H A D | e_c3d_rhs.f | 52 & shp(4,20),voldl(3,20),xl2(3,8),xsj2(3),shp2(7,8), local 150 voldl(i2,i1)=vold(i2,konl(i1)) 361 & (xl(i1,j1)+voldl(i1,j1))
|
H A D | e_c3d_cs_se.f | 76 & anisox(3,3,3,3),voldl(0:mi(2),26),vo(3,3),xloadold(2,*), local 348 voldl(i2,i1)=vold(i2,konl(i1)) 452 voldl(i2,i1)=vold(i2,konl(i1)) 472 call springstiff_n2f(xl,elas,konl,voldl,s,imat,elcon, 481 call springstiff_f2f(xl,elas,voldl,s,imat,elcon,nelcon, 510 q(i1)=xl(i1,1)+voldl(i1,1) 1020 vo(i1,j1)=vo(i1,j1)+shp(j1,k1)*voldl(i1,k1) 1207 & (xl(i1,j1)+voldl(i1,j1))
|
H A D | e_c3d_se.f | 76 & anisox(3,3,3,3),voldl(0:mi(2),26),vo(3,3),xloadold(2,*), local 387 voldl(i2,i1)=vold(i2,konl(i1)) 543 voldl(i2,i1)=vold(i2,konl(i1)) 563 call springstiff_n2f(xl,elas,konl,voldl,s,imat,elcon, 572 call springstiff_f2f(xl,elas,voldl,s,imat,elcon,nelcon, 601 q(i1)=xl(i1,1)+voldl(i1,1) 1111 vo(i1,j1)=vo(i1,j1)+shp(j1,k1)*voldl(i1,k1) 1298 & (xl(i1,j1)+voldl(i1,j1))
|
H A D | e_c3d.f | 76 & anisox(3,3,3,3),voldl(0:mi(2),26),vo(3,3),xloadold(2,*), 373 voldl(i2,i1)=vold(i2,konl(i1)) 409 voldl(i2,i1)=vold(i2,konl(i1)) 430 call springstiff_n2f(xl,elas,konl,voldl,s,imat,elcon, 440 call springstiff_f2f(xl,elas,voldl,s,imat,elcon,nelcon, 482 q(i1)=xl(i1,1)+voldl(i1,1) 1077 vo(i1,j1)=vo(i1,j1)+shp(j1,k1)*voldl(i1,k1) 1289 & (xl(i1,j1)+voldl(i1,j1))
|
H A D | e_c3d_duds.f | 79 & anisox(3,3,3,3),voldl(0:mi(2),26),vo(3,3),xloadold(2,*), 392 voldl(i2,i1)=vold(i2,konl(i1)) 567 voldl(i2,i1)=vold(i2,konl(i1)) 587 call springstiff_n2f(xl,elas,konl,voldl,s,imat,elcon, 596 call springstiff_f2f(xl,elas,voldl,s,imat,elcon,nelcon, 625 q(i1)=xl(i1,1)+voldl(i1,1) 1135 vo(i1,j1)=vo(i1,j1)+shp(j1,k1)*voldl(i1,k1) 1322 & (xl(i1,j1)+voldl(i1,j1))
|
H A D | objective_shapeener_dx.f | 61 & springarea(2,*),elconloc(21),eth(6),xkl(3,3),voldl(0:mi(2),26), 346 voldl(k,j)=vold(k,konl(j)) 668 & shp(m3,m1)*voldl(m2,m1)
|
H A D | resultsmech.f | 65 & voldl(0:mi(2),26),xikl(3,3),ener(mi(1),*), 326 voldl(k,j)=vold(k,konl(j)) 695 & shp(m3,m1)*voldl(m2,m1)
|