Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.

Resistance of full-scale beams against close-in explosions.Numerical modeling and field tests

This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.

Design and Test of Closed Hydraulic Transmission System of 4WD High-clearance Wheeled Sprayer

In conjunction with the working characteristics of the high-clearance wheeled sprayer and the benefits of the closed hydraulic system,a series of reasonable working parameters should be established,and a hydraulic system that fulfills the requisite specifications should be designed.The AMESim software model is employed to construct a closed hydraulic transmission system,and the simulation analysis is then performed according to the data of hydraulic components.According to analysis results,the prototype can be optimized and upgraded,and a verification test is further carried out.The test results demonstrate that the designed closed hydraulic transmission system meets the actual working requirements of the high-clearance wheeled sprayer and provides a stable experimental platform for intelligent control of agricultural machinery.

Rainfall-Runoff Modeling and Hydrological Responses to the Projected Climate Change for Upper Baro Basin, Ethiopia

This paper presents the results of Rainfall-Runoff modeling and simulation of hydrological responses under changing climate using HEC-HMS model. The basin spatial data was processed by HEC-GeoHMS and imported to HEC-HMS. The calibration and validation of the HEC-HMS model was done using the observed hydrometeorological data (1989-2018) and HEC-GeoHMS output data. The goodness-of-fit of the model was measured using three performance indices: Nash and Sutcliffe coefficient (NSE) = 0.8, Coefficient of Determination (R2) = 0.8, and Percent Difference (D) = 0.03, with values showing very good performance of the model. Finally, the optimized HEC-HMS model has been applied to simulate the hydrological responses of Upper Baro Basin to the projected climate change for mid-term (2040s) and long-term (2090s) A1B emission scenarios. The simulation results have shown a mean annual percent decrease of 3.6 and an increase of 8.1 for Baro River flow in the 2040s and 2090s scenarios, respectively, compared to the baseline period (2000s). A pronounced flow variation is rather observed on a seasonal basis, reaching a reduction of 50% in spring and an increase of 50% in autumn for both mid-term and long-term scenarios with respect to the base period. Generally, the rainfall-runoff model is developed to solve, in a complementary way, the two main problems in water resources management: the lack of gauged sites and future hydrological response to climate change data of the basin and the region in general. The study results imply that seasonal and time variation in the hydrologic cycle would most likely cause hydrologic extremes. And hence, the developed model and output data are of paramount importance for adaptive strategies and sustainable water resources development in the basin.

Geomechanics model test and numerical simulation of 2G-NPR bolt support effect in an active fault tunnel

Active faults are a common adverse geological phenomenon that can occur during tunnel excavation and has a very negative impact on the construction and operation of the tunnel.In this paper,the grade IV rock surrounding the cross-fault tunnel with poor geological conditions has been chosen for the study.The support capacity of 2^(nd) Generation-Negative Poisson’s Ratio(2G-NPR)bolt in an active fault tunnel has been carried out on the basis of relevant results obtained from the geomechanical model test and numerical investigations of failure model for existing unsupported fault tunnel.The investigation shows that surrounding rock of the tunnel is prone to shear deformation and crack formation along the fault,as a result,the rock mass on the upper part of the fault slips as a whole.Furthermore,small-scale deformation and loss of blocks are observed around the tunnel;however,the 2G-NPR bolt support is found to be helpful in keeping the overall tunnel intact without any damage and instability.Due to the blocking effect of fault,the stress of the surrounding rock on the upper and lower parts of the fault is significantly different,and the stress at the left shoulder of the tunnel is greater than that at the right shoulder.The asymmetrical arrangement of 2G-NPR bolts can effectively control the asymmetric deformation and instability of the surrounding rock.The present numerical scheme is in good agreement with the model test results,and can reasonably reflect the stress and displacement characteristics of the surrounding rock of the tunnel.In comparison to unsupported and ordinary PR(Poisson’s Ratio)bolt support,2G-NPR bolt can effectively limit the fault slip and control the stability of the surrounding rock of the fault tunnel.The research findings may serve as a guideline for the use of 2G-NPR bolts in fault tunnel support engineering.

