Optimal Slope of Dramix Type Fibers in Reinforced Concrete

Size: px
Start display at page:

Download "Optimal Slope of Dramix Type Fibers in Reinforced Concrete"

Transcription

1 6 th World Congresses of Structural and Multidisciplinary Optimization Rio de Janeiro, 3 May - 3 June 25, Brazil Optimal Slope of Dramix Type Fibers in Reinforced Concrete P. Prochazka 1, N. Starikov 2 1 Czech Association of Concrete Engineers, Prague, Czech Republic, petrp@fsv.cvut.cz 2 CTU Prague, Civil Engineering, Prague, Czech Republic, star - 13@mail.ru 1. Abstract Several experimental studies proved that the standard reinforcement in concrete (rebars) could be improved by short fibers to avoid local warping and cracking. In order to restrain the reinforcing rebars from corrosion, the fibers are basically at help during the curing process of concrete. Particularly in geomechanics, or foundation engineering, where one-sided moistening is expected, the reason for usage of the fibers is advantageous due to obviously irregular dehydration. The volume ratio of fibers in these cases is very low, say, about.5 percent. Another one interesting application is in construction of defense walls (bunkers), where the fiber ratio is principally higher. In all cases curvilinear fibers are preferred from then viewpoint of bearing capacity of the system. On the other hand, the financial expenses are also higher. This is why optimal shape design is required for maximum effectiveness of fiber reinforced concrete. Such a typical shape of current most produced fibers of this type is prepared by Dramix. This is why we concentrate our attention on straight fibers and also to the curvilinear fibers. The optimization of slope in Dramix type fibers is solved from pullout problem studied on unit cells of different arrangements. The steel fibers are considered and either linear or elasticplastic Mises material of concrete is taken into account. Mixed variational principle has been applied to enable one to formulate contact problem together with optimization of slopes. Transformation field analysis is employed to make easier the influence of nonlinear behavior of concrete. The critical pullout force is solved from finite element method and is identified from destruction of the aggregate. Several typical examples are solved to show the properties of the composite for various arrangements of fibers in cells. 2. Keywords: Shape optimization, fiber reinforced concrete, Dramix type fibers 3. Introduction The fibers play very important role particularly during curing process of concrete, as they suppress local cracking and warping in the composite structure and avoid a possibility of corrosion of reinforcing steel rods. Very important fact follows from plenty of experimental studies: since the steel fibers are stiffer then concrete itself, pullout problem has to be solved when requiring knowledge about the bearing capacity of the structure. In this study 3D isoparametric finite elements are used for solving the problem on unit cells. The unit impulses are applied to nodes of the 24-degree of freedom elements to create influence matrices. This approach is effective from the point of view of iterative method, which has to be applied to the solution of contact problem. Because the optimization of slope of Dramix type fibers is required, discrete problems are successively solved for different unit cells. Four typical unit cells are studied and discrete fiber slopes are solved, namely admissible pullout forces are determined from the set, 1, 3, 45, 6, and 9 degrees of angle between the axial coordinate and the slope. Interpolation of the influence of slopes on the admissible forces is available, as the results provided from these computations admit approximately continuous curve describing the influence of the slope. Natural question occurs when studying the problem of behavior of curvilinear fibers in concrete matrix. The slope of the fibers in the curvilinear part can change and the response of this change will obviously have an impact to the bearing capacity of the composite aggregate. If we leave the vertical distance between horizontal parts of the fibers constant, the slope of the middle part can change and the design parameter of the optimization is the angle of the slope. From experimental tests it has been proved that the extraction of the fibers out of the concrete are caused by penetrating of the fibers into the concrete due to attaining some admissible stress (compressive stress) or slipping of the fibers along the interfacial zones with the concrete. The optimization criteria follow this observation of experimental results. The first criterion will require that the admissible force cannot exceed the criterion of extraction of the fibers due to violation of Mohr-Coulomb hypothesis. This case prevails particularly in the case of straight or almost straight fibers. The second criterion states that near the fibers no stress attains the compressive strength. 4. Pullout problem of a steel FRC beam The pullout problem has frequently been solved in a problem of cracking of composite structures of several sorts. Previously, several numerical studies were carried out using the FEM, [1], and the BEM, [2], and the results were compared with experimental results from available literature with a good agreement. In the sections from the available literature the topics have mainly been focused on experiments and the theory has been based on various approximate hypotheses, i.e., the results from the theoretical considerations are very approximate. A remarkable feature in [1] and [2] of the mathematical solution occurred: a cracking was initiated not only at the face of the fiber-matrix system, but also inside the trial body. As the previous problems were concentrated only on straight fibers and the nature of the material used was based on classical composites with epoxy matrices and a high bearing capacity and durability, a natural question arises: what happens when the matrix is created from concrete and the aggregate behave realistic, i.e., the shape of the fibers is no more straight and the concrete matrix is taken into account. The answer to this question is the objective of this section. Some numerical experiments show on how the theoretical considerations can be applied in reality. The examples mutually differ by positioning the fibers and by volume fractions. Dramix-type fibers are used in studies of curvilinear fibers. In this study of a lag model of the system concrete matrix - steel fibers of a special shape, the geometry of the problem is considered as three-dimensional and high accurate finite isoparametric elements with 24 degrees of freedom are used together with Uzawa's algorithm, [3], for solving a contact problem fiber - matrix. A finer mesh of the finite elements is necessary for the stability of an iterative process particularly in the neighborhood of expected nonlinear behavior on the contact (interface between fibers and

