Tag Archives: strain

Jigsaw puzzling without a picture

A350 XWB passes Maximum Wing Bending test

A350 XWB passes Maximum Wing Bending test

Research sometimes feels like putting together a jigsaw puzzle without the picture or being sure you have all of the pieces.  The pieces we are trying to fit together at the moment are (i) image decomposition of strain fields [see ‘Recognising strain’ on October 28th 2015] that allows fields containing millions of data values to be represented by a feature vector with only tens of elements which is useful for comparing maps or fields of predictions from a computational model with measurements made in the real-world; (ii) evaluation of the variation in measurement uncertainty over a field of view of measured displacements or strains in a large structure [see ‘Industrial uncertainty’ on December 12th 2018] which provides information about the quality of the measurements; and (iii) a probabilistic validation metric that provides a measure of how well predictions from a computational model represent measurements made in the real world [see ‘Million to one’ on November 21st 2018].  We have found some of the missing pieces of the jigsaw, for example we have established how to represent the distribution of measurement uncertainty in the feature vector domain [see ‘From strain measurements to assessing El Niño events’ on March 17th 2021] so that it can be used to assess the significance of differences between measurements and predictions represented by their feature vectors – this connects (i) and (ii) together.  Very recently we have demonstrated a generic technique for performing image decomposition of irregularly shaped fields of data or data fields with holes [see Christian et al, 2021] which extends the applicability of our method for comparing measurements and predictions to real-world objects rather than idealised shapes.  This allows (i) to be used in industrial applications but we still have to work out how to connect this to the probabilistic metric in (iii) while also incorporating spatially-varying uncertainty.  These techniques can be used in a wide range of applications, as demonstrated in our recent work on El Niño events [see Alexiadis et al, 2021], because, by treating all fields of data as images, the techniques are agnostic about the source and format of the data.  However, at the moment, our main focus is on their application to ground tests on aircraft structures as part of the Smarter Testing project in collaboration with Airbus, Centre for Modelling & Simulation, Dassault Systèmes, GOM UK Ltd, and the National Physical Laboratory with funding from the Aerospace Technology Institute.  Together we are working towards digital continuity across virtual and physical testing of aircraft structures to provide live data fusion and enable condition-led inspections, test control and validation of computational models.  We anticipate these advances will reduce time and costs for physical tests and accelerate the development of new designs of aircraft that will contribute to global sustainability targets (the aerospace industry has committed to reduce CO2 emissions to 50% of 2005 levels by 2050).  The Smarter Testing project has an ambitious goal which reveals that our pieces of the jigsaw puzzle belong to a small section of a much larger one.

For more on the Smarter Testing project see:




Alexiadis A, Ferson S, Patterson EA. Transformation of measurement uncertainties into low-dimensional feature vector space. Royal Society open science. 8(3):201086, 2021.

Christian WJ, Dean AD, Dvurecenska K, Middleton CA, Patterson EA. Comparing full-field data from structural components with complicated geometries. Royal Society open science. 8(9):210916, 2021.

Image: http://www.airbus.com/galleries/photo-gallery

If you don’t succeed, try and try again…

