• Steel and Composite Structures
• /
• 제48권6호
• /
• pp.611-623
• /
• 2023

This work considered an optimal control formulation in the sense of Caputo derivatives. The optimality of the fractional optimal control problem. The tumor immune interaction in fractional form provides an excellent tool for the description of memory and hereditary properties of inter and intra cells. So the interaction between effector-cells, tumor cells and are modeled by using the definition of Caputo fractional order derivative that provides the system with long-time memory and gives extra degree of freedom. In addiltion, existence and local stability of fixed points are investigated for discrete model. Moreover, in order to achieve more efficient computational results of fractional-order system, a discretization process is performed to obtain its discrete counterpart. Our technique likewise allows the advancement of results, such as return time to baseline that are unrealistic with current model solvers.

• Geomechanics and Engineering
• /
• 제12권4호
• /
• pp.723-738
• /
• 2017

Rock is a heterogeneous material, which introduces complexity in the analysis of rock slopes, since both the existing discontinuities within the rock mass and the intact rock contribute to the degradation of strength. Rock failure is often catastrophic due to the brittle nature of the material, involving the sliding along structural planes and the fracturing of rock bridge. This paper proposes an advanced discretization method of rock mass based on block theory. An in-house software, GeoSMA-3D, has been developed to generate the discrete fracture network (DFN) model, considering both measured and artificial joints. Measured joints are obtained from the photogrammetry analysis on the excavation face. Statistical tools then facilitate to derive artificial joints within the rock mass. Key blocks are searched to provide guidance on potential reinforcement measures. The discretized blocky system is subsequently implemented into a discontinuous deformation analysis (DDA) code. Strength reduction technique is employed to analyze the stability of the slope, where the factor of safety can be obtained once excessive deformation of slope profile is observed. The combined analysis approach also provides the failure mode, which can be used to guide the choice of strengthening strategy if needed. Finally, an illustrated example is presented for the analysis of a rock slope of 20 m height inclined at $60^{\circ}$ using combined GeoSMA-3D and DDA calculation.

• 한국안전학회지
• /
• 제21권2호
• /
• pp.143-149
• /
• 2006

본 논문은 2차원 동적 진동문제를 공간-시간 유한요소법으로 해석하고 있다. 공간-시간 유한요소법은 공간만 분할하는 재래식 유한요소해석에 비해 보다 해를 빠르고 쉽게 얻을 수 있다. 상대적으로 큰 시간간격에 대해서 공간과 시간을 동시에 분할하는 공간-시간 유한요소 근사법을 제시한다. 가중잔차법으로 공간-시간 영역에 대해 유한요소법을 정식화하였으며 선형 사변형 공간-시간 유한요소를 선택하여 해의 안정성에 관하여 언급하였다. 일반적 동적문제에서는 상대적인 큰 시간간격으로 인하여 해의 불안정을 야기 시키고 있으나 본 연구에서는 수치의 안정성을 보여주고 있다. 비구조 공간-시간 유한요소법은 재래식 수치해석에서 흔히 발생하는 해의 불안정성에 대한 결점을 보완함은 물론 효과적인 계산방법을 지니고 있다. 이 방법의 효율성을 위해 수치예제들을 제시하였다.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2004년도 추계학술대회논문집
• /
• pp.75-78
• /
• 2004

A mesh generation algorithm adapted to the mesh density map using the Delaunay mesh generation technique is developed. In the finite element analyses of the forging processes, the numerical error increases as the process goes on because of discrete property of the finite elements or severe distortion of elements. Especially, in the region where stresses and strains are concentrated, the numerical discretization error will be highly increased. However, it is too time consuming to use a uniformly fine mesh in the whole domain to reduce the expected numerical error. Therefore, it is necessary to construct locally refined mesh at the region where the error is concentrated such as at the die corner. In this study, the point insertion algorithm is used and the mesh size is controlled by moving nodes to optimized positions according to a mesh density map constructed with a posteriori error estimation. An optimization technique is adopted to obtain a good position of nodes. And optimized smoothing techniques are also adopted to have smooth distribution of the mesh and improve the mesh element quality.

• Wind and Structures
• /
• 제14권5호
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.

• Wind and Structures
• /
• 제22권3호
• /
• pp.349-368
• /
• 2016

