首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The search for the optimum design of a given structure is a basic aspect of structural analysis. In recent years there has been a renewed interest in structural optimization, particularly in the field of composite material structures. The aim is often to obtain the best solution that is able to respect both the technical and the economic requirements. In this paper, the optimization process is pursued as the determination of the Pareto-optimal curves for composite material structures. The problem has been tackled by means of an analytical optimization methodology (trade-off method) appropriately transformed into a design procedure. A sensitivity analysis is also performed, in order to try out the measure in which the optimal solutions are affected by the variation of the design variables. Two examples are presented, concerning a hemispherical shell and a cylindrical shell.
Sommario Molti problemi della ricerca dell'ottimo progettuale presentano la difficoltà di coinvolgere nella soluzione più obiettivi da considerare simultaneamente. Il problema risulta piuttosto complesso (in particolare nel caso di strutture in materiale composito) quando gli obiettivi, spesso contrastanti, hanno natura diversa e quindi risultano non facilmente combinabili fra loro. Tale problema si pone spesso quando si vuole da un lato aumentare l'affidabilità strutturale e dall'altro diminuire i costi della struttura. Nel presente articolo, vengono affrontati alcuni aspetti dell'ottimizzazione di strutture in materiale composito: a tale scopo viene presentata una procedura di ottimizzazione per problemi multi-obiettivo che, avvalendosi del Trade-off method, è in grado di determinare la curva degli ottimi di Pareto. Successivamente viene effettuata un analisi di sensitività allo scopo di valutare la sensibilità delle soluzioni ottimali nei confronti della aleatorietà delle variabili di progetto. Vengono presentati in tal senso due esempi relativi a strutture in materiale composito.
  相似文献   

2.
Summary The problem of discontinuity wave propagation in hyperelastic shells is discussed within the framework of the nonlinear theory. A suitable linearization is proposed, such that explicit results can be obtained with direct calculations. Namely the speeds of propagation for acceleration and shock waves are evaluated, within the dependence of their values on the geometrical properties of the shell.
Sommario Si considera il problema della propagazione di onde di discontinuità in una volta sottile iperelastica nell'ambito della teoria nonlineare. Viene proposta una linearizzazione che consente di ottenere risultati espliciti con calcoli diretti. In particolare vengono determinati i valori delle velocità di propagazione delle onde d'accelerazione e d'urto e viene studiata la loro dipendenza dalle proprietà geometriche della volta.
  相似文献   

3.
In this study a numerical method for optimum synthesis of spatial function generator four-bar linkages is presented. The key features of the method are: the embedding of the assembly criterion and of the transmission ratio in the penalty function; the non-derivative minimization procedure; the parametric variation of the first-guess mechanism. The effectiveness of the method is illustrated by means of examples of optimum design. In some examples articulated mechanisms, which can replace cam systems in automatic machinery, are synthesized.
Sommario In questo studio si presenta un metodo per la sintesi ottimale di quadrilateri articolati spaziali generatori di funzione. Le principali caratteristiche del metodo proposto sono: l'inclusione del criterio di assemblabilità e dell'angolo di pressione nella funzione penalità; l'utilizzazione di un algoritmo di minimizzazione non derivativo; la variazione parametrica del meccanismo di primo tentativo. L'efficacia del metodo è messa in luce tramite una serie di esempi di progetto ottimale; in alcuni di essi sono stati sintetizzati dei quadrilateri spaziali in grado di sostituire i sistemi a camma spaziale nelle macchine automatiche.
  相似文献   

4.
Summary This paper discusses the rigid-plastic flexure of plates subjected to dynamic loading. A piecewise linear yield surface and small displacements are assumed. The plate is discretized by assuming the constant moment mixed triangular element of Herrmann. Under the above assumptions the dynamic loading of plates is turned into a problem in linear inequalities which is a priori discretized in time. An approximate technique is developed, which takes into account the influence of the strain-hardening and the strain rate sensitivity of material. Some numerical examples of plates subjected to a given velocity distribution confirm the efficiency of this approach.
Sommario Si analizza il problema della flessione rigidoplastica di una lastra soggetta a carichi dinamici. La lastra viene discretizzata assumendo l'elemento triangolare misto a stato di sollecitazione costante di Herrmann. Tale assunzione insieme alle ipotesi di piccoli spostamenti e di superficie di plasticizzazione lineare a tratti, consente di ricondurre il caricamento dinamico della lastra ad un problema retto da diseguaglianze lineari, che risulta discretizzato a priori nel tempo. Attraverso una tecnica approssimata si valuta l'influenza che hanno sulla risposta dinamica della struttura, l'incrudimento e la sensibilità alle velocità di deformazione del materiale. L'efficienza del procedimento è verificata attraverso alcuni esempi numerici di lastre di diversa geometria sottoposte ad una prefissata distribuzione di velocità.


