• 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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.

• The KIPS Transactions:PartA
• /
• v.12A no.2 s.92
• /
• pp.95-102
• /
• 2005

The Newton-Raphson iterative algorithm for finding a floating point reciprocal which is widely used for a floating point division, calculates the reciprocal by performing a fixed number of multiplications. In this paper, a variable latency Newton-Raphson's reciprocal algorithm is proposed that performs multiplications a variable number of times until the error becomes smaller than a given value. To find the reciprocal of a floating point number F, the algorithm repeats the following operations: '$'X_{i+1}=X=X_i*(2-e_r-F*X_i),\;i\in\{0,\;1,\;2,...n-1\}'$ with the initial value $'X_0=\frac{1}{F}{\pm}e_0'$. The bits to the right of p fractional bits in intermediate multiplication results are truncated, and this truncation error is less than $'e_r=2^{-p}'$. The value of p is 27 for the single precision floating point, and 57 for the double precision floating point. Let $'X_i=\frac{1}{F}+e_i{'}$, these is $'X_{i+1}=\frac{1}{F}-e_{i+1},\;where\;{'}e_{i+1}, is less than the smallest number which is representable by floating point number. So, $X_{i+1}$ is approximate to $'\frac{1}{F}{'}$. Since the number of multiplications performed by the proposed algorithm is dependent on the input values, the average number of multiplications per an operation is derived from many reciprocal tables $(X_0=\frac{1}{F}{\pm}e_0)$ with varying sizes. The superiority of this algorithm is proved by comparing this average number with the fixed number of multiplications of the conventional algorithm. Since the proposed algorithm only performs the multiplications until the error gets smaller than a given value, it can be used to improve the performance of a reciprocal unit. Also, it can be used to construct optimized approximate reciprocal tables. The results of this paper can be applied to many areas that utilize floating point numbers, such as digital signal processing, computer graphics, multimedia scientific computing, etc.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.10
• /
• pp.1427-1435
• /
• 2010

A mixed finite element model was developed using the classical plate theory to analyze the nonlinear bending of a plate. The appropriate weight functions for the constraints integrated over the domain were determined by the Lagrange multiplier method by using the principle of minimum virtual energy; which provides the constitutive relations between force-like variables and strains. All of detail terms of element wise coefficient matrices and associate tangent matrices to be used in the Newton iterative method are presented. Then, the linear solutions of the current model and those of the traditional displacement model under the SS (simple support) boundary conditions were compared with the existing analytical solution. The post-processed images of the nonlinear results of the force-like variables are presented to show the continuity of the solutions at the joint of the element boundaries. Finally, the converged nonlinear finite element solutions of the current model are compared with those of existing traditional displacement model.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• The Korean Journal of Applied Statistics
• /
• v.30 no.4
• /
• pp.539-553
• /
• 2017

We propose a multi-state model to analyze semi-competing risks data with interval-censored or missing intermediate events. This model is an extension of the three states of the illness-death model: healthy, disease, and dead. The 'diseased' state can be considered as the intermediate event. Two more states are added into the illness-death model to incorporate the missing events, which are caused by a loss of follow-up before the end of a study. One of them is a state of the lost-to-follow-up (LTF), and the other is an unobservable state that represents an intermediate event experienced after the occurrence of LTF. Given covariates, we employ the Lin and Ying additive hazards model with log-normal frailty and construct a conditional likelihood to estimate transition intensities between states in the multi-state model. A marginalization of the full likelihood is completed using adaptive importance sampling, and the optimal solution of the regression parameters is achieved through an iterative quasi-Newton algorithm. Simulation studies are performed to investigate the finite-sample performance of the proposed estimation method in terms of empirical coverage probability of true regression parameters. Our proposed method is also illustrated with a dataset adapted from Helmer et al. (2001).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.329-334
• /
• 2017

The structures, which are made up with the huge number of degrees-of-freedom and the assembly of substructures, have a great complexity. In order to increase the computational efficiency, the analysis models have to be simplified. Many substructuring techniques have been developed to simplify large-scale engineering problems. The techniques are very powerful for solving nonlinear problems which require many iterative calculations. In this paper, a modal derivatives-based model order reduction method, which is able to capture the stretching-bending coupling behavior in geometrically nonlinear systems, is adopted and investigated for its performance evaluation. The quadratic terms in nonlinear beam theory, such as Green-Lagrange strains, can be explained by the modal derivatives. They can be obtained by taking the modal directional derivatives of eigenmodes and form the second order terms of modal reduction basis. The method proposed is then applied to a co-rotational finite element formulation that is well-suited for geometrically nonlinear problems. Numerical results reveal that the end-shortening effect is very important, in which a conventional modal reduction method does not work unless the full model is used. It is demonstrated that the modal derivative approach yields the best compromised result and is very promising for substructuring large-scale geometrically nonlinear problems.