Structure Interaction Analysis Research Articles

Seismic Response of RC Frame Buildings: Theoretical Formulation and Modeling with Soil-Structure Interaction Analysis

"Seismic Response of RC Frame Buildings: Theoretical Formulation and Modeling with Soil- Structure Interaction Analysis," offers a comprehensive exploration of the seismic response of reinforced concrete (RC) framed structures, with particular emphasis on soil flexibility. Through a meticulous response spectrum analysis, the study delves into the influence of soil flexibility on buildings of varying heights, ranging from five to nine storeys, considering both fixed and flexible support conditions across diverse soil types. The research draws the following significant conclusions: Significance of Soil-Structure Interaction (SSI): Incorporating soil-structure interaction (SSI) into the analysis leads to a noteworthy extension of the natural time period of the structure. This parameter stands as a linchpin in grasping the lateral response of framed structures, underlining its pivotal role in the precision of seismic design calculations. Soil Flexibility and Amplified Base Shear: Enhanced soil flexibility is found to correlate with a proportional increase in base shear. This phenomenon is most pronounced in soft soil conditions, where base shear experiences substantial escalation, particularly in taller buildings. Roof Displacement as a Performance Indicator: Roof displacement, a crucial indicator of structural performance, exhibits an upward trajectory owing to the cumulative effects of SSI. Notably, soft soils manifest higher roof displacements, with these variations becoming more pronounced as the building's height increases. The insights garnered from this analysis underscore the profound influence of SSI on a building's seismic response. Furthermore, the study highlights that advancements in finite element methods (FEM) and computer technology have significantly facilitated the integration of SSI into structural analysis. Therefore, it is imperative to harness these advancements fully to enhance our comprehension of structural behavior and to promote the adoption of secure construction practices in seismic-prone regions.

Fluid–Structure Interaction Analysis of a Wind Turbine Blade with Passive Control by Bend–Twist Coupling

The idea of improving the energy output for small wind turbines without compromising the remaining aspects of the technology, such as costs and structural integrity, is very appealing in the context of the growing concern for global warming and the goal of providing electricity to remote and isolated regions. This work aims to contribute to the development of distributed wind generation by exploring the effects of bend–twist coupling on the performance of a wind turbine with a focus on a small rotor based on the NREL Phase VI blade geometry. After defining a structure in composite materials exhibiting the coupling behavior along with a reference counterpart, a comparative numerical analysis is performed using a Fluid–Structure Interaction (FSI) analysis. The main numerical framework is based on commercial software and consists of a finite-volume solver for fluid physics, a finite-element solver for solid physics, and a coupling interface for the interaction problem. The results, complemented with the predictions from a one-way analysis based on the blade-element momentum theory are used to define the increments in rotor torque. The analysis of the annual energy yield shows a 3% increase due to the bend–twist coupling used as a passive pitch mechanism, considering a Rayleigh distribution with an 11 m/s average wind speed. Simultaneously, the coupling causes increments of 0.2% and 0.3% for the blade root flapwise moment and the rotor thrust force, respectively, when considering parked conditions and a simplified extreme wind model.

Structural Analysis of Steel Truss and Masonry Bridge Interaction: A Case of Ali Fuat Pasha Bridge in Sakarya, Türkiye

Restoration of historical structures using new materials and techniques is widespread worldwide. In these applications, relatively new materials such as steel, concrete, reinforced concrete (RC), or composite are generally preferred. However, it is often ignored whether old materials and new materials work in compatibility. In this respect, Ali Fuat Pasha Bridge (or Bayezid II Bridge), which was built by Bayezid II in 1495 over the Sakarya River in Geyve-Sakarya, was examined. The bridge is 196.50 meters long and consists of 15 arches with different spans and three arches of the bridge were destroyed as a result of the earthquake. Then, the bridge has been restored by constructing a steel truss system in place of the destroyed arches. Within the scope of this study, the structural performance of the bridge, which is currently serving vehicle and pedestrian traffic, is examined by using finite element analyses (FEAs). Moreover, this research examines the combined behavior of steel and masonry materials and investigates the structural behavior of steel truss and masonry bridge interaction. According to the results of the analyses, there are significant behavioral differences between the masonry structure and the steel structure. The main cause for this disparity is thought to be the varying levels of stiffness and ductility in steel and masonry sections.

Analysis and Refinement of Host-Guest Interactions in Metal-Organic Frameworks.

