Comparison of reliability evaluation methods between China and Canada standards for wood structures featuring duration of load

Duration of load (DOL) is a key factor in design of wood structures, which makes the reliability analysis of wood structures more complicated. The importance of DOL is widely recognized, yet the methods and models through which it is incorporated into design codes vary substantially by country/region. Few investigations of the effect of different model assumptions of DOL and other random variables on the results of reliability analysis of wood structures can be found. In this paper, comparisons are made on the reliability analysis methods that underlie the China and the Canada standards for design of wood structures. Main characteristics of these two methods, especially the way how DOL is treated are investigated. Reliability analysis was carried out with the two methods employing the same set of material properties and load parameters. The resulted relationships between reliability index β and resistance partial factor γR* (the β–γR* curves) for four load combinations are compared to study the safety level indicated by the two methods. The comparison shows that the damage accumulation model (Foschi–Yao model) in the Canada analysis method is highly dependent on the type and duration of load, resulting in more conservative design than the China analysis method in loading cases dominated by dead load, but less conservative design in cases of high level of live loads. The characteristics of the load effect term of the performance function are also found to make considerable difference in reliability levels between the two methods. This study aims to provide references for researchers and standard developers in the field of wood structures.


Introduction
Following the general principles set out in ISO 2394 [1], most countries have adopted the reliability-based design method in their structural design codes. This method involves a concept of "limit state", which is a critical condition for a structure/member beyond which the relevant design criterion is no longer fulfilled. Two limit states, viz. the ultimate limit state and the serviceability limit state, are generally adopted, respectively, for structural safety and serviceability performance. Structural reliability represents the probability of a structure fulfilling its intended function without failure over the design reference period.
One feature of wood that distinguishes it from other structural materials (e.g., steel or concrete) is duration of load (DOL) [2], i.e., as the load on a wood member persists, the strength of the wood decreases with time. Only some of the many studies about DOL are mentioned below for a brief appreciation of the research history and the state-of-the-art progress. Ever since its discovery in 1740 [3], DOL has been widely studied by means of long-term experiments of strength of wood by researchers around the world [4][5][6][7]. Wood [8] proposed the first strength reduction curve, known as the Madison curve, by fitting a series of experimental data. This curve was purely empirical and was not reflecting the variability in DOL, but it effectively facilitated design practice and is still in use today in many design codes. Damage accumulation models [9][10][11] were later proposed which take account of the variability in DOL and allow it to be estimated given any arbitrary stress history in a member. Recently a model was also proposed [12] to better reflect the stochastic damage accumulation process in DOL. Existing models that treat DOL either as a random variable or as a stochastic process accommodate the need of covering it in reliability analysis of wood structures. However, current wood design standards in different countries differ considerably regarding whether a reliability basis is employed and the way how DOL is covered in reliability evaluation. For example, Canada and China wood standards both determine design strengths of wood on the basis of reliability analysis [13,14], in which DOL is included as a random variable. The American standard has two design methods used in parallel, one being a reliability-based load and resistance factor design (LRFD) method of a similar reliability calibration process with Canada standard [15], and the other being an empirical allowable stress design (ASD) method (including the equivalent LRFD format equations converted from it). The European design code for wood structures (EC 5) also has a reliability basis [16].
Booming of wood structures around the world [17,18] is stimulating international trade of wood products and international collaborations in timber engineering. For example, China has been importing huge volumes of wood products such as glulam, dimension lumber, woodbased structural panels and structural composite lumber from forest-rich countries like North America, Russia, EU and New Zealand, to meet the increasing demand for wood structures during the past two decades. The trend of international collaborations in wood structure engineering poses a potential necessity of comparison between design standards of different countries. Questions might arise such as which standard yields more conservative design, and how could their reliability level be compared? Since the assumptions and treatments of parameters involved in the reliability analysis of wood structures in different countries are quite different, it is of interest to investigate the influence of such assumptions and treatments as a fundamental step. Since Canadian lumber is most widely used in China and the lumber is incorporated into the national wood standard GB 50005 [19] via a newly conducted reliability analysis [14], this paper focuses on the reliability analysis methods, respectively, adopted by China and Canada standards for wood structures. Comparisons are made between the assumptions and treatments of related variables in the two reliability analysis methods, and the effect on safety level in the respective wood standard from the analysis methods investigated.