Research supported by the National Research Council (C.N.R.P.A.d.I.S.) of Italy.  相似文献   

5.
Summary This paper deals with some theoretical developments of the problem of stability of doubly curved thin shallow shells on a rectangular plan. After some critical references to results obtained in earlier studies using the linear theory of stability, the problem of elastic equilibrium in the nonlinear range is approached and solved for some types of shells using piecewise linearisation. Load-displacement curves up to post-buckling are obtained. Deflections not conforming to structural symmetry are also considered. The influence of small geometric imperfections deviates the load-displacement curve to other equilibrium positions and hence a shell with imperfections attains its snapping point at a lower load than one without. A comparison with some experimental studies by other authors is given in the appendix.
Sommario Lo studio tratta di alcuni sviluppi teorici sul problema della stabilità dell'equilibrio delle volte sottili ribassate a doppia curvatura su pianta rettangolare. Nella prima parte vengono esposte alcune considerazioni critiche circa i risultati ottenuti in parte dall'Autore in lavori precedenti con la teoria lineare della stabilità. Il problema dell'equilibrio elastico viene poi affrontato in campo non-lineare e risolto per alcuni tipi di volte mediante una linearizzazione a tratti; le curve carico-spostamento vengono ricavate fino al comportamento post-critico. Vengono altresì ricavate curve carico-spostamento corrispondenti a deformate che non conservano la simmetria di struttura e di carico.La presenza di piccole imperfezioni devia la curva caricospostamento verso altre posizioni di equilibrio e la curva relativa raggiunge il punto di instabilità progressiva per valori del carico notevolmente inferiori di quelli relativi alla stessa volta senza imperfezioni. Infine si riporta in appendice un confronto con risultati sperimentali di altri Autori.


Study supported by the C.N.R. (Consiglio Nazionale delle Ricerche).  相似文献   

6.
Gaetano Zingone 《Meccanica》1972,7(3):214-221
Summary With reference to a theory expounded in a previous note [2] on the limit state behaviour of plates resting on an elastoplastic continuum subject to distributed loads, probable collapse mechanisms are considered in order to supplying more tractable solutions for technical practice. As the case is a general one, it is shown that the results obtained are suitable for following the load carrying capacity of plates and the relative collapse mechanism as the limit resistance of the soil and the spread of the load acting on it vary.Some examples are worked out to illustrate the applicational procedure.
Sommario Con riferimento ad una teoria esposta in una nota precedente [2] sul comportamento allo stato limite delle piastre su mezzo elasto-plastico continuo soggette a carichi diffusi, si esaminano delle schematizzazioni di probabile collasso allo scopo di fornire soluzioni di più facile impiego nel campo della pratica tecnica.Trattando quindi un caso di carattere generale, si mostra come i risultati ottenuti si prestano ad essere utilizzati per seguire la capacità portante della piastra ed il relativo meccanismo di collasso al variare della resistenza limite del suolo e della estensione del carico su di essa agente.Alla fine, a chiarimento del procedimento applicativo, si sviluppano alcuni esempi.


This work was supported by CNR (Italian Research Council).  相似文献   

