Numerical Simulation of the Seismic Response of a Horizontal Storage Tank Based on a SPH-FEM Coupling Method

A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.

Dynamic response of buildings under debris flow impact

This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.

Ice-Induced Vibrational Response of Single-Pile Offshore Wind-Turbine Foundations

Important challenges must be addressed to make wind turbines sustainable renewable energy sources.A typical problem concerns the design of the foundation.If the pile diameter is larger than that of the jacket platform,traditional mechanical models cannot be used.In this study,relying on the seabed soil data of an offshore wind farm,the m-method and the equivalent embedded method are used to address the single-pile wind turbine foundation problem for different pile diameters.An approach to determine the equivalent pile length is also proposed accordingly.The results provide evidence for the effectiveness and reliability of the model based on the equivalent embedded method.

Dynamic Response of A Group of Cylindrical Storage Tanks with Baffles Considering the Effect of Soil Foundation

The sloshing in a group of rigid cylindrical tanks with baffles and on soil foundation under horizontal excitation is studied analytically.The solutions for the velocity potential are derived out by the liquid subdomain method.Equivalent models with mass-spring oscillators are established to replace continuous fluid.Combined with the least square technique,Chebyshev polynomials are employed to fit horizontal,rocking and horizontal-rocking coupling impedances of soil,respectively.A lumped parameter model for impedance is presented to describe the effects of soil on tank structures.A mechanical model for the soil-foundation-tank-liquid-baffle system with small amount of calculation and high accuracy is proposed using the substructure technique.The analytical solutions are in comparison with data from reported literature and numerical codes to validate the effectiveness and correctness of the model.Detailed dynamic properties and seismic responses of the soil-tank system are given for the baffle number,size and location as well as soil parameter.

An extended multiple-support response spectrum method incorporating fluid-structure interaction for seismic analysis of deep-water bridges

The effects of ground motion spatial variability(GMSV)or fluid-structure interaction(FSI)on the seismic responses of deep-water bridges have been extensively examined.However,there are few studies on the seismic performance of bridges considering GMSV and FSI effects simultaneously.In this study,the original multiple-support response spectrum(MSRS)method is extended to consider FSI effect for seismic analysis of deep-water bridges.The solution of hydrodynamic pressure on a pier is obtained using the radiation wave theory,and the FSI-MSRS formulation is derived according to the random vibration theory.The influence of FSI effect on the related coefficients is analyzed.A five-span steel-concrete continuous beam bridge is adopted to conduct the numerical simulations.Different load conditions are designed to investigate the variation of the bridge responses when considering the GMSV and FSI effects.The results indicate that the incoherence effect and wave passage effect decrease the bridge responses with a maximum percentage of 86%,while the FSI effect increases the responses with a maximum percentage of 26%.The GMSV and FSI effects should be included in the seismic design of deep-water bridges.

Bamboo Fiber Modified Asphalt Mixture Proportion Design and Road Performances Based on Response Surface Method

In order to comprehensively utilize the remaining bamboo residue of bamboo products,this paper presents a research on recycling the bamboo fibers from bamboo residue for improving the performance of the asphalt mixtures.First of all,the basic performance parameters of sinocalamus affinis fiber,phyllostachys pubescens fiber,green bamboo fiber were tested and analyzed,and the optimal content and length were put forward.Then,the mix ratio design of the bamboo fiber modified asphalt mixture was further designed through the response surface method,and was verified the rationality of the mix ratio.Finally,the mixture specimens were made according to the experimental design mix ratio,and the high temperature,low temperature performance and moisture susceptibility of the bamboo fiber modified mixtures asphalt were tested.The results showed that the high temperature performance,low temperature performance and moisture susceptibility of bamboo fiber modified asphalt mixtures were improved compared with the performance of SBS modified asphalt mixture.When the length of bamboo fiber is 7.25 mm and the content of 0.22%,the road performance of the asphalt mixture was optimal.Consequentially,the decomposition of bamboo residue into bamboo fiber and its application in asphalt pavement can improve the reuse of bamboo waste,with remarkable environmental benefits and great promotion value.

