首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.  相似文献   

2.
Fuji film has been widely used in studies aimed at obtaining the contact mechanics of articular joints. Once sealed for practical use in biological joints, Fuji Pressensor film has a total effective thickness of 0.30 mm, which is comparable to the cartilage thickness in the joints of many small animals. The average effective elastic modulus of Fuji film is approximately 100 MPa in compression, which is larger by a factor of 100-300 compared to that of normal articular cartilage. Therefore, inserting a Pressensor film into an articular joint will change the contact mechanics of the joint. The measurement precision of the Pressensor film has been determined systematically; however, the changes in contact mechanics associated with inserting the film into joints have not been investigated. This study was aimed at quantifying the changes in the contact mechanics associated with inserting sealed Fuji Pressensor film into joints. Spherical and cylindrical articular joint contact mechanics with and without Pressensor film and for varying degrees of surface congruency were analyzed and compared by using finite element models. The Pressensor film was taken as linearly elastic and the cartilage was assumed to be biphasic, composed of a linear elastic solid phase and an inviscid fluid phase. The present analyses showed that measurements of the joint contact pressures with Fuji Pressensor film will change the maximum true contact pressures by 10-26 percent depending on the loading, geometry of the joints, and the mechanical properties of cartilage. Considering this effect plus the measurement precision of the film (approximately 10 percent), the measured joint contact pressures in a joint may contain errors as large as 14-28 percent.  相似文献   

3.
Osteoarthritis (OA) is a multi-factor disorder of sinovial joints, which characterized by escalated degeneration and loss of articular cartilage. Treatment of OA is a critical unmet need in medicine for regeneration of damaged articular cartilage in elderly. On the other hand, lubricin, a glycoprotein specifically synthesized by chondrocytes located at the surface of articular cartilage, has been shown to provide boundary lubrication of congruent articular surfaces under conditions of high contact pressure and near zero sliding speed. Lubrication of these surfaces is critical to normal joint function, while different gene expressions of lubricin had been found in the synovium of rheumatoid arthritis (RA) and OA. Moreover, mutations or lacking of lubricin gene have been shown to link to the joint disease such as camptodactyly-arthropathy-coxa vara-pericarditis syndrome (CACP), synovial hyperplasia and failure of joint function, suggesting an important role of lubricin in the pathogenesis of these joint disease. Recent studies demonstrate that administration with recombinant lubricin in the joint cavity would be effective in the prevention of cartilage degeneration in animal OA models. Therefore, a treatment with lubricin which would protect cartilage in vivo would be desirable. This article reviews recent findings with regard to the possible role of lubricin in the progression of OA, and further discusses lubricin as a novel potential biotherapeutic approaches for the treatment of OA.  相似文献   

4.
Partial meniscectomy is believed to change the biomechanics of the knee joint through alterations in the contact of articular cartilages and menisci. Although fluid pressure plays an important role in the load support mechanism of the knee, the fluid pressurization in the cartilages and menisci has been ignored in the finite element studies of the mechanics of meniscectomy. In the present study, a 3D fibril-reinforced poromechanical model of the knee joint was used to explore the fluid flow dependent changes in articular cartilage following partial medial and lateral meniscectomies. Six partial longitudinal meniscectomies were considered under relaxation, simple creep, and combined creep loading conditions. In comparison to the intact knee, partial meniscectomy not only caused a substantial increase in the maximum fluid pressure but also shifted the location of this pressure in the femoral cartilage. Furthermore, these changes were positively correlated to the size of meniscal resection. While in the intact joint, the location of the maximum fluid pressure was dependent on the loading conditions, in the meniscectomized joint the location was predominantly determined by the site of meniscal resection. The partial meniscectomy also reduced the rate of the pressure dissipation, resulting in even larger difference between creep and relaxation times as compared to the case of the intact knee. The knee joint became stiffer after meniscectomy because of higher fluid pressure at knee compression followed by slower pressure dissipation. The present study indicated the role of fluid pressurization in the altered mechanics of meniscectomized knees.  相似文献   