7.
Summary The incremental boundary value problem for elastoplastic work-hardening continua, is discussed in the presence of large displacements but small strains, and allowing for distributed dislocations.Associated flow laws are assumed, as is the occurrence of interaction between yielding modes in singular points of the yield surface.Four extremum theorems are proved which reduce the problem to the optimization, under linear constraints, of convex quadratic forms (of a non linear form for the third theorem).Finally some results are specialized to small displacements and perfectly plastic (non hardening) cases.
Sommario Si discute della soluzione incrementale del problema elastoplastico con incrudimento nel campo dei grandi spostamenti (ma con piccole deformazioni intrinseche), ed in presenza di distorsioni distribuite.Si assumono leggi costitutive del tipo associato, ammettendo la possibilità di interazione tra i modi plastici in corrispondenza dei punti singolari della superficie di snervamento.Si dimostrano quattro proprietà di minimo della soluzione; esse riconducono il problema alla ottimizzazione di funzionali convessi (quadratici tranne che nel caso del terzo teorema) vincolata da equazioni e diseguaglianze lineari.Si specializzano infine alcuni risultati al caso dei piccoli spostamenti e della elastoplasticità senza incrudimento.


This study forms part of a research programme supported by the National (Italian) Research Council (C.N.R., Gruppo Plasticità).  相似文献   

8.
Piero Calvini 《Meccanica》1979,14(4):187-193
Summary The problem of determining stresses and displacements in a homogeneous rectangular plate rotating at constant angular speed about the barycentric axis, perpendicular to the plane of the plate, is considered. It is shown that this problem, if one assumes that the material is linear elastic with Poisson's ratio v, exhibits the same mathematical difficulties as the well-known problem of the rectangular plate subjected to parabolically-distributed tensile stress at two opposite sides.At difference from other methods employed in the past, a kind of solution based on biharmonic polynomials is proposed. This solution satisfies one boundary-condition exactly. The remaining three are satisfied in the mean. An evaluation of the error is performed.
Sommario Dal punto di vista della teoria dell'elasticità infinitesima si studia il problema della determinazione della tensione e della deformazione di una lamina rettangolare omogenea rotante con velocità angolare costante attorno all'asse baricentrico perpendicolare al piano della lamina stessa. Si dimostra che questo problema, se si vuole tenere conto del fattore di Poisson v, presenta le stesse difficoltà del ben noto problema della lamina rettangolare soggetta a forze con andamento parabolico su due lati opposti. In alternativa ai metodi già utilizzati in passato in casi come questo, si propone un tipo di soluzione basata su polinomi biarmonici che soddisfano esattamente solo una delle quattro condizioni al contorno. Le rimanenti tre sono soddisfatte in media. Si effettua una stima dell'errore.
  相似文献   

9.
Summary By means of mixing-length and turbulent Prandtl number hypothesis we solved the problem of parallel turbulent flow at constant density, both from the dynamic and thermal point of view; we then analyzed the fit with experimental data of various mixing-length formulas, and also the dependence of temperature profiles on the value of the turbulent Prandtl number.This critical analysis allowed the choice of the most suitable mixing-length formula and the value for the turbulent Prandtl number. On the basis of these results we extended the study discarding the condition of constant density; in particular we considered the case of liquids whose density was taken dependent on temperature changes across the walls, but independent of the pressure changes in flow direction.The study belongs to the case of fully developed temperature and velocity profiles.
Sommario Mediante l'introduzione del numero di Prandtl turbolento e della lunghezza di miscelamento nelle equazioni di Reynolds, viene risolto il problema della distribuzione di velocità e temperatura in un fluido a densità costante, in moto turbolento tra due piani paralleli. Le distribuzioni di velocità, ottenute con diverse espressioni della lunghezza di miscelamento, vengono poi confrontate con i dati sperimentali, allo scopo di scegliere la più opportuna di queste lunghezze; infine viene esaminata l'influenza del numero di Prandtl turbolento sulla distribuzione di temperatura.In accordo con le suddette scelte, lo studio è successivamente esteso al caso di densità dipendente dalla sola temperatura, ritenendo trascurabili le variazioni di densità per effetto del gradiente di pressione. In altri termini si limita lo studio ai liquidi.Tutti i risultati ottenuti si riferiscono a moti stabilizzati in velocità e temperatura.
  相似文献   

