Measuring and Evaluation of the Shape and Size of Initial Geometrical Imperfections of Structural Glass Members

Similar documents
Behaviour of cold bent glass plates during the shaping process

WP1 NUMERICAL BENCHMARK INVESTIGATION

PARAMETERS INFLUENCING THE PRECISION OF SLM PRODUCTION

CHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force

THE EFFECT OF THE FREE SURFACE ON THE SINGULAR STRESS FIELD AT THE FATIGUE CRACK FRONT

ksa MOS Ultra-Scan Performance Test Data

Influence of geometric imperfections on tapered roller bearings life and performance

Investigating the influence of local fiber architecture in textile composites by the help of a mapping tool

Introduction. Section 3: Structural Analysis Concepts - Review

Background CE 342. Why RISA-2D? Availability

Non-Linear Analysis of Bolted Flush End-Plate Steel Beam-to-Column Connection Nur Ashikin Latip, Redzuan Abdulla

Figure 30. Degrees of freedom of flat shell elements

Case Study - Vierendeel Frame Part of Chapter 12 from: MacLeod I A (2005) Modern Structural Analysis, ICE Publishing

Electrically tunable large aperture lens EL TC-VIS-20D

D DAVID PUBLISHING. Stability Analysis of Tubular Steel Shores. 1. Introduction

Advanced Professional Training

Determination of the aperture of the LHCb VELO RF foil

Plate Buckling ref. Hughes Chapter 12

Chapter 3 Analysis of Original Steel Post

General Physics II. Mirrors & Lenses

FE ANALYSES OF STABILITY OF SINGLE AND DOUBLE CORRUGATED BOARDS

Buckling design of axially compressed steel silos on discrete supports

Chapter 7: Geometrical Optics

EXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY

A Finite Element Solution for Lateral Track Buckling Problem.

Chapter 7: Geometrical Optics. The branch of physics which studies the properties of light using the ray model of light.

Electrically tunable large aperture lens EL TC

Mirror Example Consider a concave mirror radius -10 cm then = = Now consider a 1 cm candle s = 15 cm from the vertex Where is the image.

STRENGTH ANALYSIS OF PIN CONNECTIONS USING COMPUTER AIDED SYSTEMS

Application nr. 2 (Global Analysis) Effects of deformed geometry of the structures. Structural stability of frames. Sway frames and non-sway frames.

Electrically tunable large aperture lens EL TC

ME 475 FEA of a Composite Panel

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

Laminates can be classified according to the fiber orientation.

Finite Element Buckling Analysis Of Stiffened Plates

Nicholas J. Giordano. Chapter 24. Geometrical Optics. Marilyn Akins, PhD Broome Community College

DESIGN AND ANALYSIS OF MEMBRANE STRUCTURES IN FEM-BASED SOFTWARE MASTER THESIS

Reduction of squeal on laminated brake disc fastened with distributed contact pressure

TABLE OF CONTENTS SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 SECTION 3 WAVE REFLECTION AND TRANSMISSION IN RODS Introduction...

Using three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model

SAP 2000 Problem II Homework Problem P5.45. Recall from Lab #6 the Global and Local Reference Coordinate Systems for 2D Problems

Chapter 6. THE NORMAL DISTRIBUTION

Analysis of Composite Aerospace Structures Finite Elements Professor Kelly

Physics 1C Lecture 26A. Beginning of Chapter 26

AP Physics: Curved Mirrors and Lenses

Distance Between Two Snaked-lay of Subsea Pipeline. Yuxiao Liu1, a *

D&S Technical Note 09-2 D&S A Proposed Correction to Reflectance Measurements of Profiled Surfaces. Introduction

Example 1: Give the coordinates of the points on the graph.

SAFI Sample Projects. Design of a Steel Structure. SAFI Quality Software Inc. 3393, chemin Sainte-Foy Ste-Foy, Quebec, G1X 1S7 Canada

