• Title/Summary/Keyword: Newton-Raphson procedure

Search Result 56, Processing Time 0.026 seconds

Regression models for interval-censored semi-competing risks data with missing intermediate transition status (중간 사건이 결측되었거나 구간 중도절단된 준 경쟁 위험 자료에 대한 회귀모형)

  • Kim, Jinheum;Kim, Jayoun
    • The Korean Journal of Applied Statistics
    • /
    • v.29 no.7
    • /
    • pp.1311-1327
    • /
    • 2016
  • We propose a multi-state model for analyzing semi-competing risks data with interval-censored or missing intermediate events. This model is an extension of the 'illness-death model', which composes three states, such as 'healthy', 'diseased', and 'dead'. The state of 'diseased' can be considered as an intermediate event. Two more states are added into the illness-death model to describe missing events caused by a loss of follow-up before the end of the study. One of them is a state of 'LTF', representing a lost-to-follow-up, and the other is an unobservable state that represents the intermediate event experienced after LTF occurred. Given covariates, we employ the Cox proportional hazards model with a normal frailty and construct a full likelihood to estimate transition intensities between states in the multi-state model. Marginalization of the full likelihood is completed using the adaptive Gaussian quadrature, and the optimal solution of the regression parameters is achieved through the iterative Newton-Raphson algorithm. Simulation studies are carried out to investigate the finite-sample performance of the proposed estimation procedure in terms of the empirical coverage probability of the true regression parameter. Our proposed method is also illustrated with the dataset adapted from Helmer et al. (2001).

A Study on Buckling Behavior of Shallow Circular Arches (낮은 원호아치의 좌굴거동에 대한 연구)

  • 김연태;허택녕;오순택
    • Journal of the Earthquake Engineering Society of Korea
    • /
    • v.2 no.2
    • /
    • pp.87-94
    • /
    • 1998
  • Behavioral characteristics of shallow circular arches with dynamic loading and different end conditions are analysed. Geometric nonlinearity is modelled using Lagrangian description of the motion. The finite element analysis procedure is used to solve the dynamic equation of motion, and the Newmark method is adopted in the approximation of time integration. The behavior of arches is analysed using the buckling criterion and non-dimensional time, load and shape parameters which Humphreys suggested. But a new deflection-ratio formula including the effect of horizontal displacement plus vertical displacement is presented to apply for the non-symmetric buckling problems. Through the model analysis, it's confirmed that fix-ended arches have higher buckling stability than hinge-ended arches, and arches with the same shape parameter have the same deflection ratio at the same time parameter when loaded with the same parametric load.

  • PDF

Geometrically nonlinear dynamic analysis of FG graphene platelets-reinforced nanocomposite cylinder: MLPG method based on a modified nonlinear micromechanical model

  • Rad, Mohammad Hossein Ghadiri;Shahabian, Farzad;Hosseini, Seyed Mahmoud
    • Steel and Composite Structures
    • /
    • v.35 no.1
    • /
    • pp.77-92
    • /
    • 2020
  • The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

Analysis of Pumping Test Data and The Prediction of Drawdown for Daejong-Chun Area (대종천유역 충적대수층의 수리성 분석 및 수위강하예측에 관한 연구)

  • Choi, Jae-Jin;Sung, Won-Mo;Hahn, Jeong-Sang
    • Economic and Environmental Geology
    • /
    • v.26 no.4
    • /
    • pp.541-549
    • /
    • 1993
  • The main goal of this paper is to determine hydraulic properties and to predict drawdown for the efficient and stable development of groundwater in the Daejong-Chun area, North of Kyungsang-Do. Based on geological survey and analysis of well logging data conducted in 1991, it is found that the type of aquifer of this area is considered to be an anisotropic unconfined aquifer with saturated thickness of 19.8 m. In order to characterize this aquifer pumping test was conducted, and the resulting drawdown data were utilized for the analysis by applying both type curve matching technique and semi-log straight line method. As a result, the average specific yield of this aquifer is estimated as 32.3%, and the average ratio of $K_H$ to $K_V$ is only 2.7, which means that gravitational effect is not significant factor for this type of aquifer. For the validation of the estimated hydraulic properties, the analytical model which was developed with Newton-Raphson iteration procedure in this study, was employed to generate the drawdown. And, the resulting drawdown was compared against actual drawdown data and it shows the excellent matches. The actual drawdown data for 9 hours of pumping were used for history matching purposes and relatively satisfactory matches were achieved in this match. Then, the model was run by using the tuned parameters that are obtained during history matching stage, and the drawdown was predicted for the next 30 years of pumping with $3,000m^3/day$ of constant pumping rate. Its result indicates that the drawdown was stabilized as 1.41 m from 20 days with $3,000m^3/day$ of constant pumping rate, which is the required amount of water to be safely supplied to this area.

  • PDF