Methods
In this section, a review of the China and Canada reliability analysis methods is presented, with special attention paid to the detailed differences in the assumptions and treatments of all random variables (including the DOL factors). Parameters for reliability analysis and comparison are also discussed.

The China reliability analysis method
The design value of wood strength f d in China standard GB 50005 [19] is expressed as f d = f k K DOL /γ R . Accordingly, the performance function for ultimate limit state reliability analysis was established by Zhu et al. [14] as where f is a random variable denoting the short-term strength of wood; K A is a random variable accounting for uncertainty of the cross-sectional dimensions of wood member, defined as the ratio of actual cross-section property to the nominal value in resistance calculation; K P is a random variable accounting for uncertainty associated with resistance prediction model, defined as the ratio of true resistance to that obtained from resistance prediction model; K Q3 is a random variable accounting for DOL, defined as the ratio of long-term strength to shortterm strength of wood; f k is the characteristic strength of wood, taken as the 5% fractile of the short-term strength f; K DOL is the DOL factor specified in the design standard, taken as the mean value of random variable K Q3 ; γ R , γ G and γ Q are, respectively, the resistance partial factor, dead load partial factor and live load partial factor; ψ c is the combination factor for live load; g is a random variable, defined as the ratio of actual dead load G (random) to the characteristic value G K ; q is a random variable, defined as the ratio of actual live load Q (random) to the characteristic value Q K ; ρ is the ratio of characteristic live load to characteristic dead load, ρ = Q K /G K ; K B is a random variable accounting for uncertainty associated with load effect calculation model, defined as the ratio of true load effect to that obtained from calculation model.
The two terms on the right-hand side of Eq. (1) represent resistance and load effect of a wood member, respectively. Equation (1) is applicable only to wood members loaded parallel to grain. The short-term strength f of wood is assumed to follow a lognormal probability distribution. Its mean value f m , coefficient of variation (COV) V f and the characteristic strength f k comply with a relationship as (2) f m /f k = exp 1.645 ln 1 The mean value and COV of random variables K A , K P , K Q3 and K B are listed in Table 1, for wood members under bending, compression and tension, respectively.
In the China standard, the following load conditions are distinguished for DOL: (1) load of normal duration such as office or residential floor live load, snow load and wind load, E KQ3 = K DOL = 0.72; (2) temporary load such as construction load, E KQ3 = K DOL = 0.72 × 1.2 = 0.864; (3) load of long duration such as dead load-dominated cases when dead load contributes over 80% of the total load effect, E KQ3 = K DOL = 0.72 × 0.8 = 0.576.
With the performance function and parameters detailed above, reliability analysis for the China standard was conducted by employing first-order reliability method (FORM) for different load combinations [14].

