resume,lopatka,db

/prep7
cint,dele

Mptemp,,,,,,,,
Mptemp,1,326.85
Mptemp,2,937.85
Mpdata,ex,1,,175000
Mpdata,ex,1,,12200
Mpdata,PRXY,1,,0.3
Mpdata,PRXY,1,,0.3
Mpdata,KXX,1,,15
Mpdata,KXX,1,,25
Mpdata,C,1,,440				
Mpdata,C,1,,640				

Mptemp,,,,,,,,
Mptemp,1,0
Mpdata,dens,1,,8.4e-9			

et,1,solid87

esize,0.2
vmesh,49

esize,0.4
vmesh,17

esize,0.5
vmesh,2
vmesh,6
vmesh,16
vmesh,18

esize,1
vmesh,1
vmesh,3
vmesh,5
vmesh,7
vmesh,8
vmesh,9
vmesh,10
vmesh,11
vmesh,19

esize,1.5		!desky
vmesh,4
vmesh,12,15
vmesh,20,30
vmesh,31,40
vmesh,41,45
vmesh,46,47
vmesh,50

/solu
!! da,plocha,temp,hodnota
DA,48,TEMP,937.85
DA,308,TEMP,937.85
DA,309,TEMP,937.85
DA,44,TEMP,937.85
DA,144,TEMP,937.85
DA,302,TEMP,937.85
DA,301,TEMP,937.85
DA,165,TEMP,937.85
DA,168,TEMP,937.85
DA,311,TEMP,937.85
DA,310,TEMP,937.85

DA,46,TEMP,326.85
DA,200,TEMP,326.85
DA,42,TEMP,326.85
DA,45,TEMP,326.85
DA,142,TEMP,326.85
DA,197,TEMP,326.85
DA,196,TEMP,326.85
DA,193,TEMP,326.85
DA,190,TEMP,326.85
DA,187,TEMP,326.85
DA,184,TEMP,326.85
DA,181,TEMP,326.85
DA,178,TEMP,326.85
DA,175,TEMP,326.85
DA,172,TEMP,326.85
DA,169,TEMP,326.85

solve
finish

/prep7

ldread,temp,1,,,,file,rth

etchg


!sací strana

asel,s,,,48 $ nsla,s,1 $ F,all,Fx,0.033316
asel,s,,,43 $ nsla,s,1 $ F,all,Fx,0.018332
asel,s,,,44 $ nsla,s,1 $ F,all,Fx,0.009935
asel,s,,,144 $ nsla,s,1 $ F,all,Fx,0.009935
asel,s,,,302 $ nsla,s,1 $ F,all,Fx,0.002461


!tlaková strana

asel,s,,,310 $ nsla,s,1 $ F,all,Fx,0.030328
asel,s,,,311 $ nsla,s,1 $ F,all,Fx,0.01059
asel,s,,,168 $ nsla,s,1 $ F,all,Fx,0.004287
asel,s,,,165 $ nsla,s,1 $ F,all,Fx,0.04287
asel,s,,,301 $ nsla,s,1 $ F,all,Fx,0.00234

!okrajové podmínky
!DA,area,label -all - all dof

Da,291,all
Da,306,all

!fatigue

c=2.29e-7
m=2
tb,cgcr,2,,,paris
tbdata,1,c,m

!crack surfaces

csys,23 !lokalní soradnicovy system, na cele trhliny

!front tip
!nsel,s,loc,x,0
!nsel,r,loc,y,0 
lsel,s,,,543
nsll,s,1				
nlist
cm,trhlina,node
allsel,all,all

!surfaces

Asel,s,,,299
nsla,s,1
cm,sur_1,node
allsel,all,all

Asel,s,,,300
nsla,s,1
cm,sur_2,node
allsel,all,all

finish

/solu

antype,static
kbc,1

!sif


cint,new,1
cint,type,sifs
cint,ctnc,trhlina
cint,edir,cs,23,x
cint,norm,23,y					
cint,surf,sur_1,sur_2
cint,ncon,5

!cgrow calculation

cgrow,new,1
cgrow,cid,1
cgrow,method,smart
cgrow,fcg,meth,LC
cgrow,fcg,damx,0.5
cgrow,fcg,srat,0
cgrow,fcoption,mtab,2
nsubst,6
outres,all,all


solve
finish