Photograph of S-shaped plateYou would not think it was difficult to build a thin flat metallic plate using a digital description of the plate and a Laser Powder Bed Fusion (L-PBF) machine which can build complex components, such as hip prostheses.  But it is.  As we have discovered since we started our research project on the thermoacoustic response of additively manufactured parts (see ‘Slow start to an exciting new project on thermoacoustic response of AM metals‘ on September 9th, 2020).  L-PBF involves using a laser beam to melt selected regions of a thin layer of metal powder spread over a flat bed.  The selected regions represent a cross-section of the desired three-dimensional component and repeating the process for each successive cross-section results in the additive building of the component as each layer solidifies.  And there in those last four words lies the problem because ‘as each layer solidifies’ the temperature distribution between the layers causes different levels of thermal expansion that results in strains being locked into our thin plates.  Our plates are too thin to build with their plane surfaces horizontal or perpendicular to the laser beam so instead we build them with their plane surface parallel to the laser beam, or vertical like a street sign.  In our early attempts, the residual stresses induced by the locked-in strains caused the plate to buckle into an S-shape before it was complete (see image).  We solved this problem by building buttresses at the edges of the plate.  However, when we remove the buttresses and detach the plate from the build platform, it buckles into a dome-shape.  Actually, you can press the centre of the plate and make it snap back and forth noisily.  While we are making progress in understanding the mechanisms at work, we have some way to go before we can confidently produce flat plates using additive manufacturing that we can use in comparisons with our earlier work on the performance of conventionally, or subtractively, manufactured plates subject to the thermoacoustic loading experienced by the skin of a hypersonic vehicle [see ‘Potential dynamic buckling in hypersonic vehicle skin‘ on July 1st 2020) or the containment walls in a fusion reactor.  Sometimes research is painfully slow but no one ever talks about it.  Maybe because we quickly forget the painful parts once we have a successful outcome to brag about. But it is often precisely the painful repetitions of “try and try again” that allow us to reach the bragging stage of a successful outcome.

The research is funded jointly by the National Science Foundation (NSF) in the USA and the Engineering and Physical Sciences Research Council (EPSRC) in the UK (see Grants on the Web).


Silva AS, Sebastian CM, Lambros J and Patterson EA, 2019. High temperature modal analysis of a non-uniformly heated rectangular plate: Experiments and simulations. J. Sound & Vibration, 443, pp.397-410.

Magana-Carranza R, Sutcliffe CJ, Patterson EA, 2021, The effect of processing parameters and material properties on residual forces induced in Laser Powder Bed Fusion (L-PBF). Additive Manufacturing. 46:102192

From strain measurements to assessing El Niño events

Figure 11 from RSOS 201086One of the exciting aspects of leading a university research group is that you can never be quite sure where the research is going next.  We published a nice example of this unpredictability last week in Royal Society Open Science in a paper called ‘Transformation of measurement uncertainties into low-dimensional feature space‘ [1].  While the title is an accurate description of the contents, it does not give much away and certainly does not reveal that we proposed a new method for assessing the occurrence of El Niño events.  For some time we have been working with massive datasets of measurements from arrays of sensors and representing them by fitting polynomials in a process known as image decomposition [see ‘Recognising strain‘ on October 28th, 2015]. The relatively small number of coefficients from these polynomials can be collated into a feature vector which facilitates comparison with other datasets [see for example, ‘Out of the valley of death into a hype cycle‘ on February 24th, 2021].  Our recent paper provides a solution to the issue of representing the measurement uncertainty in the same space as the feature vector which is roughly what we set out to do.  We demonstrated our new method for representing the measurement uncertainty by calibrating and validating a computational model of a simple beam in bending using data from an earlier study in a EU-funded project called VANESSA [2] — so no surprises there.  However, then my co-author and PhD student, Antonis Alexiadis went looking for other interesting datasets with which to demonstrate the new method.  He found a set of spatially-varying uncertainties associated with a metamodel of soil moisture in a river basin in China [3] and global oceanographic temperature fields collected monthly over 11 years from 2002 to 2012 [4].  We used the latter set of data to develop a new technique for assessing the occurrence of El-Niño events in the Pacific Ocean.  Our technique is based on global ocean dynamics rather than on the small region in the Pacific Ocean which is usually used and has the added advantages of providing a confidence level on the assessment as well as enabling straightforward comparisons of predictions and measurements.  The comparison of predictions and measurements is a recurring theme in our current research but I did not expect it to lead into ocean dynamics.

Image is Figure 11 from [1] showing convex hulls fitted to the cloud of points representing the uncertainty intervals for the ocean temperature measurements for each month in 2002 using only the three most significant principal components . The lack of overlap between hulls can be interpreted as implying a significant difference in the temperature between months.


[1] Alexiadis, A. and Ferson, S. and  Patterson, E.A., , 2021. Transformation of measurement uncertainties into low-dimensional feature vector space. Royal Society Open Science, 8(3): 201086.