2 matrix), which covers the Uzawa's algorithm. A procedure, which is close to TFA, [2], [4], is used. It consists in creation of influence matrices due to unit impulses of tractions on the interface, which are applied at these points of the finite element mesh, which belong on the contact to either boundary of the fibers or the boundary of the concrete. The debonding or slipping" (jumps of responding points of the matrix and fibers boundaries in the tangential direction to the interface) can be caused due to several physical models. First, no tension tractions on the interface may occur. Also, some friction law (such as Coulomb law, Mohr-Coulomb law, shear band strength are often introduced and debond parameters are computed. Since no fiber gripping force from matrix is expected, the general Mohr-Coulomb friction law is restricted to the exclusion of shear stresses exceeding the shear strength in the tangential direction and the tension exceeding the tension strength is also excluded slang the interface between the fiber and the matrix. In our study the grip in normal direction to the fibers is neglected. The problem is obviously three dimensional, so that the most accurate results may only be obtained from very precise numerical models. The behavior of the system starts partially with observations described in [2], where one sided moistening is ensured during the curing process. A periodic system of fibers is assumed in our model and different volume tractions and positioning of the fibers are considered. In order to obtain results comparable with experiments and on the other side results from real composite structures a wide range of examples were analyzed in [1] by applying the FEM and the results from the study put forward consulted them. Our current aim is to show the results from curvilinear shape of fibers and numerical results obtained in [1] and [2] which were in a good agreement with experiments. As the solution of this problem leads to a sequence of iteration steps, an approach similar to the transformation field analysis (TFA), [4], is used for speeding substantially up the iteration. As in recent papers, [1], [2], Uzawa's algorithm is used to the pullout problem of a lag model without proof of solvability and uniqueness. The main goal of this section is to show the behavior of the pullout problem of the lag model in connection with the numerical procedural proposed. Also, Uzawa's algorithm is, in comparison to the previous papers, [1], [2], carried out in the sense of the TFA, i.e., the influence matrices needed are prepared a priori. Since more than 1 iteration steps are expected, the procedure is principally speeded up using the ideas of the TFA. Zone 1 Zone 2 Zone 3 Zone 4 Figure 1. Geometry of the model 4.1 Formulation of the problem Consider a problem of two bodies, the geometries of which are depicted in Fig. 1. The first body (fiber) occupies in undeformed state the domain Ω ' with the interfacial boundary Γ ' (rectangular parallelepiped) and the second body (concrete matrix) occupies the domain Ω' ' with the boundary Γ '', (parallelepiped). There is no external load considered but the pullout force at the face of the rectangular parallelepiped. The bodies are situated in Cartesian coordinate system Oxyz; the axis x is introduced in the radial direction of the lag. Generally, it will be assumed that the tension exceeding the tensile strength in the normal direction to the interface fiber-matrix Γ is not admitted, so that both bodies may disconnect (they mutually debond) in certain region of the common boundary Γ, which is a part of Γ ' and Γ ''. In order to hold the mathematical purity let us consider tensile strength at the contact to be zero. In practical examples, this fact will be overcome by simple introducing the tension strength but theoretical (not practical) difficulties may occur. The Mohr-Coulomb contact conditions will be taken into account. More complex study concerning the contact conditions was carried out in [1]. There a reader can find that four zones can be distinguished on the straight fiber to be pulled out of a surrounding medium. Along the first zone, additional debond can occur, if at the left face of the model described in Fig. 1 no movements in horizontal direction is admissible, for simplicity of description. The second zone is a firm bond, while in the third zone slipping in the tangential direction is possible. The zone belonging to the face (the fourth zone) is a typical opening crack. Similar situation is observed in curvilinear reinforcement. Any of the four zones can disappear but the third, which is the zone in the undeformed state. The load is considered in the following way. The pullout traction p is applied at the face of the fiber in the axial direction, prescribed by a constant value. There is symmetry about the vertical plane containing the horizontal axis x regarding both geometry and loading, so that only one half of the unit cell may be solved. Displacements are described by the vector function u = { u, v, w} of the variable x { x, y, z} =. Denote Ω Ω U Ω. The restriction,,,,,,,, of any function to Ω or Ω is denoted, respectively, by one prime or by two primes, e.g., u / Ω = u and u / Ω = u,. On the boundary displacements and tractions are prescribed in such a manner that a periodicity is assumed at all sides of a "unit cell", The loading is given in both domains due to the pullout traction p, see Fig.1, where only straight fiber is displayed. In practical examples, curvilinear fibers are assumed. Denote the set of admissible displacements u on Ω satisfying the essential boundary conditions by V. The values and the first derivatives are quadratically integrable.,,,

