首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The accurate measurement of local strain is necessary to study bone mechanics and to validate micro computed tomography (µCT) based finite element (FE) models at the tissue scale. Digital volume correlation (DVC) has been used to provide a volumetric estimation of local strain in trabecular bone sample with a reasonable accuracy. However, nothing has been reported so far for µCT based analysis of cortical bone. The goal of this study was to evaluate accuracy and precision of a deformable registration method for prediction of local zero-strains in bovine cortical and trabecular bone samples. The accuracy and precision were analyzed by comparing scans virtually displaced, repeated scans without any repositioning of the sample in the scanner and repeated scans with repositioning of the samples.  相似文献   

2.
Subchondral stiffening is a hallmark pathologic feature of osteoarthritis but its mechanical and temporal relationship to the initiation or the progression of osteoarthritis is not established. The mechanical effect of subchondral stiffening on the surrounding trabecular bone is poorly understood. This study employs a relatively new application of digital image correlation to measure strain in the trabecular region of the proximal medial tibia in normal specimens and in specimens with simulated subchondral bone stiffening. Coronal sections from eight normal human cadaveric proximal tibiae were loaded in static compression and high resolution contact radiographs were made. Repeat contact radiographs were collected after the subchondral bone near the jointline was stiffened using polymethylmethacrylate. Digital images, made from loaded and unloaded contact radiographs, were compared using in-house software to measure trabecular displacement and calculate trabecular bone strain. Overall strain was higher in the stiffened specimens suggesting experimental artifiact significantly affected our results. Consistent increases in median maximum shear strain, median maximum principal strain, median minimum principal strain, and peak shear strain were measured near the inner and outer edges of the stiffened segment. Our experiment provides direct experimental measurement of increases in trabecular bone strain caused by subchondral stiffening, however, the clinical and biologic importance of strain increases is unknown.  相似文献   

3.
It is generally accepted that the strength and stiffness of trabecular bone is strongly affected by trabecular microstructure. It has also been hypothesized that stress induced adaptation of trabecular bone is affected by trabecular tissue level stress and/or strain. At this time, however, there is no generally accepted (or easily accomplished) technique for predicting the effect of microstructure on trabecular bone apparent stiffness and strength or estimating tissue level stress or strain. In this paper, a recently developed mechanics theory specifically designed to analyze microstructured materials, called the homogenization theory, is presented and applied to analyze trabecular bone mechanics. Using the homogenization theory it is possible to perform microstructural and continuum analyses separately and then combine them in a systematic manner. Stiffness predictions from two different microstructural models of trabecular bone show reasonable agreement with experimental results, depending on metaphyseal region, (R2 greater than 0.5 for proximal humerus specimens, R2 less than 0.5 for distal femur and proximal tibia specimens). Estimates of both microstructural strain energy density (SED) and apparent SED show that there are large differences (up to 30 times) between apparent SED (as calculated by standard continuum finite element analyses) and the maximum microstructural or tissue SED. Furthermore, a strut and spherical void microstructure gave very different estimates of maximum tissue SED for the same bone volume fraction (BV/TV). The estimates from the spherical void microstructure are between 2 and 20 times greater than the strut microstructure at 10-20% BV/TV.  相似文献   

4.
Bone strain is the governing stimuli for the remodeling process necessary in the maintenance of bone's structure and mechanical strength. Strain gages are the gold standard and workhorses of human bone experimental strain analysis in vivo. The objective of this systematic literature review is to provide an overview for direct in vivo human bone strain measurement studies and place the strain results within context of current theories of bone remodeling (i.e. mechanostat theory). We employed a standardized search strategy without imposing any time restriction to find English language studies indexed in PubMed and Web of Science databases that measured human bone strain in vivo. Twenty-four studies met our final inclusion criteria. Seven human bones were subjected to strain measurements in vivo including medial tibia, second metatarsal, calcaneus, proximal femur, distal radius, lamina of vertebra and dental alveolar. Peak strain magnitude recorded was 9096 με on the medial tibia during basketball rebounding and the peak strain rate magnitude was -85,500 με/s recorded at the distal radius during forward fall from standing, landing on extended hands. The tibia was the most exposed site for in vivo strain measurements due to accessibility and being a common pathologic site of stress fracture in the lower extremity. This systematic review revealed that most of the strains measured in vivo in different bones were generally within the physiological loading zone defined by the mechanostat theory, which implies stimulation of functional adaptation necessary to maintain bone mechanical integrity.  相似文献   

