• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.173-182
• /
• 2014

The mixed convolved action principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation, convolution, and fractional calculus. In this paper, its potential in the development of numerical methods for transient problems in various dynamical systems when adopting temporally second order approximation is investigated. For this, the classical single-degree-of-freedom linear elastic dynamical systems are primarily considered to investigate computational characteristics of the developed algorithms. For the undamped system, all the developed algorithms are symplectic with respect to the time step. For the damped system, they are shown to be accurate with good convergence characteristics.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.3
• /
• pp.41-50
• /
• 1989

A basic formulation of the panel method, which is based on the potential field formulation, is reviewed for the case of two-dimensional hydrofoil problems. Numerical procedures to improve the computational efficiency of the panel method are suggested. By investigating local behavior of the flow around the trailing edge, a wedge type Kutta condition is formulated. By subdividing the trailing edge panels, where dipole strengths of the subdivided panels follow the local behavior of the potential values of the flow outside a wedge, the circulation around a hydrofoil is calculated accurately with a relatively small number of panels. The subdividing technique to improve the accuracy of the numerical Kutta condition is proved to be efficient. A local behavior of the flow around the leading edge is also investigated. By matching the flow around the leading edge with that around a parabola, a very accurate velocity distribution is obtained with relatively small number of panels. An accurate prediction of the stagnation point and the pressure distribution near the leading edge may contribute to improve the accuracy of cavity predictions and boundary layer calculations around hydrofoils.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.11 no.2
• /
• pp.15-22
• /
• 1974

As for the added mass moment of inertia of ships in torsional vibration, it seems that the works by T. Kumai[1,2] are the only systematic one available currently. The work[1] is for the calculation of the two dimensional correction factors with finitely-long elliptic cylinders as the mathematic model. In this work the authors recalculated the above factors, $J_{\tau}$, with the same mathematic model and the same problem formulation, and presented the numerical results in Fig. 1. The reason why the reinvestigation was done was that in Kumai's work he obtained the solutions of the Mathieu equations, which was derived from the problem formulation for the velocity potential, under the assumption that the dummy constant q involved in the equations was always far less than unity, whereas in fact it takes values within the region of $0<q{\leq}{\infty}$ in sequence. As a result the authors found two remarkable differences in general features of $J_{\tau}$(refer to Fg.3); one that the authors' numerical results are considerably higher than the results given in [2], and the other that for a given number of node those have properties of decreasing monotonically with increase of the beam-draft ratio while these rapidly decrease from a maximum value of near at B/T=2.00 with B/T becoming greater or less than ratio. It seems that the latter trend was resulted from the fact that the assumption of $q{\ll}1$ employed in [2] was more closely satisfied in the vicinity of B/T=2.00.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.6
• /
• pp.519-525
• /
• 2009

In this study a cohesive zone model was used to simulate the delamination phenomena which occurs by a successive crack initiation and propagation in composite laminates. The cohesive zone model was incorporated to the classical finite element method via cohesive element formulation and then implemented into the user-subroutine UEL of a commercial finite element program Abaqus. To validate the formulation and implementation of the cohesive element the finite element results were compared with the experimental data of double cantilever beam and end notched flexure tests. The numerical results well agree with the experimental load-displacement curves. Also the effect of the elastic stiffness and the size of the cohesive element on the global load-displacement curves were studied numerically. To minimize the mesh-dependency of the crack propagation path and eliminate the zig-zag patterns in the load-displacement curve, cohesive elements should be refined at the crack-tip.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.23 no.3
• /
• pp.1-1
• /
• 1987

The accurate prediction of the ship wave resistance is very important to design ships which operate satisfactorily in a wave environment. Thus, work should continue on development and validation of methods to compute ship wave patterns and wave resistance. Research efforts to improve the prediction of ship waves and wavemaking resistance are categorized in two major areas. First is the development of higher-order theories to take account of the nonlinear effect of the free surface condition and improved analytical treatment of the body boundary condition. Second is the development of direct numerical methods aimed at solving body and free-surface boundary conditions as accurately as possible. A new formulation of the slender body theory for a ship with constant speed is developed by Maruo. It is quite different from the existing slender ship theory by Vossers, Maruo and Tuck. It may be regarded as a substitute for the Neumann-Kelvin approximation. In present work, the method of asymptotic expansion of the Kelvin source is applied to obtain a new wave resistance formulation in fluid of finite depth. It takes a simple form than existing theory.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.27 no.4
• /
• pp.303-312
• /
• 1991

This paper describes the general formulation for the in-plane flexural free vibration analysis of two layered structure by the transfer influence coefficient method. The structure is regared as a distributed mass system with lumped mass and inertia moments, massless linear and rotational springs, and joints elements of releases and rolls at which the displacements are discontinuous in each layer. The results of the simple numerical examples on a personal computer demonstrate the validity of the present method, that is, the numerical high accuracy, the high speed, the flexibility for programming of the present algorithm, compared with the transfer matrix method.

• Steel and Composite Structures
• /
• v.43 no.2
• /
• pp.271-278
• /
• 2022

