/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | creeps.f | 309 elcon(3,i,nmat)=elcon(6,1,nmat) 310 elcon(4,i,nmat)=elcon(7,1,nmat) 311 elcon(5,i,nmat)=elcon(8,1,nmat) 313 elcon(3,i,nmat)=elcon(6,id,nmat) 314 elcon(4,i,nmat)=elcon(7,id,nmat) 315 elcon(5,i,nmat)=elcon(8,id,nmat) 397 elcon(13,i,nmat)=elcon(16,1,nmat) 398 elcon(14,i,nmat)=elcon(17,1,nmat) 399 elcon(15,i,nmat)=elcon(18,1,nmat) 401 elcon(13,i,nmat)=elcon(16,id,nmat) [all …]
|
H A D | surfacebehaviors.f | 117 elcon(3,1,imat)=1.5d0 153 elcon(1,1,imat)=dlog(100.d0)/elcon(1,1,imat) 179 elcon(2,1,imat)=50.d0*elcon(1,1,1) 193 elcon(1,1,imat)=elcon(1,1,1)/70000.d0 211 & elcon(2,1,imat) 232 elcon(2,1,imat)=50.d0*elcon(1,1,1) 238 & elcon(1,1,imat) 257 elcon(1,1,imat)=elcon(1,1,1)/70000.d0 265 & elcon(4,1,imat) 282 elcon(3,1,imat)=3.5d0 [all …]
|
H A D | getcontactparams.f | 42 & elcon(0:ncmat_,ntmat_,*) 53 mu=elcon(6,1,imat) 54 if(elcon(7,1,imat).le.0.0)then 57 fktauinv=1.0/elcon(7,1,imat) 60 if(elcon(8,1,imat).le.0.99)then 65 iwan=int(elcon(8,1,imat)) 77 if(elcon(3,1,imat).gt.1.4 .and. 78 & elcon(3,1,imat).lt.1.6 )then 80 p0=elcon(2,1,imat) 81 beta=1.0/(elcon(1,1,imat)) [all …]
|
H A D | hyperelastics.f | 34 real*8 elcon(0:ncmat_,ntmat_,*),um local 243 um=elcon(1,j,nmat) 245 um=2.d0*(elcon(1,j,nmat)+elcon(2,j,nmat)) 247 um=2.d0*elcon(1,j,nmat) 249 um=elcon(1,j,nmat) 251 um=elcon(1,j,nmat)+elcon(3,j,nmat) 253 um=elcon(1,j,nmat)+elcon(3,j,nmat)+elcon(5,j,nmat) 255 um=2.d0*(elcon(1,j,nmat)+elcon(2,j,nmat)) 257 um=2.d0*elcon(1,j,nmat) 259 um=elcon(1,j,nmat) [all …]
|
H A D | elastics.f | 105 & elcon(i,ntmat,nmat) 140 elcon(0,ntmat,nmat)=0.d0 200 elcon(0,ntmat,nmat)=0.d0 203 e1=elcon(1,ntmat,nmat) 204 e2=elcon(2,ntmat,nmat) 205 e3=elcon(3,ntmat,nmat) 206 un12=elcon(4,ntmat,nmat) 207 un13=elcon(5,ntmat,nmat) 208 un23=elcon(6,ntmat,nmat) 270 & elcon(8+i,ntmat,nmat) [all …]
|
H A D | plastics.f | 271 t1l=elcon(0,i,nmat) 307 elcon(11,i,nmat)=0.d0 308 elcon(12,i,nmat)=0.d0 310 elcon(14,i,nmat)=1.d0 315 elcon(12,i,nmat)=0.d0 317 elcon(14,i,nmat)=1.d0 337 t1l=elcon(0,i,nmat) 373 elcon(11,i,nmat)=0.d0 374 elcon(12,i,nmat)=0.d0 376 elcon(14,i,nmat)=1.d0 [all …]
|
H A D | materialdata_em.f | 18 subroutine materialdata_em(elcon,nelcon,alcon,nalcon, argument 29 real*8 elcon(0:ncmat_,ntmat_,*),alcon(0:6,ntmat_,*),alpha(6),t1l, local 67 call ident2(elcon(0,1,imat),t1l,nelcon(2,imat),ncmat_+1,id) 72 elconloc(k)=elcon(k,1,imat) 76 elconloc(k)=elcon(k,1,imat) 80 elconloc(k)=elcon(k,id,imat) 84 elconloc(k)=elcon(k,id,imat)+ 85 & (elcon(k,id+1,imat)-elcon(k,id,imat))* 86 & (t1l-elcon(0,id,imat))/ 87 & (elcon(0,id+1,imat)-elcon(0,id,imat))
|
H A D | materialdata_sp.f | 18 subroutine materialdata_sp(elcon,nelcon,imat,ntmat_,i,t1l, argument 28 real*8 elcon(0:ncmat_,ntmat_,*),t1l,elconloc(*), local 90 call ident2(elcon(0,1,imat),t1l,nelcon(2,imat),ncmat_+1,id) 95 elconloc(k)=elcon(k,1,imat) 99 elconloc(k)=elcon(k,1,imat) 103 elconloc(k)=elcon(k,id,imat) 107 elconloc(k)=elcon(k,id,imat)+ 108 & (elcon(k,id+1,imat)-elcon(k,id,imat))* 109 & (t1l-elcon(0,id,imat))/ 110 & (elcon(0,id+1,imat)-elcon(0,id,imat))
|
H A D | gapheatgenerations.f | 19 subroutine gapheatgenerations(inpc,textpart,elcon,nelcon, argument 36 real*8 elcon(0:ncmat_,ntmat_,*) local 84 & elcon(8+i,1,imat) 91 if(elcon(9,1,imat).le.0.d0) then 100 if((elcon(10,1,imat).lt.0.d0).or. 101 & (elcon(10,1,imat).gt.1.d0)) then 111 elcon(0,1,imat)=0.d0 126 elcon(11,1,imat)=-0.01d0 127 elseif(elcon(11,1,imat).lt.0.d0) then 128 elcon(11,1,imat)=-1.d0
|
H A D | frictions.f | 19 subroutine frictions(inpc,textpart,elcon,nelcon, argument 34 real*8 elcon(0:ncmat_,ntmat_,*) local 71 & elcon(5+i,1,imat) 78 if(elcon(6,1,imat).le.0.d0) then 85 if(elcon(7,1,imat).le.0.d0) then 90 & elcon(1,1,1)/2.d0 99 elcon(7,1,imat)=elcon(1,1,1)/2.d0 101 elcon(0,1,imat)=0.d0
|
H A D | springforc_n2f.f | 133 idof=nint(elcon(3,1,imat)) 171 idof1=nint(elcon(3,1,imat)) 172 idof2=nint(elcon(4,1,imat)) 273 if(int(elcon(3,1,imat)).eq.1) then 283 beta=elcon(1,1,imat) 292 elseif(int(elcon(3,1,imat)).eq.2) then 304 xk=elcon(2,1,imat) 312 eps=elcon(1,1,imat)*pi/xk 326 eps=elcon(1,1,imat)*pi/elcon(2,1,imat) 334 elseif(int(elcon(3,1,imat)).eq.3) then [all …]
|
H A D | contactdampings.f | 19 subroutine contactdampings(inpc,textpart,elcon,nelcon, argument 34 real*8 elcon(0:ncmat_,ntmat_,*) local 54 elcon(8,1,imat)=0.d0 59 & elcon(8,1,imat) 85 & elcon(5,1,imat) 91 elcon(0,1,imat)=0.d0
|
H A D | springstiff_n2f_th.f | 19 subroutine springstiff_n2f_th(xl,voldl,s,imat,elcon,nelcon, argument 36 & alpha,beta,elcon(0:ncmat_,ntmat_,*),xm(3),pressure, 121 if(elcon(1,1,imat).gt.0.d0) then 125 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 128 alpha=elcon(2,1,imat) 129 beta=elcon(1,1,imat) 140 eps=elcon(1,1,imat)*pi/elcon(2,1,imat) 141 pressure=-elcon(2,1,imat)*val* 158 call materialdata_sp(elcon,nelcon,imat,ntmat_,i,tmean,
|
H A D | springforc_n2f_th.f | 19 subroutine springforc_n2f_th(xl,vl,imat,elcon,nelcon, argument 37 & elcon(0:ncmat_,ntmat_,*),pproj(3),xsj2(3),xs2(3,7),val, local 115 if(elcon(1,1,imat).gt.0.d0) then 119 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 124 alpha=elcon(2,1,imat) 125 beta=elcon(1,1,imat) 136 eps=elcon(1,1,imat)*pi/elcon(2,1,imat) 137 pressure=-elcon(2,1,imat)*val* 154 call materialdata_sp(elcon,nelcon,imat,ntmat_,i,tmean,
|
H A D | magneticpermeabilitys.f | 19 subroutine magneticpermeabilitys(inpc,textpart,elcon,nelcon, argument 34 real*8 elcon(0:ncmat_,ntmat_,*) local 89 & elcon(1,ntmat,nmat) 102 elcon(2,ntmat,nmat)=idomain+0.5d0 106 & elcon(0,ntmat,nmat) 113 elcon(0,ntmat,nmat)=0.d0
|
H A D | assigndomtonodes.f | 20 & elcon,ncmat_,ntmat_,mi,ne2) argument 33 real*8 elcon(0:ncmat_,ntmat_,*) local 60 if(int(elcon(2,1,imat)).eq.2) ne2=ne2+1 65 if(inomat(node).ne.int(elcon(2,1,imat))) then 74 inomat(node)=int(elcon(2,1,imat))
|
H A D | springstiff_f2f.f | 157 if(int(elcon(3,1,imat)).eq.1) then 161 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 165 alpha=elcon(2,1,imat)*springarea(1) 166 beta=elcon(1,1,imat) 173 elseif((int(elcon(3,1,imat)).eq.2).or. 174 & (int(elcon(3,1,imat)).eq.4)) then 182 elseif(int(elcon(3,1,imat)).eq.3) then 233 if(int(elcon(3,1,imat)).eq.4) then 236 um=elcon(6,1,imat) 238 stickslope=elcon(7,1,imat)/kscale [all …]
|
H A D | springstiff_f2f_th.f | 19 subroutine springstiff_f2f_th(xl,voldl,s,imat,elcon,nelcon, argument 40 & alpha,beta,elcon(0:ncmat_,ntmat_,*),xm(3), local 162 if(int(elcon(3,1,imat)).eq.1) then 166 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 169 alpha=elcon(2,1,imat) 170 beta=elcon(1,1,imat) 176 elseif((int(elcon(3,1,imat)).eq.2).or. 177 & (int(elcon(3,1,imat)).eq.4)) then 181 pressure=-elcon(2,1,imat)*clear 182 elseif(int(elcon(3,1,imat)).eq.3) then [all …]
|
H A D | springforc_f2f_th.f | 19 subroutine springforc_f2f_th(xl,vl,imat,elcon,nelcon, argument 41 & elcon(0:ncmat_,ntmat_,*),pproj(3),clear, local 170 if(int(elcon(3,1,imat)).eq.1) then 174 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 179 alpha=elcon(2,1,imat) 180 beta=elcon(1,1,imat) 186 elseif((int(elcon(3,1,imat)).eq.2).or. 187 & (int(elcon(3,1,imat)).eq.4)) then 191 pressure=-elcon(2,1,imat)*clear 192 elseif(int(elcon(3,1,imat)).eq.3) then [all …]
|
H A D | springs.f | 21 & ncmat_,elcon,matname,irstrt,istep,istat,n,iline,ipol, argument 43 & elcon(0:ncmat_,ntmat_,*) local 152 & elcon(3,1,nmat) 166 & elcon(3,1,nmat) 168 & elcon(4,1,nmat) 183 & elcon(i,ntmat,nmat) 192 & elcon(0,ntmat,nmat) 199 elcon(0,ntmat,nmat)=0.d0 239 & elcon(3,1,nmat) 241 & elcon(4,1,nmat)
|
H A D | hyperfoams.f | 19 subroutine hyperfoams(inpc,textpart,elcon,nelcon, argument 34 real*8 elcon(0:ncmat_,ntmat_,*) local 104 & elcon(i,ntmat,nmat) 112 & elcon(0,ntmat,nmat) 133 & elcon(i,ntmat,nmat) 154 & elcon(8+i,ntmat,nmat) 162 & elcon(0,ntmat,nmat)
|
H A D | gaps.f | 21 & ncmat_,elcon,matname,irstrt,istep,istat,n,iline,ipol, argument 41 & elcon(0:ncmat_,ntmat_,*) local 98 elcon(5,ntmat,nmat)=1.d12 99 elcon(6,ntmat,nmat)=1.d-3 105 & elcon(i,ntmat,nmat) 118 & elcon(i-1,ntmat,nmat) 126 elcon(0,ntmat,nmat)=0.d0
|
H A D | springforc_f2f.f | 169 if(int(elcon(3,1,imat)).eq.1) then 173 if(dabs(elcon(2,1,imat)).lt.1.d-30) then 178 alpha=elcon(2,1,imat)*springarea(1) 179 beta=elcon(1,1,imat) 188 elseif((int(elcon(3,1,imat)).eq.2).or. 189 & (int(elcon(3,1,imat)).eq.4)) then 204 elseif(int(elcon(3,1,imat)).eq.3) then 271 if(int(elcon(3,1,imat)).eq.4) then 274 um=elcon(6,1,imat) 276 stickslope=elcon(7,1,imat)/kscale [all …]
|
H A D | materialdata_me.f | 18 subroutine materialdata_me(elcon,nelcon,rhcon,nrhcon,alcon,nalcon, argument 40 real*8 elcon(0:ncmat_,ntmat_,*),rhcon(0:1,ntmat_,*), local 181 elconloc(k)=elcon(k,1,imat) 185 elconloc(k)=elcon(k,1,imat) 392 elconloc(k)=elcon(k,1,imat) 396 elconloc(k)=elcon(k,1,imat) 400 elconloc(k)=elcon(k,id,imat) 404 elconloc(k)=elcon(k,id,imat)+ 405 & (elcon(k,id+1,imat)-elcon(k,id,imat))* 406 & (t1l-elcon(0,id,imat))/ [all …]
|
H A D | cyclichardenings.f | 20 & npmat_,plicon,nplicon,ncmat_,elcon,matname, argument 37 & elcon(0:ncmat_,ntmat_,*),plconloc(802),t1l local 159 t1l=elcon(0,i,nmat) 190 elcon(10,i,nmat)=plconloc(2) 191 elcon(11,i,nmat)=0.d0 193 elcon(10,i,nmat)=plconloc(2) 194 elcon(11,i,nmat)=(plconloc(4)-plconloc(2))/
|