• Title/Summary/Keyword: Newton-method

Search Result 1,017, Processing Time 0.026 seconds

Dynamic Modeling and Simulation of a Towing Rope using Multiple Finite Element Method (다물체 요소이론을 이용한 예인줄 동역학의 모델링 및 시뮬레이션)

  • Yoon, Hyeon-Kyu;Lee, Hong-Seok;Park, Jong-Kyu;Kim, Yeon-Gyu
    • Journal of Navigation and Port Research
    • /
    • v.36 no.5
    • /
    • pp.339-347
    • /
    • 2012
  • After towing rope connecting a barge to a tug was subdivided into multiple finite elements, then those dynamic models was established using Newton's second law and considering the external force and moment such as tension, drag, Coriolis force, gravity, buoyancy, and impact due to free surface acting on each element. While the previous research on the model of towing rope considered only translation, five-degree-of-freedom equations of motion except roll based on the body-fixed frame were established in this paper. All elements are connected by a spring and a damper, and the stiffness of the spring was set as the equivalent value of the real rope. In order to confirm the established multiple finite element model, various scenarios such as freely falling of towing rope in the air and above the free surface, accelerating of a tug which tows a barge connected by towing rope, and sinusoidal moving of a tug were set up and simulated. As the results, the trajectories of the tug, the barge, and the towing rope showed good tendencies to the ones of real expected situations.

Fracture and Hygrothermal Effects in Composite Materials (복합재의 파괴와 hygrothermal 효과에 관한 연구)

  • Kook-Chan Ahn;Nam-Kyung Kim
    • Journal of the Korean Society of Safety
    • /
    • v.11 no.4
    • /
    • pp.143-150
    • /
    • 1996
  • This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

  • PDF

A Semi-Implicit Integration for Rate-Dependent Plasticity with Nonlinear Kinematic Hardening (비선형 이동경화를 고려한 점소성 모델의 내연적 적분)

  • Yoon, Sam-Son;Lee, Soon-Bok
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.27 no.9
    • /
    • pp.1562-1570
    • /
    • 2003
  • The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

Buckling Analysis of Axisymmetric Shells by Incremental Finite Element Mothod (증분형(增分形) 유한요소법(有限要素法)에 의한 축대칭(軸對稱) Shell구조(構造)의 좌굴해석(挫屈解析))

  • J.B.,Kim;C.Y.,Kim
    • 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.

  • PDF

Independent Component Analysis for Clustering Components by Using Fixed-Point Algorithm of Secant Method and Kurtosis (할선법의 고정점 알고리즘과 첨도에 의한 군집성의 독립성분분석)

  • Cho, Yong-Hyun
    • Journal of the Korean Institute of Intelligent Systems
    • /
    • v.14 no.3
    • /
    • pp.336-341
    • /
    • 2004
  • This paper proposes an independent component analysis(ICA) of the fixed-point (FP) algorithm based on secant method and the kurtosis. The FP algorithm based on secant method is applied to improve the analysis speed and performance by simplifying the calculation process of the complex derivative in Newton method, the kurtosis is applied to cluster the components. The proposed ICA has been applied to the problems for separating the 6-mixed signals of 500 samples and 8-mixed images of $512{\times}512$ pixels, respectively. The experimental results show that the proposed ICA has always a fixed analysis sequence. The result can be solved the limit of conventional ICA based on secant method which has a variable sequence depending on the running of algorithm. Especially, the proposed ICA can be used for classifying and identifying the signals or the images.

A Study on the Elasto-Plasticity Behaviour of a Ship's Plate under Thrust According to Boundary Condition (압축력을 받는 선체판의 경계조건에 따른 탄소성거동에 관한 연구)

  • Ko Jae-Yong;Park Joo-Shin;Park Sung-Hyun
    • Journal of the Korean Society of Marine Environment & Safety
    • /
    • v.10 no.1 s.20
    • /
    • pp.29-33
    • /
    • 2004
  • Design of general steel structure had applied to achieve elastic designing concept so far. Because elastic design supposes that whole structure complies with elasticity formula so that achieve via allowable stress of material. It is concept that calculate stress distribution of construction about action external load and estimate load when the maximum stress reaches equally with allowable stress that is established by maximum safety load of the structure. But, absence that compose actuality structure by deal with external load increase small success surrender and structure hardness falls and structure in limited state finally on the whole as showing complicated process by interference between collapse and buckling under compression. Applied ANSYS (elasto-plasticity large deformation finite element method) to be mediocrity finite element program for analysis method and analysis control used in Newton-Raphson method & Arc-length method.

  • PDF