The Canada reliability analysis method
For Canada wood standard CSA O86 [20], reliability evaluation is based on the research by Foschi et al. [13]. The reliability analysis was divided into two steps, viz. short-term analysis (DOL not included) and long-term analysis (DOL included). For short-term analysis, the performance function was established as where σ s is a random variable denoting short-term strength of wood; σ 0.05 is the 5% fractile of σ s ; g is the ratio of actual dead load G (random) to nominal dead load G n ; q is the ratio of actual live load Q (random) to the nominal live load Q n ; r is the ratio of nominal dead load to nominal live load (reciprocal of ρ in Eq. (1)), r = G n /Q n ; φ, γ G and γ Q are, respectively, the resistance factor, dead load factor and live load factor.
For long-term reliability analysis, the performance function was derived as Eq. (4), where α is the damage index representing the strength reduction due to DOL. Index α accumulates from 0 up to 1, with α = 0 denoting no damage and α = 1 denoting failure of wood. Accumulation of α with time is at a rate expressed with Eq. (5), which is known as Foschi-Yao damage accumulation model [11]: where a, b, c, n and η are model parameters; σ s is the short-term strength of timber; σ(t) is the applied stress as function of time t; η is a threshold stress ratio below which damage is assumed not to be incurred by the applied stress; α(t) is the damage index at time t. Model parameters a, b, c, n and η are assumed constant for a given structural member, but vary randomly across members. Parameter a is not independent but a function of b, c, n, η and σ s [21]. Stress history σ(t) in a member can be expressed with Eq. (6). The load variables in Eq. (6) are different from those in Eq. (3). For short-term analysis the maximum load effect is concerned, therefore random variables g and q in Eq. (3) represent the maximum load over the design reference period. For long-term analysis, the load effect at any time t should be used for damage accumulation calculation, therefore, random variables g(t) and q(t) in Eq. (6) represent the load history over the design reference period. The statistical load parameters for reliability evaluation of Canadian lumber can be found in the work by Foschi et al. [13]: Instead of a close-formed integration, numerical integration is used for calculation of damage index α. The design reference period t is divided into several segments of duration Δt, and the damage accumulation is calculated segment by segment from Eq. (7). In the equation, α i is the damage index value by the end of the ith time segment; a, b, c, n and η are Foschi-Yao damage model parameters; σ i is the applied stress in the ith time segment calculated from Eq. (6), which is assumed constant over each time segment:

Table 1 Statistical parameters of random variables in China analysis method
Mean value and COV of a variable K is denoted as E K and V K , respectively. E KQ3 is multiplied with a modification factor of 1.2 and 0.8 for construction load and load of long duration, respectively Long-term reliability analysis for the Canada standard was performed via Monte Carlo simulations as follows: (1) relevant design conditions and parameters are prescribed. A resistance factor φ, load type and ratio r are selected. Statistical characteristics (the mean, COV and probability distribution) are assumed for the random variables, including load ratios g and q, short-term wood strength σ s and Foschi-Yao model parameters b, c, n and η. (2) N R sets of the parameters are generated by random sampling. Each set of the parameters include wood strength parameters (σ s , b, c, n and η) for a member and the load ratios (g i , q i ) for each time segment over the design reference period. (3) For each set of the N R samples, the performance function Z is evaluated with Eqs. (4)(5)(6)(7). The total number of the failed samples (Z ≤ 0) is counted to be N F . (4) The probability of failure corresponding to conditions prescribed in Step (1) is calculated as P f = N F /N R , and the corresponding reliability index is calculated as The relationships between reliability index β and resistance factor φ were obtained from short-term and longterm reliability analysis, respectively, as shown in Fig. 1. Corresponding to a target reliability index β 0 , the longterm resistance factor φ II is always smaller than the short-term value φ I . DOL is then defined as the ratio of resistance factor φ II to φ I , i.e., K DOL = φ II /φ I . Thus a value of K DOL = 0.8 for load of standard-term duration in Canada standard was derived from the reliability analysis (on a 30-year basis) [13]. For loads of short-term and longterm duration, the difference in value of K DOL so obtained is accounted for by introducing an adjustment factor K D , as listed in Table 2. It follows that for short-term duration .92, and for long-term duration K DOL = 0.8 × 0.65 = 0.52.