5.
The purposes of this study were to determine the in situ functional and material properties of articular cartilage in an experimental model of joint injury, and to quantify the corresponding in situ joint contact mechanics. Experiments were performed in the anterior cruciate ligament (ACL) transected knee of the cat and the corresponding, intact contralateral knee, 16 weeks following intervention. Cartilage thickness, stiffness, effective Young’s modulus, and permeability were measured and derived from six locations of the knee. The total contact area and peak pressures in the patellofemoral joint were obtained in situ using Fuji Pressensor film, and comparisons between experimental and contralateral joint were made for corresponding loading conditions. Total joint contact area and peak pressure were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint. Articular cartilage thickness and stiffness were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint in the four femoral and patellar test locations. Articular cartilage material properties (effective Young’s modulus and permeability) were the same in the ACL-transected and intact joints. These results demonstrate for the first time the effect of changes in articular cartilage properties on the load transmission across a joint. They further demonstrate a substantial change in the joint contact mechanics within 16 weeks of ACL transection. The results were corroborated by theoretical analysis of the contact mechanics in the intact and ACL-transected knee using biphasic contact analysis and direct input of cartilage properties and joint surface geometry from the experimental animals. We conclude that the joint contact mechanics in the ACL-transected cat change within 16 weeks of experimental intervention.  相似文献   

6.
Abstract

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics.  相似文献   

7.
Biphasic properties of articular cartilage allow it to be an excellent bearing material and have been studied through several simplified experiments as well as finite element modelling. However, three-dimensional biphasic finite element (FE) models of the whole joint are rare. The current study was carried out to experimentally validate FE methodology for modelling hemiarthroplasty. Material properties such as equilibrium elastic modulus and permeability of porcine acetabular cartilage were initially derived by curve-fitting an experimental deformation curve with that obtained using FE. These properties were then used in the hemiarthroplasty hip joint modelling. Each porcine acetabular cup was loaded with 400N using a 34mm diameter CoCr femoral head. A specimen-specific FE model of each acetabular cup was created using μCT and a series of software processes. Each model was analysed under conditions similar to those tested experimentally. Contact stresses and contact areas predicted by the model, immediately after loading, were then compared with the corresponding experimentally measured values. Very high peak contact stresses (maximum experimental: 14.09MPa) were recorded. A maximum difference of 12.42% was found in peak contact stresses. The corresponding error for contact area was 20.69%. Due to a fairly good agreement in predicted and measured values of contact stresses and contact areas, the integrated methodology developed in this study can be used as a basis for future work. In addition, FE predicted total fluid load support was around 80% immediately after loading. This was lower than that observed in conforming contact problems involving biphasic cartilage and was due to a smaller local contact area and variable clearance making fluid exudation easier.  相似文献   

8.
When studying joint degeneration leading to osteoarthritis (OA), it seems imperative that local joint tissue loading is known during normal everyday movement and that the adaptive/degenerative effects of this loading are quantified systematically. Philosophically, we believe the best way to approach this problem is by studying joint degeneration and osteoarthritis in long-term experimental models and by representing diarthrodial joints and the associated tissues with accurate, geometric and structural, theoretical models. Here, we present selected examples of our work representing this approach. Experimentally, we demonstrate that the local loading of joints changes continuously in experimental models of OA, not only because of the changing external and internal loading, but also because of the continuous alterations in joint contact geometry and tissue mechanical properties. Furthermore, we show that single bouts of joint loading affect gene expression, and that gene expression, as well as subsequent joint degeneration is site-specific. In fact, opposing articular surfaces that are exposed to the same loading may degenerate at completely different rates. Finally, we propose a series of theoretical models of articular cartilage and contact mechanics, demonstrating that many of the anisotropic and inhomogeneous properties can be explained by structural elements and their orientation and volumetric concentration across the tissue.  相似文献   

9.
A model of synovial fluid (SF) filtration by articular cartilage (AC) in a step-loaded spherical synovial joint at rest is presented. The effects of joint pathology (such as a depleted acetabular labrum, a depleted cartilage superficial zone consistent with early osteoarthritis and an inflammatory SF) on the squeezed synovial film are also investigated. Biphasic mixture models for AC (ideal fluid and elastic porous transversely isotropic two-layer matrix) and for SF (ideal and thixotropic fluids) are applied and the following results are obtained. If the acetabular labrum is able to seal the pressurised SF between the articular surfaces (as in the normal hip joint), the fluid in the synovial film and in the cartilage within the labral ring is homogeneously pressurised. The articular surfaces remain separated by a fluid film for minutes. If the labrum is destroyed or absent and the SF can escape across the contact edge, the fluid pressure is non-homogeneous and with a small jump at the articular surface at the very moment of load application. The ensuing synovial film filtration by porous cartilage is lower for the normal cartilage (with the intact superficial zone) than if this zone is already depleted or rubbed off as in the early stage of primary osteoarthritis. Compared with the inflammatory (Newtonian) SF, the normal (thixotropic) fluid applies favourably in the squeezed film near the contact centre only, yielding a thicker SF film there, but not affecting the minimum thickness in the fluid film profile at a fixed time. For all that, in the unsealed case for both the normal and pathological joint, the macromolecular concentration of the hyaluronic acid-protein complex in the synovial film quickly increases due to the filtration in the greater part of the contact. A stable synovial gel film, thick on the order of 10(-7)m, protecting the articular surfaces from the intimate contact, is formed within a couple of seconds. Boundary lubrication by the synovial gel is established if sliding motion follows until a fresh SF is entrained into the contact. This theoretical prediction is open for experimental verifications.  相似文献   

