Numerical tools for isogeometric analysis in the nonlinear regime

The present work deals with the development of robust numerical tools for Isogeometric Analysis suitable for problems of solid mechanics in the nonlinear regime. To that end, a new solid-shell element, based on the Assumed Natural Strain method, is proposed for the analysis of thin shell-like structures. The formulation is extensively validated using a set of well-known benchmark problems available in the literature, in both linear and nonlinear (geometric and material) regimes. It is also proposed an alternative formulation which is focused on the alleviation of the volumetric locking pathology in linear elastic problems. In addition, an introductory study in the field of contact mechanics, in the context of Isogeometric Analysis, is also presented, with special focus on the implementation of a the Point-to-Segment algorithm. All the methodologies presented in the current work were implemented in a in-house code, together with several pre- and post-processing tools. In addition, user subroutines for the commercial software Abaqus were also implemented.

A finite-strain solid–shell using local Löwdin frames and least-squares strains

A finite-strain solid–shell element is proposed. It is based on least-squares in-plane assumed strains, assumed natural transverse shear and normal strains. The singular value decomposition (SVD) is used to define local (integration-point) orthogonal frames-of-reference solely from the Jacobian matrix. The complete finite-strain formulation is derived and tested. Assumed strains obtained from least-squares fitting are an alternative to the enhanced-assumed-strain (EAS) formulations and, in contrast with these, the result is an element satisfying the Patch test. There are no additional degrees-of-freedom, as it is the case with the enhanced-assumed-strain case, even by means of static condensation. Least-squares fitting produces invariant finite strain elements which are shear-locking free and amenable to be incorporated in large-scale codes. With that goal, we use automatically generated code produced by AceGen and Mathematica. All benchmarks show excellent results, similar to the best available shell and hybrid solid elements with significantly lower computational cost.

CROSS-SECTIONAL ANALYSIS OF PRE-TWISTED THICK BEAMS USING VARIATIONAL ASYMPTOTIC METHOD

An asymptotically-exact methodology is presented for obtaining the cross-sectional stiffness matrix of a pre-twisted moderately-thick beam having rectangular cross sections and made of transversely isotropic materials. The anisotropic beam is modeled from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy of the beam is computed making use of the constitutive law and the kinematical relations derived with the inclusion of geometrical nonlinearities and initial twist. Large displacements and rotations are allowed, but small strain is assumed. The Variational Asymptotic Method is used to minimize the energy functional, thereby reducing the cross section to a point on the reference line with appropriate properties, yielding a 1-D constitutive law. In this method as applied herein, the 2-D cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged as orders of the small parameters. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that renders the 1-D strain measures well-defined. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

Analysis of Transverse Shear Strains in Pre-Twisted Thick Beams using Variational Asymptotic Method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam reference curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

A finite-strain solid–shell using local Löwdin frames and least-squares strains

A finite-strain solid–shell element is proposed. It is based on least-squares in-plane assumed strains, assumed natural transverse shear and normal strains. The singular value decomposition (SVD) is used to define local (integration-point) orthogonal frames-of- reference solely from the Jacobian matrix. The complete finite-strain formulation is derived and tested. Assumed strains obtained from least-squares fitting are an alternative to the enhanced-assumed-strain (EAS) formulations and, in contrast with these, the result is an element satisfying the Patch test. There are no additional degrees-of-freedom, as it is the case with the enhanced- assumed-strain case, even by means of static condensation. Least-squares fitting produces invariant finite strain elements which are shear-locking free and amenable to be incorporated in large-scale codes. With that goal, we use automatically generated code produced by AceGen and Mathematica. All benchmarks show excellent results, similar to the best available shell and hybrid solid elements with significantly lower computational cost.

Citizens and social media in times of natural disaster: Exploring information experience

Introduction. Social media is becoming a vital source of information in disaster or emergency situations. While a growing number of studies have explored the use of social media in natural disasters by emergency staff, military personnel, medial and other professionals, very few studies have investigated the use of social media by members of the public. The purpose of this paper is to explore citizens’ information experiences in social media during times of natural disaster. Method. A qualitative research approach was applied. Data was collected via in-depth interviews. Twenty-five people who used social media during a natural disaster in Australia participated in the study. Analysis. Audio recordings of interviews and interview transcripts provided the empirical material for data analysis. Data was analysed using structural and focussed coding methods. Results. Eight key themes depicting various aspects of participants’ information experience during a natural disaster were uncovered by the study: connected; wellbeing; coping; help; brokerage; journalism; supplementary and characteristics. Conclusion. This study contributes insights into social media’s potential for developing community disaster resilience and promotes discussion about the value of civic participation in social media when such circumstances occur. These findings also contribute to our understanding of information experiences as a new informational research object.

Some uses of the spontaneous polarization method

As a partial fulfillment of the requirements in obtaining a Professional Degree in Geophysical Engineering at the California Institute of Technology. Spontaneous Polarization method of electrical exploration was chosen as the subject of this thesis. It is also known as "self potential electrical prospecting" and "natural currents method."

The object of this thesis is to present a spontaneous polarization exploration work done by the writer, and to apply analytical interpretation methods to these field results.

The writer was confronted with the difficulty of finding the necessary information in a complete paper about this method. The available papers are all too short and repeat the usual information, giving the same examples. The decision was made to write a comprehensive paper first, including the writer's experience, and then to present the main object of the thesis.

The following paper comprises three major parts:

1 - A comprehensive treatment of the spontaneous polarization method.

2 - Report of the field work.

3 - Analytical interpretation of the field work results.

The main reason in choosing this subject is that this method is the most reliable, easiest and requires the least equipment in prospecting for sulphide orebodies on unexplored, rough terrains.