Generalized response displacement methods for seismic analysis of underground structures with complex cross section

The response displacement method(RDM)is recommended for the seismic analysis of underground structures in the transverse direction for many codes,including bases for design of structures-seismic actions for designing geotechnical works(ISO 23469)and code for seismic design of urban rail transit structures(GB 50909-2014).However,there are some obvious limitations in the application of RDM.Springs and the shear stress of the soil could be approximately evaluated for the structures having a simple cross section,such as rectangular and circular structures.It is necessary to propose simplified seismic analysis methods for structures with complex cross sections.This paper refers to the idea of RDM and proposes three generalized response displacement methods(GRDM).In GRDM1,a part of the soil surrounding a structure is selected to generate a generalized underground structure with a rectangular cross section,and the same analysis model as RDM is applied to analyze the responses of the structure.In GRDM2,a hollow soil model without a generalized structure is used to compute the equivalent load caused by the relative displacement of the soil,and the soil-structure interaction model is applied to calculate the responses of the structure.In GRDM3,a continuous soil model is applied to compute the equivalent load caused by the relative displacement and shear stress of the soil,and the soil-structure interaction model is applied to analyze the responses of the structure,which is the same as the model used in GRDM2.The time-history analysis method(THAM)is used to evaluate the accuracy of the proposed simplified methods.Results show that the error of GRDM1 is about 20%,while the error is only 5%for GRDM2 and GRDM3.Among the three proposed methods,GRDM3 has obvious advantages regarding calculation efficiency and accuracy.Therefore,it is recommended to use GRDM3 for the seismic response analysis of underground structures that have conventional simple or complex cross sections.

Coseismic site response and slope instability using periodic boundary conditions in the material point method

This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuitive particle-relocation algorithm ensuring material points always remain within the computational mesh.The explicit generalized-a time scheme was implemented in MPM to enable the damping of spurious high frequency oscillations.Firstly,the MPM was verified against finite element method(FEM).Secondly,ability of the MPM in capturing the analytical transfer function was investigated.Thirdly,a symmetric embankment was adopted to investigate the effects of ground motion arias intensity(I_(a)),geometry dimensions,and constitutive models.The results show that the larger the model size,the higher the crest runout and settlement for the same ground motion.When using a Mohr-Coulomb model,the crest runout increases with increasing I_(a).However,if the strain-softening law is activated,the results are less influenced by the ground motion.Finally,the MPM results were compared with the Newmark sliding block solution.The simplified analysis herein highlights the capabilities of MPM to capture the full deformation process for earthquake engineering applications,the importance of geometry characterization,and the selection of appropriate constitutive models when simulating coseismic site response and subsequent large deformations.

Mechanical and Permeability Analysis and Optimization of Recycled Aggregate Pervious Concrete Based on Response Surface Method

In this paper,the effects of different influencing factors and factor interaction on the compressive strength and permeability of recycled aggregate pervious concrete(RAPC)were studied based on the response surface method(RSM).By selecting the maximum aggregate size,water cement ratio and target porosity as design variables,combined with laboratory tests and numerical analysis,the influences of three factors on the compressive strength and permeability coefficient of RAPC were revealed.The regression equation of compressive strength and permeability coefficient of recycled aggregate pervious concrete were established based on RSM,and the response surface model was optimized to determine the optimal ratio of RAPC under the conditions of meeting the mechanical and permeability properties.The results show that the mismatch item of the model is not significant,the model is credible,and the accuracy and reliability of the test are high,but the degree of uncorrelation between the test data and the model is not obvious.The sensitivity of the three factors to the compressive strength is water cement ratio>maximum coarse aggregate particle size>target porosity,and the sensitivity to the permeability coefficient is target porosity>maximum coarse aggregate particle size>water cement ratio.The absolute errors of the model prediction results and the model optimization results are 1.28 MPa and 0.19 mm/s,and the relative errors are 5.06%and 4.19%,respectively.With high accuracy,RSM can match the measured results of compressive strength and permeability coefficient of RAPC.