10.
Summary We deal with the problem of time dependent stresses and, particularly, of the contact pressure in two thick-walled coaxial shrink-fitted cylinders, when neither of the two elements is rigid.We apply Rabotnov's phenomenological theory of creep in such a way that the integral equation, which gives the time dependent ratio between the contact pressure at a given moment and the starting pressure, is linear. We assume no volume deformation and a plane state of strain.We give the explicit solution in three examples both for an elastic inner cylinder and for a viscoelastic one, solving the integral equations of the problem both by quadratures and by the Laplace transformation inverting the transforms by a numerical method[21].
Sommario Si affronta il problema della ricerca dell'andamento nel tempo delle tensioni e, in particolar modo, della pressione di contatto in un accoppiamento forzato tra due cilindri cavi coassiali nel caso che nessuno dei due elementi sia rigido.Si applica la teoria fenomenologica per lo scorrimento viscoso del Rabotnov in un modo tale per cui l'equazione integrale, che dà l'andamento nel tempo del rapporto tra la pressione di contatto in un dato istante e la pressione iniziale, è lineare. Si considera nulla la variazione di volume e piano lo stato di deformazione.Si dà la soluzione esplicita in tre esempi per cilindro interno sia elastico sia viscoelastico, risolvendo le equazioni integrali del problema sia per quadrature sia con l'uso della trasformata di Laplace trovando l'antitrasformata con un particolare metodo numerico di inversione[21].

  相似文献   

11.
Summary The equilibrium configuration of a thin plate under normal pressure is studied. A variational non-linear treatment of the problem is considered, taking into account the bending stresses and allowing large deflections.Existence, uniqueness, and regularity of solutions are obtained.
Sommario Si studia il problema non lineare della configurazione di equilibrio di una piastra sottile caricata normalmente nel suo piano medio, tenendo conto delle tensioni di natura flessionale e in un regime di spostamenti moderatamente grandi.Con una tecnica variazionale si studiano l'esistenza, l'unicità, e la regolarità delle soluzioni del problema.
  相似文献   

12.
Summary Solution existence, uniqueness and regularity are analyzed for the elastic plate on unilateral elastic foundation problem. Some results are also given for the one-sided rigid foundation case.
Sommario Si studiano l'esistenza, l'unicità e la regolarità delle soluzioni del problema della piastra elastica su fondazione elastica unilaterale. Sono anche dati alcuni risultati per il caso di fondazione rigida unilaterale.
  相似文献   

13.
Modal coupling in the free nonplanar finite motion of an elastic cable   总被引:3,自引:0,他引:3  
Summary In the finite motions of a suspended elastic cable the in-plane and out-of-plane oscillations are coupled, which is in contrast with what is predicted by the theory of small oscillations. To study the phenomenon of nonlinear coupling, a simple but meaningful two degree-of-freedom model is referred here, one parameter being used to describe the in-plane motion and the other the out-of-plane motion. The solution of the dynamic equilibrium equations is accomplished by an order-three perturbational expansion, which furnishes the time solution of the two displacement parameters. The modification of the free oscillations due to the exchange of energy between the two modes in absence of internal resonance is studied for different initial conditions and the effect of modal coupling is evidenced.
Sommario Nelle oscillazioni di grande ampiezza di un cavo sospeso, il moto nel piano e fuori del piano risulta accoppiato, a differenza di quanto predetto dalla teoria delle piccole oscillazioni. Questo problema viene studiato facendo riferimento ad un modello del cavo, semplice ma significativo, a due soli gradi di libertà, dei quali uno tiene conto del moto pendolare e l'altro del moto nel piano. La soluzione delle equazioni di moto è ottenuta con una tecnica perturbativa fino al terzo ordine, adatta al problema con nonlinearità quadratiche e cubiche. Si studia la modificazione della legge del moto dovuta al trasferimento di energia tra i due modi per differenti condizioni iniziali in assenza di risonanza interna e si valutano gli effetti dell'ac-coppiamento modale nel problema nonlineare.
  相似文献   

14.
Summary In this paper we present some results from a model of Shanley type, used as a tool to investigate anomalous computed response behavior of a pin-ended beam subjected to a short pulse of transverse force. The nature of critical load values is here studied, using both exact solutions and an approximate step-by-step numerical integration of standard type. It is shown that convergence with decreasing time step may fail, near a critical value. The significance of these results for the original beam computational problem, and for other physical structures, is indicated briefly.
Sommario Il presente lavoro illustra alcuni risultati ottenuti mediante un modello alla Shanley, usato per analizzare anomalie ottenute numericamente nel comportamento di una trave incernierata soggetta a un carico trasversale, applicato nella forma di un breve impulso. Viene studiata la natura di alcuni valori critici del carico, facendo uso sia di soluzioni esatte che di metodi approssimati di integrazione numerica. Si mostra che al crescere della accuratezza può non esservi convergenza dei risultati, e che vi è anzi divergenza per alcuni valori critici del passo di integrazione. Si discute brevemente l'interesse di tali risultati nei confronti delle difficoltà numeriche ottenute studiando il problema originale di trave incernierata.
  相似文献   