ConspectusMetal-organic frameworks (MOFs) are a class of hybrid porous materials characterized by their periodic assembly using metal ions and organic ligands through coordination bonds. Their high crystallinity, extensive surface area, and adjustable pore sizes make them promising candidates for a wide array of applications. These include gas adsorption and separation, substrate binding, and catalysis, of relevance to tackling pressing global issues such as climate change, energy challenges, and pollution. In comparison to traditional porous materials such as zeolites and activated carbons, the design flexibility of organic ligands in MOFs, coupled with their orderly arrangement with associated metal centers, allows for the precise engineering of uniform pore environments. This unique feature enables a rich variety of interactions between the MOF host and adsorbed gas molecules, which are fundamental to understanding the observed uptake capacity and selectivity for target gas molecules and thus the overall performance of the material.In this Account, a data set for three-dimensional MOFs has been constructed based upon the structural analysis of host-guest interactions using the largest experimental database, the Cambridge Structural Database (CSD). A full screening was performed on structures with guest molecules of H2, C2H2, CO2, and SO2, and the relationship between the primary binding site, the isosteric heats of adsorption (Qst), and the adsorption uptake was extracted and established. We review the methodologies to refine host-guest interactions based primarily on our studies on the host-guest chemistry of MOFs. The methods include ligand functionalization, variation of metal centers, formation of defects, addition of single atom sites, and control of pore size and structure. In situ structural and dynamic investigations using diffraction and spectroscopic techniques are powerful tools to visualize the details of host-guest interactions upon the above modifications, affording key insights into functional performance at a molecular level. Finally, we give an outlook of future research priorities in the study of host-guest chemistry in MOF materials. We hope this Account will encourage the rational development and improvement of future MOF-based sorbents for applications in challenging gas adsorption, separations, and catalysis.

Interaction analysis of tobacco leaf microbial community structure and volatiles flavor compounds during cigar stacking fermentation.

Cigar stacking fermentation is a key step in tobacco aroma enhancement and miscellaneous gas reduction, which both have a great influence on increasing cigar flavor and improving industrial availability. To analyze the effect of cigar tobacco leaf (CTLs) microbial community on volatiles flavor compounds (VFCs), this study used multi-omics technology to reveal the changes in microbial community structure and VFCs of different cigar varieties during stacking fermentation, in addition to exploring the interaction mechanism of microbiome and VFCs. The results showed that the dominant microbial compositions of different CTL varieties during stacking fermentation were similar, which included Staphylococcus, Corynebacterium 1, Aerococcus, and Aspergillus. These dominant microbes mainly affected the microbial community structure and characteristic microorganisms of CTLs through microbial interactions, thereby influencing the transformation of VFCs. Characteristic microorganisms of different CTLs varieties such as Trichothecium, Trichosporon, Thioalkalicoccus and Jeotgalicoccus, were found to positively correlate with characteristic VFCs like megastigmatrienone 4, pyrazine, tetramethyl-, geranyl acetone, and 2-undecanone, 6,10-dimethyl-, respectively. This in turn affected the aroma and sensory quality of the CTLs. This study provides theoretical support for the analysis of the mechanism of microorganisms on VFCs and aroma, and development of microbial agents during cigar stacking fermentation.

Novel computational fluid dynamics-finite element analysis solution for the study of flexible material wave energy converters

The use of flexible materials for primary mover and power takeoff of wave energy converters (WECs) has attracted considerable attention in recent years, owing to their potential to enhance the reliability, survivability, and wave energy conversion efficiency. Although some reduced order models have been used to study the fluid–structure interaction (FSI) responses of flexible wave energy converters (fWECs), they are somehow inappropriate due to their limited accuracy and applicability span. To gain a deeper understanding of the physical mechanisms in fWECs, a high-fidelity approach is required. In this work, we build up a fluid–structure interaction analysis framework based on computational fluid dynamics and a finite element analysis method. The incompressible viscous flow is resolved by solving three-dimensional unsteady Navier–Stokes equations with a finite volume approach. The structure dynamics are solved by a finite element method, taking the nonlinear behavior of flexible material into consideration. A strong coupling strategy is utilized to enhance the numerical stability and convergence of the iterative process. We demonstrate the present FSI tool is able to provide rich flow field information and structural response details, such as the velocity, pressure, and structural stress distribution. This is illustrated through several case studies, including two types of fWECs. The unsteady wave–structure-interaction and the associated nonlinear phenomena are also accurately captured by this tool.