Steric Sea Level Change in Twentieth Century Historical Climate Simulation and IPCC-RCP8.5 Scenario Projection: A Comparison of Two Versions of FGOALS Model

To reveal the steric sea level change in 20th century historical climate simulations and future climate change projections under the IPCC＇s Representative Concentration Pathway 8.5 （RCP8.5） scenario, the results of two versions of LASG/IAP＇s Flexible Global Ocean-Atmosphere-Land System model （FGOALS） are analyzed. Both models reasonably reproduce the mean dynamic sea level features, with a spatial pattern correlation coefficient of 0.97 with the observation. Characteristics of steric sea level changes in the 20th century historical climate simulations and RCPS.5 scenario projections are investigated. The results show that, in the 20th century, negative trends covered most parts of the global ocean. Under the RCPS.5 scenario, global-averaged steric sea level exhibits a pronounced rising trend throughout the 21st century and the general rising trend appears in most parts of the global ocean. The magnitude of the changes in the 21st century is much larger than that in the 20th century. By the year 2100, the global-averaged steric sea level anomaly is 18 cm and 10 cm relative to the year 1850 in the second spectral version of FGOALS （FGOALS-s2） and the second grid-point version of FGOALS （FGOALS-g2）, respectively. The separate contribution of the thermosteric and halosteric components from various ocean layers is further evaluated. In the 20th century, the steric sea level changes in FGOALS-s2 （FGOALS-g2） are largely attributed to the thermosteric （halosteric） component relative to the pre-industrial control run. In contrast, in the 21st century, the thermosteric component, mainly from the upper 1000 m, dominates the steric sea level change in both models under the RCPS.5 scenario. In addition, the steric sea level change in the marginal sea of China is attributed to the thermosteric component.

Real-time multibody modeling and simulation of a scaled bogie test rig

In wheel–rail adhesion studies,most of the test rigs used are simplified designs such as a single wheel or wheelset,but the results may not be accurate.Alternatively,representing the complex system by using a full vehicle model provides accurate results but may incur complexity in design.To trade off accuracy over complexity,a bogie model can be the optimum selection.Furthermore,only a real-time model can replicate its physical counterpart in the time domain.Developing such a model requires broad expertise and appropriate software and hardware.A few published works are available which deal with real-time modeling.However,the influence of the control system has not been included in those works.To address these issues,a real-time scaled bogie test rig including the control system is essential.Therefore,a 1:4 scaled bogie roller rig is developed to study the adhesion between wheel and roller contact.To compare the performances obtained from the scaled bogie test rig and to expand the test applications,a numerical simulation model of that scaled bogie test rig is developed using Gensys multibody software.This model is the complete model of the test rig which delivers more precise results.To exactly represent the physical counterpart system in the time domain,a real-time scaled bogie test rig(RT-SBTR)is developed after four consecutive stages.Then,to simulate the RT-SBTR to solve the internal state equations and functions representing the physical counterpart system in rigs used are simplified designs such as a single wheel or wheelset,but the results may not be accurate.Alternatively,representing the complex system by using a full vehicle model provides accurate results but may incur complexity in design.To trade off accuracy over complexity,a bogie model can be the optimum selection.Furthermore,only a real-time model can replicate its physical counterpart in the time domain.Developing such a model requires broad expertise and appropriate software and hardware.A few published works are available which deal with real-time modeling.However,the influence of the control system has not been included in those works.To address these issues,a real-time scaled bogie test rig including the control system is essential.Therefore,a 1:4 scaled bogie roller rig is developed to study the adhesion between wheel and roller contact.To compare the performances obtained from the scaled bogie test rig and to expand the test applications,a numerical simulation model of that scaled bogie test rig is developed using Gensys multibody software.This model is the complete model of the test rig which delivers more precise results.To exactly represent the physical counterpart system in the time domain,a real-time scaled bogie test rig(RT-SBTR)is developed after four consecutive stages.Then,to simulate the RT-SBTR to solve the internal state equations and functions representing the physical counterpart system in equal or less than actual time,the real-time simulation environment is prepared in two stages.To such end,the computational time improved from 4 times slower than real time to 2 times faster than real time.Finally,the real-time scaled bogie model is also incorporated with the braking control system which slightly reduces the computational performances without affecting real-time capability.

