Numerical modelling of the seismic behavior of timberframed structures based on macroelements
Historical timberframed structures can be found all over the world since ancient times, being able to observe a variety of procedures and construction systems which is adapted to each zone according to the materials and available knowledge. Recent earthquakes have highlighted the good antiseismic behavior of this structural typology. Likewise, its antiquity has also evidenced a good structural capacity.
The case of Pombalino Buildings are born from the need to build an antiseismic constructive system after the earthquake of Lisbon 1755 which destroyed most part of the city. The present thesis of the autor Sonia Guerra Pinto seeks to complement the study of these buildings due to the uncertainty related to their seismic behavior since up to date there has not yet been an earthquake registered in Lisbon since 1755.
TimberFramed Structures through history
In Central and South America timber has been used since prehistoric times, there is data from the year 16.500 B.C. in Monte Verde, Chile. European archaeological excavations in the buried settlements after the eruption of Vesuvius in 79 A.C. brought into light timber frame walls in many noble buildings mainly of Herculaneum and secondarily of Pompeii.
Although there is a large amount of information regarding timber structures throughout history, the present literature review will focus on historic timber framed structures which are dated since the 17th Century in Europe (Spain) and since the 16th Century in South America, specifically in Peru.
This type of framing is assigned different names which is associated to a specific constructive system that varies according to the region in which it has been built, for example in territories from Southern Central Anatolia to the Ottoman Balkans including Black Sea Coasts of Romania, Crimea, Bulgaria, FYR Macedonia and Bosnia Herzegovina to Greece in the West we can find “hımıs” ; In Spain the timberframed construction system was called “entramado o telar”, 'Borbone' Constructive system is found in the Calabrian Region of Italy, “Dahjji Construction” in Pakistan, “Pombalino” Constructions in Lisbon – Portugal and in South America we can find “Quincha” Buildings in the coastal area of Perú, “Adobillo” Structures in the Central Zone of Chile and “Gingerbread Houses” and “Kay Peyi” in Haiti .
Behavior of Timber framed Structures under Seismic Loads
Studies show that timber structures behave very well under seismic loads, the 2010 the Haiti earthquake has revealed that traditional construction techniques have value in resisting to earthquakes.
Ambraseys & Jackson state that after the 1509 Istanbul Earthquake, the Ottoman authorities prohibited masonry and enforced the construction of timberframe houses, claiming that masonry was responsible for most of the casualties produced by the earthquake. By the end of the century the city was almost entirely built in wood. For a set of earthquakes between the beginning of the 20th century and 1980 in Turkey and Greece, it is stated that that 'the number of people killed per 100 houses destroyed by earthquakes of magnitude equal to or greater than 5,0 is only around 1 for timber constructions.
In Peru the Quincha technique reached its peak after the 1687 earthquake, when a law was passed ruling that quincha must be used for the upper storeys of any building greater than a single storey in height.
Today a number of these buildings survive, most dating from the 18th and 19th centuries, with quincha in the upper storeys, and the first storey in adobe or fired brick . In the Calabrian region of Italy after the 1783 earthquake the government ruled by Borbone dynasty enacted what is considered the first European antiseismic code; called 'Sistema costruttivo borbonico' that amongst other measures enforced the presence in the new constructions or repair of damaged buildings, of a timber skeleton as an essential requirement to ensure safety.
In Portugal the Pombalino buildings were introduced after the 1755 catastrophic earthquake as a structural solution that would provide the required seismic resistance. These were built in quarters, each block comprising an average of 10 buildings, it was arguably the first case in history of an entire town built with the purpose of providing seismic resistance to its buildings.
SEE ALSO: How to manage and mitigate disaster risk in urban areas?
MacroModel Calibration
In this thesis, the assumptions and considerations for the elaboration of the macromodel is explained and discussed.
In order to study the nonlinearities of the present model, some considerations had to be done, this has lead to a calibration process where graphs were elaborated to show the results.

