Display:

Results

Viewing 151 to 180 of 9346
2016-04-05
Technical Paper
2016-01-1680
Suresh Abasaheb Patil, Indrajit Dinkar More
Abstract This paper incorporates the on-road real time ride comfort testing and simulation analysis of the typical passenger cars. The main objective is to analyze the vibrations to the passenger’s seat through the irregularities on road surface. This analysis is carried out with the help state-of-the-art of FFT analyzer for different classes of passenger car for variety of road conditions which have been validated by simulation analysis in the 20-Sim simulation software.
2016-04-05
Journal Article
2016-01-1678
Etsuo Katsuyama, Ayana Omae
Abstract Vehicles equipped with in-wheel motors (IWMs) are capable of independent control of the driving force at each wheel. These vehicles can also control the motion of the sprung mass by driving force distribution using the suspension reaction force generated by IWM drive. However, one disadvantage of IWMs is an increase in unsprung mass. This has the effect of increasing vibrations in the 4 to 8 Hz range, which is reported to be uncomfortable to vehicle occupants, thereby reducing ride comfort. This research aimed to improve ride comfort through driving force control. Skyhook damper control is a typical ride comfort control method. Although this control is generally capable of reducing vibration around the resonance frequency of the sprung mass, it also has the trade-off effect of worsening vibration in the targeted mid-frequency 4 to 8 Hz range. This research aimed to improve mid-frequency vibration by identifying the cause of this adverse effect through the equations of motion.
2016-04-05
Technical Paper
2016-01-1616
Keiichiro Iida, Kunizo Onda, Akiyoshi Iida, Chisachi Kato, Shinobu Yoshimura, Yoshinobu Yamade, Yoshimitsu Hashizume, Yang Guo
Abstract One-way coupled simulation method that combines CFD, structural and acoustical analyses has been developed aiming at predicting the aeroacoustical interior noise for a wide range of frequency between 100 Hz and 4 kHz. Statistical Energy Analysis (SEA) has been widely used for evaluating transmission of sound through a car body and resulting interior sound field. Instead of SEA, we directly computed vibration and sound in order to investigate and understand propagation paths of vibration in a car body and sound fields. As the first step of this approach, we predicted the pressure fluctuations on the external surfaces of a car by computing the unsteady flow around the car. Secondly, the predicted pressure fluctuations were fed to the subsequent structural vibration analysis to predict vibration accelerations on the internal surfaces of the car.
2016-04-05
Technical Paper
2016-01-1617
Yoshinobu Yamade, Chisachi Kato, Shinobu Yoshimura, Akiyoshi Iida, Keiichiro Iida, Kunizo Onda, Yoshimitsu Hashizume, Yang Gou
Abstract A wall-resolving Large Eddy Simulation (LES) has been performed by using up to 40 billion grids with a minimum grid resolution of 0.1 mm for predicting the exterior hydrodynamic pressure fluctuations in the turbulent boundary layers of a test car with simplified geometry. At several sampling points on the car surface, which included a point on the side window, the door panel, and the front fender panel, the computed hydrodynamic pressure fluctuations were compared with those measured by microphones installed on the surface of the car in a wind tunnel, and effects of the grid resolution on the accuracy of the predicted frequency spectra were discussed. The power spectra of the pressure fluctuations computed with 5 billion grid LES agreed reasonably well with those measured in the wind tunnel up to around 2 kHz although they had some discrepancy with the measured ones in the low and middle frequencies.
2016-04-05
Journal Article
2016-01-1595
Haibo Wu, Jiangbin Zhou, Qian Chen, Gongwen Liu, Chaoqun Qian
Abstract In this paper we present the work which was done at Shanghai-VW for using computational aero-acoustic (CAA) simulation in the vehicle development process to assess and improve the buffeting behavior of a vehicle when the rear side window is open. In the first step, a methodology was established and validated against wind tunnel tests using a Sedan. The methodology consists of a calibration of the CAA model to represent the properties of the cabin interior of the real car in terms of damping, wall compliance and leakage followed by CAA simulations of the full vehicle at different wind speeds to obtain the transient flow field around the exterior shape and inside the passenger compartment. The interior noise spectra are directly calculated from the transient pressure inside the cabin.
2016-04-05
Technical Paper
2016-01-0411
Yosuke Akita, Kenji Abe, Yoshihiro Osawa, Yoshitsugu Goto, Yuji Nagasawa, Noboru Sugiura, Satoshi Wakamatsu, Kyoko Kosaka
Abstract If a vehicle is left in a humid environment, the coefficient of friction between the brake pads and discs increases, generating a discomforting noise during braking called brake squeal. It is assumed that this increase in the coefficient of friction in a humid environment is the effect of moisture penetrating between the brake friction surfaces. Therefore, this paper analyzes the factors causing coefficient of friction variation with moisture between the friction surfaces by dynamic observation of these surfaces. The observation was achieved by changing the disc materials from cast iron to borosilicate glass. One side of the glass brake disc was pushed onto the brake pad and the sliding surface was observed from the opposite side by a charge coupled device (CCD) camera. First, a preliminary test was carried out in a dry state using two pad materials with different wear properties to select the appropriate pad for observing the friction surfaces.
2016-04-05
Technical Paper
2016-01-0448
Yuliang Yang, Yu Yang, Ying Sun, Zhihong Dong, Yunquan Zhang
Abstract To improve the vehicle NVH performance and reduce the vibration of the exhaust system, average driving DOF displacement (ADDOFD) and dynamic analysis are used to optimize hanger locations. Based on the finite element model and rigid-flexible coupling model, exhaust system analysis model was established. According to the finite element model of the exhaust system, the free-free modal analysis is carried out, and the position of the hanging point of the exhaust system is optimized by using the ADDOFD method. Furthermore, through the dynamics analysis, the force of each hanger to the body is calculated by the dynamic analysis, then verify the rationality of the hanging position. The combination of the two methods can effectively determine the better NVH performance of the exhaust system with hanger locations in the earlier vehicle development process.
2016-04-05
Technical Paper
2016-01-0439
Tianqi Lv, Peijun Xu, Yunqing Zhang
Abstract The powertrain mounting system (PMS) plays an important role in improving the NVH (Noise, Vibration, Harshness) quality of the vehicle. In all running conditions of a vehicle, the displacements of the powertrain C.G. should be controlled in a prescribed range to avoid interference with other components in the vehicle. The conventional model of PMS is based on vibration theory, considering the rotation angles are small, ignoring the sequence of the rotations. However, the motion of PMS is in 3D space with 3 translational degrees of freedom and 3 rotational degrees of freedom, when the rotation angles are not small, the conventional model of PMS will cause errors. The errors are likely to make powertrain interfering with other components. This paper proposes a rigid body mechanics model of the powertrain mounting system. When the powertrain undergoes a large rotational motion, the rigid body mechanics model can provide more accurate calculation results.
2016-04-05
Technical Paper
2016-01-0440
Li Jie, Wang Wenzhu, Gao Xiong, Zhang Zhenwei
Abstract The ride comfort of heavy trucks is related to many factors, which include vehicle operating scenarios and vehicle structure parameters. An investigation of the influence of different factors on the ride comfort of heavy trucks was conducted. Based on the elastic theory of a uniform Euler-Bernoulli beam with both ends free, a 6 degree of freedom (DOF) half rigid-elastic vibration model of the vertical dynamic response was developed. The rigid-elastic model is more suitable to describe the actual movement of heavy trucks. The DOFs include vertical displacements of the body and each of two axles, the pitch displacement of the body, and the first and second order bending displacements of the body. The root mean square (RMS) values of body accelerations, dynamic deflections and relative dynamic loads form the evaluation index. Based on the rigid-elastic model, the influence of different factors on the ride comfort of heavy trucks is analyzed in the frequency domain.
2016-04-05
Technical Paper
2016-01-0443
Han Zhang, Gang Li, Yu Wang, Yuchuan Gu, Xiang Wang, Xuexun Guo
Abstract A vehicular hydraulic electrical energy regenerative semi-active suspension(HEERSS) was presented, and its working principle and performance were analyzed. Firstly, configuration and working principle of the HEERSS were described; Secondly, kinetic equation of HEERSS was deduced, and a skyhook controller was designed for HEERSS. The traditional skyhook control strategy should be changed for the characteristic of HEERSS, because the damping force during extension stroke could be controlled, but not in compression stroke. Thirdly, the performance of HEERSS was compared with passive suspension(PS), traditional semi-active suspension(TSS). The simulation results indicated that the performance of HEERSS would be compromise between TSS and PS, but the HEERSS could harvest vibration energy which was advanced than TSS and PS.
2016-04-05
Technical Paper
2016-01-0220
Ban Gu, Xiaohong Yuan
Abstract With great development of thermoelectric exhaust heat recovery technology, more and more attention has been paid to optimization of automotive thermoelectric generators (ATEGs). A lot of work has been done on optimization of flow field and thermal analysis. However, investigation on acoustic optimization is rather limited. In this paper, efforts have been paid to study the acoustic performance of a flat-plate TEG, and the feasibility of integration of automotive exhaust thermoelectric generator with muffler was discussed. The internal configuration of heat exchanger looks like “fishbone”. Four factors have been taken into consideration: the spacing of two fins, angle of the fins, the diameter of inlet and outlet of exchanger; and filling sound absorbing material in heat exchanger chamber. Based on these four factors, acoustic analysis was carried out.
2016-04-05
Technical Paper
2016-01-0249
Balashunmuganathan Vasanth, Kumar Sathish, Murali Govindarajalu, Mohsin Khan
In recent years reducing the automobile HVAC (Heating Ventilation and automobile conditioning) noise inside the vehicle cabin is one of the main criterions for all OEMs to provide comfort level to the passengers. The primary function of the HVAC is to deliver more air to the cabin with less noise generation for various blower speeds. Designing the optimum HVAC with less noise is one of the major challenges for all automotive manufacturers and HVAC suppliers. During the design stage, physical parts are not available and hence the simulation technique helps to evaluate the noise level of HVAC. In this study, a computational 1D (one dimensional) analysis is carried out to compute the airflow noise originated from the HVAC unit and propagated to the passenger cabin. Modeling has been done using unigraphics and the analysis is carried out using the commercial 1D software GT suite.
2016-04-05
Journal Article
2016-01-0277
Xingxing Feng, Kaimin Zhuo, Jinglai Wu, Vikas Godara, Yunqing Zhang
Abstract Interval inverse problems can be defined as problems to estimate input through given output, where the input and output are interval numbers. Many problems in engineering can be formulated as inverse problems like vehicle suspension design. Interval metrics, instead of deterministic metrics, are used for the suspension design of a vehicle vibration model with five degrees of freedom. The vibration properties of a vehicle vibration model are described by reasonable intervals and the suspension interval parameters are to be solved. A new interval inverse analysis method, which is a combination of Chebyshev inclusion function and optimization algorithm such as multi-island genetic algorithm, is presented and used for the suspension design of a vehicle vibration model with six conflicting objective functions. The interval design of suspension using such an interval inverse analysis method is shown and validated, and some useful conclusions are reached.
2016-04-05
Technical Paper
2016-01-0275
Frédéric Kihm, Andrew Halfpenny, Kurt Munson
Abstract Ground vehicle components are designed to withstand the real operational conditions they will experience during their service life. Vibration tests are performed to qualify their endurance. In order to replicate the same failure mechanism as in real conditions, the test specification must be representative of the service loads. The accelerated testing method, based on fatigue damage spectra (FDS), is a process for deriving a synthesized power spectral density (PSD) representing a random stationary Gaussian excitation and applied over a reduced duration. In real life, however, it is common that service loading includes non-Gaussian excitations. The consequences of not using a representative test signal during product validation testing are a higher field failure rate and added warranty costs. The objective of this paper is to describe a method for synthesizing a PSD test specification with a given kurtosis value, which represents a nonstationary non-Gaussian signal.
2016-04-05
Technical Paper
2016-01-0273
Richard DeJong, Se Ge Jung, John Van Baren
Abstract Methods for conducting accelerated vibration fatigue testing of structures, such as MIL-STD-810G, allow for the non-linear scaling of the test time with the inverse of the rms vibration amplitude based on the slope of the material S-N curve obtained from cyclic fatigue tests. The Fatigue Damage Spectrum (FDS) is used as a method to allow for different level scalings at different frequencies in a broadband vibration environment using the relative responses of resonances in the structure. A recent development in industry has been to mix impulses with random excitations to increase the vibration peak levels (as measured by the kurtosis), thereby accelerating the fatigue even more than would occur with a Gaussian excitation. This paper presents results from a study to determine the conditions under which high kurtosis, impulsive excitations actually produce high kurtosis responses in structural resonances thus increasing the level of the FDS.
2016-04-05
Technical Paper
2016-01-0632
Domenico Crescenzo, Viktor Olsson, Javier Arco Sola, Hongwen Wu, Andreas Cronhjort, Eric Lycke, Oskar Leufven, Ola Stenlaas
Abstract Due to demanding legislation on exhaust emissions for internal combustion engines and increasing fuel prices, automotive manufacturers have focused their efforts on optimizing turbocharging systems. Turbocharger system control optimization is difficult: Unsteady flow conditions combined with not very accurate compressor maps make the real time turbocharger rotational speed one of the most important quantities in the optimization process. This work presents a methodology designed to obtain the turbocharger rotational speed via vibration analysis. Standard knock sensors have been employed in order to achieve a robust and accurate, yet still a low-cost solution capable of being mounted on-board. Results show that the developed method gives an estimation of the turbocharger rotational speed, with errors and accuracy acceptable for the proposed application. The method has been evaluated on a heavy duty diesel engine.
2016-04-05
Technical Paper
2016-01-0643
Jian Zhang, Changwen Liu, Fengrong Bi, Yiqiang Pei, Xiaobo Bi
Abstract Knock threshold detection is the key of closed loop control of ignition in gasoline engine, and it is also the difficult point in knock measurement. In this paper, an investigation of knock detection in turbocharged gasoline engine using bispectrum slice and ensemble empirical mode decomposition (EEMD) based on the engine cylinder head vibration signals is presented. By adding some finite amplitude Gaussian white noises to the signal, EEMD keeps the signal continuous in different time span, and therefore the mode mixing inhering in the classical empirical mode decomposition (EMD) method is alleviated. Power spectrum density (PSD) estimation is used to determine the band range of the resonance frequency generated by knock component. EEMD is used to decompose the original signals, the time-frequency characteristics of the Intrinsic Mode Functions (IMF) are analyzed using Continues Wavelet Transform (CWT) due to its excellent time-frequency resolution.
2016-04-05
Journal Article
2016-01-0639
Brian C. Kaul, Benjamin Lawler, Akram Zahdeh
Abstract Engine acoustics measured by microphones near the engine have been used in controlled laboratory settings for combustion feedback and even combustion phasing control, but the use of these techniques in a vehicle where many other noise sources exist is problematic. In this study, surface-mounted acoustic emissions sensors are embedded in the block of a 2.0L turbocharged GDI engine, and the signal is analyzed to identify useful feedback features. The use of acoustic emissions sensors, which have a very high frequency response and are commonly used for detecting material failures for health monitoring, including detecting gear pitting and ring scuffing on test stands, enables detection of acoustics both within the range of human hearing and in the ultrasonic spectrum. The high-speed acoustic time-domain data are synchronized with the crank-angle-domain combustion data to investigate the acoustic emissions response caused by various engine events.
2016-04-05
Technical Paper
2016-01-1063
George Nerubenko
Abstract Up to 30% of engine noise is delivered by front end pulley combined with torsional vibration damper, and technically it is the main contributor to recorded engine noise level. So the novel solutions in terms of improving the design and performance of torsional vibration damper would help to reduce radically this component of engine noise. The results of dynamical study of patented torsional vibration damper combined with pulley are presented. Design and structure of torsional vibration damper is based on author’s US Patent 7,438,165 having the self-tuning control system for all frequencies in running engine in all operational regimes. Mathematical model has been used for the analysis of the emitting noise of engine having proposed torsional vibration damper. Attention is paid to mitigation of the sound power levels contributing by engine subsystem “end of crankshaft - torsional vibration damper - pulley”.
2016-04-05
Technical Paper
2016-01-1060
Arnab Ganguly, Niket Bhatia, Vikas Kumar Agarwal, Ulhas Mohite
Abstract With ride comfort in a motorcycle gaining significance, it is important to minimize vibration levels at the customer touch points. The reciprocating piston imparts rotary motion to the crankshaft which in turn induces unbalance forces and produces vibration in the vehicle, thus influencing the ride quality. Generally, the primary inertial forces are balanced by a combination of balancer body and crank web. However, being a commuter bike, a balancer body could not be accommodated due to cost and space constraints. In such scenario, the first order unbalance force cannot be completely eliminated but can only be redistributed by adding counterweight to the crankshaft. Proper distribution of these forces is required for optimum vibration levels at motorcycle touch sensitive points (TSP) such as handle bar, footrest etc.
2016-04-05
Technical Paper
2016-01-1061
Guiping Yue, Wenbo Niu, Jian Zhao, Dandan Kong, Yun Li, Hangsheng Hou
Abstract Gear whine noise impacts customer perception of vehicle interior quietness in general and sound quality in particular. It has been a frequently occurred annoying phenomenon during vehicle development and much discussed topic regarding transmission NVH refinement in automotive industry. This work pertains to a transmission gear whine issue encountered in prototype evaluations during a vehicle program development process. The effort centers itself on the optimization of transmission gear macroscopic and microscopic parameters to fix the issue which is deemed unacceptable for customers. Specifically, by using multi-body dynamics approach, this work carries out a transmission system whine noise simulation based on optimal gear macro parameter selection and micro tooth flank modification. The obtained results show that the proposed design changes could successfully resolve the issue, which is verified by subsequent test measurement and confirmed by subjective evaluations.
2016-04-05
Technical Paper
2016-01-1057
Masaya Miyazawa, Kei Mochizuki, Kondo Takashi
Abstract To achieve lightweight, low friction and fuel efficient engine, the crankshaft is required to be designed lightweight, small-diameter shaft, long stroke. In this case, vibration of the crankshaft is increased by reduction of shaft stiffness. The conventional way of dealing with this increased vibration used to be to add an inertia mass ring or a double mass damper. Such an approach, however, increases weight, making the balance of weight reduction and vibration reduction less readily achieved. This paper therefore reports on how the main factors causing crankshaft vibration to increase in the shaft with reduced stiffness were clarified. Based on that clarification, efforts were made to reduce crankshaft vibration without increasing the weight of the crankshaft system. Measurement and analysis were used to analyze crankshaft vibration during operation.
2016-04-05
Technical Paper
2016-01-1058
Naoki Toyama, Takashi Hoshi
Abstract The subject is technology for damping forced vibration in the multiplate wet clutches used in hybrid vehicle transmissions. As a predictive technique for forced vibration caused by the structure of the clutch, three-dimensional simulation was used in the present study to anticipate the modes of vibration that occur. Next, a one-dimensional simulation was created as a predictive technique for drivetrain torsional vibration from the engine to the driveshaft. The one-dimensional simulation created was used to extract the modes of operation that are severe with regard to forced vibration from target values for vibration anticipated from the vehicle body. The results obtained were used with three-dimensional simulation to change the clutch structure to provide greater latitude with regard to the target for forced vibration.
2016-04-05
Technical Paper
2016-01-1053
TL Balasubramanian, Vigneshwara Raja Kesavan, V Lakshminarasimhan, Vamsi Krishna Balla, Suresh Palani
Abstract This paper discusses various noise sources of cylinder head assembly and focuses on design options developed to reduce the cylinder head noise in a single cylinder, 110cc scooter engine. Various experimental procedures were used for identification and ranking of different noise sources. In case of air-cooled small engines, temperature effects are dominant and as a consequence certain noises stand out in hot condition causing severe noise discomfort. After identifying the reasons for abnormal cylinder head noise, countermeasure mechanisms for reducing unintended impacts of valve train/ rocker arm in the layout were developed. The side-effects due to introduction of these additional mechanisms are studied using performance metrics. It is essential to limit noise deterioration over time to increase customer satisfaction. Simulation cycles were developed to quantify the cylinder head noise deterioration using accelerated testing procedures.
2016-04-05
Technical Paper
2016-01-1044
Toshiaki Kobayashi
Abstract This paper describes the establishment of a new method for predicting piston skirt scuffing in the internal combustion engine of a passenger car. The authors previously constructed and reported a method that uses 3D piston motion simulation to predict piston slap noise and piston skirt friction. However, that simulation did not have a clear index for evaluation of scuffing that involves piston skirt erosion, and it impressed shortage of the predictive accuracy of a scuffing. Therefore, the authors derived a new evaluation index for piston skirt scuffing by actually operating an internal combustion engine using multiple types of pistons to reproduce the conditions under which scuffing occurs, and comparing with the results of calculating the same conditions by piston motion simulation.
2016-04-05
Technical Paper
2016-01-1040
Rostislav Hadas, Rolf Sauerstein, Radúz Zahoranský, Michael Stilgenbauer, Matthias Ruh
Abstract A new generation of gasoline engine turbochargers has been developed with a focus on high performance and excellent NVH characteristics, especially with regards to the wastegate control system. With the recent introduction of EU6 emission standards, there is a clear demand to precisely control the flow of exhaust gas through the turbine wastegate. Engine operational duty cycles measured on EU6 compatible vehicles have shown increased stresses on wastegate parts due to a higher amount of regulation strokes during operation. Recent developments in the compact design of exhaust systems together with high pressure pulsation forces acting on wastegate flaps constitute the main challenges facing turbo engineers in the effort to achieve customer durability while meeting NVH requirements. For the development of a new generation of wastegate control systems a unique load prediction model was duly developed.
2016-04-05
Technical Paper
2016-01-1094
Fang Liao, Weimin Gao, Yan Gu, Fei Kang, Yinan Li, Cheng Wang
Abstract Generally, the gear rattle noise prediction models are composed of the mass and stiffness elements. The proposals are about the gear inertia or backlash and the shaft inertia or stiffness, but there are many detailed designs in the same inertia, stiffness or backlash conditions. Therefore, these proposals can’t guide detailed designs. These models only investigate the rattle in the rotating degree, and ignore rattle contribution in the radical and axial directions. Those prediction models only consider one or several factors which affect the rattle noise performance. It is difficult to predict the influence of individual factor and multi-factors coupling on the gear rattle noise in a rattle simulation model.
2016-04-05
Technical Paper
2016-01-1069
Masayoshi Otaka, Taro Kasahara, Kenichi Komaba
Abstract As a means of further improving combustion efficiency of gasoline engine, an increase in compression ratio, which enhances the risk of knocking, is thinkable. To optimize engine combustion parameters, a technology that can precisely detect knocking is desirable. Presently skillful experts have been evaluating knocking subjectively by listening to radiation noise so far. The authors developed a device that can precisely detect knocking by means of processing sound signals, which are captured by a high-performance microphone that is sensitive in the wide frequency range. Shock waves induced by knocking cause in-cylinder gas vibrations that emits metallic hit noises from the outer engine wall. We studied how to identify the feature values of frequency characteristics when knocking occurs, under the assumption that the engine radiation noise includes more than 2nd-order harmonic components with respect to the basic frequency of the in-cylinder gas vibration mode.
2016-04-05
Journal Article
2016-01-0459
Jian Zhao, Jing Su, Bing Zhu, Jingwei Shan
Abstract Proper tire pressure is very important for multiple driving performance of a car, and it is necessary to monitor and warn the abnormal tire pressure online. Indirect Tire Pressure Monitoring System (TPMS) monitors the tire pressure based on the wheel speed signals of Anti-lock Braking System (ABS). In this paper, an indirect TPMS method is proposed to estimate the tire pressure according to its resonance frequency of circumferential vibration. Firstly, the errors of ABS wheel speed sensor system caused by the machining tolerance of the tooth ring are estimated based on the measured wheel speed using Recursive Least Squares (RLS) algorithm and the measuring errors are eliminated from the wheel speed signal. Then, the data segments with drive train torsional vibration are found out and eliminated by the methods of correlation analysis.
2016-04-05
Technical Paper
2016-01-0480
Weiguo Zhang, Mark Likich, Mac Lynch, John White
Abstract The noise radiated from the snorkel of an air induction system (AIS) can be a major noise source to the vehicle interior noise. This noise source is typically quantified as the snorkel volume velocity which is directly related to vehicle interior noise through the vehicle noise transfer function. It is important to predict the snorkel volume velocity robustly at the early design stage for the AIS development. Design For Six Sigma (DFSS) is an engineering approach that supports the new product development process. The IDDOV (Identify-Define-Develop-Optimize-Verify) method is a DFSS approach which can be used for creating innovative, low cost and trouble free products on significant short schedules. In this paper, an IDD project which is one type of DFSS project using IDDOV method is presented on developing a robust simulation process to predict the AIS snorkel volume velocity. First, the IDDOV method is overviewed and the innovative tools in each phase of IDDOV are introduced.
Viewing 151 to 180 of 9346

Filter

Subtopics