Depth Scaling Strategy Using a Flexible Damping Factor forFrequency-Domain Elastic Full Waveform Inversion

  • Oh, Ju-Won;Kim, Shin-Woong;Min, Dong-Joo;Moon, Seok-Joon;Hwang, Jong-Ha
    • Journal of the Korean earth science society
    • /
    • v.37 no.5
    • /
    • pp.277-285
    • /
    • 2016
  • We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

Development of Compressive Ultimate Strength Formulations for Ship Plating Stiffener with Cutout (선체 유공보강판의 압축최종강도에 관한 설계식 개발)

  • Ko Jae-Yong;Park Joo-Shin;Oh Dong-Kee
    • Proceedings of KOSOMES biannual meeting
    • /
    • 2004.05b
    • /
    • pp.121-125
    • /
    • 2004
  • Recently, High Tensile Steel is adapt to thin plate on the steel structure and marine structure is used widely, It is possible for buckling happens great. Specially, Initial deflection of ship structure happens in place absence necessarily by heat processing of welding or cutting etc. This Initial Deflection is exerted negative impact when thin plate absence complicated nonlinear behaviour accompanied secondary budding. As a result, must idealize initial deflection that occurrence is possible to endow stability and accuracy in the hull structure or marine structure and reflect in early structure design considering secondary buckling. Longi direction of compressive load interacts and analyzed finite element series analysis that apply various kinds initial deflection shape measured actually on occasion that is arranged simply supported condition in this research. Applied ANSYS (elasto-plasticity large deformation finite element method) to be mediocrity finite element program for analysis method and analysis control used in Newton-Raphson method & Arc-length method.

  • PDF

Enhancement of Image Reconstruction Using Region of Interest Method Based on Adaptive Threshold Value in Electrical Impedance Tomography (전기 임피던스 단층촬영법에서 적응 문턱치 기반의 관심영역 기법을 사용한 영상 복원의 개선)

  • Kim, Chang Il;Kim, Bong Seok;Kim, Kyung Youn
    • Journal of the Institute of Electronics and Information Engineers
    • /
    • v.54 no.8
    • /
    • pp.99-106
    • /
    • 2017
  • Electrical impedance tomography is a nondestructive imaging modality in which the internal resistivity distribution is reconstructed based on the injected currents and measured voltages inside a domain of interest. In this paper, an adaptive threshold value based region of interest (ROI) method is proposed to improve the spatial resolution of reconstructed images as well as to reduce the computational time of the inverse problem. Adaptive threshold value is calculated by INTERMODES method and ROI is determined from the domain based on this value. Moreover, the computational domain of image reconstruction is restricted within a ROI and iterative Gauss-Newton method is employed to estimate the resistivity distribution. To evaluate the performance of the proposed method, numerical experiments have been performed and the results are analyzed.

STK Feature Tracking Using BMA for Fast Feature Displacement Convergence (빠른 피쳐변위수렴을 위한 BMA을 이용한 STK 피쳐 추적)

  • Jin, Kyung-Chan;Cho, Jin-Ho
    • Journal of the Korean Institute of Telematics and Electronics S
    • /
    • v.36S no.8
    • /
    • pp.81-87
    • /
    • 1999
  • In general, feature detection and tracking algorithms is classified by EBGM using Garbor-jet, NNC-R and STK algorithm using pixel eigenvalue. In those algorithms, EBGM and NCC-R detect features with feature model, but STK algorithm has a characteristics of an automatic feature selection. In this paper, to solve the initial problem of NR tracking in STK algorithm, we detected features using STK algorithm in modelled feature region and tracked features with NR method. In tracking, to improve the tracking accuracy for features by NR method, we proposed BMA-NR method. We evaluated that BMA-NR method was superior to NBMA-NR in that feature tracking accuracy, since BMA-NR method was able to solve the local minimum problem due to search window size of NR.

  • PDF