Fluid–Structure Coupling and Aerodynamic Performance of a Multi-Dimensional Morphing Wing with Flexible Metastructure Skin

A multi-dimensional morphing wing skeleton mechanism is proposed with double-sided triangular pyramid units, which can realize continuous variable span-wise bend, span-wise twist, and sweep. A lockable morphing unit is designed, and its mechanism/structure characteristics, degree of freedom, and the deformable function of its deformable wing skeleton mechanism are analyzed. One kind of flexible skin is proposed to meet the performance requirements, consisting of an internal metastructure and a flexible surface bonded on both sides. The morphing wing skeleton mechanism and the equivalent treated metastructure flexible skin are then combined. Subsequently, a two-way fluid–structure interaction analysis is conducted to investigate the influence of aerodynamic loads on the flexible skin and skeleton mechanism in different deformation states, including the influence of structural passive deformation on the aerodynamic characteristics of the morphing wing. The computational fluid dynamics method is employed to analyze the aerodynamic characteristics of the morphing wing in its initial state, as well as in three deformation states, and to study its aerodynamic performance in different flight environments.

PoseEdit: enhanced ligand binding mode communication by interactive 2D diagrams.

In this article, we present PoseEdit, a new, interactive frontend of the popular pose visualization tool PoseView. PoseEdit automatically produces high-quality 2D diagrams of intermolecular interactions in 3D binding sites calculated from ligands in complex with protein, DNA, and RNA. The PoseView diagrams have been improved in several aspects, most notably in their interactivity. Thanks to the easy-to-use 2D editor of PoseEdit, the diagrams are extensively editable and extendible by the user, can be merged with other diagrams, and even be created from scratch. A large variety of graphical objects in the diagram can be moved, rotated, selected and highlighted, mirrored, removed, or even newly added. Furthermore, PoseEdit enables a synchronized 2D-3D view of macromolecule-ligand complexes simplifying the analysis of structural features and interactions. The representation of individual diagram objects regarding their visualized chemical properties, like stereochemistry, and general graphical styles, like the color of interactions, can additionally be edited. The primary objective of PoseEdit is to support scientists with an enhanced way to communicate ligand binding mode information through graphical 2D representations optimized with the scientist's input in accordance with objective criteria and individual needs. PoseEdit is freely available on the ProteinsPlus web server ( https://proteins.plus ).

A Multi-Resolution MPS/FEM Coupling Method for Three-Dimensional Fluid–Structure Interaction Analysis

This work aims to propose an efficient MPS/FEM coupling method for the simulation of fluid–structure interaction (FSI), where the MPS and FEM are respectively employed to account for fluid flows and structural deformation. The main idea of our method is to develop a multi-scale multi-resolution MPS method for efficient fluid simulations in the context of MPS/FEM coupling. In the developed multi-scale MPS method, the fluid domain is discretized into particles of different resolutions before calculation, where particles close to the interest domain will be discretized into high resolution, while the rest are discretized into low resolution. A large particle interacting with small particles is divided into several small particles virtually, and weight functions are redefined to maintain the simulation stability. A bucket-sort-based algorithm is developed for the fast search of multi-resolution neighboring particles. The capacity of a newly proposed ghost cell boundary model is further enhanced, so as to accurately treat wall boundary problems with particles of different resolutions. On this basis, the multi-resolution MPS method is coupled with the FEM for FSI simulations. Finally, several numerical examples are conducted to demonstrate the accuracy and efficiency of the development method.

Implementation of a coupled FEM-SBFEM for soil-structure interaction analysis of large-scale 3D base-isolated nuclear structures

Accurately simulating the radiation damping of an infinite domain is a crucial issue for reasonably and effectively evaluating the isolation effect of three-dimensional (3D) base-isolated nuclear structures considering soil-structure interaction (SSI). The scaled boundary finite element method (SBFEM) has excellent advantages that enable it to precisely fulfil the radiation conditions. In this study, the finite element method (FEM) and SBFEM were combined to develop a coupled FEM-SBFEM technology. The structure and finite domain were discretized using FEM, and the infinite soil was realised using SBFEM. A complete and efficient calculation procedure for SSI analysis was developed based on the coupled FEM-SBFEM, improving the calculation efficiency through parallel programming. The accuracy of the developed procedure was verified using a numerical example. Furthermore, the dynamic response of a large-scale refined 3D base-isolated reactor building considering SSI was analysed using the proposed method. The results demonstrate that the 3D base isolation technology has excellent isolation performance. Additionally, we compared the response results of the proposed coupled FEM-SBFEM and direct method. The proposed procedure has broad engineering applicability and can be used for seismic SSI analysis of large-scale nuclear structures while maintaining high computational accuracy and stability.