EULER-TYPE ALGORITHM APPLIED TO LEVELING PROCESS ANALYSIS

Final project: Design problem

Lindab We simplify construction. Table of contents

Development of Backhoe Machine By 3-D Modelling using CAD Software and Verify the Structural Design By using Finite Element Method

A new Aluminium Structural Analysis and Design Methodology

MODELING AND ANALYSIS OF LATTICE TOWERS WITH MORE ACCURATE MODELS

Simulating Reinforced Concrete Beam-Column Against Close-In Detonation using S-ALE Solver

Chapter 6. THE NORMAL DISTRIBUTION

An introduction to understanding the deflection of beams using the Push Me Pull Me models on Expedition Workshed

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =

Design and Analysis of a Single Notch Parallelogram Flexure Mechanism Based X-Y Nanopositioning Stage

APPROACHES TO THE PROCESS CAPABILITY ANALYSIS IN THE CASE OF NON- NORMALLY DISTRIBUTED PRODUCT QUALITY CHARACTERISTIC

VIRTUAL TESTING OF AIRCRAFT FUSELAGE STIFFENED PANELS

Normal Data ID1050 Quantitative & Qualitative Reasoning

INFLUENCE OF THE HIGH SPEED PART POLISHING COATING TOOLS ON THE CHANGES SURFACE TEXTURE

Element Order: Element order refers to the interpolation of an element s nodal results to the interior of the element. This determines how results can

Topic Training Load generators

Predicting the mechanical behaviour of large composite rocket motor cases

17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES

Modelling Flat Spring Performance Using FEA

Aufgabe 1: Dreipunktbiegung mit ANSYS Workbench

Exercise 1. 3-Point Bending Using the Static Structural Module of. Ansys Workbench 14.0

What makes Bolt Self-loosening Predictable?

Optimization of Steel Storage Rack Pothook Hole Base on ANSYS Workbench

Physics 1202: Lecture 17 Today s Agenda

Manual BS : 2000 steel code check

Available online at ScienceDirect. Procedia Engineering 161 (2016 ) Bohdan Stawiski a, *, Tomasz Kania a

BUCKLING COEFFICIENTS FOR SIMPLY SUPPORTED AND CLAMPED FLAT, RECTANGULAR SANDWICH PANELS UNDER EDGEWISE COMPRESSION

OPTIMIZATION OF STIFFENED LAMINATED COMPOSITE CYLINDRICAL PANELS IN THE BUCKLING AND POSTBUCKLING ANALYSIS.

Numerical Calculations of Stability of Spherical Shells

Waves & Oscillations

Application of Shell elements to buckling-analyses of thin-walled composite laminates

Multiframe Oct 2008

This document gives a detailed summary of the new features and improvements of FEM- Design version 16.

Standard dimension optimization of steel frames

13. Brewster angle measurement

THREE DIMENSIONAL ACES MODELS FOR BRIDGES

INTRODUCTION REFLECTION AND REFRACTION AT BOUNDARIES. Introduction. Reflection and refraction at boundaries. Reflection at a single surface

Optics II. Reflection and Mirrors

Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure

Advanced model of steel joints loaded by internal forces from 3D frame structures

Frequency Distributions

Module 3: Buckling of 1D Simply Supported Beam

Effect of Subdivision of Force Diagrams on the Local Buckling, Load-Path and Material Use of Founded Forms

Keywords: flows past a cylinder; detached-eddy-simulations; Spalart-Allmaras model; flow visualizations

Wave Propagation Correlations between Finite Element Simulations and Tests for Enhanced Structural Health Monitoring

Chapter 36. Image Formation

Configuration Optimization of Anchoring Devices of Frame-Supported Membrane Structures for Maximum Clamping Force

Application note Fast focusing for machine vision applications With Optotune s focus tunable lens EL C

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =

Coupling of surface roughness to the performance of computer-generated holograms

Transcription:

Measuring and Evaluation of the Shape and Size of Initial Geometrical Imperfections of Structural Glass Members ONDŘEJ PEŠEK, JINDŘICH MELCHER Institute of Metal and Timber Structures Brno University of Technology Veveří 331/95, 62 Brno CZECH REPUBLIC pesek.o@fce.vutbr.cz; melcher.j@fce.vutbr.cz http://www.kdk.fce.vutbr.cz Abstract: - Buckling resistance (including flexural buckling and lateral-torsional buckling) is major design parameter for many structural glass components. Buckling resistance is influenced by several parameters: boundary conditions, slenderness and initial imperfections. Initial imperfections, in particular initial geometrical out-of plane curvatures are very poorly presented in literature. In this paper measuring and determining of initial geometrical imperfections of fifteen specimens are presented. Measured specimens had the same dimensions but they had different composition there are made from monolithic glass or laminated glass using PVB or EVA foil bonding together annealed glass panes or fully tempered glass panes. The shape of initial overall bow was measured using CARL ZEISS testing equipment. Information about initial curvature shapes and amplitudes are major conclusions resulting from statistical evaluation of tests. Key-Words: - Structural glass, initial geometrical imperfections, flexural buckling, normal distribution, lognormal distribution 1 Introduction Glass has been established as a material of load carrying members and structures in the end of twentieth century and its importance still grows today [9]. Due to slender of glass members it is necessary to design them on stability problems flexural buckling or lateral torsional buckling. Design methods of steel and timber structures are not completely usable for glass structures because of several differences (initial imperfections, brittle behaviour and laminated glass behaviour [1]. The shape and size of initial geometrical imperfections are poorly published in recent publications. f ( w ) x = w N π x cos N N L cr Behaviour of imperfect glass columns was published in [1]. Equation (1) describe (according to the second order theory) dependency of deformation f(w ) of axially loaded (by centric force N) imperfect (initial geometrical imperfection w ) column on own size of initial imperfection considering the sinusoidal shape. Flexural deformation f(w ) increases with increasing amplitude of initial imperfection w. N is Euler's critical force, L is (1) critical length of column and x is ordinate of point of interest. Fig. 1.Flexural buckling of slender beam, loaddeflection curve of perfect and imperfect member 2 Methodology Three types of initial imperfections are distinguished according to their origin: geometrical imperfections (geometrical curvature of column), construction imperfections (actual point of load application etc.) and physical imperfections (residual stresses, inhomogeneity of the material). Geometrical imperfections are relatively easily measurable and observable before placing the member in the structure. E-ISSN: 2224-3429 253 Volume 1, 215