Optimization of Supercritical Extraction Process of Laoshan Black Tea by Response Surface Methodology

[Objectives]Laoshan black tea was subjected to supercritical CO_(2) extraction. [Methods]The extraction conditions of Laoshan black tea were studied by an orthogonal experiment and optimized by response surface methodology. [Results] The optimum extraction conditions of black tea extract by supercritical CO_(2) extraction were as follows: extraction pressure 23.53 MPa, extraction time 1.73 h, and extraction temperature 49.75 ℃, with which the extract yield could reach 5.15% theoretically. [Conclusions] Based on the traditional extraction process, a supercritical extraction method optimized by response surface methodology and a unique extraction process were formed, which enriches the extraction processes and methods of natural raw materials.

Long-term deformation analysis of Shuibuya concrete face rockfill dam based on response surface method and improved genetic algorithm

Due to the size effects of rockfill materials, the settlement difference between numerical simulation and in situ monitoring of rockfill dams is a topic of general concern.The constitutive model parameters obtained from laboratory triaxial tests often underestimate the deformation of high rockfill dams.Therefore, constitutive model parameters obtained by back analysis were used to calculate and predict the long-term deformation of rockfill dams.Instead of using artificial neural networks (ANNs), the response surface method (RSM) was employed to replace the finite element simulation used in the optimization iteration.Only 27 training samples were required for RSM, improving computational efficiency compared with ANN, which required 300 training samples.RSM can be used to describe the relationship between the constitutive model parameters and dam settlements.The inversion results of the Shuibuya concrete face rockfill dam (CFRD) show that the calculated settlements agree with the measured data, indicating the accuracy and efficiency of RSM.

Statistical Modification Analysis of Helical Planetary Gears based on Response Surface Method and Monte Carlo Simulation

Tooth modification technique is widely used in gear industry to improve the meshing performance of gearings. However, few of the present studies on tooth modification considers the influence of inevitable random errors on gear modification effects. In order to investigate the uncertainties of tooth modification amount variations on system＇s dynamic behaviors of a helical planetary gears, an analytical dynamic model including tooth modification parameters is proposed to carry out a deterministic analysis on the dynamics of a helical planetary gear. The dynamic meshing forces as well as the dynamic transmission errors of the sun-planet 1 gear pair with and without tooth modifications are computed and compared to show the effectiveness of tooth modifications on gear dynamics enhancement. By using response surface method, a fitted regression model for the dynamic transmission error（DTE） fluctuations is established to quantify the relationship between modification amounts and DTE fluctuations. By shifting the inevitable random errors arousing from manufacturing and installing process to tooth modification amount variations, a statistical tooth modification model is developed and a methodology combining Monte Carlo simulation and response surface method is presented for uncertainty analysis of tooth modifications. The uncertainly analysis reveals that the system＇s dynamic behaviors do not obey the normal distribution rule even though the design variables are normally distributed. In addition, a deterministic modification amount will not definitely achieve an optimal result for both static and dynamic transmission error fluctuation reduction simultaneously.

Distributed Collaborative Response Surface Method for Mechanical Dynamic Assembly Reliability Design