Dynamic Effect of Fluid–Structure Interaction and Modal Analysis of Crude–Pipeline Interaction System Considering the Fluctuating Pressure of Medium

This work explored the fluid–structure interaction (FSI) effect of crude oil and pipeline and the vibration characteristics of the pipe–crude oil structural system. A finite element computational analysis theory and corresponding simulation method were proposed for crude oil pipelines considering the FSI effect. Meanwhile, the spring model of the FSI interface between crude oil and pipeline was proposed. The experimental simulation and neural network were adopted to obtain the flow-excited pulsation pressure (FEPP), and the kinetic equations of the pipe–crude oil FSI were analyzed. The finite element computation (FEC) model for pipe–crude oil bi-directional FSI dynamics was established, which was composed of five support span lengths and four pipe liquid levels and can provide five different crude oil fluid velocities. One hundred orthogonal simulation test conditions were designed. In addition, the displacement response and acceleration response under FSI of the pipe wall were calculated, and the first six orders of inherent frequencies and main vibration patterns of the FEC model were analyzed. The results revealed that the dynamic effect of FSI between the crude oil transport medium and pipeline was complex and was influenced by the length of the pipeline crossing section, medium flow velocity, difference of liquid surface height, and displacement response. The acceleration response possessed prominent vibration decay characteristics, and the response peak in the same condition was about 8.2 times larger than that of the steady-flow state response. Therefore, the effect changed the inherent frequency and vibration pattern of the pipeline, to which the engineering design should pay sufficient attention.

Hypersonic Fluid–Thermal–Structural Interaction of Cone–Slice–Ramp: Computations with Experimental Validation

The paper presents an investigation of fluid–thermal–structural interaction (FTSI) on a complex geometry in hypersonic flows. The objective is to assess high-fidelity numerical tools for FTSI by validating two numerical frameworks against wind-tunnel data. We choose a cone–slice–ramp geometry with a wedge angle of 30 deg subjected to a laminar inflow at Mach 8. First, the time-averaged and the time-resolved laminar flow solutions are validated and highlight the mean flow features and their unsteady characteristics, including low-frequency bubble breathing and shear-layer flapping at a higher frequency. The numerical analysis reveals the presence of upstream and downstream legs of streamwise vortices. Next, FTSI analysis is performed at a tractable computational cost by decomposing the numerical frameworks into a quasi-steady FTSI (QS-FTSI) analysis and a transient fluid–structure interaction (FSI) analysis based on the disparity of the characteristic timescales. The QS-FTSI results highlight the shear-layer reattachment as the dominant source heating the structure, with nonuniform spanwise variations due to the streamwise vortices. The temperature increase drives the decrease in the structural modal frequencies. The FSI analysis was performed for an unheated, undeformed structure and a nominal heated and deformed structure. The average error in the predicted FSI frequency of the unheated panel is 2.65% in comparison to the experimental data; however, this value increases to 7.42% for the heated panel. The frequency peaks are associated with the structural mode shapes via dynamic mode decomposition analysis and show a close resemblance to natural modes, indicating weak FSI coupling. The FTSI study concludes that for a cone–slice–ramp configuration under the present inflow conditions, the interaction mostly occurs in a one-way coupled manner, where the flowfield drives the thermal and structural response of the panel.

Monitoring of Tidal Variation and Temperature Change-Induced Movements of an Immersed Tunnel Using Distributed Optical Fiber Sensors (DOFSs)

The short-term deformation behavior of immersed tunnels due to daily or monthly temperature changes and tidal variations is often not monitored but forms important input for a structural health assessment of the tunnel. In this study, distributed optical fiber sensors (DOFSs) are used to monitor the short-term (daily and monthly) deformation behavior of an immersed tunnel. Joint opening and the relative settlement differences between tunnel elements are monitored simultaneously at subhour intervals. Measurements show that the variation in the joint opening is strongly correlated with temperature change, and the joint gap has a tendency to open at low temperatures and to close at increasing temperatures. Simultaneously, the entire immersed section behaves more like a rigid body and moves upwards and downwards periodically due to tidal fluctuations in the river, with an observed vertical movement of slightly less than one millimeter. The tide also causes local tilting of tunnel segments, and this tilting behavior differs between winter and summer, which implies that the (seasonal) temperature-induced joint deformations affect the robustness of the tunnel to tidal loads. A soil-tunnel structure interaction analysis reveals that the cyclic vertical movement of the tunnel is driven by retardation of the tidal wave in deeper soil layers, which can be captured by a coupled flow model. This study provides new insights into the short-term deformation behavior of immersed tunnels.