Applications and verification of a computational energy dynamics model for mine climate simulations

A complete thermodynamic model is described for temperature and heat flow distribution simulation for ventilation networks in underground mines.The method is called the Computational Energy Dynamics(CED)model of the heat,mass,and energy transport.The Thermal and Humidity(TH)transport elements of the full model are described for advection,convection,and accumulation,encompassing heat capacity,radiation,latent heat for evaporation,and condensation in the airways,as well as variable heat conduction and accumulation in the rock strata.The thermal flywheel effect for time-dependent temperature field applications is included in the model solution.A CED model validation exercise is described,directly evaluating the iterated,minimized energy balance errors for the mechanical and thermal energy components for each network branch after a converged solution is determined.A simulation example relevant to mine safety and health is shown with the CED-TH model,demonstrating its capabilities in efficiency and accuracy in comparison with measurement results.

Centrifuge Model Tests and Numerical Simulations of the Impact of Underwater Explosion on an Air-Backed Steel Plate

Damage and threats to hydraulic and submarine structures by underwater explosions(UNDEXs)have raised much attention.The centrifuge model test,compared to prototype test,is a more promising way to examine the problem while reducing cost and satisfying the similitude requirements of both Mach and Froude numbers simultaneously.This study used a systematic approach employing centrifuge model tests and numerical simulations to investigate the effects of UNDEXs on an air-backed steel plate.Nineteen methodical centrifuge tests of UNDEXs were conducted.The shock wave pressure,bubble oscillation pressure,acceleration and the strain of the air-backed steel plate were recorded and compared with numerical studies using the finite element analysis(FEA)commercial software ABAQUS.By implementing empirically derived and physically measured pressures into the numerical models,the effects of the shock wave and bubble oscillation on the steel plate were investigated.Generally,the numerical results were in agreement with the experimental results.These results showed that the peak pressure of an UNDEX has a significant effect on the peak acceleration of the steel plate and that the impulse of the UNDEX pressure governs the peak strain of the steel plate.

Damage constitutive model of lunar soil simulant geopolymer under impact loading

Lunar base construction is a crucial component of the lunar exploration program,and considering the dynamic characteristics of lunar soil is important for moon construction.Therefore,investigating the dynamic properties of lunar soil by establishing a constitutive relationship is critical for providing a theoretical basis for its damage evolution.In this paper,a split Hopkinson pressure bar(SHPB)device was used to perform three sets of impact tests under different pressures on a lunar soil simulant geopolymer(LSSG)with sodium silicate(Na_(2)SiO_(3))contents of 1%,3%,5%and 7%.The dynamic stressestrain curves,failure modes,and energy variation rules of LSSG under different pressures were obtained.The equation was modified based on the ZWT viscoelastic constitutive model and was combined with the damage variable.The damage element obeys the Weibull distribution and the constitutive equation that can describe the mechanical properties of LSSG under dynamic loading was obtained.The results demonstrate that the dynamic compressive strength of LSSG has a marked strain-rate strengthening effect.Na_(2)SiO_(3) has both strengthening and deterioration effects on the dynamic compressive strength of LSSG.As Na_(2)SiO_(3) grows,the dynamic compressive strength of LSSG first increases and then decreases.At a fixed air pressure,5%Na_(2)SiO_(3) had the largest dynamic compressive strength,the largest incident energy,the smallest absorbed energy,and the lightest damage.The ZWT equation was modified according to the stress response properties of LSSG and the range of the SHPB strain rate to obtain the constitutive equation of the LSSG,and the model’s correctness was confirmed.

Quantifying major sources of uncertainty in projecting the impact of climate change on wheat grain yield in dryland environments