10.
Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error ~0.5 mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.  相似文献   

11.
Boundary lubrication of articular cartilage by conformal, molecularly thin films reduces friction and adhesion between asperities at the cartilage-cartilage contact interface when the contact conditions are not conducive to fluid film lubrication. In this study, the nanoscale friction and adhesion properties of articular cartilage from typical load-bearing and non-load-bearing joint regions were studied in the boundary lubrication regime under a range of physiological contact pressures using an atomic force microscope (AFM). Adhesion of load-bearing cartilage was found to be much lower than that of non-load-bearing cartilage. In addition, load-bearing cartilage demonstrated steady and low friction coefficient through the entire load range examined, whereas non-load-bearing cartilage showed higher friction coefficient that decreased nonlinearly with increasing normal load. AFM imaging and roughness calculations indicated that the above trends in the nanotribological properties of cartilage are not due to topographical (roughness) differences. However, immunohistochemistry revealed consistently higher surface concentration of boundary lubricant at load-bearing joint regions. The results of this study suggest that under contact conditions leading to joint starvation from fluid lubrication, the higher content of boundary lubricant at load-bearing cartilage sites preserves synovial joint function by minimizing adhesion and wear at asperity microcontacts, which are precursors for tissue degeneration.  相似文献   

12.
Patellofemoral (PF)-related pathologies, including joint laxity, patellar maltracking, cartilage degradation and anterior knee pain, affect nearly 25% of the population. Researchers have investigated the influence of articular geometry on kinematics and contact mechanics in order to gain insight into the etiology of these conditions. The purpose of the current study was to create a three-dimensional statistical shape model of the PF joint and to characterize relationships between PF shape and function (kinematics and contact mechanics). A statistical shape model of the patellar and femoral articular surfaces and their relative alignment was developed from magnetic resonance images. Using 15 shape parameters, the model characterized 97% of the variation in the training set. The first three shape modes primarily described variation in size, patella alta-baja and depth of the sulcus groove. A previously verified finite element model was used to predict kinematics and contact mechanics for each subject. Combining the shape and joint mechanics data, a statistical shape-function model was developed that established quantitative relations of how changes in the shape of the PF joint influence mechanics. The predictive capability of the shape-function model was evaluated by comparing statistical model and finite element predictions, resulting in kinematic root mean square errors of less than 3° and 2.5 mm. The key results of the study are dually in the implementation of a novel approach linking statistical shape and finite element models and the relationships elucidated between PF articular geometry and mechanics.  相似文献   

13.
Meniscectomy is a significant risk factor for osteoarthritis, involving altered cell synthesis, central fibrillation, and peripheral osteophyte formation. Though changes in articular cartilage contact pressure are known, changes in tissue-level mechanical parameters within articular cartilage are not well understood. Recent imaging research has revealed the effects of meniscectomy on the time-dependent deformation of physiologically loaded articular cartilage. To determine tissue-level cartilage mechanics that underlie observed deformation, a novel finite element modeling approach using imaging data and a contacting indenter boundary condition was developed. The indenter method reproduces observed articular surface deformation and avoids assumptions about tangential stretching. Comparison of results from an indenter model with a traditional femur-tibia model verified the method, giving errors in displacement, solid and fluid stress, and strain below 1% (RMS) and 7% (max.) of the absolute maximum of the parameters of interest. Indenter finite element models using real joint image data showed increased fluid pressure, fluid exudation, loss of fluid load support, and increased tensile strains centrally on the tibial condyle after meniscectomy-patterns corresponding to clinical observations of cartilage matrix damage and fibrillation. Peripherally there was decreased consolidation, which corresponds to reduced contact and fluid pressure in this analysis. Clinically, these areas have exhibited advance of the subchondral growth front, biological destruction of the cartilage matrix, cartilage thinning, and eventual replacement of the cartilage via endochondral ossification. Characterizing the changes in cartilage mechanics with meniscectomy and correspondence with observed tissue-level effects may help elucidate the etiology of joint-level degradation seen in osteoarthritis.  相似文献   