OpenSees
Open System for Earthquake Engineering Simulation (OpenSees), was initially developed by F. McKenna and G. L. Fenves with many other contributors at the NSF sponsored Pacific Earthquake Engineering (PEER) centre, is an objectoriented framework for simulating applications in earthquake engineering using finite element analysis . It has the capability to perform many types of analysis including static pushover, static reversedcyclic, dynamic timeseries, and uniform or multisupported excitations for inelastic timehistory analysis for both structural and geotechnical systems.
It is based in fully based on the Tcl/Tk scripting language. It is comprised of a set of four modules that perform the finite element analysis. The model builder that performs the finite element model; the analysis module specifies the analysis procedure; the domain selects quantities to be monitored during the analysis; and the recorder records the output.
The OpenSees platform originally does not contain a graphical user interface and can therefore be challenging for users unfamiliar with the Tcl/Tk scripting language. Although there are several plugins such as 'OpenSeesNavigator' for Matlab and 'OSLite V0.35' that allow us to visualize the current model, also allowing us to visualize it's stresses and modal shapes. GID program can be used in order to generate the coordinates, nodes and elements of the model that shall be developed.

Model Builder
First it is important to establish the number of degrees of freedom and the dimension of the problem (1, 2 or 3) and afterwards consider the next steps:
NonLinear Beam Column Elements
The model builder constructs as in any finite element analysis, the analyst's first step is to establish the nodes and it's coordinates, afterwards the column, beam, diagonal and girders are established as displacement beam column, truss element, linear beam column, plastic hinge, among others. In the case of the elaborated model NonLinear Beam Column Elements were considered which is an object command based on the noniterative or iterative force formulation and considers the spread of plasticity along the element.
The integration along the element is based on GaussLobatto quadrature rule (two integration points at the element ends). The element is prismatic meaning that the beam is represented by the section model identified at each integration point governing the response, in this step we must also assign the masses, in which for the present model will be considered as massnodes, although distributed mass nodes can also be considered throughout a different command.
Shell Element
This command is used to construct a plane 2D mesh object which can be represented as a Quadrilateral Element, Shell Element, Bbar Plane Strain Quadrilateral Element or Enhanced Strain Quadrilateral Element. In the case of the model elaboration Shell Element (ShellMITC4) is used, which uses a bilinear isoparametric formulation in combination with a modified shear interpolation to improve thinplate bending performance, which will be considered in order to represent masonry surfaces. In this process the properties of the section are established as a PlateFiber and Elastic Isotropic, defining the young modulus, poisson's ratio and thickness.
Section Command
The elements of the studied model will be considered as an Elastic section. The general properties of each element such as young modulus, shear modulus, torsional moment of inertia, second moment of area about the local yaxis (Iy), second moment of area about the local zaxis (Iz) and Cross Section Area of the Section will be established.
The central elements of the diagonals will be considered as uniaxial materials, where the nonlinear behavior will be concentrated as SAWS materials properties.
Linear CoOrdinate Transformation
The linear coordinate transformation object command performs a linear geometric transformation of beam stiffness and resisting force from the basic local system to the global coordinate system. In other words, this command will be established as linear and according to the movement of the model, it must be fixed at one axis and free in two axis (depending on the direction of the forces). For example in Figure 29 it is shown that according to the force application, the present model will be fixed in the z direction.
SAWS Uniaxial Material
SAWS provides the implementation of a onedimensional hysteretic model developed as part of the CUREE Caltech wood frame project. The definition and material properties are shown in Figure 30. This force deformation model is characterized by six physically identifiable parameters in developing the unloading and reloading paths. Degrading Stiffness (KP), Intercept strength of shear wall spring element (F0), Intercept strength for spring element pinching branch (FI), Spring element displacement at ult. Strength (DU), Initial stiffness of shear wall spring element (S0 or K0), Stiffness ratio of the asymptotic line (R1), Stiffness ratio of the descending branch (R2), Stiffness ratio of the unloading branch (R3), Stiffness ratio of the pinching branch (R4), Stiffness degradation parameter for the shear wall spring element (alpha), Stiffness degradation parameter for the spring element (beta).
Domain, Recorder and Analysis
The Domain object is responsible for storing the objects created by the ModelBuilder as shown in Figure 31 and for providing the Analysis and Recorder objects access to these objects.
The recorder object monitors userdefined parameters in the model during the analysis. This, for example, could be the displacement history at a node in a transient analysis, or the entire state of the model at each step of the solution procedure. Several Recorder objects are created by the analyst to monitor the analysis.
Each analysis in OpenSees consists of the following commands:
 Dimension of the problem – 1, 2 or 3
 Constraints – handle the constraints defined on the domains
 Numberer – numbers the degreesoffreedom in the domain.
 System – constructs the solving objects to store and solve the system of equations.
 Test – establish the convergence test to ensure the convergence can be achieved at the an end of iteration step.
 Algorithm – iterate from the last step to the current.
 Integrator – determines the next step for an analysis.
 Analysis – defines what type of analysis is to be performed.