Modelling the impact of climate change on cropping systems is crucial to support policy-making for farmers and stakeholders.Nevertheless,there exists inherent uncertainty in such cases.General Circulation Models(GCMs)and future climate change scenarios(different Representative Concentration Pathways(RCPs)in different future time periods)are among the major sources of uncertainty in projecting the impact of climate change on crop grain yield.This study quantified the different sources of uncertainty associated with future climate change impact on wheat grain yield in dryland environments(Shiraz,Hamedan,Sanandaj,Kermanshah and Khorramabad)in eastern and southern Iran.These five representative locations can be categorized into three climate classes:arid cold(Shiraz),semi-arid cold(Hamedan and Sanandaj)and semi-arid cool(Kermanshah and Khorramabad).Accordingly,the downscaled daily outputs of 29 GCMs under two RCPs(RCP4.5 and RCP8.5)in the near future(2030s),middle future(2050s)and far future(2080s)were used as inputs for the Agricultural Production Systems sIMulator(APSIM)-wheat model.Analysis of variance(ANOVA)was employed to quantify the sources of uncertainty in projecting the impact of climate change on wheat grain yield.Years from 1980 to 2009 were regarded as the baseline period.The projection results indicated that wheat grain yield was expected to increase by 12.30%,17.10%,and 17.70%in the near future(2030s),middle future(2050s)and far future(2080s),respectively.The increases differed under different RCPs in different future time periods,ranging from 11.70%(under RCP4.5 in the 2030s)to 20.20%(under RCP8.5 in the 2080s)by averaging all GCMs and locations,implying that future wheat grain yield depended largely upon the rising CO2 concentrations.ANOVA results revealed that more than 97.22% of the variance in future wheat grain yield was explained by locations,followed by scenarios,GCMs,and their interactions.Specifically,at the semi-arid climate locations(Hamedan,Sanandaj,Kermanshah and Khorramabad),most of the variations arose from the scenarios(77.25%),while at the arid climate location(Shiraz),GCMs(54.00%)accounted for the greatest variation.Overall,the ensemble use of a wide range of GCMs should be given priority to narrow the uncertainty when projecting wheat grain yield under changing climate conditions,particularly in dryland environments characterized by large fluctuations in rainfall and temperature.Moreover,the current research suggested some GCMs(e.g.,the IPSL-CM5B-LR,CCSM4,and BNU-ESM)that made moderate effects in projecting the impact of climate change on wheat grain yield to be used to project future climate conditions in similar environments worldwide.

Model test vs virtual simulation of a VLCC FPSO hookup

This paper describes the model test and the virtual simulation respectively for the VLCC class FPSO hookup, as well as addresses their different applications to the mating operation between the FPSO and the soft yoke mooring system （SYMS） in extremely shallow water. The scope of the model test and the virtual simulation covers various installation stages including a series of positioning trials, positioning keeping and temporary mooring to the pre-installed SYMS mooring tower, pendulum mating, and yoke ballasting to storm-safe. The model test is to accurately verify bollard pull capacity to keep the FPSO in position and assess motion responses and mooring loads for the FPSO and installation vessels during various installation stages. The virtual simulation is to provide a virtual-reality environment, thus realistically replicating the hookup operation at the Simulation Test Center （STC） facility and identifying any deficiencies in key installation personnel, execution plan, or operation procedures. The methodologies of the model test and the virtual simulation addressed here can be easily extended to the deepwater applications such as positioning and installation operations of various floating systems.

Multi-Year Simulations and Experimental Seasonal Predictions for Rainy Seasons in China by Using a Nested Regional Climate Model (RegCM_NCC) Part Ⅱ:The Experimental Seasonal Prediction