15.
An analytical approach to the dynamics of rotating shafts   总被引:1,自引:0,他引:1  
An analytical procedure, based on the dynamic stiffness method, is proposed for the study of rotor dynamics problems. On the grounds of the governing differential equations of a continuous beam, the dynamic stiffness matrix of the rotating Timoshenko beam is derived, including the effects of translational and rotational inertia, gyroscopic moments and shear deformation of the shaft. Concentrated disks and isotropic, elastic bearings are taken into account in the element formulation. The results obtained by the proposed method are compared with both classical closed form solutions and finite element analyses taken from the literature.
Sommario Nel presente lavoro viene esposto un procedimento analitico, basato sul metodo della rigidezza dinamica, per lo studio del comportamento dinamico dei rotori. In particolare, è ricavata la matrice di rigidezza dinamica della trave di Timoshenko rotante, che consente di tenere conto degli effetti dell'inerzia alla traslazione e alla rotazione, del momento giroscopico distribuito e della deformazione a taglio della trave. Nell'ambito della formulazione proposta, è possibile considerare, mediante inserimento diretto nella matrice di rigidezza dinamica, gli effetti dovuti alla presenza di dischi e di supporti elastici isotropi. Il metodo proposto è applicato alla determinazione delle velocità critiche di sistemi rotanti, con riferimento sia a casi semplici, per i quali è nota la soluzione analitica, sia al caso piu' complesso di un rotore già studiato in letteratura; per quest'ultimo esempio, in particolare, è presentato un confronto con i risultati delle analisi agli elementi finiti svolte da altri autori.
  相似文献   

16.
Summary A partial differential equation is considered whose dominant linear term involves third order time derivatives and which contains a nonlinear term multiplied by a small parameter. The solutions of initial-boundary-value problems for the equation are investigated using an extension of the Krylov-Bogoliubov-Mitropolskii asymptotic method. The vibrations of a viscoelastic rod lead to an equation of the type considered when the material constitutive relation is that of a standard linear solid together with additional small nonlinearities of behaviour, and this particular physical problem is investigated in detail. It is shown that under appropriate initial conditions only a single mode of vibration is excited to lowest order, and some numerical results are given which illustrate the effects of the nonlinearity on these single mode solutions.
Sommario Viene considerata una equazione differenziale alle derivate parziali il cui termine lineare dominante implica derivate del terzo ordine del tempo e che contiene un termine non lineare moltiplicato per un parametro di piccola entità. Vengono studiate le soluzioni dell'equazione con le condizioni iniziali e al contorno usando una estensione del metodo asintotico di Krylov-Bogoliubov-Mitropolskii. Le vibrazioni in una sbarra viscoelastica portano ad una equazione del tipo considerato quando la soluzione costitutiva è quella di un normale solido lineare con l'aggiunta di piccole non-linearità di comportamento. Viene studiato nel dettaglio questo particolare problema fisico. Viene mostrato che, con appropriate condizioni iniziali, viene eccitato solo il primo modo di vibrare. Vengono inoltre forniti alcuni risultati numerici che illustrano gli effetti della non linearità per questo particolare modo.
  相似文献   

17.
Milvio Capovani 《Meccanica》1968,3(4):252-258
Summary Amplitudes of linear vibrations of a slightly unbalanced shaft, rotating on lubricated bearings near one of its critical speeds, are determined, taking into account the damping action of lubricant forces and couples on the journal movement. The analysis (carried out on the lines of earlier work by Sternlicht, Warner and others) shows that the effect of lubricant couples is never relevant in the normal range of variability of parameters. The analysis puts instead in evidence the effect of the total damping action of the lubricant and hence the importance in design of forecasts of amplitude of vibration at running speed. The attention s often limited to forecasts of critical speeds; such procedure may be misleading particularly near higher critical speeds.
Sommario Si studiano le vibrazioni lineari d'un asse leggermente sbilanciato che ruota su cuscinetti lubrificati, tenendo conto dell'azione delle forze e momenti generati dalla presenza del lubrificante. L'analisi (sviluppata sulla traccia di precedenti lavori di Sternlicht, Warner ed altri) mentre mostra che l'effetto dei momenti è trascurabile nel normale campo di variabilità dei parametri, mette invece in risalto l'effetto ammortizzante totale del lubrificante. Appare così importante, in sede di progetto di dispositivi ruotanti, la determinazione delle ampiezze di vibrazione alle varie velocità di rotazione piuttosto che delle sole velocità critiche.

  相似文献   