5.
Dependence of yield strain of human trabecular bone on anatomic site   总被引:17,自引:0,他引:17  
Understanding the dependence of human trabecular bone strength behavior on anatomic site provides insight into structure-function relationships and is essential to the increased success of site-specific finite element models of whole bones. To investigate the hypothesis that the yield strains of human trabecular bone depend on anatomic site, the uniaxial tensile and compressive yield properties were compared for cylindrical specimens from the vertebra (n=61), proximal tibia (n=31), femoral greater trochanter (n=23), and femoral neck (n=27) taken from 61 donors (67+/-15years). Test protocols were used that minimized end artifacts and loaded specimens along the main trabecular orientation. Yield strains by site (mean+/-S.D.) ranged from 0.70+/-0.05% for the trochanter to 0.85+/-0.10% for the femoral neck in compression, from 0.61+/-0.05% for the trochanter to 0.70+/-0.05% for the vertebra in tension, and were always higher in compression than tension (p<0.001). The compressive yield strain was higher for the femoral neck than for all other sites (p<0.001), as was the tensile yield strain for the vertebra (p<0.007). Analysis of covariance, with apparent density as the covariate, indicated that inter-site differences existed in yield stress even after adjusting statistically for density (p<0.035). Coefficients of variation in yield strain within each site ranged from only 5-12%, consistent with the strong linear correlations (r(2)=0.94-0.98) found between yield stress and modulus. These results establish that the yield strains of human trabecular bone can differ across sites, but that yield strain may be considered uniform within a given site despite substantial variation in elastic modulus and yield stress.  相似文献   

6.
Rat and mouse femur and tibia fracture calluses were collected over various time increments of healing. Serial sections were produced at spatial segments across the fracture callus. Standard histological methods and in situ hybridization to col1a1 and col2a1 mRNAs were used to define areas of cartilage and bone formation as well as tissue areas undergoing remodeling. Computer-assisted reconstructions of histological sections were used to generate three-dimensional images of the spatial morphogenesis of the fracture calluses. Endochondral bone formation occurred in an asymmetrical manner in both the femur and tibia, with cartilage tissues seen primarily proximal or distal to the fractures in the respective calluses of these bones. Remodeling of the calcified cartilage proceeded from the edges of the callus inward toward the fracture producing an inner-supporting trabecular structure over which a thin outer cortical shell forms. These data suggest that the specific developmental mechanisms that control the asymmetrical pattern of endochondral bone formation in fracture healing recapitulated the original asymmetry of development of a given bone because femur and tibia grow predominantly from their respective distal and proximal physis. These data further show that remodeling of the calcified cartilage produces a trabecular bone structure unique to fracture healing that provides the rapid regain in weight-bearing capacity to the injured bone.  相似文献   

7.
The combination of three-dimensional (3-D) models with dual fluoroscopy is increasingly popular for evaluating joint function in vivo. Applying these modalities to study knee motion with high accuracy requires reliable anatomical coordinate systems (ACSs) for the femur and tibia. Therefore, a robust method for creating ACSs from 3-D models of the femur and tibia is required. We present and evaluate an automated method for constructing ACSs for the distal femur and proximal tibia based solely on 3-D bone models. The algorithm requires no observer interactions and uses model cross-sectional area, center of mass, principal axes of inertia, and cylindrical surface fitting to construct the ACSs. The algorithm was applied to the femur and tibia of 10 (unpaired) human cadaveric knees. Due to the automated nature of the algorithm, the within specimen variability is zero for a given bone model. The algorithm’s repeatability was evaluated by calculating variability in ACS location and orientation across specimens. Differences in ACS location and orientation between specimens were low (<1.5 mm and <2.5°). Variability arose primarily from natural anatomical and morphological differences between specimens. The presented algorithm provides an alternative method for automatically determining subject-specific ACSs from the distal femur and proximal tibia.  相似文献   