Pombalino Case Study – Lots 210 to 220 of Rua (Street) da Prata.
Pombalino construction was standardized to be used for the reconstruction of Lisbon after the 1755 earthquake, but its use was not limited to Lisbon. Mass production of building components, a basic principle of Pombaline architecture in Portugal was also used in Vila Real de Santo Antonio, and in many other Portuguese cities.
The form construction of these buildings was very similar to that of Lisbon, the party and exterior walls being in stone, while inside there were timber framed partition walls with St. Andrew crosses incorporated to them. Arches in bricks were used to tie the foundations together and for some walls.
The plan of Santo Antonio, as in Lisbon, consisted of a rectangle, with one of the long sides facing the river, to the east. The rectangle was cut by 5 streets in a NorthSouth direction and orthogonally in an EastWest direction. All the streets were the same width and they contained 43 blocks; 32 of which were identical in size, being 240 palms by 100.
In Vila Real, there are four quite distinct architectural types: the river front buildings, the buildings in the square, the single storey houses with towers and the single storey houses without towers. The walls of the residential floor had an antiseismic wooden structure incorporated into them, similar to that seen in the Pombalino quarter in Lisbon. The roof structure was very simple and repetitive, covered by wooden boards on which the tiles were laid.
The properties of the macromodel were used as a basis to elaborate a numerical model of a global pombalino building.
The plan that created the downtown Pombalino was defined by a regular scheme of streets and squares. Each square was divided into lots of buildings with different front widths, but maintaining the same depth and height as shown in Figure 57, which also shows the studies building.
The total area of the square was of 2.000 m2, varying the areas of the lots between one hundred and three hundred square meters. The studied lots have a very wide facade (Figure 57).
In terms of the initial plan, the width of seventeen meters corresponds to six modules of facade, each having 2.8 meters. In the Floor Plan (Figure 58) there are two shops and the entrance of the building. The upper floors (2nd to 4th) are composed of a floorplan type shown as shown in Figure 59, in which we can find: on the first floor, a row of rods and a commercial store; In the second, an office; In the third, housing; In the fourth, a house with a balcony; The fifth floor is composed of two attics (Figure 60).
Conclusions
In the studied case the rigidity of the masonry walls consequently increases the major stress concentration on timberframe (nonbraced) elements. Although it is shown that gaiola (mixed timber masonry) structure has an important contribution in absorbing the stresses.
Overall, the studied pombalino building has a good seismic behavior in which the high stresses are concentrated on the top section of the building. This is a good indication because the structure will not fail from the middle floors, although it is important to consider the large stress concentrations on the 3rd floor.
Seismic assessment can be performed from the obtained results in order to reduce the stress in this section (it could be recommended to include some longitudinal gaiola bracing element).
Would you like to know more? Leave your email below to receive the full academic thesis behind this article:
Subscribe to our Newsletter
Receive every month our top articles about rehabilitation and preservation of historical constructions, information about job opportunities and much more.