3 Consider Hooke s law in the form:,,,,,,,,, ε ij ( u' ) = L ijkl σ kl ε ij ( u' ') = L ijkl σ kl (1),,, where L and L are the material stiffness matrices of the fiber and the concrete, respectively. When assuming "small deformation" theory, it may be satisfactory to formulate the essential boundary condition on the contact boundary as follows (Signorini's conditions): Denote,,, [ ] u u < u n = n n a.e. on Γ (2) ;[ u] H { u... in... V n < a.e on Γ } (3) The set H is a cone of admissible displacements with respect to the essential boundary and contact conditions, V are displacements from the space of continuous functions. Suppose that we disconnect both bodies under consideration, but keep the stress and deformation state in them "frozen". Then the p = p, p, p must be introduced and their equilibrium (action and reaction law) says that: vector of contact tractions { } x y z p = p + p,,, = (4) Let us write the total energy J of both bodies assuming them separately: J(u, p ) = Π (u ) I( u, p ), (5) where 1,,, ( u ) = a( u,u) p, u, I ( p,u) = ( p' u' + p'' u'' )dγ, ( u,u) = ( ' ε ' + Π dγ 2 Γ where σ and ε are the stresses and the strains in both bodies. The Mohr-Coulomb contact conditions read as: [] < u n, n <, p [ ] = Γ a σ σ '' ε '')dω (6) Ω u n p n, p t < k ( p n ) c, if pt > k( pt ) c then p t = c (7) and k is Heaviside's function, c is the cohesion (shear strength): Formulation of the contact problem of pull out can be declared as: Find minimum u from H and maximum p of J obeying Eqs. (4) and (7). Such a posted formulation leads to the well-known Uzawa's algorithm; see [3], for example. In this way all examples were solved. Although the Coulomb influence is absent in the law (7), a slip part can occur, Heaviside's function is applied in the formulas given by Eq. (7) because of assuring a natural requirement: once the tensile strength is violated, no shear forces occur either at the trial points. 5 Straight fibers Four examples were tested to show the behavior of different shapes of fibers, which are pulled out of the concrete. In each example the stainless wire fibers are considered, for which (concrete) E' = 179 GPa, v' =.3, (steel) E'' = 41 GPa, v'' =.16 for pseudo-linear computation. The results from the examples are displayed in compact form and described as Experiment-1 to Experiment-4. First, the vertical cut crossing the fibers in the ''unit cell illustrates the geometry of the aggregate in Fig. 2. The vertical size is always identical with the horizontal depth. From the geometry one can calculate the fiber volume fraction. The fiber volume fractions range from.8 to 1.5 percent. First an approximate case is considered for straight fibers imbedded in concrete matrix. We suppose the linear behavior of both fibers and concrete matrices. Moreover, ideal bond along the interfacial zone between fibers and concrete is also assumed. Such an assumption is too rough and this is why additional assumption is introduced to improve the model, which then behaves more versatile: The flaw of various lengths appears at the face of the aggregate. In order to hold the linear behavior in the formulation of pullout problem, the flaws are considered to open successively at nodal points on the interface and the responses of unit pullout force are observed. The unit force is selected as.15 kn. The results are illustrated in the pictures, which follow in this text. In Fig. 1 vertical cut of a typical straight fiber is depicted together with the numbering of nodes, the flaws in which are introduced. The vertical cut of unit cells together with finite element meshing shows the fiber position in the cell in Fig. 2. From this picture it is seen that the position of fibers is symmetric in each cell and one symmetric part (one quarter) is solved in numerical interpretation with appropriate boundary conditions expressing the symmetry. The unit cells are selected in such a manner that typical periodic structures are acquired. Periodicity is introduced along the entire boundary but the face. It means that along both horizontal

4 boundaries and the left vertical boundary the periodicity conditions are considered in the classical sense, not in the sense of homogenized composites. The right boundary (face) is free but the fiber where the pullout force is applied. Zero in parentheses means the slope of curvilinear fibers, which are generally considered in this paper Figure 1. Shape and numbering of nodes of trial fibers Experiment-1() Experiment-2() 6,5 6 12, Experiment-3(),5 16,5 6,5 12, ,5 8,5 Experiment-4() Figure 2 Vertical cut of four cases of fiber positions in unit cells with FEM meshing In the following four pictures the responses in transversal tractions on unit.15 kn pullout force are shown. Fig. 3 depicts the distribution of transversal tractions for Experiment 1, Fig. 4 for Experiment 2, Fig. 5 for Experiment 3, and Fig. 6 for Experiment 4. On the vertical axis values of tractions are depicted, on the horizontal axis the interface is described and the orientation is enforced by numbers of nodes, see Fig. 1. The curves describe the distribution of transversal tractions on the interface from unit pullout force and various flaws. Series 1 belongs to the flaw at node 2, i.e., the fiber is debonded from points Similarly, Series i = 2,,19 belong to the flaws at nodes i + 1. From the pictures one can derive the bearing capacity of the aggregate in the sense, which is explained in what follows: If the tensile strength is, say,.3 MPa, Experiments 1 3 fail, as the tensile strength is exceeded at all nodal points. As for the Experiment 4, the strength is overcome from nodes 9 21, where debond can also be expected. On the other hand, at nodes 1 to 8 the strength is not exceeded and fiber presence increases the bearing capacity of concrete but the crack starting from nodal point 9. If, for example, the tensile strength is.4 MPa, Experiment 1 is undeceived as for the admissible pullout force, in Experiment 2 the real flaw is between nodes 17 and 18, in Experiment 3 the force is pulled out of the aggregate, and the force in Experiment 4 stiffens the aggregate, of course with respect to the previous case (tensile strength is.3). In this manner we can estimate the bearing capacity for other particular cases considered in this study. Similarly, assuming the unit force introduced as.15 kn, and, say, in Experiment 1 we require the fiber is somewhere bonded and interacts with concrete for tensile strength to be equal to.2 MPa, the force must be less then.75 kn (one half of the unit force). In this way we can discuss various cases under study, as the problem is linear, although debond (flaws) is considered.

5 MPa,7,6,5,4,3,2,1 -,1 -,2 Experiment Number of nodes Series1 Series2 Series3 Series4 Series5 Series6 Series7 Series8 Series9 Series1 Series11 Series12 Series13 Series14 Series15 Series16 Series17 Series18 Series19 Figure 3. Distribution of transversal tractions - various flaws in Experiment 1 MPa,6,5,4,3,2,1 -,1 Experiment Number of nodes Series1 Series2 Series3 Series4 Series5 Series6 Series7 Series8 Series9 Series1 Series11 Series12 Series13 Series14 Series15 Series16 Series17 Series18 Series19 Figure 4. Distribution of transversal tractions - various flaws in Experiment 2 MPa,8,7,6,5,4,3,2,1 -,1 -,2 Experiment Number of nodes Series1 Series2 Series3 Series4 Series5 Series6 Series7 Series8 Series9 Series1 Series11 Series12 Series13 Series14 Series15 Series16 Series17 Series18 Series19 Figure 5. Distribution of transversal tractions - various flaws in Experiment 3