A nested regional climate model has been experimentally used in the seasonal prediction at the China National Climate Center （NCC） since 2001. The NCC/IAP （Institute of Atmospheric Physics） T63 coupled GCM （CGCM） provides the boundary and initial conditions for driving the regional climate model （RegCM_NCC）. The latter has a 60-km horizontal resolution and improved physical parameterization schemes including the mass flux cumulus parameterization scheme, the turbulent kinetic energy closure scheme （TKE） and an improved land process model （LPM）. The large-scale terrain features such as the Tibetan Plateau are included in the larger domain to produce the topographic forcing on the rain-producing systems. A sensitivity study of the East Asian climate with regard to the above physical processes has been presented in the first part of the present paper. This is the second part, as a continuation of Part Ⅰ. In order to verify the performance of the nested regional climate model, a ten-year simulation driven by NCEP reanalysis datasets has been made to explore the performance of the East Asian climate simulation and to identify the model＇s systematic errors. At the same time, comparative simulation experiments for 5 years between the RegCM2 and RegCM_NCC have been done to further understand their differences in simulation performance. Also, a ten-year hindcast （1991-2000） for summer （June-August）, the rainy season in China, has been undertaken. The preliminary results have shown that the RegCM_NCC is capable of predicting the major seasonal rain belts. The best predicted regions with high anomaly correlation coefficient （ACC） are located in the eastern part of West China, in Northeast China and in North China, where the CGCM has maximum prediction skill as well. This fact may reflect the importance of the largescale forcing. One significant improvement of the prediction derived from RegCM_NCC is the increase of ACC in the Yangtze River valley where the CGCM has a very low, even a negative, ACC. The reason behind this improvement is likely to be related to the more realistic representation of the large-scale terrain features of the Tibetan Plateau. Presumably, many rain-producing systems may be generated over or near the Tibetan Plateau and may then move eastward along the Yangtze River basin steered by upper-level westerly airflow, thus leading to enhancement of rainfalls in the mid and lower basins of the Yangtze River. The real-time experimental predictions for summer in 2001, 2002, 2003 and 2004 by using this nested RegCM-NCC were made. The results are basically reasonable compared with the observations.

Simulations of a Hydrological Model as Coupled to a Regional Climate Model

Considering a detailed hydrologic model in the land surface scheme helps to improve the simulation of regional hydro-climatology. A hydrologic model, which includes spatial heterogeneities in precipitation and infiltration, is constructed and incorporated into the land surface scheme BATS. Via the coupled-model (i.e., a regional climate model) simulations, the following major conclusions are obtained: the simulation of surface hydrology is sensitive to the inclusion of heterogeneities in precipitation and infiltration; the runoff ratio is increased after considering the infiltration heterogeneity, a result which is more consistent with the observations of surface moisture balance over humid areas; the introduction of the parameterization of infiltration heterogeneity can have a greater influence on the regional hydro-climatology than the precipitation heterogeneity; and the consideration of the impermeable fraction for the region reveals some features that are closer to the trend of aridification over northern China.

SIMULATION OF PRESENT CLIMATE OVER EAST ASIA BY A REGIONAL CLIMATE MODEL

A 15-year simulation of climate over East Asia is conducted with the latest version of a regional climate model RegCM3 nested in one-way mode to the ERA40 Re-analysis data. The performance of the model in simulating present climate over East Asia and China is investigated. Results show that RegCM3 can reproduce well the atmospheric circulation over East Asia. The simulation of the main distribution patterns of surface air temperature and precipitation over China and their seasonal cycle/evolution, are basically agree with that of the observation. Meanwhile a general cold bias is found in the simulation. As for the precipitation, the model tends to overestimate the precipitation in northern China while underestimate it in southern China, particularly in winter. In general, the model has better performance in simulating temperature than precipitation.

Centrifugal model tests and numerical simulations for barrier dam break due to overtopping

The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.

On modeling approach for embedded real-time software simulation testing

Modeling technology has been introduced into software testing field. However, how to carry through the testing modeling effectively is still a difficulty. Based on combination of simulation modeling technology and embedded real-time software testing method, the process of simulation testing modeling is studied first. And then, the supporting environment of simulation testing modeling is put forward. Furthermore, an approach of embedded real-time software simulation testing modeling including modeling of cross-linked equipments of system under testing （SUT）, test case, testing scheduling, and testing system service is brought forward. Finally, the formalized description and execution system of testing models are given, with which we can realize real-time, closed loop, mad automated system testing for embedded real-time software.

