• Structural Engineering and Mechanics
• /
• 제10권2호
• /
• pp.181-195
• /
• 2000

A two-step procedure for the application of non linear constrained programming to the limit analysis of rigid brick-block systems with no-tension and frictional interface is implemented and applied to various masonry structures. In the first step, a linear problem of programming, obtained by applying the upper bound theorem of limit analysis to systems of blocks interacting through no-tension and dilatant interfaces, is solved. The solution of this linear program is then employed as initial guess for a non linear and non convex problem of programming, obtained applying both the 'mechanism' and the 'equilibrium' approaches to the same block system with no-tension and frictional interfaces. The optimiser used is based on the sequential quadratic programming. The gradients of the constraints required are provided directly in symbolic form. In this way the program easily converges to the optimal solution even for systems with many degrees of freedom. Various numerical analyses showed that the procedure allows a reliable investigation of the ultimate behaviour of jointed structures, such as stone masonry structures, under statical load conditions.

• Structural Engineering and Mechanics
• /
• 제22권1호
• /
• pp.107-130
• /
• 2006

The paper deals with the numerical solution of the dynamic problem of masonry structures. Masonry is modelled as a non-linear elastic material with zero tensile strength and infinite compressive strength. Due to the non-linearity of the adopted constitutive equation, the equations of the motion must be integrated directly. In particular, we apply the Newmark or the Hilber-Hughes-Taylor methods implemented in code NOSA to perform the time integration of the system of ordinary differential equations obtained from discretising the structure into finite elements. Moreover, with the aim of evaluating the effectiveness of these two methods, some dynamic problems, whose explicit solutions are known, have been solved numerically. Comparisons between the exact solutions and the corresponding approximate solutions obtained via the Newmark and Hilber-Hughes-Taylor methods show that in the cases under consideration both numerical methods yield satisfactory results.

• Computers and Concrete
• /
• 제16권1호
• /
• pp.17-35
• /
• 2015

In this study, linear and non-linear response of a masonry wall that includes an opening was presented. The masonry wall was modeled with two-dimensional finite elements. Smeared crack model that includes the strain softening behavior was selected to the masonry wall material. For the numerical application, linear and non-linear analyses of the masonry wall were carried out using east-west and vertical components of the 1992 Erzincan and 2003 $Bing{\ddot{o}}l$ earthquake acceleration records. Linear and non-linear solutions were compared each other. The displacement and stress results at the selected points of the masonry wall and crack propagation in the masonry wall were presented for both earthquake acceleration records.

• Structural Engineering and Mechanics
• /
• 제17권6호
• /
• pp.751-766
• /
• 2004

This paper extends the use of the hierarchic degenerated shell element to geometric non-linear analysis of composite laminated skew plates by the p-version of the finite element method. For the geometric non-linear analysis, the total Lagrangian formulation is adopted with moderately large displacement and small strain being accounted for in the sense of von Karman hypothesis. The present model is based on equivalent-single layer laminate theory with the first order shear deformation including a shear correction factor of 5/6. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. A wide variety of linear and non-linear results obtained by the p-version finite element model are presented for the laminated skew plates as well as laminated square plates. A numerical analysis is made to illustrate the influence of the geometric non-linear effect on the transverse deflections and the stresses with respect to width/depth ratio (a/h), skew angle (${\beta}$), and stacking sequence of layers. The present results are in good agreement with the results in literatures.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2004년도 추계학술대회 논문집
• /
• pp.411-416
• /
• 2004

This paper introduces a fast Fourier transform (FFT)-based spectral analysis method for the transient responses as well as the steady-state responses of linear dynamic systems with non-proportional damping. The force vibration of a non-proportionally damped three-DOF system is considered as the illustrative numerical example. The proposed spectral analysis method is evaluated by comparing with the numerical solution obtained by the Runge-Kutta method

• Structural Engineering and Mechanics
• /
• 제1권1호
• /
• pp.87-102
• /
• 1993