A temperature-driven approach for quantitative assessment of strengthening effect of continuous bridges using structural health monitoring data

The stiffness degeneration of small to medium span bridges has been increasingly observed in recent years, and it has become a major concern of government and bridge owners. A quantitative evaluation method for the bridge performance with strengthening measures is highly desired. Due to the advantages of uninterrupted traffic and long-term tracking capability, a temperature-driven approach for characterization of the correlation pattern between bridge temperature-induced strains and bridge status was proposed in the present study by using structural health monitoring data. First, a theoretical solution of the simplified bridge model was derived to establish the correlation between the stress and deterioration extent under temperature gradient load. After a numerical simulation that combines the thermal–structural interaction analysis and the vehicle–bridge interaction analysis, the strain range was proposed as an assessment index to ensure the stability and effectiveness of the evaluation results. Next, the Generalized Extreme Studentized Deviate method was used for detecting the outliers. The statistical results of the assessment index for different strengthening methods were compared to evaluate the associated strengthening efficiency, and the associated equivalent section height was calculated for visualizing the bridge condition after strengthening measures were taken. The results demonstrated that the proposed temperature-driven method was able to quantitatively evaluate the bridge strengthening effects with a high efficiency.

A suggested dynamic soil – structure interaction analysis

In this paper, a new methodology for time domain analysis of buildings on raft foundations considering soil – structure interaction is proposed. Sub-structuring technique is used to separate the building as super-structure and the underneath soil as sub-structure. The super-structure can be modeled using any numerical method. However, in this paper the super-structure is modeled via the BEM to consider the real interaction area between column and slabs. The dynamic load is considered as an earthquake acceleration record that can be transformed to equivalent dynamic loads acting on the super-structure floors. The sub-structure is analyzed using the dual reciprocity boundary element method as closed domain. New iterative coupling technique is proposed between the super and sub-structures to reduce the computational effort and required storage. An example is presented to demonstrate the strength and the practicality of the proposed methodology.

Loading–shape interaction analysis of membrane structures subjected to deformation dependent hydrostatic loadings

Due to their flexibility, membrane structures loaded and/or supported by gas and/or fluid undergo large deformations when subjected to external loads, for example when a relatively flat membrane deflects under accumulative ponding. The most noteworthy feature during the deformation process is the loading–shape interaction since changes in the form of the structure induce changes in the loading profiles, and vice versa. Description of the hydrostatic loading is complex since the discontinuity of structural domain and fluid domain occurs with the deformation of membrane. In this paper, an interactive model dealing with accurate hydrostatic load description and nonlinear finite element analysis is employed to analysis the structural response. The hydrostatic loading is described in vectors by fluid level and loading direction in current mesh of the membrane. In case to conserve a specific volume, the free surface height is calculated with a Newton algorithm. Equivalent nodal forces evaluation on different elements are deduced especially when they are partially or fully loaded by judging the spatial relation between the element and the fluid domain, and integrating with shape functions and area coordinates in a local coordinate system. In case of multiple fluid domains, the hydrostatic load description procedure is proposed to be conducted for each fluid domain with care in identifying respective loading direction and then sum them up to get the total hydrostatic loading forces. Two numerical examples, (a) single fluid domain case: a flat membrane subjected to accumulative water loading, and (b) multiple fluid domains case: an air-inflated dam partially filled with water and loaded by lateral pressures of water, are presented to demonstrate correctness, performance and effectiveness of the proposed method in large displacement analysis of flexible membrane structures subjected to deformation dependent hydrostatic loading.

Determination of Optimum Design Parameters for Gyroid Scaffolds to Mimic a Real Bone-Like Condition In Vitro: A Fluid Structure Interaction Study

