基于Plant Simulation仿真技术的装配生产线优化研究

【目的】优化装备生产线,缩短产品交付周期。【方法】基于Plant Simulation仿真技术,对装配生产线进行建模、编程、仿真、分析和优化,有效计算产品产量和成本,识别并优化装配生产线的瓶颈工位。【结果】优化后的三维产线仿真模型产量增幅接近10%,生产效率明显提升。【结论】通过将智能制造技能竞赛和科研教学活动相结合,能够凝练总结竞赛内容,促使教师紧盯前沿知识,创新改革教学内容,实现以赛促教、以赛促学、以赛促改、以赛促建的多重目标。

Analysis of Micromechanical Properties at the Interface of Pre-wet SBS Modified Asphalt Mixture Based on Molecular Simulation Technology

The pre-wetting of aggregate surface is a means to improve the interface performance of SBS modified asphalt and aggregate.The effect of pre-wetting technology on the interaction between SBS modified asphalt and aggregate was analyzed by molecular dynamics simulation.The diffusion coefficient and concentration distribution of SBS modified asphalt on aggregate surface are included.The simulation results show that the diffusion coefficient of the aggregate surface of SBS modified asphalt is increased by 47.6%and 70.5%respectively after 110#asphalt and 130#asphalt are pre-wetted.The concentration distribution of SBS modified asphalt on the aggregate surface after pre-wetting is more uniform.According to the results of interface energy calculation,the interface energy of SBS modified bitumen and aggregate can be increased by about 5%after pre-wetting.According to the results of molecular dynamics simulation,the pre-wetting technology can effectively improve the interface workability of SBS modified bitumen and aggregate,so as to improve the interface performance.

Molecular simulation study of the microstructures and properties of pyridinium ionic liquid[HPy][BF_(4)]mixed with acetonitrile

The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.

An Algorithm for Cloud-based Web Service Combination Optimization Through Plant Growth Simulation

In order to improve the efficiency of cloud-based web services,an improved plant growth simulation algorithm scheduling model.This model first used mathematical methods to describe the relationships between cloud-based web services and the constraints of system resources.Then,a light-induced plant growth simulation algorithm was established.The performance of the algorithm was compared through several plant types,and the best plant model was selected as the setting for the system.Experimental results show that when the number of test cloud-based web services reaches 2048,the model being 2.14 times faster than PSO,2.8 times faster than the ant colony algorithm,2.9 times faster than the bee colony algorithm,and a remarkable 8.38 times faster than the genetic algorithm.

Plastic flow and interfacial bonding behaviors of embedded linear friction welding process:Numerical simulation combined with thermophysical experiment

In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.

Numerical Simulation of Blood Flow Dynamics in a Stenosed Artery Enhanced by Copper and Alumina Nanoparticles

Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.

Omnidirectional simulation analysis of thermomechanical coupling mechanism in inertia friction welding of Ni-based superalloy

The coupling between heat and pressure is the kernel of inertia friction welding(IFW)and is still not fully understood.A novel 3D fully coupled finite element model based on a plastic friction pair was developed to simulate the IFW process of a Ni-based superalloy and reveal the omnidirectional thermo-mechanical coupling mechanism of the friction interface.The numerical model successfully simulated the deceleration,deformation processes,and peak torsional moments in IFW and captured the evolution of temperature,contact pressure,and stress.The simulated results were validated through measured thermal history,optical macrography,and axial shortening.The results indicated that interfacial friction heat was the primary heat source,and plastic deformation energy only accounted for 4%of the total.The increase in initial rotational speed and friction pressure elevated the peak temperature,reaching a maximum of 1525.5K at an initial rotational speed of 2000 r/min and friction pressure of 400 MPa.The interface heat generation could form an axial temperature gradient exceeding 320K/mm.The radial inhomogeneities of heat generation and temperature were manifested in a concentric ring distribution with maximum heat flux and temperature ranging from 2/5 to 2/3 radius.The radial inhomogeneities were caused by increasing linear velocity along the radius and an opposite distribution of contact pressure,which could reach 1.7 times the set pressure at the center.The circumferential inhomogeneity of thermomechanical distribution during rotary friction welding was revealed for the first time,benefiting from the 3D model.The deflection and transformation of distribution in contact pressure and Mises stress were indicators of plastic deformation and transition of quasi-steady state welding.The critical Mises stress was 0.5 times the friction pressure in this study.The presented modeling provides a reliable insight into the thermo-mechanical coupling mechanism of IFW and lays a solid foundation for predicting the microstructures and mechanical properties of inertia friction welded joints.