6 MPa,7,6,5,4,3,2,1 -,1 -,2 Experiment Number of nodes Series1 Series2 Series3 Series4 Series5 Series6 Series7 Series8 Series9 Series1 Series11 Series12 Series13 Series14 Series15 Series16 Series17 Series18 Series19 Figure 6. Distribution of transversal tractions for various flaws in Experiment 4 6. Curvilinear fibers - plastic matrix Now we concentrate our attention on nonlinear behavior of concrete according to [3]. Also, the interfacial zones behave nonlinearly according to the law described in chapter 4. Fig. 7 shows the possible geometry and numbering used for denotation of values on the interface Figure 7. Numbering of nodes of fiber for the slope of 45 Since the optimization problem strongly depends on the meshing, some types of meshing in the neighborhood of fibers are shown in Fig. 8. The meshings for slopes of successively, 1, 3, 45, 6, and 9 degrees are illustrated in these pictures in the vicinity of fibers. () Degree (1) Degree (3) Degree (45) Degree (6) Degree (9) Degree Figure 8. Meshing in the neighborhood of fibers Although the mesh changes sometimes irregularly, the results seem to be relatively reliable and reasonable. In each example the stainless wire fibers are considered. The material properties are described as follows: (concrete) E' = 179 GPa, E' p = 154 GPa, E' r =

7 1 GPa, v' =.33 v' p =.45, and (steel) E'' = 33 GPa, E'' p = 22 GPa, E'' r = 5 GPa, v'' =.16, v'' =.38. The primed quantities belong to steel and two primed belong to concrete. The subscript p denotes peak and r residual values. Along the interfacial boundary the brittle behavior is considered. It means that ideally elasto-plastic law is obeyed on this boundary. Two cases are studied in details: 1. In the first case, the shear bond strength is c = 4.35 MPa. The tensile strength is.42 Mpa 2. The second case is characterized by interfacial properties given by the shear bond strength c = 2.3 Mpa and the tensile strength is.2 MPa. In both cases the same properties of fibers and concrete matrix are considered, but the interfacial material properties differ basically: In the first case relatively high parameters are implied and in the second case the interface is not as stiff as in the first case. The optimization procedure leads to the admissible pullout forces acting at the face of the samples, the distribution of which is seen in Figs. 9 to 12. In Fig. 9 results from Experiment 1 are depicted, in Fig. 1 the admissible forces from Experiment 2 are described, in Fig. 11 from Experiment 3 and in Fig. 12 forces from Experiment 4 are described, all in dependence of slope angles, which are considered as, 1, 3, 45, 6, and 9 degrees. Assuming the dependence of pullout force slope of the curvilinear fibers to be continuous, in all experiments the straight fibers have the optimal shapes. Similarly, for the second case the distributions of admissible forces are displayed in Figs. 13 to 16. Experiment-1,6,5,4,3,2, Figure 9. Distribution of admissible forces in Experiment 1 the first case Experiment-2,7,6,5,4,3,2, Figure 1. Distribution of admissible forces in Experiment 2 the first case Experiment 1 considers the volume fraction ratio about.7 percent and larger distances between fibers in vertical direction. Long debond zone, but relatively short sliding zones are recorded in this example. Due to the bending of the fiber the axial contact forces are lower in the upper part and vice versa for lower part. This observation is valid for all experiments. Sometimes an additional debond is seen at the points in from of the slope for slopes starting at 4 degrees. Sudden drop of transversal tractions indicates occurrence of a crack.

8 Experiment-3,5,48,46,44,42,4, Figure 11. Distribution of admissible forces in Experiment 3 the first case Experiment-4,14,12,1,8,6,4, Figure 12. Distribution of admissible forces in Experiment 4 the first case Experiment-1,25,2,15,1, Figure 13. Distribution of admissible forces in Experiment 1 the second case In Experiment 2 we consider relatively higher volume fraction ratio, about 1.5. On the other hand dense distribution of fibers in vertical direction is introduced. The unit cell is the smallest from all four cases. Basically higher stresses are observed at the decisive points. The debonding region is relatively short.

9 ,25,2,15,1,5 Experiment Figure 14. Distribution of admissible forces in Experiment 2 the second case Experiment-3,1,8,6,4,2, Figure 15. Distribution of admissible forces in Experiment 1 the second case Experiment-4,14,12,1,8,6,4,2, Figure 16. Distribution of admissible forces in Experiment 4 the second case Experiment 3 has extremely low volume fraction ratio, namely.4, and very long distances between fibers in vertical direction. It leads to decreasing of stresses concentrated at the point of observation. The admissible force appears to be much higher then in the previous case.