Main differences in the performance functions
The China analysis method takes account of more random variables (K A , K P and K B ) that may affect the resistance and load effect in the performance function. These random variables, however, have limited effect on the reliability analysis results as all their mean is 1 and COV is no higher than 0.05 (see Table 1). A more significant difference lies in the way of dealing DOL. In the performance function of China analysis method, DOL of wood is reflected by K Q3 and K DOL , where K DOL is a constant discounting the design strength of wood, and K Q3 is a variable reflecting the random value of DOL. In Eq. (1), K Q3 in the resistance term is randomly valued, with a mean being equal to the design parameter K DOL in the load effect term. When evaluating the reliability, mean value of K Q3 and value of K DOL in Eq. (1) actually offset each other, what left and really affects the reliability results is the variation of DOL. Therefore, influence of DOL on reliability analysis is actually influence of COV of K Q3 (V KQ3 = 0.12 regardless of the load level and its duration over the design reference period, see Table 1). In contrast, in the Canada analysis method, reliability is directly affected by the DOL that is reflected by damage index α. Variation of DOL is introduced by the model parameters (see Table 3). The COV of these model parameters is considerably high (0.075-0.574) compared to that of K Q3 in the China analysis method.

Material and load parameters in reliability analysis
For assessing the relative safety level achieved from the China and Canada reliability analysis methods, the input parameters for wood and load profiles are carefully chosen such that they are as closely comparable as possible.

Material parameters
Western Hemlock dimension lumber is used to conduct the reliability analyses. To explore the influence of different variations of strength, members under bending and members under compression are included in the study. Specifically, the bending strength of lumber (2″ × 6″, No. 2 and better) [13] and the compression strength of lumber (2″ × 4″, No. 2) [22] are used, as listed in Table 3. The COV of strength is 0.410 for bending and 0.246 for compression, respectively. The lumber strength is assumed to follow a lognormal probability distribution, and the relationship in Eq. (2) applies. In the analysis the characteristic strength f k in the China method is equivalent to the nominal strength σ 0.05 in the Canada method.
The estimated Foschi-Yao model parameters [13] are adopted in the current analysis, e.g., mean (b) = 37.161, COV (b) = 0.281 (see Table 3). It is noted that uncertainties of parameter estimates for damage accumulation models do exist, viz. different estimates of the model parameters might be obtained if model fitting were carried out for different materials, different samples of the same material, or even the same material sample employing different fitting methods. Influence of such uncertainties was initially examined by Foschi et al. [13] and more recently Chun-Hao Yang et al. [23] proposed a Bayesian framework for model fitting with better likelihood. The uncertainties of model parameter estimates are a rather complicated problem beyond the scope of this study.

Load parameters
Statistical data from the China load code GB 50009 [24] are adopted to derive the load parameters for reliability analysis, as tabulated in Tables 4, 5, 6. The design reference period is 50 years. For the China analysis method, the statistical characteristics of load are listed in Table 4.
For the Canada analysis method, the damage accumulation model requires division of the design reference period into proper number of segments as well as determination of load ratio g and q over each segment. The load statistical parameters are listed in Table 5. Loads at different points of time are randomly sampled according to the statistical characteristics. Specifically, for a given member, the value of dead load is sampled only once, and is assumed to last all through the whole design reference period. Floor live load is superposition of a sustained part and an extraordinary part, both alter once every 10 years, and therefore, one random sampling is made for every 10-year period. The sustained part of floor live load is assumed to last over the whole period, while the extraordinary floor live load is assumed to be intermittent, lasting 4 h per week [25].
Reference snow load specified in load code GB 50009 [24] is derived on the basis of statistical data of annual maximum snow pressure S 0 , which is a random variable following an Extreme Type I probability distribution. The cumulative probability distribution of S 0 is expressed in Eq. (8), and the maximum snow pressure S 0,R of return period R is expressed in Eq. (9): where u and α are probability distribution parameters; R is the return period, in years. S 0,10 , S 0,50 and S 0,100 for some cities in China are tabulated in load code GB 50009 [24], from which the parameters u and α can be backcalculated via Eq. (9). For example, the snow statistical parameters for Beijing are listed in Table 6. S 0,50 is taken as the reference snow pressure for a 50-year design reference period, and the nominal snow pressure on roof is calculated as μ r S 0,50 , where μ r is a snow distribution factor dependent on roof shape, wind direction and velocity.
For the China reliability analysis method, the maximum snow load (over a 50-year design reference period) is used, the statistical parameters of which are included in Table 4.
For the Canada reliability analysis method, the snow pressure history (instead of the maximum pressure) within the design reference period is needed for calculation of damage accumulation. The snow load history is derived from the annual maximum snow pressure S 0 referring to research by Foschi et al. [13]. Assume that winter is 5-month-long each year and it is divided into N segments, each of duration Δt. Snow pressure in each segment Δt is assumed to remain constant. The cumulative probability distribution of snow pressure in each segment is assumed the same but mutually independent, which is denoted F S (x). The cumulative probability distribution of annual maximum snow pressure is denoted F(x) (see Eq. (8)). It's evident that the probability distribution of the maximum value among the N segments should equal the probability distribution of the annual maximum snow pressure. Assume that p s0 and p 0 are, respectively, the probability of not snowing in a segment Δt and in a year (N·Δt), the expressions in Eqs. (10,11) stand. The cumulative distribution function F S (x) can be derived as Eq. (12): For Monte Carlo simulations in the Canada analysis method, the snow load S Δt = x in a time segment can be sampled with Eq. (12) by assigning a random number (uniformly distributed between 0 and 1) to F s (x). Load ratio q(t) for use in Eq. (6) is determined as