A $C^{\circ}$ continuous finite element formulation of a higher order displacement theory is presented for predicting linear and geometrically non-linear in the sense of von Karman transient responses of composite and sandwich plates. The displacement model accounts for non-linear cubic variation of tangential displacement components through the thickness of the laminate and the theory requires no shear correction coefficients. In the time domain, the explicit central difference integrator is used in conjunction with the special mass matrix diagonalization scheme which conserves the total mass of the element and included effects due to rotary inertia terms. The parametric effects of the time step, finite element mesh, lamination scheme and orthotropy on the linear and geometrically non-linear responses are investigated. Numerical results for central transverse deflection, stresses and stress resultants are presented for square/rectangular composite and sandwich plates under various boundary conditions and loadings and these are compared with the results from other sources. Some new results are also tabulated for future reference.

• 한국해양공학회지
• /
• 제17권6호
• /
• pp.32-37
• /
• 2003

In this paper, finite element method is applied for the numerical analysis of wave height distribution. The mild-slope equation is used as the basic equation. The key of this model is to impose the effect of nonlinear waves. Numerical results are presented and agreed well with the results from experimental measurements and other numerical analysis. The present method to determine wave height distribution can be broadly utilized for the analysis of new harbor and port designs in the future.

• Earthquakes and Structures
• /
• 제4권3호
• /
• pp.241-264
• /
• 2013

Since the beginning of the twentieth century, the progressive and rapid spread of reinforced concrete (RC) has led to the adoption of mixed masonry-RC solutions, such as the confined masonry. However, together with structures conceived with a definite role for earthquake behaviour, the spreading of RC technology has caused the birth of mixed solutions inspired more by functional aspects than by structural ones, such as: internal masonry walls replaced by RC frames, RC walls inserted to build staircases or raising made from RC frames. Usually, since these interventions rise from a spontaneous build-up, any capacity design or ductility concepts are neglected being designed only to bear vertical loads: thus, the vulnerability assessment of this class becomes crucial. To investigate the non-linear seismic response of these structures, suitable models and effective numerical tools are needed. Among the various modelling approaches proposed in the literature and codes, the authors focus their attention on the equivalent frame model. After a brief description of the adopted model and its numerical validation, the authors aim to point out some specific peculiarities of the seismic response of mixed masonry-RC structures and their repercussions on safety verification procedures (referring in particular way to the non-linear static ones). In particular, the results of non-linear static analyses performed parametrically to various configurations representative of different interventions are discussed.

• 한국철도학회논문집
• /
• 제5권1호
• /
• pp.40-47
• /
• 2002

In this study, and analytical method for ride sensitivity analysis of a train with non-linear suspension elements are proposed. Non-linear characteristics of springs and dampers for primary and secondary suspensions of a train are parameterized using polynomial interpolation. Vertical dynamic model of a three-body train running on straight rail with the predetermined roughness expressed in terms of spectral density function is set up and its equations of motion for ride analysis are derived. Using the direct differentiation method, sensitivity equations of the vertical dynamic model with respect to design parameters associated with non-linearity of suspensions are obtained. Based on the sensitivity analysis, improvement of ride is achieved by varying appropriate suspension parameters.

• Geomechanics and Engineering
• /
• 제27권5호
• /
• pp.465-479
• /
• 2021

One of the important causes of building and infrastructure failure, such as bridges on pile foundations, is the placement of the piles in liquefiable soil that can become unstable under seismic loads. Therefore, the overarching aim of this study is to investigate the seismic behavior of a soil-pile system in liquefiable soil using three-dimensional numerical FEM analysis, including soil-pile interaction. Effective parameters on concrete pile response, involving the pile diameter, pile length, soil type, and base acceleration, were considered in the framework of finite element non-linear dynamic analysis. The constitutive model of soil was considered as elasto-plastic kinematic-isotropic hardening. First, the finite element model was verified by comparing the variations on the pile response with the measured data from the centrifuge tests, and there was a strong agreement between the numerical and experimental results. Totally 64 non-linear time-history analyses were conducted, and the responses were investigated in terms of the lateral displacement of the pile, the effect of the base acceleration in the pile behavior, the bending moment distribution in the pile body, and the pore pressure. The numerical analysis results demonstrated that the relationship between the pile lateral displacement and the maximum base acceleration is non-linear. Furthermore, increasing the pile diameter results in an increase in the passive pressure of the soil. Also, piles with small and big diameters are subjected to yielding under bending and shear states, respectively. It is concluded that an effective stress-based ground response analysis should be conducted when there is a liquefaction condition in order to determine the maximum bending moment and shear force generated within the pile.