A preconditioning technique is presented for a simultaneous solution to wind-membrane interaction. In the simultaneous equations, a linear elastic model was employed to deal with the fluid-structure data transfer at the interface. A Lagrange multiplier was introduced to impose the specified boundary conditions at the interface and strongly coupled simultaneous equations are derived after space and time discretization. An initial linear elastic model preconditioner and modified one were derived by treating the linearized elastic model equation as a saddle point problem, respectively. Accordingly, initial and modified fluid-structure interaction (FSI) preconditioner for the simultaneous equations were derived based on the initial and modified linear elastic model preconditioners, respectively. Wind-membrane interaction analysis by the proposed preconditioners, for two and three dimensional membranous structures respectively, was performed. Comparison was made between the performance of initial and modified preconditioners by comparing parameters such as iteration numbers, relative residuals and convergence in FSI computation. The results show that the proposed preconditioning technique greatly improves calculation accuracy and efficiency. The priority of the modified FSI preconditioner is verified. The proposed preconditioning technique provides an efficient solution procedure and paves the way for practical application of simultaneous solution for wind-structure interaction computation.

• Structural Engineering and Mechanics
• /
• 제45권4호
• /
• pp.521-542
• /
• 2013

In this study a truss model is used for the geometrically nonlinear static and dynamic analysis of a thin shallow arch subject to snap-through. Thanks to the very simple geometry of a truss, the equilibrium conditions can be easily written and the global stiffness matrix can be easily updated with respect to the deformed structure, within each step of the analysis. A very coarse discretization is applied; so, in a very simple way, the high frequency modes are suppressed from the beginning and there is no need to develop a complicated reduced-order technique. Two short computer programs have been developed for the geometrically nonlinear static analysis by displacement control of a plane truss model of a structure as well as for its dynamic analysis by the step-by-step time integration algorithm of trapezoidal rule, combined with a predictor-corrector technique. These two short, fully documented computer programs are applied on the geometrically nonlinear static and dynamic analysis of a specific thin shallow arch subject to snap-through.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1997년도 하계학술대회 논문집 B
• /
• pp.600-602
• /
• 1997

The kalman filter which has good estimating capabilities by means of the recursive computation from the previously known of obtained data is usually used for the system estimation in the case of not being directly measurable. The best estimating technique is still open issues on the PWR reactor control system to increase operating contingencies and to predict the safety margins for safer reactor operation. This paper addressed its estimating technique using kalman filter for the more flexible reactor control and showed the reasonable approach for discretization of the continuos-time system for reduction of computation errors.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2000년도 추계학술대회논문집
• /
• pp.674-679
• /
• 2000

This research is concerned with the modeling technique and active vibration controller design for slewing smart structures. When cantilever beam rotates about axes perpendicular to the undeformed beam's longitudinal axis, it experiences inertial loading. Hence, the beam vibrates from the initial stage of slewing, In this paper, the analytical model for a single slewing flexible beam with surface bonded piezoelectric sensor and actuator is developed using the Hamilton's principle with discretization by the assumed mode method. It is found from experiments that the theoretical model lacks the frictional effect. The frictional effect is incorporated into the equations of motion by employing the coupling factor. Theoretical and experimental results show problems arising in modeling and controller design.

• 전산구조공학
• /
• 제5권1호
• /
• pp.75-82
• /
• 1992

h-version 유한요소에서 평활 곡선경계는 충분한 갯수의 직선경계에 의해 근사될 수 있다. 그러나, 일반적으로 곡선경계가 충분하지 않은 갯수의 직선변을 갖는 다각형요소, 또는 곡선요소등에 의한 사상이 정확하지 않을 경우 해가 수렴되지 않을 뿐만아니라 특히, 곡면에 수직방향의 응력은 다른 방향의 응력요소에 비해 수렴속도가 늦거나 틀린 해를 보여준다. 한편, p-version 유한요소는 사용되는 요소의 크기가 클 뿐아니라 변형되는 정도가 크므로 이러한 이산오차를 피하기 위해 초유한 보간기법에 제안되어 정확한 사상을 하게 된다. 본 연구에서는 직선경계는 물론 곡선경계에 초유한 사상을 h-version과 p-version에 적용하는 방법과 이에 필요한 초유한 보간자를 유도하여 세 문제의 예제를 통해 그 적용성과 우월성을 보이고자 한다.