The intention of the writer in compiling the theoretical and analytical information has been mainly to prepare a reference paper about this method.

The writer wishes to express his appreciation to Dr. G. W. Potapenko, Associate Professor of Physics at California Institute of Technology, for his generous help.

Effective Length K-Factors For Flexural Buckling Strengths Of Web Members In Open Web Steel Joists

Open web steel joists are designed in the United States following the governing specification published by the Steel Joist Institute. For compression members in joists, this specification employs an effective length factor, or K-factor, in confirming their adequacy. In most cases, these K-factors have been conservatively assumed equal to 1.0 for compression web members, regardless of the fact that intuition and limited experimental work indicate that smaller values could be justified. Given that smaller K-factors could result in more economical designs without a loss in safety, the research presented in this thesis aims to suggest procedures for obtaining more rational values. Three different methods for computing in-plane and out-of-plane K-factors are investigated, including (1) a hand calculation method based on the use of alignment charts, (2) computational critical load (eigenvalue) analyses using uniformly distributed loads, and (3) computational analyses using a compressive strain approach. The latter method is novel and allows for computing the individual buckling load of a specific member within a system, such as a joist. Four different joist configurations are investigated, including an 18K3, 28K10, and two variations of a 32LH06. Based on these methods and the very limited number of joists studied, it appears promising that in-plane and out-of-plane K-factors of 0.75 and 0.85, respectively, could be used in computing the flexural buckling strength of web members in routine steel joist design. Recommendations for future work, which include systematically investigating a wider range of joist configurations and connection restraint, are provided.

固有频率法评估损伤的阈值研究

通过分析损伤可测性的影响因素,对固有频率法判断损伤阈值的估算方法进行了研究,并以悬臂梁模型为例,应用该方法研究了悬臂梁的损伤阈值.数值及理论计算的结果表明:在已知频率测量精度及方法偏差时,该方法可以判别同有频率法损伤识别的阈值,在早期损伤阶段,该方法判别的阈值与数值模拟结果吻合良好.

Increasing allowable deformation criteria through application of level II LTSM approach

In the central part of the Delft railway tunnel project, an underground railway station is being built at very close distance to the existing station building, which is still in operation. Although elaborate sensitivity analyses were made, some unforeseen deformations were encountered during the first phases of the execution process. Especially the installation of temporary sheet pile walls as well as the installation of a huge amount of grout anchor piles resulted in deformations exceeding the predicted final deformations as well as the boundary values defined by a level I limiting tensile strain method (LTSM) approach. In order to ensure the execution process, supplementary analyses were made to predict future deformations, and this for multiple cross sections. These deformations were implemented into a finite element model of the masonry of the building in order to define probable crack formation. This Level II LTSM approach made it possible to increase the initially foreseen deformation criteria and the continuation of the works. Design steps, design models and monitoring results will be explained within this paper.

Assessment of the settlement vulnerability of masonry buildings

Excavation works in urban areas require a preliminary risk damage assessment. In historical cities, the prediction of building response to settlements is necessary to reduce the risk of damage of the architectural heritage. The current method used to predict the building damage due to ground deformations is the Limiting Tensile Strain Method (LTSM). In this approach the building is modelled as an elastic beam subjected to imposed Greenfield settlements and the induced tensile strains are compared with a limit value for the material. These assumptions can lead to a non realistic evaluation of the damage. In this paper, the possibility to apply a settlement risk assessment derived from the seismic vulnerability approach is considered. The parameters that influence the structural response to settlements can be defined through numerical analyses which take into account the nonlinear behaviour of masonry and the soil-structure interaction. The effects of factors like material quality, geometry of the structure, amount of openings, type of foundation or the actual state of preservation can be included in a global vulnerability index, which should indicate the building susceptibility to damage by differential settlements of a given magnitude. Vulnerability curves will represent the expected damage of each vulnerability class of building as a function of the settlement.

Computational modelling of building damage due to tunnelling

Underground constructions in soft ground may lead to settlement damage to existing buildings. In The Netherlands the situation is particularly complex, because of the combination of soft soil, fragile pile foundations and brittle, unreinforced masonry façades. The tunnelling design process in urban areas requires a reliable risk damage assessment. In the engineering practice the current preliminary damage assessment is based on the limiting tensile strain method (LTSM). Essentially this is an uncoupled analysis, in which the building is modelled as an elastic beam subject to imposed Greenfield settlements and the induced tensile strains are compared with a limit value for the material. The soil-structure interaction is included only as a ratio between the soil and the building stiffness. In this paper, a coupled approach is evaluated. The soil-structure interaction in terms of normal and shear behaviour is represented by interface elements and a cracking model for masonry is included. This project aims to improve the existing damage classification system for masonry buildings subjected to tunnel-induced settlement, in order to evaluate the necessity of strengthening techniques or mitigation measures.

Coupled analysis of building damage due to tunnelling

Excavation works in urban areas require a preliminary risk damage assessment. In historical cities, the prediction of building response to settlements is necessary to reduce the risk of damage of the architectural heritage. The current method used to predict the building damage due to ground deformations is the Limiting Tensile Strain Method (LTSM). This method is based on an uncoupled soil-structure analysis, in which the building is modelled as an elastic beam subject to imposed greenfield settlements and the induced tensile strains are compared with a limit value for the material. This approach neglects many factors which play an important rule in the response of the structure to tunneling induced settlements. In this paper, the possibility to apply a settlement risk assessment derived from the seismic vulnerability approach is considered. The parameters that influence the structural response to settlements can be defined through numerical coupled analyses which take into account the nonlinear behaviour of masonry and the soil-structure interaction.