14.
Abnormal, excessive stresses acting on articular joint surfaces are speculated to be one of the causes for joint degeneration. However, articular surface stresses have not been studied systematically, since it is technically difficult to measure in vivo contact areas and pressures in dynamic situations. Therefore, we implemented a numerical model of articular surface contact using accurate surface geometries. The model was developed for the cat patellofemoral joint. We demonstrated that small misalignments of the patella relative to the femur change the joint contact mechanics substantially for a given external load. These results suggest that misalignment might be studied as one of the factors causing articular cartilage disorder and joint degeneration.  相似文献   

15.
Non-invasive techniques for quantifying early biochemical and biomechanical changes in articular cartilage may provide a means of more precisely assessing osteoarthritis (OA) progression. The goals of this study were to determine the relationship between T1rho magnetic resonance (MR) imaging relaxation times and changes in cartilage composition, cartilage mechanical properties, and synovial fluid biomarker levels and to demonstrate the application of T1rho imaging to evaluate cartilage composition in human subjects in vivo. Femoral condyles and synovial fluid were harvested from healthy and OA porcine knee joints. Sagittal T1rho relaxation MR images of the condyles were acquired. OA regions of OA joints exhibited an increase in T1rho relaxation times as compared to non-OA regions. Furthermore in these regions, cartilage sGAG content and aggregate modulus decreased, while percent degraded collagen and water content increased. In OA joints, synovial fluid concentrations of sGAG decreased and C2C concentrations increased compared to healthy joints. T1rho relaxation times were negatively correlated with cartilage and synovial fluid sGAG concentrations and aggregate modulus and positively correlated with water content and permeability. Additionally, we demonstrated the application of these in vitro findings to the study of human subjects. Specifically, we demonstrated that walking results in decreased T1rho relaxation times, consistent with water exudation and an increase in proteoglycan concentration with in vivo loading. Together, these findings demonstrate that cartilage MR imaging and synovial fluid biomarkers provide powerful non-invasive tools for characterizing changes in the biochemical and biomechanical environments of the joint.  相似文献   

16.
The Phan-Thien and Tanner (PTT) model is one of the most widely used rheological models. It can properly describe the common characteristics of viscoelastic non-Newtonian fluids. There is evidence that synovial fluid in human joints, which also lubricates artificial joints, is viscoelastic. Modeling the geometry of the total hip replacement, the PTT model is applied in spherical coordinates for a thin confined fluid film. A modified Reynolds equation is developed for this geometry. Several simplified illustrative problems are solved. The effect of the edge boundary condition on load-carrying normal stress is discussed. Solutions are also obtained for a simple squeezing flow. The effect of both the relaxation time and the PTT shear parameter is to reduce the load relative to a Newtonian fluid with the same viscosity. This implies that the Newtonian model is not conservative and may overpredict the load capacity. The PTT theory is a good candidate model to use for joint replacement lubrication. It is well regarded and derivable from molecular considerations. The most important non-Newtonian characteristics can be described with only three primary material parameters.  相似文献   

17.
The objective of this study is to establish and verify the set of boundary conditions at the interface between a biphasic mixture (articular cartilage) and a Newtonian or non-Newtonian fluid (synovial fluid) such that a set of well-posed mathematical problems may be formulated to investigate joint lubrication problems. A "pseudo-no-slip" kinematic boundary condition is proposed based upon the principle that the conditions at the interface between mixtures or mixtures and fluids must reduce to those boundary conditions in single phase continuum mechanics. From this proposed kinematic boundary condition, and balances of mass, momentum and energy, the boundary conditions at the interface between a biphasic mixture and a Newtonian or non-Newtonian fluid are mathematically derived. Based upon these general results, the appropriate boundary conditions needed in modeling the cartilage-synovial fluid-cartilage lubrication problem are deduced. For two simple cases where a Newtonian viscous fluid is forced to flow (with imposed Couette or Poiseuille flow conditions) over a porous-permeable biphasic material of relatively low permeability, the well known empirical Taylor slip condition may be derived using matched asymptotic analysis of the boundary layer at the interface.  相似文献   