2.1 Geometrical Imperfections of Glass Members and Normative Documents The differences of nominal dimensions are generated during manufacturing processes and they should be lower than values given by product standards. Following tolerances should be checked: (i) glass pane thickness, (ii) glass pane length, wide and rectangularity, (iii) edge deformations due to vertical production (does not apply on float glass) and (iv) planarity (flatness). Geometrical deformations (curvature) are the results of glass tempering (manufacturing process of fully tempered glass or heat strengthened glass). The size of deformations depends on type of glass (coated glass, patterned glass etc.), on glass dimensions and aspect ratio, on nominal thickness and on type of tempering process (vertical or horizontal). horizontally tempered glass (another types) and (iii) 5. mm/m for vertically tempered glass (all types). 2.2 Specimens Initial geometrical imperfections (global bow) were measured on specimens listed in Table 2. Specimen Designation Description Glass Foil Pcs T1 T3 ESG 12 Safety glass FTG - 3 T4 T6 VG 66.2 T7 T12 VSG 66.3 T13 T15 VSG 444.33 Laminated glass Laminated safety glass Laminated safety glass Table 2 ANG PVB 3 FTG EVASAFE 6 FTG EVASAFE 3 Every specimens made of soda lime silica glass has the same dimensions length was 15 mm and width was 15 mm. Specimens was from monolithic glass or double laminated glass or triple laminated glass (made from annealed glass or fully tempered glass) and sum of the glass panes thicknesses was 12 mm for all specimens. Fig. 2 Glass pane deformations types Four types of deformations are follows: (a) overall bow, (b) roller wave (only for horizontally tempered glass), (c) curvature of edges (only for horizontally tempered glass) and (d) local bow (only for horizontally tempered glass) see Fig. 2. 2.3 Imperfections Measurement Method Overall bow shape of glass specimens was analyzed used mechanical measuring system Carl Zeiss. Measured specimen was fixed in vertical position (due to dead load deflections elimination) and supported by two timber blocks (in according to EN 1863-1 and EN 1215-1) and stability was ensured by steel stand. Glass type Float glass without coating Position at tempering Fully tempered glass and heat strengthened glass Global bow mm/m Local bow mm/3 mm Horizontal 3..3 Others Horizontal 4..5 All Vertical 5. 1. Table 1 The maximum allowable values of the overall bow deflection according to the EN 1863-1 and EN 1215-1 are: (i) 3. mm/m for float glass horizontally tempered, (ii) 4. mm/m for 2.4 Measurement Device Description Measuring device Carl Zeiss was composed from height-adjustable guiding rail situated on two supporting threaded columns, which was independently height-adjustable by two locking screws. Wagon with mechanical deflection sensor (Carl Zeiss 3 19 with accuracy.1 mm) was moving on guiding rail (x axis). Wagon was able to move with mechanical indicator in vertical (z axis) and lateral (y axis). Measuring set up was tuned so that the tip of indicator watches was situated 1 mm under measured specimen edge. Measuring set up is plotted in Fig. 3. E-ISSN: 2224-3429 254 Volume 1, 215