Load combinations
Four load combinations are considered, these include dead load alone and three other combinations of dead load and a single live load, i.e., office floor live, residential The load partial factors specified in load code GB 50009 [24] are used. For the case of dead load alone, the load partial factor γ G = 1.35; for the combinations of dead load and a single live load, the dead load partial factor γ G = 1.2, the live load partial factor γ Q = 1.4, and the combination factor ψ c = 1.0.

Unified format of result presentation
The China analysis method introduces DOL explicitly with a factor K DOL (see "The China reliability analysis method" section), while the Canada analysis method has the DOL factor implicitly included in the resistance factor φ (see "The Canada reliability analysis method" section, φ II = φ I K DOL ). Term γ R /K DOL in the China analysis method is thus equivalent to term 1/φ in the Canada analysis method. To enable comparison between the reliability analysis results, a unified resistance partial factor γ R * is defined and used in place of the two equivalent terms, i.e., γ R * = 1/φ for the Canada analysis method, and γ R * = γ R /K DOL for the China analysis method.

Results and discussion
Reliability analysis of hemlock dimension lumber under bending and under compression was conducted using the China analysis method and the Canada analysis method, respectively. The results are discussed as follows. Note that only long-term analysis was carried out for the Canada analysis method, the results from which are comparable to those from the China analysis method.

Dead load alone
The β-γ R * curves for the case of dead load alone are plotted in Fig. 2. It is observed that the China curves lie above the corresponding Canada curves. Given the same reliability index β, the resistance partial factor γ R * resulted  from the China reliability analysis method is smaller than that from the Canada analysis method. This means that for dead load alone, the China analysis method results in less conservative design than the Canada analysis method. Specifically, at a reliability index of 3.2, the target reliability index by the China code requirement [26], the γ R * value for bending by the China analysis method and the Canada analysis method is 2.63 and 3.19, respectively, resulting in a ratio of about 0.82. This higher conservatism of the Canada analysis method is due to the damage accumulation model it assumes.
The variations of DOL model parameters are high in the Canada analysis method, as described in "Main differences in the performance functions" section, this leads to severer DOL, hence the higher conservatism. E. g., for dead load alone, the DOL factor in the China wood standard GB 50005 [19] is 0.576 for 50-year duration ("The China reliability analysis method" section), while that in the Canada standard CSA O86 [20] is 0.52 for 30-year duration ("The Canada reliability analysis method" section), as calculated by the Canada reliability analysis. Although variables K A , K P and K B introduce extra variations in the China analysis method, they are overweighed by the smaller DOL value in the Canada analysis method, and the China method still results less conservative design.

