afgrow.rar_afgrow_pairchv_washaw5_应力强度因子_经验公式

function f_0=funf0(m,n,s) x=[0.2 0.2 0.2 0.2 0.2 0.2 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.4 0.4 0.4 0.4 0.4 0.4 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 1 1 1 1 1 1 ... 1 1 1 1 1 1 ... 1 1 1 1 1 1 ... 1 1 1 1 1 1 ... 1 1 1 1 1 1 ... 1 1 1 1 1 1]'; y=[0 0 0 0 0 0 ... 0.1 0.1 0.1 0.1 0.1 0.1 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.8 0.8 0.8 0.8 0.8 0.8 ... 1 1 1 1 1 1 ... 0 0 0 0 0 0 ... 0.1 0.1 0.1 0.1 0.1 0.1 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.8 0.8 0.8 0.8 0.8 0.8 ... 1 1 1 1 1 1 ... 0 0 0 0 0 0 ... 0.1 0.1 0.1 0.1 0.1 0.1 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.8 0.8 0.8 0.8 0.8 0.8 ... 1 1 1 1 1 1 ... 0 0 0 0 0 0 ... 0.1 0.1 0.1 0.1 0.1 0.1 ... 0.2 0.2 0.2 0.2 0.2 0.2 ... 0.5 0.5 0.5 0.5 0.5 0.5 ... 0.8 0.8 0.8 0.8 0.8 0.8 ... 1 1 1 1 1 1]'; z=[0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1 ... 0 0.1 0.2 0.5 0.8 1]'; v=[1.28 1.285 1.291 1.33 1.39 1.441 ... 1.311 1.318 1.322 1.355 1.408 1.452 ... 1.374 1.38 1.388 1.42 1.47 1.513 ... 1.752 1.787 1.833 1.993 2.243 2.448 ... 2.498 2.663 2.832 3.528 4.603 5.491 ... 3.286 3.585 3.931 5.34 7.514 9.311 ... 1.173 1.177 1.183 1.244 1.314 1.365 ... 1.198 1.201 1.206 1.267 1.343 1.398 ... 1.241 1.246 1.257 1.337 1.417 1.488 ... 1.488 1.521 1.567 1.749 2.072 2.349 ... 1.985 2.069 2.198 2.781 3.816 4.723 ... 2.504 2.638 2.861 3.972 5.875 7.559 ... 1.158 1.16 1.165 1.22 1.302 1.381 ... 1.179 1.18 1.185 1.245 1.339 1.424 ... 1.211 1.217 1.228 1.31 1.417 1.511 ... 1.414 1.446 1.492 1.684 2.045 2.355 ... 1.827 1.899 2.018 2.574 3.609 4.516 ... 2.26 2.368 2.558 3.568 5.38 7.059 ... 1.138 1.141 1.144 1.192 1.343 1.523 ... 1.142 1.144 1.152 1.22 1.399 1.609 ... 1.145 1.154 1.172 1.267 1.486 1.685 ... 1.232 1.251 1.309 1.547 2.056 2.514 ... 1.413 1.47 1.565 2.075 3.171 4.162 ... 1.615 1.686 1.801 2.555 4.19 5.977]'; F=scatteredInterpolant(x,y,z,v,'linear'); f_0=F(m,n,s); end