18.
Quantitative computed tomography-based finite element models of proximal femora must be validated with cadaveric experiments before using them to assess fracture risk in osteoporotic patients. During validation, it is essential to carefully assess whether the boundary condition (BC) modeling matches the experimental conditions. This study evaluated proximal femur stiffness results predicted by six different BC methods on a sample of 30 cadaveric femora and compared the predictions with experimental data. The average stiffness varied by 280% among the six BCs. Compared with experimental data, the predictions ranged from overestimating the average stiffness by 65% to underestimating it by 41%. In addition, we found that the BC that distributed the load to the contact surfaces similar to the expected contact mechanics predictions had the best agreement with experimental stiffness. We concluded that BC modeling introduced large variations in proximal femora stiffness predictions.  相似文献   

19.
The stiffness and hydraulic permeability of soft contact lenses may influence its clinical performance, e.g., on-eye movement, fitting, and wettability, and may be related to the occurrence of complications; e.g., lesions. It is therefore important to determine these properties in the design of comfortable contact lenses. Micro-indentation provides a nondestructive means of measuring mechanical properties of soft, hydrated contact lenses. However, certain geometrical and material considerations must be taken into account when analyzing output force-displacement (F-D) data. Rather than solely having a solid response, mechanical behavior of hydrogel contact lenses can be described as the coupled interaction between fluid transport through pores and solid matrix deformation. In addition, indentation of thin membranes ( approximately 100 microm) requires special consideration of boundary conditions at lens surfaces and at the indenter contact region. In this study, a biphasic finite element model was developed to simulate the micro-indentation of a hydrogel contact lens. The model accounts for a curved, thin hydrogel membrane supported on an impermeable mold. A time-varying boundary condition was implemented to model the contact interface between the impermeable spherical indenter and the lens. Parametric studies varying the indentation velocities and hydraulic permeability show F-D curves have a sensitive region outside of which the force response reaches asymptotic limits governed by either the solid matrix (slow indentation velocity, large permeability) or the fluid transport (high indentation velocity, low permeability). Using these results, biphasic properties (Young's modulus and hydraulic permeability) were estimated by fitting model results to F-D curves obtained at multiple indentation velocities (1.2 and 20 microm/s). Fitting to micro-indentation tests of Etafilcon A resulted in an estimated permeability range of 1.0 x 10(-15) to 5.0 x 10(-15) m(4)N s and Young's modulus range of 130 to 170 kPa.  相似文献   

20.
The transmission of load through the human hip joint   总被引:8,自引:0,他引:8  
This paper describes the results of loading experiments carried out on human hip joints. The unloaded surfaces of the femoral head and the acetabulum are slightly incongruous. The location and magnitude of the contact areas between the surfaces therefore depend on the magnitude and direction of the applied load. The contact areas were determined experimentally for a variety of loads typical of normal walking. Two distinct contact areas were found on the anterior and posterior aspects of the acetabulum at light loads, gradually merging with increasing load until, at a certain transition load, the dome of the acetabulum comes into contact and contact is then complete. The value of the transition load depends on the rate of loading, due to creep of the cartilage, and was found to vary from 50 per cent of body weight in young specimens to 25 per cent of body weight for elderly specimens for rates of loading typical of normal walking. Thus, the dome of the acetabulum is out of contact for a substantial portion of the swing phase of normal walking.

The analysis of a much simplified model of the hip joint is presented. The dependence of contact area on load is demonstrated, but also a method of determining the transition load for complete contact from the load/deflection relation for the hip is suggested. The values of the transition load quoted above were obtained by this method. The analysis further indicates that the distribution of pressure between the articular surfaces depends critically on the distribution of cartilage thickness throughout the joint. It is suggested that the distribution of cartilage thickness is such as to lead to a state of uniform pressure at the upper end of the physiological load range. Some experimental evidence is presented in support of this suggestion.

It is concluded that the function of joint incongruity is to allow the articular surfaces to come out of contact at light loads so that the cartilage may be exposed to synovial fluid for the purposes of nutrition and lubrication. At large loads, the distribution of cartilage thickness ensures that a state of hydrostatic pressure is achieved in order that cartilage, with a large fluid content, may transmit large pressures without flow and consequent loss of its integrity.  相似文献   


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

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