Dead load combined with a live load
The results for combination of dead and residential floor live load are similar to that for combination of dead and office floor live load, the results for the latter combination are demonstrated in Figs. 3 and 4. Comparison of the β-γ R * curves at load ratios of ρ = 0.25, 1.0, 2.0 and 4.0 is shown in Fig. 3. It is observed that at a lower ratio of ρ = 0.25, the China curves lie considerably higher, i.e., less conservative, than the corresponding Canada curves. As load ratio ρ increases, the Canada curves gradually approach (when ρ ≤ 2.0) and finally go slightly higher than (ρ = 4.0) the China curves. Influence of load ratio ρ on the β-γ R * curves is shown in Fig. 4. The four clusters of curves all demonstrate a growing trend of β with load ratio ρ. The growth rate gradually decreases, and the curves in each cluster asymptotically approach a final location. The change with respect to load ratio ρ between the Canada curves is more significant than that between the China curves, indicating a higher sensitivity to load ratio ρ.
For combination of dead and snow load, snow load characteristics of Beijing are used as an example. The βγ R * curves for this load combination are shown in Figs. 5 and 6. As reflected in Fig. 5, at a lower ratio of ρ = 0.25, the China curves still lie higher, i.e., less conservative, than the Canada curves, but the Canada curves overtake quickly when load ratio ρ increases and the Canada curves remains a bit higher than China curves when ρ ≥ 2.0. It can be seen from Fig. 6 that for this load combination, the change of the β-γ R * curves with respect to load ratio ρ is obviously less significant than combination of dead and office floor live load as described above. The clusters of the China curves in Fig. 6a, c demonstrate hardly any change with load ratio ρ, and the change of the Canada curves in Fig. 6b, d is also minor compared to the counterparts in Fig. 4.
In general, for combinations of dead load and a single live load, the Canada analysis method gives not as conservative results as it does for the case of dead load alone. When load ratio ρ increases, i.e., when live load makes up an increasing portion of the total load, the Foschi-Yao model in the Canada analysis method captures the changes in stress level and duration of load, resulting in increasingly alleviated DOL. Contrarily, DOL is not altered in the China analysis method in this process. This results in a higher sensitivity of the Canada analysis method to load ratio ρ in comparison with the China analysis method.
It's further observed that differences in the analysis results are not only made by the load ratio ρ but also by the specific type of live load. The mechanism how characteristics of the load effect term affect the reliability index values is further discussed in "Influence from characteristics of the load effect term" section.

Influence from characteristics of the load effect term
As can be seen from the analysis above, the load effect term (more specifically the load ratio ρ as well as the load type) is of significant influence on the relative safety level Influence from statistical characteristics of the load effect term on reliability index is investigated as follows.
In addition to the unified resistance partial factor γ R *, the load effect term in the Canada and the China analysis methods is mainly dependent on items [g(t)r + q(t)]/ (γ G r + γ Q ) and (g + qρ)/(γ G+ γ Q ρ), respectively (see Eqs. (1) and (6)). Note that little difference is made by variable K B (mean = 1.0, COV = 0.05) in the China load effect term. Item (g + qρ)/(γ G+ γ Q ρ) in the China analysis method is analyzed here as an example, its mean and COV are expressed, respectively, as (14) g + qρ where (X) m and V (X) denote the mean and COV of a random variable X, respectively. As load ratio ρ increases, the mean of item (g + qρ)/ (γ G+ γ Q ρ) should start from g m /γ G (when ρ = 0) and asymptotically approach a final value of q m/ γ Q (when ρ → ∞). If g m /γ G = q m/ γ Q , no change is incurred by change of ρ, for this case a critical value of q m can be identified as q m,cr = g m γ Q /γ G . In a similar way, the COV of item (g + qρ)/(γ G+ γ Q ρ) changes from V g (when ρ = 0) to a final value of V q (when ρ → ∞). Note that the mean of item (g + qρ)/(γ G+ γ Q ρ) is affected to some extent by the γ G and γ Q values, while the COV of it is not affected. Assuming γ G = 1.2 and γ Q = 1.4, and taking the statistical parameters of g and q from Table 4, the trend lines of mean and COV of item (g + qρ)/(γ G+ γ Q ρ) with ρ are plotted in Fig. 7a, b. It's observed that the change rate of mean and COV of load effect is large when ρ is small, and the rate diminishes gradually as ρ goes higher. Among the three live loads, snow load is of a mean load ratio (q m = 1.04, see Table 4) closest to the critical value (q m,cr = g m γ Q /γ G = 1.06 × 1.4/1.2 = 1.24), therefore, it triggers a minimum extent of change in the mean load effect as shown in Fig. 7a. A change in the load partial factors, for example, to γ G = 1.25 and γ Q = 1.5 as specified in the Canada standard, makes minor difference which is not altering the general trend of curves in Fig. 7a.
Moreover, for the cases of residential or office floor live load, the beneficial decrease of mean value would overweigh the detrimental increase of COV of the load effect term. As a result, the reliability index resulted from the China analysis method increases with load ratio ρ for these two load combinations, see Fig. 4a, c. The variation rate of the β-γ R * curves with load ratio ρ is gradually diminishing, which matches the trend of variation in Fig. 7. For combination of snow load, the beneficial decrease of the mean would approximately counterbalance the detrimental increase of COV of the load effect term. Therefore, the β-γ R * curves vary little between different load ratios ρ, see Fig. 6a, c.
Similar trends exist in item [g(t)r + q(t)]/(γ G r + γ Q ) of load effect term in the Canada analysis method, except for that variables g(t) and q(t) have some different statistical parameters as listed in Table 5. The beneficial decrease in the mean value of load effect is more significant than that for the China analysis method.