18.
Michele Capurso 《Meccanica》1972,7(4):236-247
Summary This paper discusses the dynamics of rigid-viscoplastic bodies, bearing in mind large displacement effects. The extremal properties of the incremental solution of the dynamic problem are demonstrated. It is shown that the true instantaneous acceleration field minimizes some quadratic functionals in a class of functions constrained by linear equations and inequalities.It is also shown that, by idealizing the body as an assemblage of finite elements, the above principles reduce the problem to the solution of a classical quadratic programming problem that can be solved by several well-known efficient algorithms. Thus a general purpose finite element program for the dynamic analysis of rigid-viscoplastic structures can be derived by using the concepts shown in this paper.
Sommario Viene discusso il problema della dinamica delle strutture rigido plastiche nel campo dei grandi spostamenti. Vengono dimostrate le proprietà estremali della soluzione incrementale dimostrando che il campo di accelerazioni minimizza alcuni funzionali quadratici in una classe di funzioni vincolata dal rispetto di opportune condizioni di eguaglianza e disegnaglianza lineare. Discretizzando il corpo col metodo degli elementi finiti, i suddetti principi riportano il problema nel campo della programmazione quadratica e consentono quindi di utilizzare i ben noti ed efficienti algoritmi numerici di soluzione dei detti problemi. Si sottolinea quindi come, avvalendosi dei concetti esposti in questo lavoro, si possa pervenire a programmi di calcolo numerico del tutto generali per la soluzione del problema in questione.


The results presented in this paper were obtained in the course of research sponsored by the National (Italian) Research Council (C.N.R.).  相似文献   

19.
Robert S. D. Thomas 《Meccanica》1991,26(2-3):111-115
In this paper, we obtain the modes and velocities of acceleration waves on a thin hyperelastic shell in terms of the second fundamental form, which represents the geometrical properties of the shell, and of seven elastic moduli derived from the velocities in a plate of the same material. Some examples are studied, and approximations obtained in the case of a shallow shell.
Sommario In questo lavoro si ottengono i modi e le velocità delle onde di accelerazione in una volta sottile iperelastica, con riferimento alla seconda forma fondamentale che rappresenta le proprietà geometriche della volta e a sette moduli elastici derivati dalle velocità in una piastra dello stesso materiale. Si studiano alcuni esempi e si presentano soluzioni approssimate nel caso di una volta ribassata.
  相似文献   

20.
Summary Here we study impact-initiated disturbances in fluid-filled elastic tubes. The undeformed tube diameter, wall thickness, and elastic modulus of the tube material are assumed to be functions of the distance along the centre line of the tube. The linearized version of the governing equations are solved by the Laplace transform, which is inverted by means of an approximate method. The original non-linear system of governing equations is solved numerically by the method of characteristics. Relationships between the axial fluid velocity and axial coordinate as well as between velocity and time are displayed for fixed values of time and axial coordinate respectively for the linear and nonlinear theory for ease of comparison.
Sommario Si studiano le perturbazioni prodotte da impatto in tubi elastici riempiti di fluidi. Si assume che il diametro del tubo indeformato, lo spessore della parete e il modulo elastico del materiale del tubo siano funzioni della distanza misurata lungo la linea centrale del tubo. La versione linearizzata delle equazioni che governano il fenomeno è risolta con la trasformata di Laplace, invertita con un metodo approssimato. Il sistema originale non lineare delle equazioni è risolto numericamente con il metodo delle caratteristiche. Vengono rappresentate, per valori fissati del tempo e della coordinata assiale, le relazioni tra la velocità assiale del fluido e la coordinata assiale oltre alle relazioni fra velocità e tempo, sia per la teoria lineare sia per quella non lineare.

  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号