期刊文献+
共找到24篇文章
< 1 2 >
每页显示 20 50 100
捆扎线缆束电磁耦合效应限值模型研究 被引量:1
1
作者 侯兴勃 史贤俊 廖剑 《仪表技术》 2014年第7期27-29,共3页
捆扎线缆束中线缆间的串扰研究是复杂电子系统电磁兼容预测分析的关键。提出了一种预测随机捆扎线缆束中线线耦合的限值模型,并设计了针对性试验进行验证。仿真结果与试验结果吻合良好,证明了限值模型的正确性和有效性。
关键词 电磁兼容 随机捆扎线缆束 线线耦合 限值模型
下载PDF
离散和限值因变量模型在微观消费研究中的价值和作用
2
作者 柳鹏程 《理论月刊》 2004年第8期53-56,共4页
传统的微观消费模型往往假设因变量是可以连续变化的,但是随着对消费问题的深入探讨,调查中出现了部分消费为零的情况,有时甚至要面对在几个不同性质的选项中进行选择,这些都极大的违反了“连续性”的假设。而因变量为离散和限值因变量... 传统的微观消费模型往往假设因变量是可以连续变化的,但是随着对消费问题的深入探讨,调查中出现了部分消费为零的情况,有时甚至要面对在几个不同性质的选项中进行选择,这些都极大的违反了“连续性”的假设。而因变量为离散和限值因变量的计量模型则是针对这类研究的模型。本文较为系统地介绍了三类单方程模型:托宾模型、二分变量模型和多段式因变量模型,并对以上模型在微观消费研究中的应用进行了介绍和评述。 展开更多
关键词 离散和因变量模型 消费 托宾模型 二分变量 分段式因变量模型
下载PDF
CISPR 36:2020中磁场强度限值探究
3
作者 崔强 张强 付君 《安全与电磁兼容》 2022年第1期55-59,76,共6页
为了理解CISPR 36:2020中规定的磁场强度限值,依据30 MHz以下辐射骚扰场强CISPR限值模型,根据电动车辆和混合动力车辆的实际使用和需要保护的无线电业务,对模型中的每个分量进行了详细的分析和研究,推导了CISPR 36:2020中30 MHz以下辐... 为了理解CISPR 36:2020中规定的磁场强度限值,依据30 MHz以下辐射骚扰场强CISPR限值模型,根据电动车辆和混合动力车辆的实际使用和需要保护的无线电业务,对模型中的每个分量进行了详细的分析和研究,推导了CISPR 36:2020中30 MHz以下辐射骚扰磁场强度限值,帮助读者理解限值的来源和限值的推导。 展开更多
关键词 限值模型 辐射骚扰 磁场强度
下载PDF
燃煤电厂脱硫废水排放指标限值的计算方法研究 被引量:1
4
作者 赵保华 范紫瑄 +2 位作者 温佳琪 张金柱 马双忱 《化工进展》 EI CAS CSCD 北大核心 2018年第A01期213-218,共6页
目前我国燃煤电厂脱硫废水标准DL/T997—2006的排放指标与限制内容已不符合社会发展需要,为此,本文提出了燃煤电厂脱硫废水排放指标限值相关计算方法。论文调研了美国和国内的相关规范,对排放指标确定范围的具体数值和算法模型展开深入... 目前我国燃煤电厂脱硫废水标准DL/T997—2006的排放指标与限制内容已不符合社会发展需要,为此,本文提出了燃煤电厂脱硫废水排放指标限值相关计算方法。论文调研了美国和国内的相关规范,对排放指标确定范围的具体数值和算法模型展开深入研究,结合我国行业发展状况和国情给出了具体的修订建议,通过计算模型得出脱硫废水污染物控制参数的直接排放限值,氯化物日最大排放限值≤500mg/L,总溶解固体(TDS)日最大排放限值≤215mg/L,硒≤1.5mg/L,汞≤0.005mg/L等。 展开更多
关键词 脱硫废水 标准修订 污染物排放指标 浓度计算模型
下载PDF
Numerical study of fatigue damage of asphalt concrete using cohesive zone model 被引量:5
5
作者 金光来 黄晓明 +1 位作者 张苏龙 梁彦龙 《Journal of Southeast University(English Edition)》 EI CAS 2013年第4期431-435,共5页
In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is es... In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model (CZM) is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete. 展开更多
关键词 fatigue damage indirect tensile fatigue test asphalt concrete cohesive zone model numerical simulation finite element method
下载PDF
Numerical simulation of residual stress field in green power metallurgy compacts by modified Drucker-Prager Cap model 被引量:6
6
作者 周蕊 张连洪 +1 位作者 何柏岩 刘玉红 《Transactions of Nonferrous Metals Society of China》 SCIE EI CAS CSCD 2013年第8期2374-2382,共9页
Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the invest... Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces. 展开更多
关键词 residual stress green compacts numerical simulation modified Drucker-Prager Cap model
下载PDF
重组后的CISPR 16系列标准
7
作者 陈世钢 《安全与电磁兼容》 2003年第6期3-3,共1页
IEC(国际电工委员会)/CISPR(国际无线电干扰特别委员会)2003年11月发布了重组后的CISPR16(无线电骚扰和抗扰度测量设备和方法规范》(Specification for radio distur-bance and immunity measuring apparatus and methods)系列标准。该... IEC(国际电工委员会)/CISPR(国际无线电干扰特别委员会)2003年11月发布了重组后的CISPR16(无线电骚扰和抗扰度测量设备和方法规范》(Specification for radio distur-bance and immunity measuring apparatus and methods)系列标准。该系列标准目前包括以下14个子标准,它们取代以前出版的CISPR16-1、CISPR16-2、CISPR16-3、CISPR16-4及其各修正案。 第一部分:无线电骚扰和抗扰度测量设备(Radio distur-bance and immunity measuring apparatus) CISPR 16-1-1:测量设备(Measuring apparatus) 展开更多
关键词 CISPR16 标准 无线电骚扰 抗扰度测量设备 不确定度 限值模型
下载PDF
3D modeling and deformation analysis for electromagnetic sheet forming process 被引量:5
8
作者 CUI Xiao-hui MO Jian-hua ZHU Ying 《Transactions of Nonferrous Metals Society of China》 SCIE EI CAS CSCD 2012年第1期164-169,共6页
Electromagnetic forming (EMF) is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling a... Electromagnetic forming (EMF) is a high-speed forming method which can be quite effective in increasing the forming limits of metal sheet. However, the EMF process is complicated due to magnetic-structure coupling analysis. Numerical simulation offers an opportunity to overcome the problem. Nevertheless, most present models for EMF process are limited to 2D axisymmetric model. So, a three-dimensional (3D) finite element model was established to analyze the electromagnetic sheet bulging. The contact between the sheet and the die and the effect of sheet deformation on the magnetic field analysis were both taken into consideration during the forming process. The simulation results of deflection at the sheet center and 20 mm away from the center were in agreement with the experimental ones. The plastic strain energy and plastic strain were analyzed. 展开更多
关键词 electromagnetic forming high velocity forming numerical simulation finite element mode magnetic force
下载PDF
GIS-based Numerical Modelling of Debris Flow Motion across Three-dimensional Terrain 被引量:2
9
作者 WU Jian CHEN Guang-qi +1 位作者 ZHENG Lu ZHANG Ying-bin 《Journal of Mountain Science》 SCIE CSCD 2013年第4期522-531,共10页
The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of eros... The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective. 展开更多
关键词 Debris flow Numerical simulation GIS MOVEMENT EROSION DEPOSITION Equilibriumconcentration
下载PDF
Mechanical Behavior of a Glass-fiber Reinforced Composite to Steel Joint for Ships 被引量:3
10
作者 Xiaowen Li Ping Li Zhuang Lin Dongmei Yang 《Journal of Marine Science and Application》 CSCD 2015年第1期39-45,共7页
The use of a glass-fiber reinforced composite in marine structures is becoming more common, particularly due to the potential weight savings. The mechanical response of the joint between a glass-fiber reinforced polym... The use of a glass-fiber reinforced composite in marine structures is becoming more common, particularly due to the potential weight savings. The mechanical response of the joint between a glass-fiber reinforced polymer (GRP) superstructure and a steel hull formed is examined and subsequently modified to improve performance through a combined program of modeling and testing. A finite-element model is developed to predict the response of the joint. The model takes into account the contact at the interface between different materials, progressive damage, large deformation theory, and a non-linear stress-strain relationship. To predict the progressive failure, the analysis combines Hashin failure criteria and maximum stress failure criteria. The results show stress response has a great influence on the strength and bearing of the joint. The Balsawood-steel interface is proved to be critical to the mechanical behavior of the joint. Good agreement between experimental results and numerical predictions is observed. 展开更多
关键词 glass-fiber reinforced composite marine structure mechanical behavior steel joint finite-element model progressivefailure Hashin failure criteria
下载PDF
Numerical simulation of tire/soil interaction using a verified 3D finite element model 被引量:6
11
作者 Namjoo Moslem Golbakhshi Hossein 《Journal of Central South University》 SCIE EI CAS 2014年第2期817-821,共5页
The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at dif... The compaction and stress generation on terrain were always investigated based on empirical approaches or testing methods for tire/soil interaction.However,the analysis should be performed for various tires and at different soil strengths.With the increasing capacity of numerical computers and simulation software,finite element modeling of tire/terrain interaction seems a good approach for predicting the effect of change on the parameters.In this work,an elaborated 3D model fully complianning with the geometry of radial tire 115/60R13 was established,using commercial code Solidwork Simulation.The hyper-elastic and incompressible rubber as tire main material was analyzed by Moony-Rivlin model.The Drucker-Prager yield criterion was used to model the soil compaction.Results show that the model realistically predicts the laboratory tests outputs of the modeled tire on the soft soil. 展开更多
关键词 tire/soil interaction finite element method(FEM) soil compaction stress distribution inflation pressure
下载PDF
Numerical analysis on Joule heating of double-loop channel induction furnaces 被引量:1
12
作者 赵涛 周孑民 +1 位作者 熊家政 李勇 《Journal of Central South University》 SCIE EI CAS 2009年第5期851-856,共6页
In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with ... In order to investigate Joule heating power,a three-dimensional finite element model(FEM) was developed to predict Joule heating power in the channels of double-loop inductor. The simulated results were compared with experimental data from low load trials for a 400 kW inductor. The results,such as power factor and Joule heating power,show reasonable correlation with experimental data,and Joule heating rate reaches the maximum at the corners and the minimum at the centre of the cross-section area. With increasing relative permeability of iron core,length of coils,current frequency and resistivity of metal melt,the power factor and Joule heating power change. It is concluded that current frequency,the resistivity and length of the coil play a critical role in determining the power factor and Joule heating power,whereas relative permeability of the magnetic core shows no significant influence on them. 展开更多
关键词 numerical simulation channel induction furnace Joule heating power factor
下载PDF
Study on Collision Between Two Ships Using Selected Parameters in Collision Simulation 被引量:2
13
作者 Dong-Myung Bae Aditya Rio Prabowo +2 位作者 Bo Cao Ahmad Fauzan Zakki Gunawan Dwi Haryadi 《Journal of Marine Science and Application》 CSCD 2016年第1期63-72,共10页
In the present analysis, several parameters used in a numerical simulation are investigated in an integrated study to obtain their influence on the process and results of this simulation. The parameters studied are el... In the present analysis, several parameters used in a numerical simulation are investigated in an integrated study to obtain their influence on the process and results of this simulation. The parameters studied are element formulation, friction coefficient, and material model. Numerical simulations using the non-linear finite element method are conducted to produce virtual experimental data for several collision scenarios. Pattern and size damages caused by collision in a real accident case are assumed as real experimental data, and these are used to validate the method. The element model study performed indicates that the Belytschko-Tsay element formulation should be recommended for use in virtual experiments. It is recommended that the real value of the friction coefficient for materials involved is applied in simulations. For the study of the material model, the application of materials with high yield strength is recommended for use in the side hull structure. 展开更多
关键词 ship collision collision accident non-linear finite element collision parameter hull structure
下载PDF
Numerical Simulation of Nonlinear Wave Force on a Quasi-ellipse Caisson 被引量:4
14
作者 YongxueWang Xiaozhong Ren GuoyuWang 《Journal of Marine Science and Application》 2011年第3期265-271,共7页
A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume... A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force. 展开更多
关键词 quasi-ellipse caisson wave force VOF method partial cell method
下载PDF
Flux Limiter Lattice Boltzmann Scheme Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number
15
作者 甘延标 许爱国 +1 位作者 张广财 李英骏 《Communications in Theoretical Physics》 SCIE CAS CSCD 2011年第9期490-498,共9页
We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which... We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems. 展开更多
关键词 lattice Boltzmann method flux limiter compressible flows Prandtl number
下载PDF
Study of numerical simulation on dual-frequency IP method with FEM 被引量:2
16
作者 LIN Jiayong DIN Maobin +1 位作者 TANG Jingtian YAN Hong 《Global Geology》 2010年第2期99-103,共5页
Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary valu... Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary value problem and variational equations, then replaced the complex resistivity of the model with the Cole-Cole model's parameters under ignoring the EM effect. Through solving the last linear equations, electric potentials of all the model's points were obtained. With changing model's parameters, the authors got different curves of the Fs and phases. According to the results of the simulation, the algorithm is proved to be correct and adaptable. 展开更多
关键词 dual-frequency IP numerical simulation frequency domain FEM
下载PDF
Numerical simulations of the coastal effects of tsunami waves caused by the 1993 Hokkaido-Nansei-Oki earthquake
17
作者 潘文亮 王盛安 蔡树群 《Chinese Journal of Oceanology and Limnology》 SCIE CAS CSCD 2010年第5期1029-1039,共11页
Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,w... Generally,tsunami waves become hazardous only when approaching the coast.Studying the runup and inundation of tsunami waves is important for understanding the tsunami evolution and for tsunami hazard assessment.Here,we simulated the 1993 Hokkaido-Nansei-Oki(HNO) tsunami using a finite-difference model based on nonlinear shallow-water equations.We focused on the runup and inundation of tsunami wave propagation onto coastal area of Okushiri Island near Hokkaido,Japan,and investigate the relationship of different runup heights with the morphology and bathymetry of the seashore.In the simulation,a nested 4-layer grid system and moving boundary technique are adopted to study runup and inundation.The calculated tsunami heights and inundations in the region agreed well with field measurements.The local bathymetric and topographic characteristics had a first-order effect on the runup.Numerical experiments show that the focusing of certain local bathymetric features would amplify both wave height and current velocity remarkably.The results show that computation on dense grids is necessary to reproduce the observed runup heights,and inundation velocity is an important factor preventing tsunami devastation.In addition,we discussed the potential capability of sediment transport to illustrate the impact of tsunami waves on coastal geomorphology. 展开更多
关键词 TSUNAMI coastal effects RUNUP INUNDATION numerical modeling
下载PDF
Relaxation Property and Stability Analysis of the Quasispecies Models
18
作者 FEGN Xiao-Li LI Yu-Xiao +1 位作者 GU Jian-Zhong ZHUO Yi-Zhong 《Communications in Theoretical Physics》 SCIE CAS CSCD 2009年第10期726-734,共9页
The relaxation property of both Eigen model and Crow-Kimura model with a single peak fitness landscape is studied from phase transition point of view. We first analyze the eigenvalue spectra of the replication mutatio... The relaxation property of both Eigen model and Crow-Kimura model with a single peak fitness landscape is studied from phase transition point of view. We first analyze the eigenvalue spectra of the replication mutation matrices. For sufficiently long sequences, the almost crossing point between the largest and seeond-largest eigenvalues locates the error threshold at which critical slowing down behavior appears. We calculate the critical exponent in the limit of infinite sequence lengths and compare it with the result from numerical curve fittings at sufficiently long sequences. We find that for both models the relaxation time diverges with exponent 1 at the error (mutation) threshold point. Results obtained from both methods agree quite well. From the unlimited correlation length feature, the first order phase transition is further confirmed. Finally with linear stability theory, we show that the two model systems are stable for all ranges of mutation rate. The Igigen model is asymptotically stable in terms of mutant classes, and the Crow-Kimura model is completely stable. 展开更多
关键词 relaxation time critical exponent error threshold
下载PDF
Mechanical analysis of fixed geosynthetic technique of GRPS embankment
19
作者 张军 郑俊杰 马强 《Journal of Central South University》 SCIE EI CAS 2013年第5期1368-1375,共8页
To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed a... To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data. 展开更多
关键词 mechanical model load transfer mechanism fixed geosynthetic technique GRPS embankment field test
下载PDF
Selective laser melting 3D printing of Ni-based superalloy: understanding thermodynamic mechanisms 被引量:14
20
作者 Mujian Xia Dongdong Gu +3 位作者 Guanqun Yu Donghua Dai Hongyu Chen Qimin Shi 《Science Bulletin》 SCIE EI CAS CSCD 2016年第13期1013-1022,共10页
A mesoscopic model has been established to investigate the thermodynamic mechanisms and densification behavior of nickel-based superalloy during additive manufacturing/three-dimensional (3D) printing (AM/3DP) by n... A mesoscopic model has been established to investigate the thermodynamic mechanisms and densification behavior of nickel-based superalloy during additive manufacturing/three-dimensional (3D) printing (AM/3DP) by numerical simulation, using a finite volume method (FVM). The influence of the applied linear energy density (LED) on dimensions of the molten pool, thermodynamic mechanisms within the pool, bubbles migration and resultant densification behavior of AM/3DP-processed superalloy has been discussed. It reveals that the center of the molten pool slightly shifts with a lagging of 4 ktm towards the center of the moving laser beam. The Mar- angoni convection, which has various flow patterns, plays a crucial role in intensifying the convective heat and mass transfer, which is responsible for the bubbles migration and densification behavior of AM/3DP-processed parts. At an optimized LED of 221.5 J/m, the outward convection favors the numerous bubbles to escape from the molten pool easily and the resultant considerably high relative density of 98.9 % is achieved. However, as the applied LED further increases over 249.5 J/m, the convection pattern is apparently intensified with the formation of vortexes and the bubbles tend to be entrapped by the rotating flow within the molten pool, resulting in a large amount of residual porosity and a sharp reduction in densification of the superalloy. The change rules of the relative density and the corresponding distribution of porosity obtained by experiments are in accordance with the simulation results. 展开更多
关键词 Selective laser melting 3D printing Mesoscopic simulation - Thermodynamics DENSIFICATION POROSITY
原文传递
上一页 1 2 下一页 到第
使用帮助 返回顶部