SolarDesign:An online photovoltaic device simulation and design platform

Solar Design(https://solardesign.cn/)is an online photovoltaic device simulation and design platform that provides engineering modeling analysis for crystalline silicon solar cells,as well as emerging high-efficiency solar cells such as organic,perovskite,and tandem cells.The platform offers user-updatable libraries of basic photovoltaic materials and devices,device-level multi-physics simulations involving optical–electrical–thermal interactions,and circuit-level compact model simulations based on detailed balance theory.Employing internationally advanced numerical methods,the platform accurately,rapidly,and efficiently solves optical absorption,electrical transport,and compact circuit models.It achieves multi-level photovoltaic simulation technology from“materials to devices to circuits”with fully independent intellectual property rights.Compared to commercial softwares,the platform achieves high accuracy and improves speed by more than an order of magnitude.Additionally,it can simulate unique electrical transport processes in emerging solar cells,such as quantum tunneling,exciton dissociation,and ion migration.

Numerical Simulation of Flow and Temperature Distribution in a Bottom-Blown Copper Bath

Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.

Boundary fluid constraints during electrochemical jet machining of large size emerging titanium alloy aerospace parts in gas–liquid flows:Experimental and numerical simulation

Large size titanium alloy parts are widely used in aerospace.However,they are difficult to manufacture using mechanical cutting technology because of severe tool wear.Electrochemical jet machining is a promising technology to achieve high efficiency,because it has high machining flexibility and no machining tool wear.However,reports on the macro electrochemical jet machining of large size titanium alloy parts are very scarce,because it is difficult to achieve effective constraint of the flow field in macro electrochemical jet machining.In addition,titanium alloy is very sensitive to fluctuation of the flow field,and a turbulent flow field would lead to serious stray corrosion.This paper reports a series of investigations of the electrochemical jet machining of titanium alloy parts.Based on the flow analysis and experiments,the machining flow field was effectively constrained.TB6 titanium alloy part with a perimeter of one meter was machined.The machined surface was smooth with no obvious machining defects.The machining process was particularly stable with no obvious spark discharge.The research provides a reference for the application of electrochemical jet machining technology to achieve large allowance material removal in the machining of large titanium alloy parts.

Plastic deformation mechanism of γ-phase U–Mo alloy studied by molecular dynamics simulations

Uranium–molybdenum(U–Mo) alloys are critical for nuclear power generation and propulsion because of their superior thermal conductivity, irradiation stability, and anti-swelling properties. This study explores the plastic deformation mechanisms of γ-phase U–Mo alloys using molecular dynamics(MD) simulations. In the slip model, the generalized stacking fault energy(GSFE) and the modified Peierls–Nabarro(P–N) model are used to determine the competitive relationships among different slip systems. In the twinning model, the generalized plane fault energy(GPFE) is assessed to evaluate the competition between slip and twinning. The findings reveal that among the three slip systems, the {110}<111>slip system is preferentially activated, while in the {112}<111> system, twinning is favored over slip, as confirmed by MD tensile simulations conducted in various directions. Additionally, the impact of Mo content on deformation behavior is emphasized. Insights are provided for optimizing process conditions to avoid γ → α′′ transitions, thereby maintaining a higher proportion of γ-phase U–Mo alloys for practical applications.

Land use/cover change and ecological network in Gansu Province,China during 2000-2020 and their simulations in 2050

Land use/cover change(LUCC)constitutes the spatial and temporal patterns of ecological security,and the construction of ecological networks is an effective way to ensure ecological security.Exploring the spatial and temporal change characteristics of ecological network and analyzing the integrated relationship between LUCC and ecological security are crucial for ensuring regional ecological security.Gansu is one of the provinces with fragile ecological environment in China,and rapid changes in land use patterns in recent decades have threatened ecological security.Therefore,taking Gansu Province as the study area,this study simulated its land use pattern in 2050 using patch-generating land use simulation(PLUS)model based on the LUCC trend from 2000 to 2020 and integrated the LUCC into morphological spatial pattern analysis(MSPA)to identify ecological sources and extract the ecological corridors to construct ecological network using circuit theory.The results revealed that,according to the prediction results in 2050,the areas of cultivated land,forest land,grassland,water body,construction land,and unused land would be 63,447.52,39,510.80,148,115.18,4605.21,8368.89,and 161,752.40 km^(2),respectively.The number of ecological sources in Gansu Province would increase to 80,with a total area of 99,927.18 km^(2).The number of ecological corridors would increase to 191,with an estimated total length of 6120.66 km.Both ecological sources and ecological corridors showed a sparse distribution in the northwest and dense distribution in the southeast of the province at the spatial scale.The number of ecological pinch points would reach 312 and the total area would expect to increase to 842.84 km^(2),with the most pronounced increase in the Longdong region.Compared with 2020,the number and area of ecological barriers in 2050 would decrease significantly by 63 and 370.71 km^(2),respectively.In general,based on the prediction results,the connectivity of ecological network of Gansu Province would increase in 2050.To achieve the predicted ecological network in 2050,emphasis should be placed on the protection of cultivated land and ecological land,the establishment of ecological sources in desert areas,the reinforcement of the protection for existing ecological sources,and the construction of ecological corridors to enhance the stability of ecological network.This study provides valuable theoretical support and references for the future construction of ecological networks and regional land resource management decision-making.

Synthesis and high frequency structure simulator electromagnetic simulation of hollow NC@CeO_(2)nanospheres for broad absorption bandwidth

Recent progress in microwave absorption materials stimulates the extensive exploration of rare earth oxide materials.Herein,we report the synthesis of a hollow sphere-based carbon material compounded with rare earth oxides.Hollow N-doped carbon nano-spheres loaded ceria composites(H-NC@CeO_(2))were designed and prepared by the template method,combined with in-situ coating,pyrolysis and chemical etching.By controlling the loading content of H-NC@CeO_(2)and adjusting the impedance matching of the material,the H-NC@CeO_(2)/PS(polystyrene)composite exhibited a minimum reflection loss(RL)of-50.8 dB and an effective absorption band-width(EAB)of 4.64 GHz at a filler ratio of 20wt%and a thickness of 2 mm.In accordance with measured electromagnetic parameters,simulations using the high frequency structure simulator(HFSS)software were conducted to investigate the impact of the honeycomb structure on the electromagnetic wave performance of H-NC@CeO_(2)/PS.By calculating the surface electric field and the material’s bulk loss density,the mechanism of electromagnetic loss for the honeycomb structure was elaborated.A method for structural design and man-ufacturing of broadband absorbing devices was proposed and a broadband absorber with an EAB of 11.9 GHz was prepared.This study presents an innovative approach to designing advanced electromagnetic(EM)wave absorbing materials with broad absorption band-widths.

Molecular dynamics simulations of helium transport through inorganic mineral nanopores

Helium transport through nanoscale inorganic mineral pores and pore throats is essential for its overall migration.To elucidate helium's transport dynamics within nanopores,we employed equilibrium and non-equilibrium molecular dynamics simulations to investigate helium's static self-diffusion and pressure-driven flow in quartz slit-shaped nanopores.We also introduced water and various gases,including hydrogen,methane,ethane,nitrogen,and carbon dioxide,into the nanopores to assess their influence on helium transport.Our findings indicate minimal helium adsorption on quartz pore surfaces.Under conditions where the pore size is less than 5 nm and the pressure under 10 MPa,environmental factors markedly influence helium diffusion.Large pore sizes,high temperatures,and low gas pressures enhance helium desorption and facilitate faster diffusion.We observed a positive correlation between helium flow velocity and factors such as pore size,pressure gradient,and surface smoothness of the pores.Notably,the presence of pore water and carrier gases in quartz nanopores,which diffuse more slowly than helium,tends to reduce helium surface adsorption and slow its diffusion.Among the carrier gases studied,nitrogen showed similar adsorption capacity,diffusivity,and stability to helium,while carbon dioxide displayed the highest adsorption capacity and the slowest diffusion rate,markedly differing from helium.Based on the simulation results,we concluded that water and carrier gases primarily function as transport mediums in helium migration,moving together with helium.Nitrogen,which shares similar properties with helium,effectively assists in this co-migration process.Conversely,carbon dioxide,due to its high adsorption capacity and slow diffusion,tends to be lost during co-migration.As a result,gas reservoirs with high nitrogen levels and low carbon dioxide levels are more likely to have higher helium concentrations.Additionally,the smaller pore sizes and higher gas pressures in caprocks can impede helium's diffusion,favoring its preservation in reservoirs.Moreover,the presence of water and carrier gases significantly obstructs these pores,further hindering helium's escape.

基于Simulation X的多轴转向系统功能安全仿真分析

基于ISO 26262《道路车辆功能安全》的标准要求,对特种车的多轴电液转向系统进行分析,以提高系统的安全性和可靠性。运用Simulation X软件建立详细的多轴特种车仿真模型,通过模拟故障模式注入进行仿真试验。对仿真结果和数据进行分析,评估得到故障的严重性、暴露度和可控性,从而确定相应的汽车安全完整性等级。基于故障注入仿真的汽车功能安全分析方法,在系统早期设计阶段可以作为评估架构安全性的重要手段。

基于Plant Simulation的化纤自动落丝系统仿真实验分析

基于Plant Simulation软件,构建一种化纤自动落丝系统的3D数字化模型,对系统效率关键输入因子进行参数化设置,将系统各工位三维数字模型导入软件中建立层次化运动机构图形结构,并运用Simltalk语言实现三维动作仿真,通过基础物理参数设置,保证了仿真数字化模型与现实系统更具一致性。通过多级实验设计分析了系统影响因子对于系统效率的影响特性曲线,并进一步通过动态参数化实验方法,计算出双输入因子对系统效能影响的敏感度。该实验结果可为化纤自动落丝系统的建设成本控制与可行性分析提供指导,具有较好的工程应用价值。

The Probability Density Function Related to Shallow Cumulus Entrainment Rate and Its Influencing Factors in a Large-Eddy Simulation

The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.

Global hybrid simulations of soft X-ray emissions in the Earth’s magnetosheath

Earth’s magnetopause is a thin boundary separating the shocked solar wind plasma from the magnetospheric plasmas,and it is also the boundary of the solar wind energy transport to the magnetosphere.Soft X-ray imaging allows investigation of the large-scale magnetopause by providing a two-dimensional(2-D)global view from a satellite.By performing 3-D global hybrid-particle-in-cell(hybrid-PIC)simulations,we obtain soft X-ray images of Earth’s magnetopause under different solar wind conditions,such as different plasma densities and directions of the southward interplanetary magnetic field.In all cases,magnetic reconnection occurs at low latitude magnetopause.The soft X-ray images observed by a hypothetical satellite are shown,with all of the following identified:the boundary of the magnetopause,the cusps,and the magnetosheath.Local X-ray emissivity in the magnetosheath is characterized by large amplitude fluctuations(up to 160%);however,the maximum line-of-sight-integrated X-ray intensity matches the tangent directions of the magnetopause well,indicating that these fluctuations have limited impact on identifying the magnetopause boundary in the X-ray images.Moreover,the magnetopause boundary can be identified using multiple viewing geometries.We also find that solar wind conditions have little effect on the magnetopause identification.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission will provide X-ray images of the magnetopause for the first time,and our global hybrid-PIC simulation results can help better understand the 2-D X-ray images of the magnetopause from a 3-D perspective,with particle kinetic effects considered.

Hybrid-Vlasov simulation of soft X-ray emissions at the Earth’s dayside magnetospheric boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.

Mixed D-vine copula-based conditional quantile model for stochastic monthly streamflow simulation

Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.