Fig. 4 Evaluation of guiding rail initial imperfections Fig. 3 Measurement set-up Because of elimination of measurement errors due to measuring set up geometrical imperfections (guiding rail) every specimen and both its edges were measured in two locations (positive location and negative location) and in two positions (position 1 and position 2. Lateral horizontal deflection (y coordinate) was measured in steps of 5 mm. Both edges were measured two times for both locations and positions of specimen. In the first step of evaluation the measured coordinates was transferred to new coordinate system, so that coordinate y = at both ends of specimen. In the second step the average coordinates was made from two measuring to received coordinates of deflections (imperfections) u (x). In the last step the evaluated deflections was corrected by the value of guiding rail imperfections Using this principle, the shape imperfections of the guiding rail u,rail (x) was determined with (2). Repeating this process several times, the imperfections of the guiding rail could be reproduced with high accuracy. The resulting corrected geometrical imperfection u (x) was obtained using (3). u u, rail ( x) ( x) u ( x) u, uncorr,1, uncorr,2 = (2) 2 = u x u x (3) ( x) ( ) ( ), uncorr ±, rail Shape imperfections of the guiding rail are plotted in Fig. 5. To evaluation of initial imperfections of all specimens the mean value of guiding rail imperfections was taken into account. The shape of initial imperfections from all measuring is similar, but on the other hand the amplitudes are relatively different. Mean value in the point of maximum deflection (close to the midspan) was.46 mm with accuracy ±.11 mm. 2.5 Measuring Device Imperfection Measuring set up itself was subjected to implicit errors precision of mechanical deflection device, deviations of the theoretical parallel placement of specimen and test set up and geometrical imperfections of the guiding rail. The first two errors are not significant and were not taken into account. Errors due to shape imperfections of the guiding rail were significant and they were taken into account by compensating method described in Fig. 3. The geometrical imperfection of the guiding rail was deducted from measuring the initial shape imperfections of the same glass specimen twice: once in the conventional position (u,uncorr,1 ) and once in the mirrored position (u,uncorr,2 ). Fig. 5 Guiding rail imperfections 3 RESULTS AND DISCUSSION Evaluation of the measurement was carried out according to the approach first presented by Belis et al in [2]. The shapes of initial imperfections of all specimens are plotted in graphs according to the glass composition (see Tab. 1) in Appendix. E-ISSN: 2224-3429 255 Volume 1, 215

WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS Imperfections of all specimens are evaluated with positive sign. Fig. 7 Typical and untypical specimen: Shapes of imperfection and their approximation Fig. 6 Amplitudes of initial imperfections of tested specimens To quantitatively describe the approximation quality of the initial imperfection shape, a fitting error (fepar for the parabola and fesin for the sinusoid) is defined as the ratio of the maximal error ( par for the parabola and sin for the sinusoid) and the maximal initial imperfection u: The value of initial imperfection u/l is representing maximum relative amplitude, which is not generally situated in the half of length of specimen. Imperfections are plotted in Fig. 6. 3.1 Overall Bow Shape From imperfection shapes plotted in Fig. 7 is clear that initial imperfection shape should be approximated by second order parabola (4) or half sinusoidal wave (5): a par (x ) = u u (,5 L) 2 π x a sin ( x ) = u cos L x2 fe par = par u 1 (6) fe sin = sin u 1 (7) (4) (5) Overall bow shape of typical (most of the samples) and untypical (only unique case) specimen is plotted in Fig. 7. In the graph there approximation curves (sinusoid and parabola) are plotted. On secondary vertical axis there is plotted error of approximations. Fig. 8 Fitting errors and mean values; approximation by sinusoid and parabola E-ISSN: 2224-3429 256 Volume 1, 215

The error of approximation at position x is defined by (8) and (9). ( x) = u par ( x) a ( x) ( x) = u sin ( x) a ( x) par (8) sin (9) The values of fitting errors fe par and fe sin are plotted in graphs in Fig. 8. The mean value of the positive fitting error is higher for a sinusoidal approximation (+15.8 %) than for parabolic approximation (+11.8 %). In the case of negative fitting error it is oppositely: parabola (-23.4 %) and sinusoid (-2. %). Adding and subtracting of absolute values of fitting errors for parabolic approximation (11.8 + 23.4 = 35.2 % or 23.4 11.8 = 11.6 %) and sinusoid approximation (15.8 + 2. = 35.8 % and 2. 15.8 = 4.2 %) leads to the conclusion, that the sinusoidal approximation is more stable. Both functions can be used for implementing initial imperfections in numerical models. However, preference is given to the sinusoidal function because it is corresponding to the eigen-mode, which is often adopted as initial geometrical imperfection in buckling analysis. 3.2 Amplitude of Imperfection Influence of glass type is illustrated on Fig. 9 upper. Initial imperfection of annealed glass (mean.242 mm/m) is significantly less than for fully tempered glass (mean.599 mm/m).influence of lamination is plotted in Fig. 9 bottom. From comparison of mean values of imperfections for monolithic glass (.468 mm/m) and laminated glass (.543 mm/m) it should be noted that this influence is negligible. Fig. 9 Influences on imperfections; glass type effect (upper); lamination effect (bottom) 3.3 Characteristic value of initial imperfection For practice design it is necessary to know the characteristic value of initial imperfection (u /L) k, which entering into buckling analysis. Generally, characteristic value is considered as 5% quantile. 3.3.1 Normal Distribution The results of initial imperfections measuring ware evaluated such that all imperfection values ware positive. Actually the curvature might be convex or concave and mean value of large population is theoretically equal to zero. Available data set of 15 absolute values is assumed as truncated normal distribution with truncation at zero, in agreement with the Probabilistic Model Code [3] by JCSS (The Joint Committee on Structural Safety). This truncated normal distribution is equivalent to a normal distribution which corresponds to a modified population. More specifically, the available data set was considered twice, whereas one time it was given a negative sign and the other it kept a positive sign. Consequently, a new population of 3 imperfections values was obtained with a mean equal to zero and with normal distribution, Fig. 1. For new population the characteristic value was calculated. In general, 5% quantile is calculated by (1), where µ is mean value and σ is standard deviation. 5% quantile = µ + 1. 645 σ (1) Because of a doubled population (to get correct data set), 2.5% quantile on each tail of the normal distribution was calculated using (11): 5% quantile = 2.5% quantile truncated dist. untruncate = d distr. = µ + 1.96 σ (11) E-ISSN: 2224-3429 257 Volume 1, 215

The characteristic value of initial imperfection based on this method is (u /L) k = 1.276 mm/m. Appendix: Shapes of glass specimens sorted according to the glass composition are plotted in graphs. Fig. 1 Normal distribution of doubled population 3.3.2 Log-normal Distribution An alternative approach is to analyze directly the asymmetric probability density function based on original 15 imperfection measurement values only. Histogram of original data set with log-normal distribution function is illustrated in Fig 11. Characteristic value - 5% quantile calculated by STATISTICA software [4] is (u /L) k = 1.186 mm/m. Fig. 11 Log-normal distribution of original population 4 Conclusions Characteristic value of initial imperfection 1.276 mm/m is corresponding to curvature L/8 (for both normal and log-normal distribution). This value is less than allowable values according to the product standards: 3 mm/m for horizontally tempered float glass. From study performed by Belis et al [2] resulted characteristic value (u /L) k = 2.5 mm/m (L/4). E-ISSN: 2224-3429 258 Volume 1, 215

Acknowledgment: This paper has been supported by the Czech Ministry of Education, Youth and Sports within specific research FAST-J-14-2374. Engineering Series. Vol. x, No. 5, 212, pp. 234-239. ISSN: 2227-4588. References: [1] O. Pešek, J. Melcher: Numerical analysis of behaviour of compression members made of laminated structural glass. In Transactions of the VSB - Technical University of Ostrava. Construction Series. Volume 23, Issue 2, Pages 117 126, ISSN (Online) 184-4824, ISSN (Print) 1213-1962, DOI: 1.2478/tvsb-213-18. [2] J. Belis, D. Mocibob, A. Luible, M. Vandebroek: On the size and shape of initial out-of-plane curvatures in structural glass components. Elsevier Ltd., Construction and Building Materials, Volume 25, Issue 5/211, str. 27-2712. ISSN: 95-618. [3] JCSS probabilistic model code. The joint committee on structural safety [online]. 21 [cit. 214-12-15]. Available from: https://www.jcss.byg.dtu.dk. [4] StatSoft Inc. STATISTICA Cz 12, [software]. Available from: https://www.statsoft.com. [5] EN 1863-1. Glass in building Heat strengthened soda lime silicate glass Part 1: Definition and description. European Committee for Standardization, Brussel. ref. number EN 1863-1:211 E. [6] EN 1215-1. Glass in building Thermally toughened soda lime silicate glass Part 1: Definition and description. European Committee for Standardization, Brussel. ref. number EN 1215-1:2 E. [7] EN 572-1. Glass in building Basic soda lime silicate glass products Part 1: Definitions and general physical and mechanical properties. European Committee for Standardization, Brussel. ref. number EN 572-1:212 E. [8] EN 572-2. Glass in building Basic soda lime silicate glass products Part 2: Float glass. European Committee for Standardization, Brussel. ref. number EN 572-2:212 E. [9] O. Pešek, J. Melcher: Influence of Shear Forces on Deformation of Structural Glass Beams. Mathematics and Computers in Science and Engineering Series. Vol. x, No. 5, 212, pp. 24-245. ISSN: 2227-4588. [1] O. Pešek, J. Melcher, M. Pilgr: Numerical Analysis of Hybrid Steel-Glass Beam. Mathematics and Computers in Science and E-ISSN: 2224-3429 259 Volume 1, 215