Influence from COV of material strength
In the current analysis, two variation levels of dimension lumber strength are included, viz. COV = 0.410 for bending strength and COV = 0.246 for compression strength of lumber (Table 3). It's straightforward that as COV of a b c d Fig. 4 Influence of load ratio ρ for combination of dead and office floor live load lumber strength decreases, the reliability index β tends to increase in both analysis methods, as reflected in Figs. 2, 3, 4, 5, and 6. However, COV level of lumber strength is not found to make significant difference in the safety level comparison between the two analysis methods, as reflected by comparison between the 'bending' and 'compression' results in Figs. 2, 3 and 5.

Conclusions
The reliability analysis methods underlying China and Canada standards for wood structures are investigated, with special attention paid to the way how DOL is treated. To compare the safety level achieved by the two methods, reliability analysis was carried out on the basis of the same set of material properties and load parameters. Examples in this study reflect that the influence on structural safety of the China and the Canada analysis methods varies from case to case, depending on the variation level of variables in their performance functions. If the same reliability level were given, the Canada analysis method would tend to result in more conservative design than the China analysis method for dead load-dominated cases. This is a result of more significant influence from DOL in the Canada analysis method for this loading condition, which considerably overweighs the influence of extra random variables (K A , K P and K B ) in the China analysis method. As live loads make up an increasing portion of the total load, DOL resulted from the Canada analysis method is considerably alleviated while DOL in the China analysis method is virtually not affected. Therefore, the Canada analysis method results in not as conservative design as it does in the dead load-dominated case. In a b c d In a more general sense, though target reliability index β is an important measure for controlling and unifying the level of risk in construction engineering in a country/ region, it should be carefully used for comparison when different reliability analysis methods are involved (which is typically the case when international design codes are evaluated). Safety level is not only dependent on the target reliability index β adopted by a design standard, but also related to the assumptions and treatments of variables or parameters and specific analysis method employed. Safety level of wood structure design is significantly affected by the inclusion of different random variables, and more importantly, different DOL models in the reliability analysis.
Finally, we would like to conclude that, this study is not to address the real reliability issues in either the China wood standard GB 50005 or the Canada wood standard CSA O86, but to make comparisons of the influence on the safety issues by different reliability analysis methods employing different assumptions and different treatments of some variables. Understanding of the similarities and differences between national wood standards would help enhance international exchanges and collaborations in wood industry.
Abbreviations DOL: Duration of load; COV: Coefficient of variation.