8.
This study investigated the effects of green tea polyphenols (GTP) and alfacalcidol on bone microstructure and strength along with possible mechanisms in rats with chronic inflammation. A 12-week study using a 2 (no GTP vs. 0.5%, w/v GTP in drinking water)×2 (no alfacalcidol vs. 0.05 μg/kg alfacalcidol orally, 5×/week) factorial design was employed in lipopolysaccharide (LPS)-administered female rats. A group receiving placebo administration was used to compare with a group receiving LPS administration only to evaluate the effect of LPS. Changes in tibial and femoral microarchitecture and strength of femur were evaluated. Difference in expression of tumor necrosis factor-α (TNF-α) in proximal tibia using immunohistochemistry was examined. Compared to the placebo group, the LPS-administered-only group had significantly lower femoral mass, trabecular volume, thickness and number in proximal tibia and femur, and lower periosteal bone formation rate in tibial shafts but had significantly higher trabecular separation and osteoclast number in proximal tibia and eroded surface in endocortical tibial shafts. Both GTP and alfacalcidol reversed these LPS-induced detrimental changes in femur, proximal tibia and endocortical tibial shaft. Both GTP and alfacalcidol also significantly improved femoral strength, while significantly suppressed TNF-α expression in proximal tibia. There were significant interactions in femoral mass and strength, trabecular separation, osteoclast number and TNF-α expression in proximal tibia. A combination of both showed to sustain bone microarchitecture and strength. We conclude that a protective impact of GTP and alfacalcidol in bone microarchitecture during chronic inflammation may be due to a suppression of TNF-α.  相似文献   

9.
In this study, we developed a numerical framework that computationally determines simultaneous and interactive structural changes of cortical and trabecular bone types during bone remodeling, and we investigated the structural correlation between the two bone types in human proximal femur. We implemented a surface remodeling technique that performs bone remodeling in the exterior layer of the cortical bone while keeping its interior area unchanged. A micro-finite element (μFE) model was constructed that represents the entire cortical bone and full trabecular architecture in human proximal femur. This study simulated and compared the bone adaptation processes of two different structures: (1) femoral bone that has normal cortical bone shape and (2) perturbed femoral bone that has an artificial bone lump in the inferomedial cortex. Using the proposed numerical method in conjunction with design space optimization, we successfully obtained numerical results that resemble actual human proximal femur. The results revealed that actual cortical bone, as well as the trabecular bone, in human proximal femur has structurally optimal shapes, and it was also shown that a bone abnormality that has little contribution to bone structural integrity tends to disappear. This study also quantitatively determined the structural contribution of each bone: when the trabecular adaptation was complete, the trabecular bone supported 54% of the total load in the human proximal femur while the cortical bone carried 46%.  相似文献   

10.
Digital Volume Correlation (DVC) is used to measure internal displacements and strains in bone. Recent studies have shown that Synchrotron radiation micro-computed tomography (SR-microCT) can improve the accuracy and precision of DVC. However, only zero-strain or virtually-moved test have been used to quantify the DVC uncertainties, leading to potential underestimation of the measurement errors.In this study, for the first time, the uncertainties of a global DVC approach have been evaluated on repeated SR-microCT scans of bovine cortical bone (voxel size: 1.6 μm), which were virtually deformed for different magnitudes and along different directions.The results showed that systematic and random errors of the normal strain components along the deformation direction were higher than the errors along unstrained directions. The systematic percentage errors were smaller for larger virtual deformations. The random percentage error was in the order of 10% of the virtual deformation. However, higher errors were localized at the boundary of the volumes of interest, perpendicular to the deformation direction. When only the central region of the samples was considered (100 µm layers removed from the borders where the deformation was applied), the errors in the direction of virtual deformation were comparable to the errors in the unstrained directions.In conclusion, the method presented to estimate the uncertainties of DVC is suitable for testing anisotropic specimens as cortical bone. The good agreement between the uncertainties in measurements of strain components obtained with this approach and with the simpler zero-strain-test suggests that the latter is adequate in the tested deformation scenarios.  相似文献   