10 Experiment 4 has similar geometry as the first one, and the behavior of the admissible force is also similar. The fibers are almost continuous in horizontal direction. In the latter two cases it is necessary to count with long debonding region in the front part of the fibers. It was shown that in the case of high volume fiber ratio the admissible pullout force is basically lower then in the case of course distribution of fibers. On the other hand, it is necessary to note that the pullout force in this case bears distributed loading, i.e., the tensile force being applied to a cross section of the beam is distributed in such a way that smaller part of the overall loading is applied to this force. Very interesting is a comparison of the previous pseudo-linear case (PL) and the non-linear case (NL) from the last study. The results are summed up in Table 1. Table 1. Comparison of admissible pullout forces for two cases studied in this paper Experiment 1 Experiment 2 Experiment 3 Experiment 4 Pseudo-linear Non-linear Concerning the data from this section, very similar results are obtained for Experiment 2, namely the admissible force in the PL is.79 kn and in the NL.6 kn. The differences in other cases are also closed each other. The higher critical forces in the PL are, of course, due to too stiff concrete matrix consideration. 5. Observations While in the most papers dealing with pullout problem the fibers are straight, in this paper also different shapes of fibers are taken into account. The cohesion of fibers and the concrete matrix is considered and the optimal pullout force is studied. Large test examples were treated and several interfacial conditions were applied. From these examples four experiments are selected to show, first, the material behavior of the interface in the case of the shape envisaged, and second, a possibility of the application of the numerical method introduced at the beginning of this paper. In comparison with the papers [1] and [2] and the papers cited in [1] and [2] it uniquely follows that for the cement-based matrix with steel fibers the shape plays very important role for the interfacial mechanical behavior and, consequently, for the bearing capacity of the aggregate. The optimal slope of Dramix type fibers was the principal goal of chapter 4. It has been shown that the shape of the fibers has very important role in the bearing capacity of fiber-concrete composite. From typical examples of unit cells numerical studies showed important impacts on the behavior of the aggregate. Two cases of interfacial conditions have been studied. In the first case (relatively stiff interface) the best solution is found to be for straight fibers for all arrangements of fibers in unit cells. In the second case (the interfacial zone is chosen in such a way that a relatively lower stiffness is considered), the slope ranging from 45 to 6 degrees seems to be optimal under condition that has been taken into consideration, i.e., the pullout force damages the composite either because of failure of contact conditions or due to attaining the admissible stress (compressive strength). Moreover, the maximum value of admissible force is the higher the course the density of distribution of fibers is. But it is worth noting that the force on the overall cross section is approximately regularly distributed on fibers. The mechanical behavior of fibers is considered as linear elastic, but the material of concrete obeys von Mises Huber Hencky elastic-plastic criterion. This will desire more sophisticated treatment then in the case of straight fibers with linear concrete matrix. It will start with generalized Transformation field analysis, which involves an influence of eigenparameters, see, e.g., [4]. 6. Conclusions For four typical unit cells pullout problem of generally curvilinear fibers out of concrete matrix is formulated and solved in this paper. Starting with modified Uzawa s algorithm and tricks taken from the transformation field analysis extreme pullout force is described for discrete set of slopes of fibers and from this optimal slope is attained. Simplified version is considered for straight fibers as a special case of the general approach. It is seen that the deviations of values from the more accurate computations are not decisive and can be estimated by some constant. Acknowledgment: This paper has been financially supported by the Grand Agency of the Czech Republic, projects No. 13/3/183 and 13/4/ References 1. Prochazka P and Sejnoha, M. Development of Debond Region of Lag Model. Int. J. Computers & Structures 1995, 55(2): Procházka P and Sejnoha M. Pull out problem of lag model. Proc. of BETECH 96, 1996, Honolulu, Hawaii, Duvaut G and Lions J L. Les inéquations en méchanique et en physique. Paris: Dunod, Dvorak G J Prochazka P. Thick-walled Composite Cylinders with Optimal Fiber Prestress, Composites, Part B, 27B, 1996:

Debonding process in composites using BEM

Debonding process in composites using BEM Boundary Elements XXVII 331 Debonding process in composites using BEM P. Prochazka & M. Valek Czech Technical University, Prague, Czech Republic Abstract The paper deals with the debonding fiber-matrix

More information

Homogenization of debonding composites

Homogenization of debonding composites Homogenization of debonding composites P. Prochazka, A. Bilek-Larrea Department of Structural Mechanics, Faculty ofmechanical Engineering, Czech Technical University, Thakurova 7, 166 29 Prague 6, Czech

More information

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in Sabah Shawkat Cabinet of Structural Engineering 17 3.6 Shear walls Walls carrying vertical loads should be designed as columns. Basically walls are designed in the same manner as columns, but there are

More information

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS 9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING

More information

5 ADVANCED FRACTURE MODELS

5 ADVANCED FRACTURE MODELS Essentially, all models are wrong, but some are useful George E.P. Box, (Box and Draper, 1987) 5 ADVANCED FRACTURE MODELS In the previous chapter it was shown that the MOR parameter cannot be relied upon

More information

Abstract. 1 Introduction

Abstract. 1 Introduction Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut

More information

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations W. Van Paepegem *, I. De Baere and J. Degrieck Ghent

More information

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic

More information

Interface properties between a steel pre-stressing strand and an epoxy matrix.

Interface properties between a steel pre-stressing strand and an epoxy matrix. Interface properties between a steel pre-stressing strand and an epoxy matrix. J. Van Vooren 1, B. Van Vooren 2, D. Van Hemelrijck 1, 1 Free University of Brussel, Dept. of Mechanics of Materials and Constructions,

More information

University of Sheffield The development of finite elements for 3D structural analysis in fire

University of Sheffield The development of finite elements for 3D structural analysis in fire The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures

More information

A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES

A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES Konuralp Girgin (Ph.D. Thesis, Institute of Science and Technology,

More information

Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas

Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas San Miguel de Tucuman, Argentina September 14 th, 2011 Seminary on Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas Antonio Caggiano 1, Guillermo Etse 2, Enzo

More information

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric

More information

Classical fracture and failure hypotheses

Classical fracture and failure hypotheses : Chapter 2 Classical fracture and failure hypotheses In this chapter, a brief outline on classical fracture and failure hypotheses for materials under static loading will be given. The word classical

More information

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,

More information

Chapter 3. Load and Stress Analysis

Chapter 3. Load and Stress Analysis Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3

More information

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd

More information

Table of Contents. Preface...xvii. Part 1. Level

Table of Contents. Preface...xvii. Part 1. Level Preface...xvii Part 1. Level 1... 1 Chapter 1. The Basics of Linear Elastic Behavior... 3 1.1. Cohesion forces... 4 1.2. The notion of stress... 6 1.2.1. Definition... 6 1.2.2. Graphical representation...

More information

THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.

THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE CASE OF A TERRACE UNIT. John N Karadelis 1. INTRODUCTION. Aim to replicate the behaviour of reinforced concrete in a multi-scale

More information

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment 7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment à It is more difficult to obtain an exact solution to this problem since the presence of the shear force means that

More information

Basic Energy Principles in Stiffness Analysis

Basic Energy Principles in Stiffness Analysis Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting

More information

Discrete Element Modelling of a Reinforced Concrete Structure

Discrete Element Modelling of a Reinforced Concrete Structure Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr

More information

Practice Final Examination. Please initial the statement below to show that you have read it

Practice Final Examination. Please initial the statement below to show that you have read it EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

More information

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global

More information

Chapter 8. Shear and Diagonal Tension

Chapter 8. Shear and Diagonal Tension Chapter 8. and Diagonal Tension 8.1. READING ASSIGNMENT Text Chapter 4; Sections 4.1-4.5 Code Chapter 11; Sections 11.1.1, 11.3, 11.5.1, 11.5.3, 11.5.4, 11.5.5.1, and 11.5.6 8.2. INTRODUCTION OF SHEAR

More information

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi

More information

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS Abdul Qadir Bhatti 1, Norimitsu Kishi 2 and Khaliq U Rehman Shad 3 1 Assistant Professor, Dept. of Structural Engineering,

More information

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,

More information

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having

More information

SHEAR RESISTANCE BETWEEN CONCRETE-CONCRETE SURFACES

SHEAR RESISTANCE BETWEEN CONCRETE-CONCRETE SURFACES , DOI: 10.2478/sjce-2013-0018 M. KOVAČOVIC SHEAR RESISTANCE BETWEEN CONCRETE-CONCRETE SURFACES Marek KOVAČOVIC Email: marek.kovacovic@gmail.com Research field: composite structures Address: Department

More information

Flexure: Behavior and Nominal Strength of Beam Sections

Flexure: Behavior and Nominal Strength of Beam Sections 4 5000 4000 (increased d ) (increased f (increased A s or f y ) c or b) Flexure: Behavior and Nominal Strength of Beam Sections Moment (kip-in.) 3000 2000 1000 0 0 (basic) (A s 0.5A s ) 0.0005 0.001 0.0015

More information

Structural Analysis I Chapter 4 - Torsion TORSION

Structural Analysis I Chapter 4 - Torsion TORSION ORSION orsional stress results from the action of torsional or twisting moments acting about the longitudinal axis of a shaft. he effect of the application of a torsional moment, combined with appropriate

More information

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield

More information

Numerical analysis of the mechanical response of wood glulam beams reinforced through the thickness by FRP rods.

Numerical analysis of the mechanical response of wood glulam beams reinforced through the thickness by FRP rods. Numerical analysis of the mechanical response of wood glulam beams reinforced through the thickness by FRP rods. Giuseppe Giambanco 1, Tiziana Turetta 1, Alessia Cottone 1 1 Department of Structural, Aerospace

More information

THE COLLAPSE LOAD IN SUBMARINE PIPELINES UNDER COMPRESSIVE LOAD AND INTERNAL PRESSURE

THE COLLAPSE LOAD IN SUBMARINE PIPELINES UNDER COMPRESSIVE LOAD AND INTERNAL PRESSURE SDSS Rio 010 STABILITY AND DUCTILITY OF STEEL STRUCTURES E. Batista,. Vellasco, L. de Lima (Eds.) Rio de Janeiro, Brazil, September 8-10, 010 THE COLLASE LOAD IN SUBMARINE IELINES UNDER COMRESSIVE LOAD

More information

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 Samantha Ramirez TENSION AND COMPRESSION TESTS Tension and compression tests are used primarily to determine the relationship between σ avg and ε avg in any material.

More information

DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)

DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS) DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types

More information

NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS

NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.

More information

Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING )

Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress

More information

Finite Element Modelling with Plastic Hinges

Finite Element Modelling with Plastic Hinges 01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated post-yield behaviour in one or more degrees of freedom. Hinges only

More information

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS 5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.

More information

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro

More information

Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP

Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP Proceedings of the World Congress on Engineering 21 Vol II WCE 21, June 2 - July 1, 21, London, U.K. Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP Lihua Huang,

More information

ME 243. Mechanics of Solids

ME 243. Mechanics of Solids ME 243 Mechanics of Solids Lecture 2: Stress and Strain Ahmad Shahedi Shakil Lecturer, Dept. of Mechanical Engg, BUET E-mail: sshakil@me.buet.ac.bd, shakil6791@gmail.com Website: teacher.buet.ac.bd/sshakil

More information

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS 6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement

More information

Deflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis

Deflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis Rotating Machinery, 10(4): 283 291, 2004 Copyright c Taylor & Francis Inc. ISSN: 1023-621X print / 1542-3034 online DOI: 10.1080/10236210490447728 Deflections and Strains in Cracked Shafts due to Rotating

More information

IMPACT RESPONSE ANALYSIS OF LARGE SCALE RC GIRDER WITH SAND CUSHION

IMPACT RESPONSE ANALYSIS OF LARGE SCALE RC GIRDER WITH SAND CUSHION -Technical Paper- IMPACT RESPONSE ANALYSIS OF LARGE SCALE RC GIRDER WITH SAND CUSHION Abdul Qadir BHATTI *1, Norimitsu KISHI *2, Shin-ya OKADA *3 and Hisashi KONNO *4 ABSTRACT In order to establish a proper

More information

A NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES

A NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDFCATO Publishers, D-79104 Freiburg, Germany A NEW SMPLFED AND EFFCENT TECHNQUE FOR FRACTURE BEHAVOR ANALYSS OF CONCRETE STRUCTURES K.

More information

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.

More information

Static Failure (pg 206)

Static Failure (pg 206) Static Failure (pg 06) All material followed Hookeʹs law which states that strain is proportional to stress applied, until it exceed the proportional limits. It will reach and exceed the elastic limit

More information

REGRESSION MODELING FOR STRENGTH AND TOUGHNESS EVALUATION OF HYBRID FIBRE REINFORCED CONCRETE

REGRESSION MODELING FOR STRENGTH AND TOUGHNESS EVALUATION OF HYBRID FIBRE REINFORCED CONCRETE REGRESSION MODELING FOR STRENGTH AND TOUGHNESS EVALUATION OF HYBRID FIBRE REINFORCED CONCRETE S. Eswari 1, P. N. Raghunath and S. Kothandaraman 1 1 Department of Civil Engineering, Pondicherry Engineering

More information

A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES

A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES Roque Luiz Pitangueira 1 and Raul Rosas e Silva 2 1 Department of Structural Engineering -Federal University of Minas Gerais

More information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These

More information

MODELING OF THE BEHAVIOR OF WOVEN LAMINATED COMPOSITES UNTIL RUPTURE

