seismic analysis using robot millennium (m nedelcu)
TRANSCRIPT
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
1/10
33
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Seismic analysis using Robot Millennium
Mihai NedelcuDepartment of Structural Mechanics , Technical University of Cluj-Napoca, 400020, Romania
Summary
The paper presents the modelling and structural analysis of a complex structureunder seismic loads using the Robot Millennium software package. In order to testthe output results, the author chose a building already designed in thedocumentation [1].
The above mentioned documentation was created as a guideline for applying thedesign standard P100-1/2006. Following the given structure architecture andloads, the author of this paper remodelled the structure using Robot Millenniumv.20.1.
The objective of the paper is to show:1. the Robot Millennium instruments used for modelling this type of
structures
2. various modelling ways for the same structure highlighting the mostconvenient one (in the authors opinion)
3. the results accuracy of the analysis in comparison with the output given bythe work of Professor Tudor Postelnicu
KEYWORDS: seismic analysis, modal analysis, stories, rigid links, panel cut,reinforced concrete walls.
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
http://www.utcluj.ro/english/civil_engineering/catedra_mecanica_constructiilor.phphttp://www.utcluj.ro/english/civil_engineering/catedra_mecanica_constructiilor.php -
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
2/10
Mihai Nedelcu
34
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
1. INTRODUCTION
The chosen building is fully described in the 3rd
example of a very professional anduseful documentation: P100-1/ Building seismic design, volume 2-B. Commentsand calculus examples, responsible author: Tudor Postelnicu, Ph.D., Prof.,U.T.C.B. For the seismic analysis, the authors of the above mentioneddocumentation used the lateral force method and 3D linear-elastic computationgenerated by means of the ETABS program. The author of this paper wanted to seeif the Robot Millennium software package could be a good alternative in this kindof structures analysis.
The analyzed building is located in Bucharest and has:3 underground levels (h=3m) + ground floor (h=6m) and 10 stories (h=3m)5 longitudinal spans x 8m and 5 transversal spans 2x7+1x4+2x7m
The structural characteristics:R.C. walls (both uncoupled and coupled by spandrel walls), columns and
beamsconcrete class C25/30steel PC52
ag=0.24g, Tc=1.6sec, ductility class H, importance coefficient g1=1.2
A current floor plan is shown in Figure 1.
2. MODELLING PRESENTATION
For the seismic analysis, the structure is considered fixed at the ground floor base.
Using Robot Millennium v.20.1, the author modeled the structure using bars forcolumns/beams and panels for slabs/ R.C. walls (see Figure 2). Of course after themesh generation (even with large finite elements) the great number of equations (oforder 10^5) led to a significant slowing down of the analysis. As a good alternativethe Rigid Links additional attribute can be used instead of panels for the slabs
modeling. By introducing the Membrane rigid link, the user can connect the nodesof each floor according to any DOF, in this case the X,Y displacements and the RZrotation. And so the slab effect as a rigid body is fully covered. At the same timethe beams have to be modeled as T-section, taking into account the correspondingslab rigidity. A slab width of 3xhp (the slab thickness) was taken on each side for
the interior beams and of 2xhp on each side for the marginal beams. Having no slabfinite elements reduces considerably the analysis time with absolutely no damage
to the results. The model is presented in Figure 3.
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
3/10
Seismic analysis using Robot Millennium
35
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Figure 1. The building current floor plan
Figure 2. The structure model showing: a) section shapes, b) rigid links
Following the instructions of P100-2006, the rigidity of the elements is takendifferently, depending on the type of analysis. First, the goal is to make the modal
analysis, in which the rigidity of walls, columns and beams is taken EI=0.5EcIc. Asfor the spandrel walls, their Young modulus is Espandrel =0.4* Ewalls and the other
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
4/10
Mihai Nedelcu
36
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
sectional characteristics: A=0.2Ac, I=0.2Ic. Reducing the rigidity of the elements,the elements are considered in a plastic stage and so the building behavior duringan earthquake is well approximated according to P100-2006. To achieve that theuser can define a new material with a modified Young modulus, in this case named
C25/30_0.5 (Job Preferences/Materials/Modification). The characteristics of alongitudinal wall are set as presented in Figure 3. As for the spandrel walls the
double reduction is achieved by multiplying the thickness of 50cm with the 0.2factor.
Figure 3. Defining the section properties of: a) longitudinal walls, b) spandrel walls
As for the bar elements, there is an alternative. The program allows the reductionof moment of inertia according to local axes x, y, z (see Figure 4).
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
5/10
Seismic analysis using Robot Millennium
37
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Figure 4. Defining the properties of: a) columns 70x70cm, b) beams 30x70cm
Another problem in the modal analysis is the modeling of the additionaleccentricities: 5% in each direction. The code explains in detail how their effecthas to be taken into account. The level seismic forces, computed by hand or by
some user defined computer program have to be applied to the resistance elementsdepending on their rigidity. Because there are 2 main directions, for each one 2signs and the mass eccentricities can also be positive or negative, there are 8combinations to be made. This work is time consuming and the probability of usererrors is highly increased. A good alternative seems to be the Modal AnalysisParameters/Definition of mass eccentricities that the program offers (see Figure 5).The sign of eccentricity can be changed but is impossible to have both signs in thesame model. So, different models have to be made. Of course, if the structure issomeway symmetrical, the work is much reduced.
Figure 5. Introducing the additional mass eccentricities
In the 3rd
example the loads are generally given. Because of this generality theresults of this paper are not identical to the original ones, but the differences areacceptable. Given the loads, the load to mass conversion is made automatically.The size of FE-s should be set as big as possible, since from the modal analysis
point of view the mesh refinement will not change the results. Doing that the modalanalysis was made in a few seconds. The results are compared to the original ones(see Figures 6 and 7). Also the shape of the eigenvectors for the first 3 modes isshown in XY view (see Figure 8). The additional eccentricities effect for the firsttwo modes can be seen.
Figure 6. The Robot Millennium modal analysis results
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
6/10
Mihai Nedelcu
38
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Figure 7. The original modal analysis results
Figure 8. The shape of the eigenvectors for the first 3 modes in XY view
The next step is to define the seismic analysis (see Figure 9). The behavior factor(q) has different values according direction X and Y. For the ductility class H, on
X-longitudinal direction, the walls are considered as cantilevers ( )
while on Y-transversal direction, the walls are coupled, and so results a higher q
( ).
1/4 uq =
1/4 uq =
Figure 9. Introducing the seismic analysis parameters
Since the seismic analysis is concerned, only the special load combinations areintroduced. To assess the deformations, no changes are made to the FE mesh or to
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
7/10
Seismic analysis using Robot Millennium
39
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
the sectional rigidities. The allowed displacements are given by P100-2006 inrelation to the elastic drift.
The program offers another useful tool by allowing the user to define the stories ofthe building: Geometry/Code Parameters/Stories. The elements of each storey are
presented in a different color in the Figure 10. The walls are not colored because ofthe large FE mesh.
Figure 10. The stories display
The characteristics of each storey can be seen by the Story Table/Values option. Inthe Figure 11, are presented the mass, the gravity and torsion centers, eccentricitye0 and the applied additional eccentricities e2.
Figure 11. Story Table: Values
Using Story Table/Displacements option, the stories maximum absolute andrelative displacements (dre) can be seen (Figure 12).
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
8/10
Mihai Nedelcu
40
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Figure 11. Story Table: Values
And so the elastic and inelastic drift can be computed. Table 1 is showing acomparison between the results of this paper and the original ones for X direction.
Table 1. Name of the table
Robot Millennium Example 3
Level dreX driftSLS driftULS driftSLS driftULS(cm)
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
9/10
Seismic analysis using Robot Millennium
41
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Figure 12. Sign convention
Figure 13. Choosing a) the cut position, b) the reduced forces to be displayed
There are two possibilities for modeling the wall columns. Until now they weremodeled as bars. In this case, for a RC wall the internal forces of both wall andwall columns must be combined in order to have the result for the entire wall
section. A second possibility is to eliminate the wall columns and to use anequivalent wall section from the primary moment of inertia point of view. Thebuilding stiffness will slightly decrease due to the lack of compression/tensionabsorbed by the wall columns. Comparison of the 2 modeling methods with theoriginal results is given for a longitudinal wall (Table 2). The forces are calculatedat the ground floor level from one seismic combination on X direction. The
difference between the axial forces has to come from the lack of knowledge in theloads distribution over each storey since the evaluation of the building total weightwas checked and almost perfect similarity was achieved.
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis
-
7/27/2019 Seismic Analysis Using Robot Millennium (M Nedelcu)
10/10
Mihai Nedelcu
42
Seismic Analysis
SSN 1582-3024
http://www.intersections.ro
Table 2. Final internal forces for a longitudinal wall
Wall PL1 N (kN) M (kNm) T (kN)
Wall -7741.22 -32648.3 2986.15
Wall column 1 -2760.13 -178.81 324
Wall column 2 5175.29 -251.61 389Final forces -10156.4 64389.96 3699.15
Wall without columns -9939.21 -63955.5 3468.02
Example 3 -11456 65463 3778
3. CONCLUSIONS
The Robot Millennium software package gives us a reliable and elegant alternativefor seismic analysis. Its features fully cover the demands of P100-2006. However,
careful examination of its tools has to be done in order to choose the optimal
modeling strategy.
References
1. Postelnicu Tudor, P100-1/Proiectarea Seismica a cladirilor vol.2-B. Comentarii si exemple decalcul ctr. U.T.C.B. nr. 158/02.08.2005 M.T.C.T.
2. Robot Millennium v.20.1.Autodesk Inc, 20083. CR 2-1-1.1, Cod de proiectare a constructiilor cu pereti structurali de beton armat4. P100-1/2006, Cod de proiectare seismica
Article No.3, Intersections/Intersecii, Vol.5, 2008, No.3, Seismic Analysis