12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 1 Multivariate Modelling of Extreme Load Combinations for Wind Turbines Nikolay Dimitrov Research Scientist, Dept. of Wind Energy, Technical University of Denmark, Denmark ABSTRACT: We demonstrate a model for estimating the joint probability distribution of two load components acting on a wind turbine blade cross section. The model addresses the problem of modelling the probability distribution of load time histories with large periodic components by dividing the signal into a periodic part and a perturbation term, where each part has a known probability distribution. The proposed model shows good agreement with simulated data under stationary conditions, and a design load envelope based on this model is comparable to the load envelope estimated using the standard procedure for determining contemporaneous loads. By defining a joint probability distribution and full return-period contours for multiple load components, the suggested procedure gives the possibility for determining the most critical loading direction in a blade cross section, or for carrying out reliability analysis on an entire cross section. 1. INTRODUCTION Wind turbines as moving machinery subjected to random environmental influence experience loads and deformations along multiple degrees of freedom in the structure. A fully three-dimensional stress analysis of the entire structure under dynamic loading is extremely expensive computationally. Therefore, the design loads on wind turbines are usually estimated by dynamic finite-element or modal simulations using a simplified representation of the structure in terms of beam elements. Since each node in a three-dimensional beam element has six degrees of freedom, any given load condition at a beam node, i.e., a cross section of a turbine component, will be a combination of six load components. Due to the non-symmetric geometry of some of the turbine components such as the blades, using only the vector magnitude of the loads is not sufficient as the maximum allowable loads will vary for different load directions. In current wind turbine design guidelines (i.e., IEC 61400-1, ed.3 from 2005) the problem of finding the contemporaneous values of different load components is addressed in a relatively simple way and methods for estimating the actual joint probabilities are not available. The present study demonstrates a multivariate model for the joint probability distribution of contemporaneous extreme loads on wind turbines. The model considers the load components as a combination of a cyclic, gravity-driven component and a randomly distributed perturbation. The subtraction of the sinusoidal signal from the random data results in better possibilities for assessing the correlations between loads in different directions, and thus leads to an improved model for the joint distribution of loads. The parameters of the marginal distributions are estimated from dynamic simulations of the wind turbine behavior under turbulent incoming wind. The joint distribution of the perturbations in different directions is obtained from the marginals using the Nataf multivariate transformation (Liu and Der Kiureghian, 1986), but can also be done with a Rosenblatt transformation (Rosenblatt, 1952), or principal component analysis (Natarajan and Velerst, 2011). A joint distribution for the loads is then obtained by convolution of the sinusoidal signal with the perturbations and integration over the rotor azimuth angles. The obvious 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 2 application of this distribution model is to define the return period contours for extreme load combinations and thus the design load envelope, but it also gives the possibility for carrying out reliability analysis of wind turbine components. 2. EXTREME LOADS FROM AEROELASTIC SIMULATIONS 2.1. Input description The limit state considered in this study is the multi-axial extreme load acting on a cross section from a wind turbine blade. The wind turbine used as an example is the DTU 10MW reference wind turbine (Bak et al., 2013) with rotor diameter of 179m and hub height of 119m. The loads are determined from aeroelastic simulations using the Hawc2 code (Larsen and Hansen, 2012). Fitting a representative probabilistic model of the loads is only feasible under stationary conditions, i.e., constant mean wind speed, turbulence, etc., and if long-term load distributions under varying conditions are required, they are estimated by fitting a number of short-term distributions under stationary conditions, followed by integration over the respective distribution of long-term variables (see e.g. the IEC61400-1, ed.3 standard, 2005). For simplicity we consider stationary time series with a fixed mean wind speed of 11m/s, turbulence intensity 0.176 (IEC turbulence class B), and wind shear exponent of 0.2. The loads on a blade cross section are also influenced by the blade pitch angle, where an increase in the pitch angle normally leads to a reduction in the aerodynamic thrust acting on the rotor, and hence to reduction in the flapwise loads. We choose to take into account the time periods in the simulations when the blade pitch angle is nearly zero (less than 0.5β), which can be considered as a worst-case scenario since the flapwise loads are highest at zero pitch angle and 11m/s is the wind speed at which the highest extreme loads are typically attained. This set of conditions is sufficient for demonstrating the methods for modelling of the joint distribution of extreme loads. However, for design purposes it will be relevant to fit a number of short-term distributions under different wind conditions and taking account of different blade pitch angles, and then estimate the design load envelope by integrating the short-term distributions into a long-term distribution. 2.2. Extraction of load extremes The ultimate limit state design loads for a wind turbine are specified as the extreme loads with a given long period of recurrence (i.e. 1 year or 50 years). Estimating the distribution of the extremes requires extracting independent peaks from the time series under consideration. In the multivariate case, the definition of an extreme event is more complicated as all related load components have to be taken into consideration. We adopt an approach where a point in time π‘ = π is considered a load peak if it satisfies the following two conditions: 1) There is a load reversal of the flapwise load ππ₯ at π‘ = π, i.e, |ππ₯(π)| > |ππ₯(π β Ξπ‘)| and |ππ₯(π)| > |ππ₯(π + Ξπ‘)| 2) There is no larger peak within 1s before and after π‘ = π: maxMx(π‘)(π‘ β [π β 1π , π + 1π ]) = ππ₯(π). Figure 1: Example of load peaks extracted from simulated time history for blade bending moments ππ₯ and ππ¦. Figure 1 shows an example of the time series for the resultant of the flapwise and 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 3 edgewise bending moments ππ = οΏ½ππ₯2 +ππ¦2 and the peaks identified using the criterion given above. When a probability distribution is fit to the load peaks extracted with the above procedure, the probability values will be relative to the number of peaks used to fit the distribution. This is a so-called local distribution of peaks. In order to determine the actual short-term distribution of peaks with respect to a specific reference period, π, the local distribution has to be raised to the power π(π), where π is the expected number of load peaks within one time period π: πΉπ βπππ‘βπ‘π‘ππ‘(π|π,π) = οΏ½πΉπππππ(π|π)οΏ½π(π) (1) For example, if the desired reference time period is π = 10πππ (the typical duration of a single aeroelastic simulation) , π(π) equals the number of load peaks extracted per time series. 3. JOINT PROBABILITY DISTRIBUTION OF MULTI-AXIAL LOADS For the purpose of detailed reliability analysis or structural design of a wind turbine blade section in general, we need the probability of observing a given combination of load components acting on the cross section. The components are typically given in six degrees of freedom; that is moments and forces acting in three axes of a Cartesian coordinate system. The most critical of these six components are the transverse bending moments ππ₯ and ππ¦ (flapwise and edgewise bending moments). The other four components can either be sufficiently represented as functions of ππ₯ and ππ¦, or their magnitude is sufficiently small as to not influence the design significantly. Therefore, the major part of the problem of determining the probability of given load combination amounts to finding the joint distribution of the extreme flapwise and edgewise loads, πππ₯ππ¦: πππ₯ππ¦ = πππ₯(ππ₯)πππ¦οΏ½ππ¦οΏ½ππ₯οΏ½= πππ¦οΏ½ππ¦οΏ½πππ₯οΏ½ππ₯οΏ½ππ¦οΏ½ (2) Figure 2 shows an example of the joint probability density of the extremes of ππ₯ and ππ¦ for a typical cross section in a wind turbine blade. Figure 2: Joint pdf of flapwise bending moment ππ₯ and edgewise bending moment ππ¦ for a wind turbine blade cross section. Based on 100 Hawc2 simulations for a 10MW turbine at 11m/s wind speed, section located 46.6m from the blade root. 3.1. Marginal distributions of load components The edgewise loads are mainly gravity-driven and as such have a strong sinusoidal component caused by the rotation of the turbine. As a result, the marginal distribution of ππ¦ is bimodal (Figure 3). The presence of the sinusoidal signal means that ππ¦ can be considered as the sum of a sinusoidal function ππ¦ with amplitude π and phase π, and a random βnoiseβ component ππ¦: ππ¦ = ππ¦ + π β sin(ππ‘ + π) + ππ¦= ππ¦ + ππ¦ (3) where ππ¦ is the mean moment, π is the angular velocity, π‘ is the time, and ππ‘ = πΎ is the rotor azimuth. The sinusoidal signal considered as a random variable can be represented either as a sinewave-function of a uniformly-distributed azimuth variable, or as a standalone random variable using the arcsine distribution: 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 4 πππ¦(ππ¦) =1ποΏ½π2 β οΏ½ππ¦ β ππ¦οΏ½2 (4) The perturbation term ππ¦ can be represented by any suitable statistical distribution, and in the case when the sinusoidal component is much bigger than the noise component, a normal distribution for ππ¦ seems to be a good choice. The probability distribution of the edgewise moment ππ¦ is then considered as the distribution of the sum of two independent random variables, given by the convolution of the marginal distribution functions of ππ¦ and ππ¦: πππ¦+ππ¦(ππ¦)= οΏ½ πππ¦(π)πππ¦οΏ½ππ¦ β ποΏ½ππβββ (5) An example fit of the distribution model from equation (4) to data from aeroelastic simulations is shown on Figure 3. Figure 3: Marginal probability density function of the extreme edgewise bending moment ππ¦: comparison between simulations and a distribution model based on convolution of a sinusoidal function and a normal distribution. The flapwise bending moment is mainly driven by turbulence, however there is also a small sinusoidal component due to the blade pitch axis not being exactly in the rotor plane. As a result, the distribution of ππ₯ can also be modelled as the convolution of a sinusoidal distribution and a perturbation, where the perturbation term will have the dominating contribution, and in this case is modelled by the Generalized Extreme Value (GEV) distribution. Figure 4 compares the fitted distribution for ππ₯ to the empirical distribution from observations. Since for the flapwise load we expect that the distribution tails have higher importance than the values close to the mean, the comparison is shown in terms of the cumulative distribution function πΉππ₯(ππ₯) and its complement: πππ₯(ππ₯) = 1 β πΉππ₯(ππ₯). Figure 4: Marginal distribution function of the extreme flapwise bending moment ππ₯: comparison between simulations and a distribution model based on convolution of a sinusoidal function and Generalized Extreme Value distribution. In the present case the GEV distribution fits reasonably well both the higher and lower tails of the marginal distribution of the flapwise moment, ππ₯ . However, this is not necessarily always true, and in some other cases more advanced fitting may be required, e.g. separate fits to the upper and lower tails of the distribution, or other approaches. 3.2. Modelling of the joint distribution of flapwise and edgewise blade moments Several methods for obtaining design load combinations based on the time series for multiple load components have been employed in previous studies and in design codes: 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 5 β’ The IEC 61400-1 standard suggests a simple technique where the load combination case combines the highest extreme flapwise load with the mean edgewise load, or in a more conservative manner the combination can also consist of the extremes for both flapwise and edgewise loads. β’ The joint distribution of ππ₯ and ππ¦ is expressed in terms of the marginal distributions of uncorrelated standard variables by the use of a multivariate transformation such as the Nataf transform, the Rosenblatt transformation, or Principal Component Analysis (see e.g. Natarajan et al.). Obviously, the procedure recommended in the IEC 61400-1 standard is simple, however not necessarily providing the designer with the right information. The methods employing multivariate transformations also have limitations: the Rosenblatt transformation requires fitting of a number of conditional distributions, which limits its usability for very high load values where the number of observations is small; and the Nataf transform and Principal Component Analysis require the use of the correlation matrix with Pearsonβs correlation coefficients, which are only representative for data with linear correlation. As seen on Figure 1, this is not the case for ππ₯ and ππ¦. Figure 5 demonstrates the limitations of this approach, showing a comparison of the joint PDF of ππ₯ and ππ¦ observed from simulations to the joint PDF obtained by using a Nataf transformation (Liu & Der Kiureghian, 1998): πππ₯ππ¦(ππ₯,ππ¦)= πππ₯(ππ₯)πππ¦(ππ¦)ππ§(π³,πβ²)π(π§1)π(π§2) (6) where π§π = Ξ¦β1οΏ½πΉππ(ππ)οΏ½ , π = π₯,π¦ , and ππ§(π³,πβ²) is the n-dimensional normal PDF for variables with zero mean, unit variance and correlation matrix πβ² . The correlation coefficients πππβ² in πβ² are function of the correlation coefficients of the data, πππ, and can be determined by an iterative solution or, for some distributions, calculated from tables, see Liu & Der Kiureghian for details. The marginal distributions of ππ₯ and ππ¦ are calculated based on the convolution rule as defined in Equation (5): πππ₯ = πππ₯+ππ₯ and πππ¦ = πππ¦+ππ¦. Figure 5: Comparison of the joint probability density for blade flapwise and edgewise extreme bending moments, estimated from data and by directly using a multivariate transformation on the marginal distributions of ππ₯ and ππ¦. The transformation used to create Figure 5 requires the use of the correlation coefficient between ππ₯ and ππ¦, which according to the data is πππ₯ππ¦ = β0.24 . However, as mentioned earlier Pearsonβs correlation coefficient describes well only linearly dependent variables, which is not the case for the present problem β and subsequently the multivariate distribution based on this approach does not fit the data well. The nonlinearity in our case is caused by the bi-modal behavior of the loads due to the cyclic components in them, and once the sinusoidal part is eliminated, the dependence between the perturbation terms ππ₯ and ππ¦ is well described by a linear correlation coefficient. As a result, the joint distribution πππ₯ππ¦ can be modelled adequately by with Equation 6, as Figure 6 shows. For comparison, the correlation coefficient πππ₯ππ¦ is β0.69 , significantly larger than πππ₯ππ¦. 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 6 Figure 6: Joint probability density function of the noise terms ππ₯ and ππ¦. Comparison between simulations and a multivariate distribution fit. Based on this observation and knowing that the bending moments ππ₯ and ππ¦ can be expressed as function of the perturbation terms ππ₯, ππ¦ and the rotor azimuth πΎ, we propose that the joint distribution πππ₯ππ¦ can be modelled using the joint distribution of ππ₯, ππ¦ and πΎ: πππ₯ππ¦οΏ½ππ₯,ππ¦οΏ½ == οΏ½ ππΎ(πΎ)πππ¦οΏ½ππ¦οΏ½πΎοΏ½πππ₯(ππ₯|πΎ, ππ¦)2π0ππΎ= οΏ½ ππΎ(πΎ)πππ¦ οΏ½ππ¦οΏ½ππ¦οΏ½πΎ,ππ¦οΏ½οΏ½2π0β πππ₯οΏ½ππ₯οΏ½ππ₯(πΎ,ππ₯), ππ¦οΏ½ππΎ (7) where πππ₯ππ¦ = οΏ½ ππΎππ₯ππ¦οΏ½πΎ, ππ₯, ππ¦οΏ½ππΎ2π0= οΏ½ ππΎ(πΎ)πππ¦οΏ½ππ¦οΏ½πΎοΏ½πππ₯(ππ₯|πΎ, ππ¦)ππΎ2π0 (8) This approach results in a significantly improved match to the joint probability distribution of ππ₯ and ππ¦, which is shown in Figure 7. Figure 7: Comparison of the joint probability density for blade flapwise and edgewise extreme bending moments, estimated from data and by using a multivariate model based on the marginal distributions of perturbation terms ππ₯ and ππ¦ and integrated over the rotor azimuth πΎ. 3.3. Return period contours for the joint distribution of ππ₯ and ππ¦ Return period contours are evaluated using the inverse-FORM (i-FORM) procedure, Winterstein et al. (1993). The i-FORM calculations are carried out in terms of standard normal, uncorrelated variables, meaning that for using the method, the distributions of ππ₯ and ππ¦ have to undergo two transformations: 1) transformation to a pair of dependent standard normal variables π1 and π2 , and 2) to account for dependency and transform π1 and π2 into a pair of independent standard variables, π1 and π2. The latter transformation is done using Cholesky decomposition of the correlation matrix πβ² (derived from the correlation matrix of the data π via the Nataf transformation), obtaining a transformation matrix π where ππβ² = πβ² , and the final result is found by π = ππ. In order to verify the above calculation, a large number of time series under the same conditions described in Section 2 were simulated for obtaining a reference. Based on a total of 30000 time series (approximately equal to 7 months) under stationary conditions, load combination contours corresponding to 1-day and 1-month return periods are estimated. Figure 8 shows the 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 7 comparison between the reference simulations and the joint distribution model. The agreement between the model and the reference data is good in general, with the exception of a slight deviation in the region of maximal flapwise load. Considering that the model is based on just 100 time series compared to the 30000 used for the reference calculation, at least part of the disagreement can be attributed to statistical uncertainty. Figure 8: Probability contours corresponding to 1-day and 1-month return periods, estimated with the multivariate extreme load model, compared to contours estimated from 30000 time series. 3.4. Comparison to contemporaneous loads calculation The IEC 61400-1 standard, ed.3 first amendment (2010) recommends a simplified procedure for obtaining contemporaneous loads, where the characteristic extreme value for a given load component is found as the mean of the maxima from each realization (e.g. turbulence seed), and the contemporaneous load for other components is found as the mean of the contemporaneous values from each realization. The design load is then calculated by multiplying the characteristic loads by a safety factor of 1.35. Considering that extreme design load cases in the IEC61400-1 are normally tuned to represent events with 50-years recurrence periods, on Figure 9 we compare the contemporaneous design loads obtained using the IEC-procedure to the 50-year return period contour estimated using the proposed joint distribution model. As the figure shows, the contemporaneous design loads based on the standard procedure agree with the 50-year joint probability contour, but with a reduced safety factor of 1.25 instead of the standard 1.35. This observation however can only be verified for the simulation conditions used for the present study. For other conditions (e.g. load cases with extreme turbulence as DLC1.3 which are often design-driving for blades), the standard contemporaneous loads and the joint probability distribution model might compare differently. Figure 9: Comparison of the 1-year and 50-year return period contours estimated by the multivariate model to the design envelope determined using the standard contemporaneous loads procedure. The value of the safety factor used is tuned so the loads multiplied with it coincide with the 50-year contour. 4. DISCUSSION AND CONCLUSIONS In the present paper, we demonstrated a model for estimating the joint probability distribution of two wind turbine blade load components. The study addressed the problem of modelling the probability distribution of load time histories with large periodic components by dividing the signal into a periodic part and a perturbation term, where each part has a known probability distribution. The proposed model shows good agreement with simulated data under stationary conditions, and a design load envelope based on 12th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP12 Vancouver, Canada, July 12-15, 2015 8 this model is comparable to the load envelope estimated using the standard procedure for determining contemporaneous loads. By defining a joint probability distribution and full return-period contours for multiple load components, the suggested procedure gives the possibility for determining the most critical loading direction in a blade cross section, or for carrying out reliability analysis on an entire cross section. This makes the model a potentially valuable tool to wind turbine designers as it provides additional information compared to what is learned by assessing just a few selected loading directions. ACKNOWLEDGEMENTS The work presented herein is a part of the Danish Energy Technology Development and Demonstration (EUDP) project titled, βDemonstration of a basis for tall wind turbine designβ, Project no 64011-0352. The financial support is greatly appreciated. REFERENCES Liu, P. L and Der Kiureghian, A. (1986) βMultivariate distribution models with prescribed marginals and covariancesβ Probabilistic Engineering Mechanics, 1(2), 105-112 Rosenblatt, M. (1952) βRemarks on a multivariate transformationβ Annals of Mathematical Statistics, 23, 470-472 Natarajan, A. and Velerst, D. R. (2011) βOutlier robustness for wind turbine extrapolated extreme loadsβ, Wind Energy, DOI 10.1002/we.497 Bak, C., Zahle, F., Bitsche, R., Kim, T., Yde, A., Henriksen, L.C., Natarajan, A., Hansen, M. (2013) βDescription of the DTU 10 MW reference wind turbineβ, DTU Wind Energy Report-I-0092 Larsen, T. J., Hansen, A. M. (2012) βHow to HAWC2, the userβs manual.β Tech. Rep. RisΓΈ-R-1597(ver.4-3) (EN), DTU Wind Energy, Roskilde, Denmark, April 2012 IEC. (2005) βInternational Standard IEC 61400-1: Wind Turbines β Part 1: Design Guidelinesβ Winterstein, S. R., Ude, T. C., Cornell, C. A., Bjerager, P., Haver, S. (1993) βEnvironmental parameters for extreme response: inverse FORM with omission factorsβ, in Proceedings of the ICOSSAR-93 conference, Innsbruck, Austria. IEC. (2010) βInternational Standard IEC61400-1, Ammendment 1. Wind Turbines - Part 1: Design Guidelinesβ
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015) /
- Multivariate modelling of extreme load combinations...
Open Collections
International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP) (12th : 2015)
Multivariate modelling of extreme load combinations for wind turbines Dimitrov, Nikolay Jul 31, 2015
pdf
Page Metadata
Item Metadata
Title | Multivariate modelling of extreme load combinations for wind turbines |
Creator |
Dimitrov, Nikolay |
Contributor | International Conference on Applications of Statistics and Probability (12th : 2015 : Vancouver, B.C.) |
Date Issued | 2015-07 |
Description | We demonstrate a model for estimating the joint probability distribution of two load components acting on a wind turbine blade cross section. The model addresses the problem of modelling the probability distribution of load time histories with large periodic components by dividing the signal into a periodic part and a perturbation term, where each part has a known probability distribution. The proposed model shows good agreement with simulated data under stationary conditions, and a design load envelope based on this model is comparable to the load envelope estimated using the standard procedure for determining contemporaneous loads. By defining a joint probability distribution and full return-period contours for multiple load components, the suggested procedure gives the possibility for determining the most critical loading direction in a blade cross section, or for carrying out reliability analysis on an entire cross section. |
Genre |
Conference Paper |
Type |
Text |
Language | eng |
Notes | This collection contains the proceedings of ICASP12, the 12th International Conference on Applications of Statistics and Probability in Civil Engineering held in Vancouver, Canada on July 12-15, 2015. Abstracts were peer-reviewed and authors of accepted abstracts were invited to submit full papers. Also full papers were peer reviewed. The editor for this collection is Professor Terje Haukaas, Department of Civil Engineering, UBC Vancouver. |
Date Available | 2015-05-21 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0076153 |
URI | http://hdl.handle.net/2429/53321 |
Affiliation |
Non UBC |
Citation | Haukaas, T. (Ed.) (2015). Proceedings of the 12th International Conference on Applications of Statistics and Probability in Civil Engineering (ICASP12), Vancouver, Canada, July 12-15. |
Peer Review Status | Unreviewed |
Scholarly Level | Faculty |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 53032-Paper_372_Dimitrov.pdf [ 370.82kB ]
- Metadata
- JSON: 53032-1.0076153.json
- JSON-LD: 53032-1.0076153-ld.json
- RDF/XML (Pretty): 53032-1.0076153-rdf.xml
- RDF/JSON: 53032-1.0076153-rdf.json
- Turtle: 53032-1.0076153-turtle.txt
- N-Triples: 53032-1.0076153-rdf-ntriples.txt
- Original Record: 53032-1.0076153-source.json
- Full Text
- 53032-1.0076153-fulltext.txt
- Citation
- 53032-1.0076153.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.53032.1-0076153/manifest