Buckling is one of the most critical phoneme in the design of steel structures. Lateral torsional buckling (LTB) is particularly significant for slender beams generally subjected to loading in plane. The web distortion effects on LTB are not addressed explicitly in standards for flexural design of steel I-section members. Hence, the present study is focused to predict the influence of the web distortion on the elastic (L_r) limiting lengths given in American Institute of Steel Construction (AISC) code for the lateral torsional buckling (LTB) behavior of steel beams due to no provision in the code for consideration of web distortion. For this aim, the W44x335 beam is adopted in the buckling analysis carried out by the ABAQUS finite element (FE) program since it is one of the most critical sections in terms of lateral torsional buckling (LTB). The strain results at mid-height of the web at mid-span of the beam are taken into account as the monitoring parameters. The web strain results are found to be relatively greater than the yield strain value when L/L_r is equal to 1.0. In other words, the ratio of L/L_r is estimated from the numerical analysis to be about 1.5 when the beam reaches its first yielding at mid-span of the beam at mid-height of the section. Due to the effect of web distortion, the elastic limiting length (L_r) from the numerical analysis is obtained to be considered as greater than the calculated length from the code formulation. It is suggested that the formulations of the limiting length proposed in the code can be corrected considering the influence of the web distortion. This correction can be a modification factor or a shape factor that reduces sectional slenderness for the LTB formulation in the code.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.437-454
• /
• 2022

Elliptic-braced simple resisting frame as a new lateral bracing system installed in the middle bay of frame in building facades has been recently introduced. This system not only creates a problem for opening space from the architectural viewpoint but also improves the structural behavior. Despite the researches on the seismic performance of lateral bracing systems, there are few studies performed on the effect of the stiffness parameters on the elastic story drift and calculation of period in simple braced steel frames. To overcome this shortcoming, in this paper, for the first time, an analytical solution is presented for calculating elastic lateral stiffness in a simple steel frame equipped with elliptic brace subjected to lateral load. In addition, for the first time, in this study, a precise formulation has been developed to evaluate the elastic stiffness variation in a steel frame equipped with a two-dimensional single-story single-span elliptic brace using strain energy and Castigliano's theorem. Thus, all the effective factors, including axial and shear loads as well as bending moments of elliptic brace could be considered. At the end of the analysis, the lateral stiffness can be calculated by an improved and innovative relation through the energy method based on the geometrical properties of the employed sections and specification of the used material. Also, an equivalent element of an elliptic brace was presented for the ease of modeling and use in linear designs. Application of the proposed relation have been verified through a variety of examples in OpenSees software. Based on the results, the error percentage between the elastic stiffness derived from the developed equations and the numerical analyses of finite element models was very low and negligible.

• Geomechanics and Engineering
• /
• v.36 no.2
• /
• pp.183-204
• /
• 2024

In current investigation, a novel beam finite element model is formulated to analyze the buckling and free vibration responses of functionally graded porous beams resting on Winkler-Pasternak elastic foundations. The novelty lies in the formulation of a simplified finite element model with only three degrees of freedom per node, integrating both C⁰ and C¹ continuity requirements according to Lagrange and Hermite interpolations, respectively, in isoparametric coordinate while emphasizing the impact of z-coordinate-dependent porosity on vibration and buckling responses. The proposed model has been validated and demonstrating high accuracy when compared to previously published solutions. A detailed parametric examination is performed, highlighting the influence of porosity distribution, foundation parameters, slenderness ratio, and boundary conditions. Unlike existing numerical techniques, the proposed element achieves a high rate of convergence with reduced computational complexity. Additionally, the model's adaptability to various mechanical problems and structural geometries is showcased through the numerical evaluation of elastic foundations, with results in strong agreement with the theoretical formulation. In light of the findings, porosity significantly affects the mechanical integrity of FGP beams on elastic foundations, with the advanced beam element offering a stable, efficient model for future research and this in-depth investigation enriches porous structure simulations in a field with limited current research, necessitating additional exploration and investigation.

• Structural Engineering and Mechanics
• /
• v.16 no.3
• /
• pp.259-274
• /
• 2003

In this paper, a novel numerical solution technique, the differential cubature method is employed to study the buckling problems of thick plates with arbitrary quadrilateral planforms and non-uniform boundary constraints based on the first order shear deformation theory. By using this method, the governing differential equations at each discrete point are transformed into sets of linear homogeneous algebraic equations. Boundary conditions are implemented through discrete grid points by constraining displacements, bending moments and rotations of the plate. Detailed formulation and implementation of this method are presented. The buckling parameters are calculated through solving a standard eigenvalue problem by subspace iterative method. Convergence and comparison studies are carried out to verify the reliability and accuracy of the numerical solutions. The applicability, efficiency, and simplicity of the present method are demonstrated through solving several sample plate buckling problems with various mixed boundary constraints. It is shown that the differential cubature method yields comparable numerical solutions with 2.77-times less degrees of freedom than the differential quadrature element method and 2-times less degrees of freedom than the energy method. Due to the lack of published solutions for buckling of thick rectangular plates with mixed edge conditions, the present solutions may serve as benchmark values for further studies in the future.