MODELING OF THE BEHAVIOR OF WOVEN LAMINATED COMPOSITES UNTIL RUPTURE MODELING OF THE BEHAVIOR OF WOVEN LAMINATED COMPOSITES UNTIL RUPTURE Jean Paul Charles, Christian Hochard,3, Pierre Antoine Aubourg,3 Eurocopter, 375 Marignane cedex, France Unimeca, 6 rue J. Curie, 3453

More information

Special edition paper

Special edition paper Development of New Aseismatic Structure Using Escalators Kazunori Sasaki* Atsushi Hayashi* Hajime Yoshida** Toru Masuda* Aseismatic reinforcement work is often carried out in parallel with improvement

More information

Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis

Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis R.E. Hammah, T. Yacoub, B. Corkum & F. Wibowo Rocscience Inc., Toronto, Canada J.H. Curran Department of Civil Engineering

More information

FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING WEB DEPTH

FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING WEB DEPTH Journal of Engineering Science and Technology Vol. 12, No. 11 (2017) 2839-2854 School of Engineering, Taylor s University FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING

More information

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT

More information

Continuous Modeling Technique of Fiber Pullout from a Cement Matrix with Different Interface Mechanical Properties Using Finite Element Program

Continuous Modeling Technique of Fiber Pullout from a Cement Matrix with Different Interface Mechanical Properties Using Finite Element Program 1937 Continuous Modeling Technique of Fiber Pullout from a Cement Matrix with Different Interface Mechanical Properties Using Finite Element Program Abstract Fiber-matrix interface performance has a great

More information

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the

More information

3D dynamic crack propagation analysis with PDS-FEM

3D dynamic crack propagation analysis with PDS-FEM JAMSTEC-R IFREE Special Issue, November 2009 3D dynamic crack propagation analysis with PDS-FEM 1*, Hide Sakaguchi 1, Kenji Oguni 2, Muneo Hori 2 discretizations. Keywords +81-45-778-5972 lalith@jamstec.go.jp

More information

NUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL

NUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL Numerical simulations of corners in RC frames using Strut-and-Tie Method and CDP model XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J.

More information

Adhesive Joints Theory (and use of innovative joints) ERIK SERRANO STRUCTURAL MECHANICS, LUND UNIVERSITY

Adhesive Joints Theory (and use of innovative joints) ERIK SERRANO STRUCTURAL MECHANICS, LUND UNIVERSITY Adhesive Joints Theory (and use of innovative joints) ERIK SERRANO STRUCTURAL MECHANICS, LUND UNIVERSITY Wood and Timber Why I m intrigued From this to this! via this Fibre deviation close to knots and

More information

CFRP. FRP FRP. Abaqus

CFRP. FRP FRP. Abaqus C Epoxy C baqus C FE baqus C amene.kia@gmail.com EROG 1 Debonding : / C (mm) E c (MPa) f' c (MPa) f' t (MPa) a S L d h b 1 2 3 1 2 3 baqus C Leong Kok Leong Kok CDP E p (GPa) (MPa) E a (MPa) G a (MPa)

More information

Chapter Two: Mechanical Properties of materials

Chapter Two: Mechanical Properties of materials Chapter Two: Mechanical Properties of materials Time : 16 Hours An important consideration in the choice of a material is the way it behave when subjected to force. The mechanical properties of a material

More information

ME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft.

ME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft. ME 323 - Final Exam Name December 15, 2015 Instructor (circle) PROEM NO. 4 Part A (2 points max.) Krousgrill 11:30AM-12:20PM Ghosh 2:30-3:20PM Gonzalez 12:30-1:20PM Zhao 4:30-5:20PM M (x) y 20 kip ft 0.2

More information

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens. a). Cohesive Failure b). Interfacial Failure c). Oscillatory Failure d). Alternating Failure Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double

More information

Forces on piles preventing debris slope slips

Forces on piles preventing debris slope slips Risk Analysis VII PI-637 Forces on piles preventing debris slope slips J. Vacek & S. Hrachová Czech Technical University, Klokner Institute, Czech Republic Abstract Failure of rock mass is not a static

More information

Identification of interface properties using Fibre Bragg Grating sensors in a fibre pull-out test Gabriel Dunkel, Laurent Humbert and John Botsis

Identification of interface properties using Fibre Bragg Grating sensors in a fibre pull-out test Gabriel Dunkel, Laurent Humbert and John Botsis Identification of interface properties using Fibre Bragg Grating sensors in a fibre pull-out test Gabriel Dunkel, Laurent Humbert and John Botsis Laboratory of Applied Mechanics and Reliability Analysis

More information

THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT

THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT Víctor Martínez 1, Alfredo Güemes 2, Norbert Blanco 1, Josep Costa 1 1 Escola Politècnica Superior. Universitat de Girona. Girona, Spain (17071) 2

More information

1. Background. is usually significantly lower than it is in uniaxial tension

1. Background. is usually significantly lower than it is in uniaxial tension NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and

More information

6. Bending CHAPTER OBJECTIVES

6. Bending CHAPTER OBJECTIVES CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where

More information

Mechanics of Materials Primer

Mechanics of Materials Primer Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus

More information

Pressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials

Pressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials Pressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials Pressure Vessels: In the previous lectures we have discussed elements subjected

More information

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e

More information

Role of Force Resultant Interaction on Ultra-High Performance Concrete

Role of Force Resultant Interaction on Ultra-High Performance Concrete First International Interactive Symposium on UHPC 216 Role of Force Resultant Interaction on Ultra-High Performance Concrete Author(s) & Affiliation: Titchenda Chan (1), Kevin R. Mackie (1), and Jun Xia

More information

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components Salim Mirza Element Materials Technology Hitchin, UK Introduction Fibre reinforced elastomers are used in many applications,

More information

MODELLING NON-LINEAR BEHAVIOUR OF STEEL FIBRE REINFORCED CONCRETE