11.
Dual energy X-ray absorptiometry (DXA) is the standard for assessing fragility fracture risk using areal bone mineral density (aBMD), but only explains 60–70% of the variation in bone strength. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides 3D-measures of bone microarchitecture and volumetric bone mineral density (vBMD), but only at the wrist and ankle. Finite element (FE) models can estimate bone strength with 86–95% precision. The purpose of this study is to determine how well vBMD and FE bone strength at the wrist and ankle relate to fracture strength at the hip and spine, and to compare these relationships with DXA measured directly at those axial sites. Cadaveric samples (radius, tibia, femur and L4 vertebra) were compared within the same body. The radius and tibia specimens were assessed using HR-pQCT to determine vBMD and FE failure load. aBMD from DXA was measured at the femur and L4 vertebra. The femur and L4 vertebra specimens were biomechanically tested to determine failure load. aBMD measures of the axial skeletal sites strongly correlated with the biomechanical strength for the L4 vertebra (r = 0.77) and proximal femur (r = 0.89). The radius correlated significantly with biomechanical strength of the L4 vertebra for vBMD (r = 0.85) and FE-derived strength (r = 0.72), but not with femur strength. vBMD at the tibia correlated significantly with femoral biomechanical strength (r = 0.74) and FE-estimated strength (r = 0.83), and vertebral biomechanical strength for vBMD (r = 0.97) and FE-estimated strength (r = 0.91). The higher correlations at the tibia compared to radius are likely due to the tibia’s weight-bearing function.  相似文献   

12.
Trabecular bone modulus-density relationships depend on anatomic site   总被引:6,自引:0,他引:6  
One outstanding issue regarding the relationship between elastic modulus and density for trabecular bone is whether the relationship depends on anatomic site. To address this, on-axis elastic moduli and apparent densities were measured for 142 specimens of human trabecular bone from the vertebra (n=61), proximal tibia (n=31), femoral greater trochanter (n=23), and femoral neck (n=27). Specimens were obtained from 61 cadavers (mean+/-SD age=67+/-15 years). Experimental protocols were used that minimized end-artifact errors and controlled for specimen orientation. Tissue moduli were computed for a subset of 18 specimens using high-resolution linear finite element analyses and also using two previously developed theoretical relationships (Bone 25 (1999) 481; J. Elasticity 53 (1999) 125). Resultant power law regressions between modulus and density did depend on anatomic site, as determined via an analysis of covariance. The inter-site differences were among the leading coefficients (p<0.02), but not the exponents (p>0.08), which ranged 1.49-2.18. At a given density, specimens from the tibia had higher moduli than those from the vertebra (p=0.01) and femoral neck (p=0.002); those from the trochanter had higher moduli than the vertebra (p=0.02). These differences could be as large as almost 50%, and errors in predicted values of modulus increased by up to 65% when site-dependence was ignored. These results indicate that there is no universal modulus-density relationship for on-axis loading. Tissue moduli computed using methods that account for inter-site architectural variations did not differ across site (p>0.15), suggesting that the site-specificity in apparent modulus-density relationships may be attributed to differences in architecture.  相似文献   

13.
In 1961, Evans and King documented the mechanical properties of trabecular bone from multiple locations in the proximal human femur. Since this time, many investigators have cataloged the distribution of trabecular bone material properties from multiple locations within the human skeleton to include femur, tibia, humerus, radius, vertebral bodies, and iliac crest. The results of these studies have revealed tremendous variations in material properties and anisotropy. These variations have been attributed to functional remodeling as dictated by Wolff's Law. Both linear and power functions have been found to explain the relationship between trabecular bone density and material properties. Recent studies have re-emphasized the need to accurately quantify trabecular bone architecture proposing several algorithms capable of determining the anisotropy, connectivity and morphology of the bone. These past studies, as well as continuing work, have significantly increased the accuracy of analytical and experimental models investigating bone, and bone/implant interfaces as well as enhanced our perspective towards understanding the factors which may influence bone formation or resorption.  相似文献   