Study on the Effects of Land Surface Heterogeneities in Temperature and Moisture on Annual Scale Regional Climate Simulation

The ＂combined approach＂, which is suitable to represent subgrid land surface heterogeneity in both interpatch and intra-patch variabilities, is employed in the BiOsphere/Atmosphere Transfer Scheme （BATS） as a land surface component of the regional climate model RegCM3 to consider the heterogeneities in temperature and moisture at the land surface, and then annual-scale simulations for 5 years （1988-1992） were conducted. Results showed that on the annual scale, the model＇s response to the heterogeneities is quite sensitive, and that the effect of the temperature heterogeneity （TH） is more pronounced than the moisture heterogeneity （MH）. On the intraannual scale, TH may lead to more （less） precipitation in warm （cold） seasons, and hence lead to larger intraannual variability in precipitation; the major MH effects may be lagged by about 1 month during the warm, rainy seasons, inducing -6% more precipitation for some sub-regions. Additionally, the modeled climate for the northern sub-regions shows larger sensitivities to the land surface heterogeneities than those for the southern sub-regions. Since state-of-art land surface models seldom account for surface intra-patch variabilities, this study emphasizes the importance of including this kind of variability in the land surface models.

Numerical Simulation for the Impact of Deforestation on Climate in China and Its Neighboring Regions

In this paper, the CCMOB model is used to study the effect of the deforestation on the climate of China and its neighboring regions. On the assumption that the forest in China would be replaced by the vegetation(such as grassland), the distribution of the albedo changed was calculated. The initial fields used were taken from the FGGE zonal mean data on 16 July, 1979. In the control simulation, the observed albedo data were used to modify the physical parameters of the original model. The control and sensitive experiments were run each for 210 days, in which the external forcing fields were fixed in July. As a result, we find that the East Asian Monsoon, Hadley cell and troposphere easterly jet are weakened for the deforestation in China. The precipitation and cloud amount over China are also decreased. The changes in evaporation and surface temperature are small. The results also show that the deforestation in China exerts a remarkable effect on the climate in the neighboring regions of China.

Feasibility and simulation model of a pilot scale membrane bioreactor for wastewater treatment and reuse from Chinese traditional medicine

The lack and pollution of water resource make wastewater reuse necessary. The pilot scale long-term tests for submerged membrane bioreactor were conducted to treat the effluents of anaerobic or aerobic treatment process for the high-strength Chinese traditional medicine wastewater. This article was focused on the feasibility of the wastewater treatment and reuse at shorter hydraulic retention time （HRT） of 5.0, 3.2 and 2.13 h. MLSS growth, membrane flux, vacuum values and chemical cleaning periods were also investigated. The experimental results of treating two-phase anaerobic treatment effluent demonstrated that the CODfilt was less than 100 mg/L when the influent COD was between 500-10000 mg/L at HRT of 5.0 h, which could satisfy the normal discharged standard in China. The experimental results to treat cross flow aerobic reactor effluent demonstrated that the average value of CODfilt was 17.28 mg/L when the average value of influent COD was 192.84 mg/L at HRT of 2.13 h during 106 d, which could completely meet the normal standard for water reuse. The maximum MLSS and MLVSS reached 24000 and 14500 mg/L at HRT of 3.2 h respectively. Membrane flux had maximal resume degrees of 94.7% at vacuum value of 0.02 MPa after cleaning. Chemical cleaning periods of membrane module were 150 d. A simulation model of operational parameters was also established based on the theory of back propagation neural network and linear regression of traditional mathematical model. The simulation model showed that the optimum operational parameters were suggested as follows： HRT was 5.0 h, SRT was 100 d, the range of COD loading rate was between 10.664-20.451 kg/（m3.d）, the range of MLSS was between 7543-13694 mg/L.