鋼板剪力牆折減係數0﹒743,其舊結構可靠度C語言原始碼如下所示。fi是Try&Error的數值。函數不一定是正態分布,應根據實際狀況使用函數。以結構可靠度攻讀博士學位,其不成文規定為懂折減係數,就可以逕行發出博士學位。Rosenblatt求逆,我在清華學習手算方法,並不是數學軟體計算方法,我據其過程提供數學軟體方法。
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_rng.h>
int main(void) {
const gsl_rng_type *T;
gsl_rng *r;
int i, x;
double Fy, B, t, fi, Vn, Z, times;
gsl_rng_env_setup();
T = gsl_rng_default;
r = gsl_rng_alloc(T);
x = 0;
fi = 0.743;
times = pow(10, 8);
for (i = 0; i < times; i++) {
Fy = 193.5 + gsl_rng_uniform(r) * 43;
B = 2490 + gsl_rng_uniform(r) * 20;
t = 49 + gsl_rng_uniform(r) * 2;
Vn = 4137676.8 + gsl_rng_uniform(r) * 2068838.4;
Z = sqrt(3) / 9 * Fy * B * t - fi * Vn;
if (Z < 0)
x++;
}
printf("%10.8f\n", x / times);
gsl_rng_free(r);
return 0;
}
C:\>main
0.00007983
Mathematica
#include <stdio.h>
#include <math.h>
#include <gsl/gsl_rng.h>
int main(void) {
const gsl_rng_type *T;
gsl_rng *r;
int i, x;
double Fy, B, t, fi, Vn, Z, times;
gsl_rng_env_setup();
T = gsl_rng_default;
r = gsl_rng_alloc(T);
x = 0;
fi = 0.743;
times = pow(10, 8);
for (i = 0; i < times; i++) {
Fy = 193.5 + gsl_rng_uniform(r) * 43;
B = 2490 + gsl_rng_uniform(r) * 20;
t = 49 + gsl_rng_uniform(r) * 2;
Vn = 4137676.8 + gsl_rng_uniform(r) * 2068838.4;
Z = sqrt(3) / 9 * Fy * B * t - fi * Vn;
if (Z < 0)
x++;
}
printf("%10.8f\n", x / times);
gsl_rng_free(r);
return 0;
}
C:\>main
0.00007983
Mathematica
結構可靠度結構功能函數 |