14.
With the prevalent use of DXA-measured BMD to assess pathologic hip fractures and its recently reported lack of reliability to predict fracture or account for efficacy of anti-resorptive therapy, it is reasonable to assess whether variations in the primary and secondary tensile and compressive trabecular microstructure can account for variations in proximal femur strength in comparison to DXA-measured BMD. To that end, microstructural and densitometric measures of trabecular bone specimens, from discrete sites within the proximal femur, were correlated with their mechanical properties. We hypothesize that accounting for regional variations in trabecular microstructure will improve predictions of proximal femur strength and stiffness compared to bone density measured by DXA. Forty-seven samples (seven donors) from seven distinct sites of human proximal femur underwent DXA and muCT imaging and mechanical testing. The results revealed significant variations in BMC, morphometric indices and mechanical properties within the proximal femur. This work has demonstrated that the mechanical performance of each sub-region is highly dependent on the corresponding trabecular microstructure. BMD measured by DXA at standard regions of interest cannot resolve the variations in trabecular density and microstructure that govern the mechanical behavior of the proximal femur. This work suggests that a quantitative Singh index that uses high resolution QCT to monitor the trabecular microstructure at specific sub-regions of the proximal femur may allow better predictions of hip fracture risk in individual patients and an improved assessment of changing bone structure in response to pharmacological interventions.  相似文献   

15.
The macroscopic and microscopic distribution of intramuscularly injected, essentially monomeric, 239Pu was studied in the skeleton of the adult tree shrew (Tupaia belangeri). Data for the period between 15 and 50 months after injection are presented and compared with the data from earlier time points. Between 83 and 500 days after injection the nuclide content and the wet weight of the skeleton decreased to a constant level at about 55 per cent of the maximum values. The microscopic distribution has been analysed in distal femora, proximal humerus, proximal tibia and lumbar vertebra over the whole observation time; additionally at some selected time points proximal femur, femur shaft, distal humerus and distal tibia were analysed. The initial endosteal surface activity ranged from 3.8 to 5.3 Bq/cm2 and decreased to a minimum at about 1000 days after injection and increased thereafter. A similar behaviour was found for the dose rate near bone surfaces which was initially about 0.075 Gy/day on endosteal surfaces. In the deep bone and the deep marrow the dose rate was negligible, about 0.008 Gy/day and 0.001 Gy/day, respectively. The average cumulative dose 1500 days after injection was about 67 Gy on the endosteum, six times greater than the cumulative dose calculated from the mean concentration of plutonium in the whole skeleton. All values are normalized to an injected activity of 37 kBq/kg body weight. The tupaia data are discussed in relation to the available data from monkeys, dogs and rats.  相似文献   

16.
For a better understanding of traumatic bone fractures, knowledge of the mechanical behavior of bone at high strain rates is required. Importantly, it needs to be clarified how quasistatic mechanical testing experiments relate to real bone fracture. This merits investigating the mechanical behavior of bone with an increase in strain rate. Various studies examined how cortical and trabecular bone behave at varying strain rates, but no one has yet looked at this question using individual trabeculae. In this study, three-point bending tests were carried out on bovine single trabeculae excised from a proximal femur to test the trabecular material's strain rate sensitivity. An experimental setup was designed, capable of measuring local strains at the surface of such small specimens, using digital image correlation. Microdamage was detected using the bone whitening effect. Samples were tested through two orders of magnitude, at strain rates varying between 0.01 and 3.39 s(-1). No linear relationship was observed between the strain rate and the Young's modulus (1.13-16.46 GPa), the amount of microdamage, the maximum tensile strain at failure (14.22-61.65%) and at microdamage initiation (1.95-12.29%). The results obtained in this study conflict with previous studies reporting various trends for macroscopic cortical and trabecular bone samples with an increase in strain rate. This discrepancy might be explained by the bone type, the small sample geometry, the limited range of strain rates tested here, the type of loading and the method of microdamage detection. Based on the results of this study, the strain rate can be ignored when modeling trabecular bone.  相似文献   