Because of the randomness of many impact factors influencing the dynamic assembly relationship of complex machinery, the reliability analysis of dynamic assembly relationship needs to be accomplished considering the randomness from a probabilistic perspective. To improve the accuracy and efficiency of dynamic assembly relationship reliability analysis, the mechanical dynamic assembly reliability（MDAR） theory and a distributed collaborative response surface method（DCRSM） are proposed. The mathematic model of DCRSM is established based on the quadratic response surface function, and verified by the assembly relationship reliability analysis of aeroengine high pressure turbine（HPT） blade-tip radial running clearance（BTRRC）. Through the comparison of the DCRSM, traditional response surface method（RSM） and Monte Carlo Method（MCM）, the results show that the DCRSM is not able to accomplish the computational task which is impossible for the other methods when the number of simulation is more than 100 000 times, but also the computational precision for the DCRSM is basically consistent with the MCM and improved by 0.40-4.63% to the RSM, furthermore, the computational efficiency of DCRSM is up to about 188 times of the MCM and 55 times of the RSM under 10000 times simulations. The DCRSM is demonstrated to be a feasible and effective approach for markedly improving the computational efficiency and accuracy of MDAR analysis. Thus, the proposed research provides the promising theory and method for the MDAR design and optimization, and opens a novel research direction of probabilistic analysis for developing the high-performance and high-reliability of aeroengine.

Optimization of High-Protein Glutinous Rice Flour Production Using Response Surface Method

A response surface method was employed to study the effect of α-amylase concentration, hydrolysis temperature and time on the production of high protein glutinous rice flour(HPGRF). The suspension of glutinous rice flour(15%) that contained 6.52% protein was gelatinized and subsequently hydrolyzed by thermostable α-amylase. The hydrolysis yielded 0.144–0.222 g/g HPGRF with 29.4%–45.4% protein content. Hydrolysis time exerted a significant effect, while enzyme concentration and hydrolysis temperature showed insignificant effect on the protein content and production yield of HPGRF. The result of response surface method showed that the optimum condition for the production of HPGRF that contained at least 36% protein was treating gelatinized 15% glutinous rice flour suspension with 0.90 Kilo Novo α-amylase Unit(KNU)/g α-amylase at 80 oC for 99 min. By carrying out the predicted hydrolysis condition, HPGRF with 35.9% protein and 61.8% carbohydrates was resulted. The process yielded 0.172 g/g HPGRF. HPGRF contained higher amount of essential amino acids compared to glutinous rice flour. HPGRF had higher solubility and lower swelling power, and also showed no pasting peak compared with glutinous rice flour.

Interval analysis method and convex models for impulsive response of structures with uncertain-but-bounded external loads

Two non-probabilistic, set-theoretical methods for determining the maximum and minimum impulsive responses of structures to uncertain-but-bounded impulses are presented. They are, respectively, based on the theories of interval mathematics and convex models. The uncertain-but-bounded impulses are assumed to be a convex set, hyper-rectangle or ellipsoid. For the two non-probabilistic methods, less prior information is required about the uncertain nature of impulses than the probabilistic model. Comparisons between the interval analysis method and the convex model, which are developed as an anti-optimization problem of finding the least favorable impulsive response and the most favorable impulsive response, are made through mathematical analyses and numerical calculations. The results of this study indicate that under the condition of the interval vector being determined from an ellipsoid containing the uncertain impulses, the width of the impulsive responses predicted by the interval analysis method is larger than that by the convex model; under the condition of the ellipsoid being determined from an interval vector containing the uncertain impulses, the width of the interval impulsive responses obtained by the interval analysis method is smaller than that by the convex model.

Grouping response method for equivalent static wind loads based on a modified LRC method

Wind loading is one of the most important loads for controlling the design of large-span roof structures. Equivalent static wind loads, which can generally aim at determining a specific response, are widely used by structural designers. A method for equivalent static wind loads applicable to multi-responses is proposed in this paper. A modified load- response-correlation （LRC） method corresponding to a particular peak response is presented, and the similarity algorithm implemented for the group response is described. The main idea of the algorithm is that two responses can be put into one group if the value of one response is close to that of the other response, when the structure is subjected to equivalent static wind loads aiming at the other response. Based on the modified LRC, the grouping response method is put forward to construct equivalent static wind loading. This technique can simultaneously reproduce peak responses for some grouped responses. To verify its computational accuracy, the method is applied to an actual large-span roof structure. Calculation results show that when the similarity of responses in the same group is high, equivalent static wind loads with high accuracy and reasonable magnitude of equivalent static wind distribution can be achieved.

