PlateR6: R6 analysis for Surface Cracks in Plates

In ShinsukeSakai0321/FFS:

Description

• K-solution: Raju_Newman method

• L-solution: Dillstrom, P.and Sattari-Far

Usage

PlateR6(E = 1.92e+11, nu = 0.3, Su = 4.9e+08, Sy = 2.7e+08, J1c = 1e+05, a = 0.02, cc = 0.015, b = 0.1, t = 0.04, P = 2e+06, M = 0)

Arguments

E
nu
Su
Sy
J1c
a
cc
b
t
P
M

Author(s)

Shinsuke Sakai

Examples

#Kouzou Kenzensei Hyouka Handbook
#subject 2.5
a<-0.02; cc<-0.015;
res <- PlateR6(Su = 4.9e+08, Sy = 2.7e+08,a = a, cc = cc, t = 0.04, b = 0.1, P=2.7e6,M=0)
cat("Kr=",res[1],",Lr=",res[2])
cat("Safety Margin=",SafetyMargin(Kr=res[1],Lr=res[2]))
DrawR6(res[1],res[2],Su = 4.9e+08, Sy = 2.7e+08)
# solution
#Kr= 0.3266297 ,Lr= 0.9459494

#subject 2.6
a<-0.02; cc<-0.015;
for(i in 1:300000){
  Kmax<-Raju_Newman(a = a, cc = cc, t = 0.04, b = 0.1, P=8e5,M=0)
  da<-dadN1(R=0.0,dK=Kmax$KA/1e6)/1e3
  dc<-dadN1(R=0.0,dK=Kmax$KB/1e6)/1e3
  a<-a+da
  cc<-cc+dc
  #cat(i,"a=",a,"2c=",cc*2,"\n")
}
cat("a=",a,"2c=",cc*2)
#expected solution
#a= 0.03044039 2c= 0.07068406
res<-PlateR6(Su = 4.9e+08, Sy = 2.7e+08, a = a, cc = cc,  t = 0.04, b = 0.1, P=2.7e6, M = 0)
DrawR6(Kr=res[1],Lr=res[2],Su = 4.9e+08, Sy = 2.7e+08)
cat("Safety Margin=",SafetyMargin(Kr=res[1],Lr=res[2]))

# Example of probabilistic fracture mechanics
n<-500
a<-rnorm(n,mean=0.02,sd=0.02*0.2)
cc<-rnorm(n,mean=0.015,sd=0.015*0.2)
P<-rnorm(n,mean=2.7e6,sd=2.7e6*0.1)
Sy<-rnorm(n,mean=2.7e8,sd=2.7e8*0.05)
for(i in 1:n){
  res <- PlateR6(Su = 4.9e+08, Sy = Sy[i],a = a[i], cc = cc[i], t = 0.04, b = 0.1, P=2.7e6,M=0)
  DrawR6(res[1],res[2],Su = 4.9e+08, Sy = 2.7e+08) ; par(new=T)
}

ShinsukeSakai0321/FFS documentation built on May 20, 2019, 5:09 p.m.