[2] Lampeas G, Pasialis V, Lin X, Patterson EA. 2015.  On the validation of solid mechanics models using optical measurements and data decomposition. Simulation Modelling Practice and Theory 52, 92-107.

[3] Kang J, Jin R, Li X, Zhang Y. 2017, Block Kriging with measurement errors: a case study of the spatial prediction of soil moisture in the middle reaches of Heihe River Basin. IEEE Geoscience and Remote Sensing Letters, 14, 87-91.

[4] Gaillard F, Reynaud T, Thierry V, Kolodziejczyk N, von Schuckmann K. 2016. In situ-based reanalysis of the global ocean temperature and salinity with ISAS: variability of the heat content and steric height. J. Climate. 29, 1305-1323.

Poleidoscope (=polariscope + kaleidoscope)

A section from a photoelastic model of turbine disc with a single blade viewed in polarised light to reveal the stress distribution.Last month I wrote about the tedium of collecting data 35 years ago without digital instrumentation and how it led me to work on automation and digitalisation in experimental mechanics [see ‘35 years later and still working on a PhD thesis‘ on September 16th, 2020].  Thirty years ago, one of the leading methods for determining stresses in components was photoelasticity, which uses polarised light to generate fringe patterns in transparent components or models that correspond to the distribution of stress.  The photoelastic fringes can be analysed in a polariscope, of which the basic principles are explained in a note at the end of this post.  During my PhD, I took hundreds of black and white photographs in a polariscope using sheets of 4×5 film, which came in boxes of 25 sheets that you can still buy, and then scanned these negatives using a microdensitometer to digitise the position of the fringes.  About 15 years after my PhD, together with my collaborators, I patented the poleidoscope which is a combination of a polariscope and a kaleidoscope [US patents 6441972 & 5978087] that removes all of that tedium.  It uses the concept of the multi-faceted lens in a child’s kaleidoscope to create several polariscopes within a compound lens attached to a digital camera.  Each polariscope has different polarising elements such that photoelastic fringes are phase-shifted between the set of images generated by the multi-faceted lens.  The phase-shifted fringe patterns can be digitally processed to yield maps of stress much faster and more reliably than any other method.  Photoelastic stress analysis is no longer popular in mainstream engineering or experimental mechanics due to the simplicity and power of digital image correlation [see ‘256 shades of grey‘ on January 22nd, 2014]; however, the poleidoscope has found a market as an inspection device that provides real-time information on residual stresses in glass sheets and silicon wafers during their production.  In 2003, I took study leave for the summer to work with Jon Lesniak at Glass Photonics in Madison, Wisconsin on the commercialisation of the poleidoscope.  Subsequently, Glass Photonics have  sold more than 250 instruments worldwide.

For more information on the poleidoscope see: Lesniak JR, Zhang SJ & Patterson EA, The design and evaluation of the poleidoscope: a novel digital polariscope, Experimental Mechanics, 44(2):128-135, 2004

Note on the Basic principles of photoelasticity: At any point in a loaded component there is a stress acting in every direction. The directions in which the stresses have the maximum and minimum values for the point are known as principal directions. The corresponding stresses are known as maximum and minimum principal stresses. When polarised light enters a loaded transparent component, it is split into two beams. Both beams travel along the same path, but each vibrates along a principal direction and travels at a speed proportional to the associated principal stress. Consequently, the light emerges as two beams vibrating out of phase with one another which when combined produce an interference pattern.   The polarised light is produced by the polariser in the polariscope and the analyser performs the combination. The interference pattern is observed in the polariscope, and the fringes are contours of principal stress difference which are known as isochromatics. When plane polarised light is used black fringes known as isoclinics are superimposed on the isochromatic pattern. Isoclinics indicate points at which the principal directions are aligned to the polarising axes of the polariser and analyser.

Image: a section from a photoelastic model of turbine disc with a single blade viewed in polarised light to reveal the stress distribution.