17.
The present study compared the bone anabolic effects of graded doses of alfacalcidol in proximal femurs (hematopoietic, red marrow skeletal site) and distal tibiae (fatty, yellow marrow skeletal site). One group of 8.5-month-old female Sprague-Dawley rats were killed at baseline and 4 groups were treated 5 days on/2 days off/week for 12 weeks with 0, 0.025, 0.05 and 0.1 microg alfacalcidol/kg by oral gavage. The proximal femur, bone site with hematopoietic marrow, as well as the distal tibia bone site with fatty marrow, were processed undecalcified for quantitative bone histomorphometry. In the red marrow site of the proximal femoral metaphysis (PFM), 0.1 microg alfacalcidol/kg induced increased cancellous bone mass, improved architecture (decreased trabecular separation, increased connectivity), and stimulated local bone formation of bone 'boutons' (localized bone formation) on trabecular surfaces. There was an imbalance in bone resorption and formation, in which the magnitude of depressed bone resorption greater than depressed bone formation resulted in a positive bone balance. In addition, bone 'bouton' formation contributed to an increase in bone mass. In contrast, the yellow marrow site of the distal tibial metaphysis (DTM), the 0.1 microg alfacalcidol/kg dose induced a non-significant increased cancellous bone mass. The treatment decreased bone resorption equal to the magnitude of decreased bone formation. No bone 'bouton' formation was observed. These findings indicate that the highest dose of 0.1 microg alfacalcidol/kg for 12 weeks increased bone mass (anabolic effect) at the skeletal site with hematopoietic marrow of the proximal femoral metaphysis, but the increased bone mass was greatly attenuated at the fatty marrow site of the distal tibial metaphysis. In addition, the magnitude of the bone gain induced by alfacalcidol treatment in red marrow cancellous bone sites of the proximal femoral metaphysis was half that of the lumbar vertebral body. The latter data were from a previous report from the same animal and protocol. These findings indicated that alfacalcidol as an osteoporosis therapy is less efficacious as a positive bone balance agent that increased trabecular bone mass in a non-vertebral skeletal site where bone marrow is less hematopoietic.  相似文献   

18.
This study aims to evaluate the effects of Remifemin (isopropanolic extract of Cimicifuga Racemosa) on postmenopausal osteoporosis. 120 female Sprague-Dawley rats were randomly assigned to four groups: sham surgery with vehicle, ovariectomy with vehicle, ovariectomy with estradiol valerate, or ovariectomy with Remifemin. Daily oral administrations of the vehicle, estradiol valerate, or Remifemin began 2 weeks after surgery and lasted to 4, 8, or 12 weeks. Ten rats in each group were sacrificed at each timestep with assessment of bone mineral density, trabecular bone structure, and biomechanical parameters of the femur and lumbar vertebra. Bone turnover markers were evaluated 12 weeks after surgery. Both drugs prevented bone density loss in the distal end of the femur and preserved the trabecular bone structure in both the lumbar vertebra and distal end of the femur following ovariectomy. Both drugs protected bone stiffness at the tested regions and reduced bone reabsorption in ovariectomized rats. The preventive effects of Remifemin against bone-loss can rival those of estradiol valerate if treatment duration is adequately extended. In conclusion, Remifemin may demonstrate equivalent effects to estradiol valerate in terms of preventing postmenopausal osteoporosis.  相似文献   

19.
The safety and effectiveness were examined of the spirulina alga on bone metabolism in ovariectomized estrogen-deficient rats and hindlimb-unloaded mice. The dosage range was from an amount equal to that recommended in so-called health foods for humans (0.08 g/kg BW/day) to a 100-fold higher dose. The bone mineral density (BMD) of the whole femur and tibia of ovariectomized rats in the any spirulina-treated groups was not significantly different from that of the ovariectomized group, although BMD of the distal femur and proximal tibia was significantly lower in the spirulina-treated groups than in the ovariectomized group after a 6 week-experimental period. BMD of the femur and tibia was not affected by treatment with any dose of spirulina in hindlimb-unloaded mice. These results suggest that the intake of spirulina decreased BMD in the trabecular bone of rodents under estrogen-deficient conditions.  相似文献   

20.
Two-dimensional, finite element studies were conducted of the proximal tibia before and after joint arthroplasty. Equivalent-thickness models projected onto the mid-frontal plane were created for the natural, proximal tibia and for the proximal tibia with four different types of tibial plateau components. All components simulated bony ingrowth fixation, i.e. no cement layer existed between component and bone. In addition, the interface between component and bone was assumed to be intimately connected, representing complete bony ingrowth and a rigid state of fixation. Loads consisted of bi-condylar and uni-condylar forces. Results indicated that conventional plateau designs with central posts or multiple pegs led to higher stress magnitudes in the trabecular bone near the distal ends of the post/pegs and stress shielding at more proximal locations. A design without posts or pegs whose interface geometry mimics the epiphyseal plate minimizes bone stress shielding. An implant consisting of separate components covering each condyle was found effective in limiting component tilting and the consequent tensile stresses caused by non-symmetrical, uni-condylar loading.  相似文献   

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

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