Abstract A suitable scaffold architecture is always desirable to get a biomimetic scaffold for bone tissue engineering. In this regard, a fluid structure interaction analysis was carried out on different Micro-CTs (μCTs) and gyroids to observe the in vitro mechanical responses due to fluid flow. Computational fluid dynamics method was used to evaluate the permeability and wall shear stress (WSS), followed by a finite element method to obtain the mechanical stress within scaffolds. Different types of gyroids were designed based on the number of unit cells and porosity, where porosity of gyroids was kept same as μCTs. The main objective of the study is to examine the variations of permeability, WSS and mechanical stress with respect to the number of unit cells and porosity for different gyroids and μCTs. Mechanical responses were also compared between gyroids and μCTs. The results of this study highlighted that permeability and WSS of μCTs came close to the gyroids with eight unit cells but had significant differences in mechanical stress. The permeability of gyroids increased with the increase of porosity but decreased with the increase in number of unit cells. The opposite trend was shown in case of WSS within gyroids. This study will guide us in predicting an ideal scaffold for trabecular bone replacement.

Strongly coupled fluid–structure interaction analysis of aquatic flapping wings based on flexible multibody dynamics and the modified unsteady vortex lattice method

Because of their superior efficiency and low detectability compared to conventional submarine-shaped vehicles, bionic underwater vehicles such as Mantabot are playing increasingly important roles in ocean exploration. However, current research on manta-inspired robots is limited to structural design and pure hydrodynamics analysis of pectoral fin performance, with other aspects such as hydroelasticity and optimal design of flapping wings rarely addressed. This work presents a novel method for analyzing aquatic flapping wings by coupling flexible multibody dynamics with a modified unsteady vortex lattice method. The flexible multibody system is solved by using the absolute nodal coordinate formulation, and the conventional unsteady vortex lattice method is modified to accommodate dynamic stall. Compactly supported radial basis functions are used to transfer the hydroelastic forces and structural displacements across the interface meshes while satisfying global energy conservation, and a predictor–corrector method is used to stabilize the iteration procedure. Three different experiments are used to validate the flexible multibody dynamics solver, the modified unsteady vortex lattice solver, and the proposed fluid–structure interaction framework, respectively. Finally, the hydroelastic framework is demonstrated for a fin-like wing, and how flexibility and movement affect the flapping wing dynamics is examined. Numerical examples show that reducing the structural stiffness of the aquatic flapping wing within a certain range improves the propulsive efficiency but also reduces the thrust coefficient, and how the stiffness distribution affects the thrust and propulsive efficiency remains consistent over a range of the Strouhal number.

CFD-FEA based model to predict leak-points in a 90-degree pipe elbow

The aim of this paper is to numerically investigate Vibration-Based Leak Detection (VBLD) method in pipeline systems based on Fluid–Structure Interaction (FSI) analysis to predict leakages. In previous investigations, laboratory tests were widely used to study the VBLD technique in small-diameter water loop system pipes. The current project uses Ansys Workbench to extend these findings by integrating Computational Fluid Dynamics (CFD) with Finite Element Analysis (FEA). The study outlines a numerical method for VBLD to identify leakages in a 90-degree pipe elbow by predicting variations in vibration signals, with applications in the oil and gas industry. Firstly, changes in fluid behaviour (centrifugal force, pressure drop, secondary flow, and frictional force) experienced in the internal pipe wall resulting from a probable leakage (modelled as an additional outlet) are determined using CFD. Subsequently, the CFD results are coupled with FEA to model structural responses of the pipe walls subjected to different forces. This in turn allows the variations in vibration signals to be measured. The numerical approach presented in this paper based on FSI and incorporating the VBLD method provides a practical and convenient early detection tool that can complement physical vibration monitoring equipment in the field.

An efficient numerical algorithm for solving the dynamic impedance function of arbitrary-shaped foundations in layered soil

In soil–structure interaction analysis, the dynamic impedance solution can be obtained using Green's function based on frequency dependence. To satisfy the mixed boundary conditions of the ground–ground interface, the interface should be frequently divided into several evenly spaced sub-elements instead of finite elements. Generally, the accuracy increases and the computational efficiency decreases with an increase in the number of sub-elements. To obtain reliable solutions and improve efficiency, a numerical algorithm for solving the dynamic impedance function of arbitrary-shaped foundations in layered soil was investigated. Based on the wave propagation equation of a horizontally layered half space, the relationship between the number and radius of sub-disk elements and their dynamic impedance matrix was derived using the high-precision method of solving Green's function for the individual sub-disk elements. The derived equation provides an efficient and accurate method to calculate the dynamic impedance of rigid layered foundations and was verified based on numerical examples. The numerical examples revealed that the proposed method reproduced the results in the literature. With the given accuracy, this method reduces the calculation time and improves the calculation efficiency. This proposed high-precision method for solving the foundation dynamic response is well suited to practical engineering applications.