MODELLING NON-LINEAR BEHAVIOUR OF STEEL FIBRE REINFORCED CONCRETE 6th RILEM Symposium on Fibre-Reinforced Concretes (FRC) - BEFIB - September, Varenna, Italy MODELLING NON-LINEAR BEHAVIOUR OF STEEL FIBRE REINFORCED CONCRETE W. A. Elsaigh, J. M. Robberts and E.P. Kearsley

More information

CHAPTER 3 EXPERIMENTAL STUDY

CHAPTER 3 EXPERIMENTAL STUDY Experimental Study 42 CHAPTER 3 EXPERIMENTAL STUDY 3.1. INTRODUCTION The experimental study that has been carried out in this thesis has two main objectives: 1. Characterise the concrete behaviour in mode

More information

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE 1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for

More information

Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation

Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,

More information

Loading capacity of yielding connections used in steel arch roadway supports

Loading capacity of yielding connections used in steel arch roadway supports Ground Support 2013 Y. Potvin and B. Brady (eds) 2013 Australian Centre for Geomechanics, Perth, ISBN 978-0-9806154-7-0 https://papers.acg.uwa.edu.au/p/1304_31_horyl/ Loading capacity of yielding connections

More information

Advanced numerical modelling methods of rock bolt performance in underground mines

Advanced numerical modelling methods of rock bolt performance in underground mines University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2010 Advanced numerical modelling methods of rock bolt performance in underground mines

More information

FRACTURE MECHANICS APPROACHES STRENGTHENING USING FRP MATERIALS

FRACTURE MECHANICS APPROACHES STRENGTHENING USING FRP MATERIALS Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDIFICATIO Publishers, D-79104 Freiburg, Germany FRACTURE MECHANICS APPROACHES STRENGTHENING USING FRP MATERIALS Triantafillou Department

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

More information

five Mechanics of Materials 1 ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2017 lecture

five Mechanics of Materials 1 ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2017 lecture ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2017 lecture five mechanics www.carttalk.com of materials Mechanics of Materials 1 Mechanics of Materials MECHANICS MATERIALS

More information

Finite element analysis of diagonal tension failure in RC beams

Finite element analysis of diagonal tension failure in RC beams Finite element analysis of diagonal tension failure in RC beams T. Hasegawa Institute of Technology, Shimizu Corporation, Tokyo, Japan ABSTRACT: Finite element analysis of diagonal tension failure in a

More information

ME 2570 MECHANICS OF MATERIALS

ME 2570 MECHANICS OF MATERIALS ME 2570 MECHANICS OF MATERIALS Chapter III. Mechanical Properties of Materials 1 Tension and Compression Test The strength of a material depends on its ability to sustain a load without undue deformation

More information

Numerical simulation of delamination onset and growth in laminated composites

Numerical simulation of delamination onset and growth in laminated composites Numerical simulation of delamination onset and growth in laminated composites G. Wimmer, C. Schuecker, H.E. Pettermann Austrian Aeronautics Research (AAR) / Network for Materials and Engineering at the

More information

U.S. South America Workshop. Mechanics and Advanced Materials Research and Education. Rio de Janeiro, Brazil. August 2 6, Steven L.

U.S. South America Workshop. Mechanics and Advanced Materials Research and Education. Rio de Janeiro, Brazil. August 2 6, Steven L. Computational Modeling of Composite and Functionally Graded Materials U.S. South America Workshop Mechanics and Advanced Materials Research and Education Rio de Janeiro, Brazil August 2 6, 2002 Steven

More information

MECHANICS OF MATERIALS

MECHANICS OF MATERIALS Third E CHAPTER 2 Stress MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech University and Strain Axial Loading Contents Stress & Strain:

More information

MATERIALS FOR CIVIL AND CONSTRUCTION ENGINEERS

MATERIALS FOR CIVIL AND CONSTRUCTION ENGINEERS MATERIALS FOR CIVIL AND CONSTRUCTION ENGINEERS 3 rd Edition Michael S. Mamlouk Arizona State University John P. Zaniewski West Virginia University Solution Manual FOREWORD This solution manual includes

More information

Mechanical analysis of timber connection using 3D finite element model

Mechanical analysis of timber connection using 3D finite element model Mechanical analysis of timber connection using 3D finite element model Bohan XU Ph.D Student Civil Engineering Laboratory (CUST) Clermont-Ferrand, France Mustapha TAAZOUNT Dr-Ing Civil Engineering Laboratory

More information

Nonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads

Nonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads Transactions of the 7 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 7) Prague, Czech Republic, August 7, 00 Paper # J0- Nonlinear Analysis of Reinforced Concrete Shells

More information

The Frictional Regime

The Frictional Regime The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation

More information

Appendix J. Example of Proposed Changes

Appendix J. Example of Proposed Changes Appendix J Example of Proposed Changes J.1 Introduction The proposed changes are illustrated with reference to a 200-ft, single span, Washington DOT WF bridge girder with debonded strands and no skew.

More information

Numerical Characterization of Concrete Heterogeneity

Numerical Characterization of Concrete Heterogeneity Vol. Materials 5, No. Research, 3, 2002Vol. 5, No. 3, Statistical 309-314, 2002. Characterization of the Concrete Numerical Modeling of Size Effect In Heterogeneity 2002 309 Numerical Characterization

More information

SERVICEABILITY OF BEAMS AND ONE-WAY SLABS

SERVICEABILITY OF BEAMS AND ONE-WAY SLABS CHAPTER REINFORCED CONCRETE Reinforced Concrete Design A Fundamental Approach - Fifth Edition Fifth Edition SERVICEABILITY OF BEAMS AND ONE-WAY SLABS A. J. Clark School of Engineering Department of Civil

More information

3D Finite Element analysis of stud anchors with large head and embedment depth

3D Finite Element analysis of stud anchors with large head and embedment depth 3D Finite Element analysis of stud anchors with large head and embedment depth G. Periškić, J. Ožbolt & R. Eligehausen Institute for Construction Materials, University of Stuttgart, Stuttgart, Germany

More information