Large deflection behavior and stability of slender bars under self weight

  • Goncalves, Paulo B.;Jurjo, Daniel Leonardo B.R.;Magluta, Carlos;Roitman, Ney;Pamplona, Djenane
    • Structural Engineering and Mechanics
    • /
    • v.24 no.6
    • /
    • pp.709-725
    • /
    • 2006
  • In this paper the buckling and post-buckling behavior of slender bars under self-weight are studied. In order to study the post-buckling behavior of the bar, a geometrically exact formulation for the non-linear analysis of uni-directional structural elements is presented, considering arbitrary load distribution and boundary conditions. From this formulation one obtains a set of first-order coupled nonlinear equations which, together with the boundary conditions at the bar ends, form a two-point boundary value problem. This problem is solved by the simultaneous use of the Runge-Kutta integration scheme and the Newton-Raphson method. By virtue of a continuation algorithm, accurate solutions can be obtained for a variety of stability problems exhibiting either limit point or bifurcational-type buckling. Using this formulation, a detailed parametric analysis is conducted in order to study the buckling and post-buckling behavior of slender bars under self-weight, including the influence of boundary conditions on the stability and large deflection behavior of the bar. In order to evaluate the quality and accuracy of the results, an experimental analysis was conducted considering a clamped-free thin-walled metal bar. As this kind of structure presents a high index of slenderness, its answers could be affected by the introduction of conventional sensors. In this paper, an experimental methodology was developed, allowing the measurement of static or dynamic displacements without making contact with the structure, using digital image processing techniques. The proposed experimental procedure can be used to a wide class of problems involving large deflections and deformations. The experimental buckling and post-buckling behavior compared favorably with the theoretical and numerical results.

Performance-based structural fire design of steel frames using conventional computer software

  • Chan, Y.K.;Iu, C.K.;Chan, S.L.;Albermani, F.G.
    • Steel and Composite Structures
    • /
    • v.10 no.3
    • /
    • pp.207-222
    • /
    • 2010
  • Fire incident in buildings is common, so the fire safety design of the framed structure is imperative, especially for the unprotected or partly protected bare steel frames. However, software for structural fire analysis is not widely available. As a result, the performance-based structural fire design is urged on the basis of using user-friendly and conventional nonlinear computer analysis programs so that engineers do not need to acquire new structural analysis software for structural fire analysis and design. The tool is desired to have the capacity of simulating the different fire scenarios and associated detrimental effects efficiently, which includes second-order P-D and P-d effects and material yielding. Also the nonlinear behaviour of large-scale structure becomes complicated when under fire, and thus its simulation relies on an efficient and effective numerical analysis to cope with intricate nonlinear effects due to fire. To this end, the present fire study utilizes a second-order elastic/plastic analysis software NIDA to predict structural behaviour of bare steel framed structures at elevated temperatures. This fire study considers thermal expansion and material degradation due to heating. Degradation of material strength with increasing temperature is included by a set of temperature-stress-strain curves according to BS5950 Part 8 mainly, which implicitly allows for creep deformation. This finite element stiffness formulation of beam-column elements is derived from the fifth-order PEP element which facilitates the computer modeling by one member per element. The Newton-Raphson method is used in the nonlinear solution procedure in order to trace the nonlinear equilibrium path at specified elevated temperatures. Several numerical and experimental verifications of framed structures are presented and compared against solutions in literature. The proposed method permits engineers to adopt the performance-based structural fire analysis and design using typical second-order nonlinear structural analysis software.