An Intelligent Method for Structural Reliability Analysis Based on Response Surface

As water depth increases, the structural safety and reliability of a system become more and more important and challenging. Therefore, the structural reliability method must be applied in ocean engineering design such as offshore platform design. If the performance function is known in structural reliability analysis, the first-order second-moment method is often used. If the performance function could not be definitely expressed, the response surface method is always used because it has a very clear train of thought and simple programming. However, the traditional response surface method fits the response surface of quadratic polynomials where the problem of accuracy could not be solved, because the true limit state surface can be fitted well only in the area near the checking point. In this paper, an intelligent computing method based on the whole response surface is proposed, which can be used for the situation where the performance function could not be definitely expressed in structural reliability analysis. In this method, a response surface of the fuzzy neural network for the whole area should be constructed first, and then the structural reliability can be calculated by the genetic algorithm. In the proposed method, all the sample points for the training network come from the whole area, so the true limit state surface in the whole area can be fitted. Through calculational examples and comparative analysis, it can be known that the proposed method is much better than the traditional response surface method of quadratic polynomials, because, the amount of calculation of finite element analysis is largely reduced, the accuracy of calculation is improved, and the true limit state surface can be fitted very well in the whole area. So, the method proposed in this paper is suitable for engineering application.

Response surface method using grey relational analysis for decision making in weapon system selection

A proper weapon system is very important for a na- tional defense system. Generally, it means selecting the optimal weapon system among many alternatives, which is a multiple- attribute decision making （MADM） problem. This paper proposes a new mathematical model based on the response surface method （RSM） and the grey relational analysis （GRA）. RSM is used to obtain the experimental points and analyze the factors that have a significant impact on the selection results. GRA is used to an- alyze the trend relationship between alternatives and reference series. And then an RSM model is obtained, which can be used to calculate all alternatives and obtain ranking results. A real world application is introduced to illustrate the utilization of the model for the weapon selection problem. The results show that this model can be used to help decision-makers to make a quick comparison of alternatives and select a proper weapon system from multiple alternatives, which is an effective and adaptable method for solving the weapon system selection problem.

Advanced response surface method for mechanical reliability analysis

Based on the classical response surface method （RSM）, a novel RSM using improved experimental points （EPs） is presented for reliability analysis. Two novel points are included in the presented method. One is the use of linear interpolation, from which the total EPs for determining the RS are selected to be closer to the actual failure surface; the other is the application of sequential linear interpolation to control the distance be- tween the surrounding EPs and the center EP, by which the presented method can ensure that the RS fits the actual failure surface in the region of maximum likelihood as the center EPs converge to the actual most probable point （MPP）. Since the fitting precision of the RS to the actual failure surface in the vicinity of the MPP, which has significant contribution to the probability of the failure surface being exceeded, is increased by the presented method, the precision of the failure probability calculated by RS is increased as well. Numerical examples illustrate the accuracy and efficiency of the presented method.

An Improved Simplified Method for Predicting the Hydroelastic Response of Mat-Like VLFS

Very Large Floating Structures (VLFS) have received considerable attention recently. Efficient and accurate estimation of their hydroelastic responses in waves is very important for the design. The most efficient approach would obviously be the analytical one, Within the category of analytical approaches, the simplified method proposed by Ohkusu and his colleague are of special characteristics. However, when one studies their methods, several questions arise. The purpose of this paper is to critically study the simplified methods proposed by Ohkusu and his colleague in order to answer these questions. Some problems in their original methods have been found and possible improvements are suggested. It is concluded that the improved simplified method using the same idea of Ohkusu and his colleague could provide a reasonable estimate of the hydroelastic response of mat-like VLFS in a certain range of incident angles of waves.