期刊文献+
共找到1,390篇文章
< 1 2 70 >
每页显示 20 50 100
Improving convergence of incremental harmonic balance method using homotopy analysis method 被引量:3
1
作者 Yanmao Chen Jike Liu 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2009年第5期707-712,共6页
We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial va... We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence. 展开更多
关键词 incremental harmonic balance method Homotopy analysis method Initial value CONVERGENCE
下载PDF
DEFORMATION ANALYSIS OF SHEET METAL SINGLE-POINT INCREMENTAL FORMING BY FINITE ELEMENT METHOD SIMULATION 被引量:3
2
作者 MA Linwei MO Jianhua 《Chinese Journal of Mechanical Engineering》 SCIE EI CAS CSCD 2008年第1期31-35,共5页
Single-point incremental forming (SPIF) is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation a... Single-point incremental forming (SPIF) is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method(FEM) model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters' impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger. 展开更多
关键词 Sheet metal incremental forming DEFORMATION Finite element method(FEM) Numerical simulation
下载PDF
INCREMENTAL HARMONIC BALANCE METHOD FOR AIRFOIL FLUTTER WITH MULTIPLE STRONG NONLINEARITIES 被引量:1
3
作者 蔡铭 刘济科 李军 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 2006年第7期953-958,共6页
The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered... The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities. 展开更多
关键词 strongly nonlinear flutter incremental harmonic balance method BIFURCATION limit cycle
下载PDF
Incremental harmonic balance method for periodic forced oscillation of a dielectric elastomer balloon 被引量:1
4
作者 Yin WANG Ling ZHANG Jinxiong ZHOU 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI CSCD 2020年第3期459-470,共12页
Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices cal... Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy. 展开更多
关键词 dielectric elastomer nonlinear oscillation incremental harmonic method arc-length continuation
下载PDF
The Study of Force and the Mechanical Characteristic of Incremental Launching Construction Method on a Steel-Concrete Continuous Beam Bridge 被引量:2
5
作者 Xu Luo 《Journal of Architectural Research and Development》 2019年第6期46-50,共5页
The usage of steel-mixed composite beams is quite extensive today.During an event of constructing steel-mixed composite bridges,the incremental launching construction method is generally adopted.This paper mainly anal... The usage of steel-mixed composite beams is quite extensive today.During an event of constructing steel-mixed composite bridges,the incremental launching construction method is generally adopted.This paper mainly analyzes the force of incremental launching construction on a steel-concrete continuous beam bridge,the classification of incremental launching construction,the application of incremental launching construction in steel-mixed composite beams,the temporary facilities existing in incremental launching construction as well as their existing problems.Lastly,the analysis of the stress of composite beams in incremental launching construction is described by using the reference for the construction of mixed composite continuous beam bridges provided. 展开更多
关键词 STEEL-CONCRETE Composite BEAM incremental LAUNCHING Construction method Force DISPLACEMENT
下载PDF
THE ADAPTIVE PARAMETER INCREMENTAL METHOD FOR THE ANALYSIS OF SNAPPING PROBLEMS
6
作者 赵琪 叶天麒 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 1995年第9期851-858,共8页
By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effect... By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effectiveness of this method isdemonstrated by solving the snapping response of spherical caps subjected to centrallydistributed pressures. 展开更多
关键词 snapping problems adaptive parameter incremental method
下载PDF
A note on a family of proximal gradient methods for quasi-static incremental problems in elastoplastic analysis
7
作者 Yoshihiro Kanno 《Theoretical & Applied Mechanics Letters》 CAS CSCD 2020年第5期315-320,共6页
Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria.It has been demonstrated through numerical ... Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria.It has been demonstrated through numerical experiments that these methods can outperform conventional optimization-based approaches in computational plasticity.However,in literature these algorithms are described individually for specific yield criteria,and hence there exists no guide for application of the algorithms to other yield criteria.This short paper presents a general form of algorithm design,independent of specific forms of yield criteria,that unifies the existing proximal gradient methods.Clear interpretation is also given to each step of the presented general algorithm so that each update rule is linked to the underlying physical laws in terms of mechanical quantities. 展开更多
关键词 Elastoplastic analysis incremental problem Nonsmooth convex optimization First-order optimization method Proximal gradient method
下载PDF
Springback prediction for incremental sheet forming based on FEM-PSONN technology 被引量:6
8
作者 韩飞 莫健华 +3 位作者 祁宏伟 龙睿芬 崔晓辉 李中伟 《Transactions of Nonferrous Metals Society of China》 SCIE EI CAS CSCD 2013年第4期1061-1071,共11页
In the incremental sheet forming (ISF) process, springback is a very important factor that affects the quality of parts. Predicting and controlling springback accurately is essential for the design of the toolpath f... In the incremental sheet forming (ISF) process, springback is a very important factor that affects the quality of parts. Predicting and controlling springback accurately is essential for the design of the toolpath for ISF. A three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those from the experiment. The springback angle was found to be in accordance with the experimental result, proving the FEM to be effective. A coupled artificial neural networks (ANN) and finite element method technique was developed to simulate and predict springback responses to changes in the processing parameters. A particle swarm optimization (PSO) algorithm was used to optimize the weights and thresholds of the neural network model. The neural network was trained using available FEM simulation data. The results showed that a more accurate prediction of s!oringback can be acquired using the FEM-PSONN model. 展开更多
关键词 incremental sheet forming (ISF) springback prediction finite element method (FEM) artificial neural network (ANN) particle swarm optimization (PSO) algorithm
下载PDF
Numerical simulation and experimental investigation of incremental sheet forming process 被引量:4
9
作者 韩飞 莫健华 《Journal of Central South University of Technology》 EI 2008年第5期581-587,共7页
In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the pr... In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming. 展开更多
关键词 incremental sheet forming (ISF) sheet metal forming numerical simulation finite element method
下载PDF
Incremental-based Nonlinear Detection Algorithm for MIMO System 被引量:2
10
作者 李会勇 何子述 +1 位作者 刘本永 刘钊 《Journal of Electronic Science and Technology of China》 2006年第3期253-256,共4页
A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then sum... A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then summed with weights. The weight coefficients are incrementally computed to avoid direct computation of the inverse of a matrix, which greatly reduce the computational complexity. Simulation and comparison show that the proposed algorithm can obtain better performance of Bit Error Rate (BER) than linear Minimum Mean Square Error (MMSE). 展开更多
关键词 nonlinear detection incremental method wireless communication Multiple-Input Multiple-Output (MIMO)
下载PDF
INCREMENTAL MICRO-MECHANICAL MODEL OF PLAIN WOVEN FABRIC
11
作者 ZhangYitong HaoYongjiang LiCuiyu 《Acta Mechanica Solida Sinica》 SCIE EI 2004年第2期131-139,共9页
Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. I... Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. In this paper the expressions of incremental components of strain tensor are derived, the non-linear model of woven fabric is linearized physically and its geometric non-linearity survives. The convenience of determining the total deformation is shown by the choice of the coordinate system of the principal axes of the material, with the convergence of the incremental methods illustrated by examples. This incremental model furnishes a basis for numerical simulations of fabric draping and wrinkling based on the micro-mechanical model of fabric. 展开更多
关键词 micro-weaving structure MICRO-MECHANICS incremental constitutive model prin- cipal axes of material woven fabric incremental methods
下载PDF
Force-Based Quadrilateral Plate Bending Element for Plate Using Large Increment Method
12
作者 贾红学 刘西拉 《Journal of Donghua University(English Edition)》 EI CAS 2015年第3期345-350,共6页
A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successful... A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis. 展开更多
关键词 large increment method(LIM) displacement-based finite element method(FEM) Mindlin-Reissner plate theory spurious zero energy modes shear locking
下载PDF
Static and Dynamic Analysis of Mooring Lines by Nonlinear Finite Element Method 被引量:10
13
作者 杨敏冬 滕斌 《China Ocean Engineering》 SCIE EI 2010年第3期417-430,共14页
This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two node... This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters. 展开更多
关键词 geometrically nonlinear finite element mooring line incremental iteration method
下载PDF
ELASTO-PLASTICITY ANALYSIS BASED ON COLLOCATION WITH THE MOVING LEAST SQUARE METHOD 被引量:2
14
作者 Song Kangzu Zhang Xiong Lu Mingwan (Department of Engineering Mechanics,Tsinghua University,Beijing 100084,China) 《Acta Mechanica Solida Sinica》 SCIE EI 2003年第2期162-170,共9页
A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed ... A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis. 展开更多
关键词 meshless method ELASTO-PLASTICITY collocation method moving least square incremental formulation
下载PDF
Decoupling Conditions for Elasto-plastic Consolidation Question Based on Numerical Modeling Method 被引量:1
15
作者 Cheng Tao Wang Jingtao Dong Bichang 《Journal of China University of Geosciences》 SCIE CSCD 2005年第4期363-368,共6页
Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applie... Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot's consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that (1) the effect of the stress path should be taken into account in the simulation of the soil consolidation question; (2) this decoupling method can predict the evolvement of pore water pressure; (3) the settlement using cam-clay model is less than that using numerical model because of dilatancy. 展开更多
关键词 numerical modeling method Blot's consolidation stress path constitutive model liquisolid coupling decouple incremental diffusion equation dilatancy.
下载PDF
NONLINEAR QUASI-CONFORMING FINITE ELEMENT METHOD
16
作者 关玉璞 唐立民 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 1993年第3期269-276,共8页
The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation ... The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given. 展开更多
关键词 incremental principle of potential energy quasi-conforming finite element method NONLINEARITY
下载PDF
Geometrically Nonlinear Analysis of Structures Using Various Higher Order Solution Methods: A Comparative Analysis for Large Deformation
17
作者 Ali Maghami Farzad Shahabian Seyed Mahmoud Hosseini 《Computer Modeling in Engineering & Sciences》 SCIE EI 2019年第12期877-907,共31页
The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quad... The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures. 展开更多
关键词 Geometrically nonlinear analysis higher order methods predictor-corrector algorithms convergence rate sensitivity to the increment size
下载PDF
Analysis of Office-Teaching Comprehensive Buildings Using a Modified Seismic Performance Evaluation Method
18
作者 Hanbo Zhu Changqing Miao Meiling Zhuang 《Computer Modeling in Engineering & Sciences》 SCIE EI 2019年第3期471-491,共21页
Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is u... Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is urgent.The evaluation in this paper is based on a 11-story reinforced concrete office building,incremental dynamic analysis(IDA)is conducted in Perform 3D program using models capable to simulate all possible limit states up to collapse.Next,the probability distribution of post-earthquake casualties,rebuild costs repair costs and business downtime loss are calculated in PACT software for the studied building considering the modified component vulnerability groups and population flow models.The evaluation procedure can also shed light on other types of buildings in China.For non-typical functional building structures,this article proposes to build a finite element model of structural components and to classify the vulnerability groups based on the construction drawings,and to supply and improve the vulnerability library of appendages in FEMA P-58 according to the actual situation.In this way,the application scope of building seismic performance evaluation can be expanded. 展开更多
关键词 A MODIFIED SEISMIC capacity ASSESSMENT method frame-shear wall structure incremental dynamic analysis(IDA) SEISMIC performance ASSESSMENT ANALYSIS software(PACT)
下载PDF
基于分布式潮流控制器的海上风电系统谐波治理方法和控制策略 被引量:2
19
作者 唐爱红 宋幸 +3 位作者 尚宇菲 郭国伟 余梦琪 詹细妹 《电力系统自动化》 EI CSCD 北大核心 2024年第2期20-28,共9页
由于电力电缆的电容效应,海上风电经电缆汇集系统极易出现谐波谐振放大的现象,造成电能质量的下降。分布式潮流控制器属于基于电压源换流器的装置,在进行潮流调节的同时也能进行谐波治理。文中首先构建了海上风电系统的频域相关模型,基... 由于电力电缆的电容效应,海上风电经电缆汇集系统极易出现谐波谐振放大的现象,造成电能质量的下降。分布式潮流控制器属于基于电压源换流器的装置,在进行潮流调节的同时也能进行谐波治理。文中首先构建了海上风电系统的频域相关模型,基于该模型分析了谐波谐振放大的原因;随后,采用了将分布式潮流控制器串入海上风电系统的谐波治理方式,推导并得到了含分布式潮流控制器的海上风电系统的谐波特性。基于该谐波特性,设计了一种控制策略。该策略通过控制分布式潮流控制器实时跟踪使并网点谐波电压幅值为零的谐波补偿电压,从而降低并网点的谐波电压含量。仿真结果表明,所提出的基于分布式潮流控制器的海上风电系统谐波治理方法和控制策略能够有效地降低并网点的谐波电压,改善电能质量。 展开更多
关键词 海上风电 电能质量 谐波治理 分布式潮流控制器 变增量电导增量法
下载PDF
海外已开发油田新井项目经济评价方法及影响参数研究 被引量:1
20
作者 董翠 干卫星 《国际石油经济》 2024年第5期97-105,共9页
中国油气企业对中国境内已开发油气田开展新钻井经济评价工作时,主要采用直接增量法。受诸多因素制约,海外已开发油气田新钻井的经济评价不能简单应用直接增量法。国内与国外的已开发油气田在合同模式、财税制度、开发生产期限、弃置费... 中国油气企业对中国境内已开发油气田开展新钻井经济评价工作时,主要采用直接增量法。受诸多因素制约,海外已开发油气田新钻井的经济评价不能简单应用直接增量法。国内与国外的已开发油气田在合同模式、财税制度、开发生产期限、弃置费等方面存在明显差异。对于海外已开发油气田新钻井进行经济评价,需要考虑合同模式和财税制度的国别差异、评价期的选取、投资成本估算、弃置作业及弃置费处理这些关键影响因素,采用“有无对比法”评价新钻井项目的经济性。应以合同有效期和油田经济寿命期二者中的较短值为评价期。对投资成本的处理应遵循“相关性”原则。新钻井延期支付油田弃置费的时间价值依资源国、合同条款的具体规定而异。 展开更多
关键词 海外已开发油田 新钻井 经济评价 有无对比法 简单增量法 经济可行性
下载PDF
上一页 1 2 70 下一页 到第
使用帮助 返回顶部