• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1992.10a
• /
• pp.39-44
• /
• 1992

This paper is concerned with formulations of the hierarchical $C^{o}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element is that it is still difficult to construct elements based on h-version concepts which are accurate and stable against the shear locking effects. An adaptive mesh refinement and selective p-distribution of the polynomial degree using hp-version of the finite element method we proposed to verify the superior convergence and algorithmic efficiency with the help of the clamped L-shaped plate problems.s.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.118-127
• /
• 2010

Resistivity logging instruments are designed to measure the electrical resistivity of a formation, and this can be directly interpreted to provide a water-saturation profile. However, resistivity logs are sensitive to borehole and shoulder-bed effects, which often result in misinterpretation of the results. These effects are emphasised more in the presence of tool eccentricity. For precise interpretation of short- and long-normal logging measurements in the presence of tool eccentricity, we simulate and analyse eccentricity effects by combining the use of a Fourier series expansion in a new system of coordinates with a 2D goal-oriented high-order self-adaptive hp finite-element refinement strategy, where h denotes the element size and p the polynomial order of approximation within each element. The algorithm automatically performs local mesh refinement to construct an optimal grid for the problem under consideration. In addition, the proper combination of h and p refinements produces highly accurate simulations even in the presence of high electrical resistivity contrasts. Numerical results demonstrate that our algorithm provides highly accurate and reliable simulation results. Eccentricity effects are more noticeable when the borehole is large or resistive, or when the formation is highly conductive.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.95-105
• /
• 2009

The flow fields around the HARTII rotor were numerically investigated using a viscous flow solver on adaptive unstructured meshes. An overset mesh and a deforming mesh technique were used to handle the blade motion including blade deflection, which was obtain from the HARTII experimental data. A solution-adaptive mesh refinement technique was also used to capture the rotor wake effectively. Comparison of the sectional normal force and pitching moment at 87% radial station between the two cases, with and without the blade deflection, showed that the blade loading is significantly affected by blade torsion. It was found that as the mesh was refined, the strength of tip vortex is better preserved, and the magnitude of high frequency blade loading, caused by blade-vortex interaction (BVI), is further magnified. It was also found that a proper time step size, which corresponds to the cell size, should be used to predict unsteady solutions accurately. In general, the numerical results in terms of the unsteady blade loading and the rotor wake show good agreement with the experimental data.

• Steel and Composite Structures
• /
• v.26 no.2
• /
• pp.171-182
• /
• 2018

This paper presents an isogeometric discretization of Kirchhoff-Love thin shells using truncated hierarchical B-splines (THB-splines). It is demonstrated that the underlying basis functions are ideally appropriate for adaptive refinement of the so-called thin shell structures in the framework of isogeometric analysis. The proposed approach provides sufficient flexibility for refining basis functions independent of their order. The main advantage of local THB-spline evaluation is that it provides higher degree analysis on tight meshes of arbitrary geometry which makes it well suited for discretizing the Kirchhoff-Love shell formulation. Numerical results show the versatility and high accuracy of the present method. This study is a part of the efforts by the authors to bridge the gap between CAD and CAE.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.151-160
• /
• 1993

This paper is concerned with formulations of the hierarchical $C^{\circ}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element based on Integrals of Legendre shape functions is that it is still difficult to construct elements based on h-version concepts which are accurate and stable against the shear locking effects. An adaptive mesh refinement and selective p-distribution of the polynomial degree using hp-version of the finite element method are proposed to verify the superior convergence and algorithmic efficiency with the help of the simply supported L-shaped plate problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.49-56
• /
• 2024

An automated dynamic structural analysis module stands as a crucial element within a structural integrated mitigation system. This module must deliver prompt real-time responses to enable timely actions, such as evacuation or warnings, in response to the severity posed by the structural system. The finite element method, a widely adopted approximate structural analysis approach globally, owes its popularity in part to its user-friendly nature. However, the computational efficiency and accuracy of results depend on the user-provided finite element mesh, with the number of elements and their quality playing pivotal roles. This paper introduces a computationally efficient adaptive mesh generation scheme that optimally combines the h-method of node movement and the r-method of element division for mesh refinement. Adaptive mesh generation schemes automatically create finite element meshes, and in this case, representative strain values for a given mesh are employed for error estimates. When applied to dynamic problems analyzed in the time domain, meshes need to be modified at each time step, considering a few hundred or thousand steps. The algorithm's specifics are demonstrated through a standard cantilever beam example subjected to a concentrated load at the free end. Additionally, a portal frame example showcases the generation of various robust meshes. These examples illustrate the adaptive algorithm's capability to produce robust meshes, ensuring reasonable accuracy and efficient computing time. Moreover, the study highlights the potential for the scheme's effective application in complex structural dynamic problems, such as those subjected to seismic or erratic wind loads. It also emphasizes its suitability for general nonlinear analysis problems, establishing the versatility and reliability of the proposed adaptive mesh generation scheme.

• Economic and Environmental Geology
• /
• v.43 no.3
• /
• pp.259-267
• /
• 2010

Resistivity logging instruments were designed to measure electrical resistivity of formation, which can be directly interpreted to provide water-saturation profile. Short and long normal logging measurements are made under groundwater level. In some investigation sites, groundwater level reaches to a depth of a few meters. It has come to attention that the proximity of groundwater level might distort short and long normal logging readings, when the measurements are made near groundwater level, owing to the proximity of an insulating air. This study investigates the effects of the proximity of groundwater level (and also the proximity of earth surface) on the normal by simulating normal logging measurements near groundwater level. In the simulation, we consider all the details of real logging situation, i.e., the presence of wellbore, the tool mandrel with current and potential electrodes, and currentreturn and reference-potential electrodes. We also model the air to include the earth’'s surface in the simulation rather than the customary choice of imposing a boundary condition. To obtain apparent resistivity, we compute the voltage, i.e., potential difference between monitoring and reference electrodes. For the simulation, we use a twodimensional, goal-oriented and high-order self-adaptive hp finite element refinement strategy (h denotes the element size and p the polynomial order of approximation within each element) to obtain accurate simulation results. Numerical results indicate that distortion on the normal logging is greater when the reference potential electrode is closer to the borehole and distortions on long normal logging are larger than those on short normal logging.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.217-228
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of the p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity using the constitutive equation for work-hardening materials, and the associated flow rule. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the circular plate with uniformly distributed load. Those results are compared with the theoretical solutions and the numerical solutions of ADINA

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.2
• /
• pp.77-82
• /
• 2023

The structural analysis module is an essential part of any integrated structural system. Diverse integrated systems today require, from the analysis module, efficient real-time responses to real-time input such as earthquake signals, extreme weather-related forces, and man-made accidents. An integrated system may also be for the entire life span of a civil structure conceived during the initial conception, developed throughout various design stages, effectively used in construction, and utilized during usage and maintenance. All these integrated systems' essential part is the structural analysis module, which must be automated and computationally efficient so that responses may be almost immediate. The finite element method is often used for structural analysis, and for automation, many effective finite element meshes must be automatically generated for a given analysis. A computationally efficient finite element mesh generation scheme based on the r-h method of mesh refinement using strain deviations from the values at the Gauss points as error estimates from the previous mesh is described. Shape factors are used to sort out overly distorted elements. A standard cantilever beam analyzed by four-node plane stress elements is used as an example to show the effectiveness of the automated algorithm for a time-domain dynamic analysis. Although recent developments in computer hardware and software have made many new applications in integrated structural systems possible, structural analysis still needs to be executed efficiently in real-time. The algorithm applies to diverse integrated systems, including nonlinear analyses and general dynamic problems in earthquake engineering.