期刊文献+
共找到158篇文章
< 1 2 8 >
每页显示 20 50 100
Mg-9Al合金铸造凝固模型 被引量:13
1
作者 王业双 张咏波 +3 位作者 王渠东 马春江 丁文江 朱燕萍 《金属学报》 SCIE EI CAS CSCD 北大核心 2002年第5期539-543,共5页
在热力学基础上建立了Mg-9Al合金的铸造凝固模型,考察了Mg-9Al合金在铸造凝固条件下的固相分数、体积收缩和微观偏析行为,模拟结果表明.在本凝固条件下.合金实际凝固结束温度在共晶温度附近而不是固相线温度.凝固过程中的体积变化也... 在热力学基础上建立了Mg-9Al合金的铸造凝固模型,考察了Mg-9Al合金在铸造凝固条件下的固相分数、体积收缩和微观偏析行为,模拟结果表明.在本凝固条件下.合金实际凝固结束温度在共晶温度附近而不是固相线温度.凝固过程中的体积变化也是非线性的,在凝固初始阶段收缩较大:当剩余液相到达共晶成分时收缩较小,在凝固后期由于发生共晶转变.产生剧烈的凝固收缩,模拟结果与实验测试相符合. 展开更多
关键词 Mg-9Al合金 铸造 凝固模型
下载PDF
机械振动对纯Al晶粒细化及凝固收缩的影响 被引量:20
2
作者 王红霞 张国平 +1 位作者 许春香 梁伟 《铸造设备研究》 2007年第1期28-31,共4页
研究了振动场对纯Al晶粒细化和凝固收缩的影响,并着重进行理论分析。实验结果表明:振动能明显细化晶粒,改善收缩。随着振动的加强,纯Al中心的等轴晶区逐渐增大,边缘部位的柱状晶区逐渐减小,当振击力达到81.87N时整个试样端面全部变为等... 研究了振动场对纯Al晶粒细化和凝固收缩的影响,并着重进行理论分析。实验结果表明:振动能明显细化晶粒,改善收缩。随着振动的加强,纯Al中心的等轴晶区逐渐增大,边缘部位的柱状晶区逐渐减小,当振击力达到81.87N时整个试样端面全部变为等轴晶。理论分析认为:振动引起的强烈冲击和搅拌作用使金属液产生晶粒游离、增殖、快速冷却是晶粒细化和收缩改善的根本原因。 展开更多
关键词 机械振动 纯Al 晶粒细化 凝固收缩
下载PDF
水泥固化高含水率淤泥的收缩性质 被引量:17
3
作者 张春雷 朱伟 范公俊 《河海大学学报(自然科学版)》 CAS CSCD 北大核心 2010年第3期295-299,共5页
针对高含水率疏浚淤泥水泥固化后的失水收缩性质,通过收缩试验对水泥添加量、养护龄期和淤泥初始含水率对水泥固化后的收缩规律进行了研究.结果表明:对于不同水泥添加量、养护龄期和初始含水率的固化淤泥,收缩曲线的形式与天然黏土相似... 针对高含水率疏浚淤泥水泥固化后的失水收缩性质,通过收缩试验对水泥添加量、养护龄期和淤泥初始含水率对水泥固化后的收缩规律进行了研究.结果表明:对于不同水泥添加量、养护龄期和初始含水率的固化淤泥,收缩曲线的形式与天然黏土相似;固化淤泥的收缩主要是由于水分蒸发引起,水泥添加量越高、养护龄期越长、淤泥的初始含水率越低对固化淤泥的收缩性控制越有利. 展开更多
关键词 淤泥 水泥 固化 高含水率 收缩
下载PDF
大型铸钢件凝固过程数值模拟 被引量:17
4
作者 刘晶峰 李洪友 江开勇 《铸造技术》 CAS 北大核心 2011年第4期443-445,共3页
铸钢件凝固收缩大,铸造工艺复杂,进行凝固过程计算机模拟十分必要。本文基于有限差分法建立铸钢件凝固过程传热计算的数学模型,开发铸钢件凝固过程三维温度场数值模拟分析软件。利用该软件计算了大齿轮毛坯铸钢件凝固过程的温度场,结果... 铸钢件凝固收缩大,铸造工艺复杂,进行凝固过程计算机模拟十分必要。本文基于有限差分法建立铸钢件凝固过程传热计算的数学模型,开发铸钢件凝固过程三维温度场数值模拟分析软件。利用该软件计算了大齿轮毛坯铸钢件凝固过程的温度场,结果与实际情况吻合较好。 展开更多
关键词 凝固收缩 温度场 有限差分法 数值模拟
下载PDF
铸钢件凝固过程计算机模拟技术应用研究 被引量:4
5
作者 范柳光 《大型铸锻件》 2002年第2期1-7,共7页
介绍清华大学开发的计算机三维数值凝固模拟软件FT—Star的功能及其在二重集团公司的工艺性试验和应用研究。解决了软件实际应用于铸钢件凝固过程中出现的一些问题,摸索出一些计算经验和规律。对不同铸钢件如何恰当地剖分网格、确定缩... 介绍清华大学开发的计算机三维数值凝固模拟软件FT—Star的功能及其在二重集团公司的工艺性试验和应用研究。解决了软件实际应用于铸钢件凝固过程中出现的一些问题,摸索出一些计算经验和规律。对不同铸钢件如何恰当地剖分网格、确定缩孔缩松判据参数、优化补贴工艺设计等都进行了深入研究,保证了计算机数值模拟结果的精确、可靠。实现了铸造缩孔缩松缺陷的预测,避免过去靠经验来进行铸造工艺设计,使工艺设计更具有科学性,大大降低了铸造生产成本。 展开更多
关键词 铸钢件 凝固过程 计算机模拟 铸造工艺 缩孔 缩松
下载PDF
铸钢件缩松判据的改进 被引量:6
6
作者 田学雷 李光友 +1 位作者 李成栋 王桂进 《铸造》 CAS CSCD 北大核心 2001年第6期346-348,共3页
对铸钢件凝固过程数值模拟中通常用的缩松判据G/R进行了改进 ,将金属液的静压力P0 引入到缩松的判据中 ,新的缩松判据为P0 G/R。通过试验证明 ,当P0 G/R≤ 3 0× 10 4Pa (s℃ ) 1/2 /cm时 ,铸钢件中将产生缩松。同时说明以P0 G/R作... 对铸钢件凝固过程数值模拟中通常用的缩松判据G/R进行了改进 ,将金属液的静压力P0 引入到缩松的判据中 ,新的缩松判据为P0 G/R。通过试验证明 ,当P0 G/R≤ 3 0× 10 4Pa (s℃ ) 1/2 /cm时 ,铸钢件中将产生缩松。同时说明以P0 G/R作为铸钢件凝固过程数值模拟中的缩松判据比G/R更加准确。 展开更多
关键词 铸钢件 凝固模拟 缩松 预测 判据
下载PDF
土壤固化剂稳定粉质黏土性能的试验研究 被引量:16
7
作者 徐平 乐金朝 +1 位作者 马清文 闫东明 《建筑材料学报》 EI CAS CSCD 北大核心 2010年第1期57-61,115,共6页
以河南省境内某段高速公路附近的粉质黏土为原材料,通过室内试验,首先对比了掺或不掺固化剂的情况下,水泥质量分数分别为5%,10%,15%和20%时集合料的7 d无侧限抗压强度,结果表明掺固化剂能明显提高集合料的无侧限抗压强度;其次制备了不... 以河南省境内某段高速公路附近的粉质黏土为原材料,通过室内试验,首先对比了掺或不掺固化剂的情况下,水泥质量分数分别为5%,10%,15%和20%时集合料的7 d无侧限抗压强度,结果表明掺固化剂能明显提高集合料的无侧限抗压强度;其次制备了不同配合比的集合料试件,分析了水泥和固化剂含量对集合料最大干密度、最佳含水量、7 d和28 d无侧抗压强度2、8 d间接抗拉强度的影响,综合试验结果与成本等因素,选定了最优配合比为m(水泥)∶m(固化剂)∶m(粉质黏土)=4∶8∶88;最后,对该最优配合比进行了延迟和干缩等路用性能试验及SEM微观分析,并在某高速公路选取100 m的底基层进行了现场试验和追踪调查,结果表明上述最优配合比集合料完整密实且强度较高,未出现病害,路用性能较好. 展开更多
关键词 粉质黏土 固化剂 无侧限抗压强度 间接抗拉强度 延迟 干缩
下载PDF
铸造工艺的数值模拟优化 被引量:15
8
作者 胡红军 杨明波 +1 位作者 龚喜兵 李国瑞 《兵器材料科学与工程》 CAS CSCD 北大核心 2006年第6期51-53,共3页
为了研究和预测铸造工艺对铸件质量的影响,设置合理的军用汽车转向臂的铸造浇冒口系统和工艺参数。应用铸造模拟软件对转向臂的三种不同工艺方案进行凝固模拟,根据凝固模拟结果显示的缺陷及内部缩松情况,提出改进工艺方案并对其进行凝... 为了研究和预测铸造工艺对铸件质量的影响,设置合理的军用汽车转向臂的铸造浇冒口系统和工艺参数。应用铸造模拟软件对转向臂的三种不同工艺方案进行凝固模拟,根据凝固模拟结果显示的缺陷及内部缩松情况,提出改进工艺方案并对其进行凝固模拟,选择最佳方案应用于生产。研究表明,3#是最合理的浇冒口布置方式,最优的浇注温度825℃,浇注时间15 s,采用水平分型。应用表明,铸造模拟软件能够准确地预测充型凝固过程中可能产生的缺陷,从而辅助工艺人员进行工艺优化。 展开更多
关键词 凝固模拟 军用汽车转向臂 铸造工艺优化 浇冒口系统 缩孔 铸造模拟软件
下载PDF
真空自耗电弧熔炼TC4铸锭的凝固组织和缩松缩孔的模拟 被引量:14
9
作者 张颖娟 寇宏超 +4 位作者 李鹏飞 钟宏 胡锐 李金山 周廉 《特种铸造及有色合金》 CAS CSCD 北大核心 2012年第5期418-421,共4页
建立三维轴对称有限元模型,采用移动边界法模拟了真空自耗电弧熔炼(VAR)TC4合金铸锭凝固过程的温度场分布,并基于元胞自动机和有限元耦合法(CAFE法)计算了VAR过程中TC4铸锭的凝固组织、缩松和缩孔的形成。结果表明,模拟结果与试验观察... 建立三维轴对称有限元模型,采用移动边界法模拟了真空自耗电弧熔炼(VAR)TC4合金铸锭凝固过程的温度场分布,并基于元胞自动机和有限元耦合法(CAFE法)计算了VAR过程中TC4铸锭的凝固组织、缩松和缩孔的形成。结果表明,模拟结果与试验观察在晶粒结构、柱状晶生长方向,柱状晶-等轴晶转变及缩松和缩孔等方面相符合。 展开更多
关键词 VAR CAFE模型 凝固组织 缩松缩孔 钛合金
原文传递
无凝固收缩铝硅合金的研制 被引量:9
10
作者 常芳娥 坚增运 《铸造技术》 EI CAS 北大核心 2005年第4期300-302,共3页
研究了铝硅合金凝固收缩率与硅含量的关系和无凝固收缩铝硅合金的力学性能。结果发现:铝硅合金的凝固收缩率随硅含量增加而减小,在硅含量为25%时,其凝固收缩率为零。熔体温度处理可大大提高无凝固收缩铝硅合金的力学性能,熔体温度处理... 研究了铝硅合金凝固收缩率与硅含量的关系和无凝固收缩铝硅合金的力学性能。结果发现:铝硅合金的凝固收缩率随硅含量增加而减小,在硅含量为25%时,其凝固收缩率为零。熔体温度处理可大大提高无凝固收缩铝硅合金的力学性能,熔体温度处理后无凝固收缩铝硅合金在铸态和T6处理下的硬度分别达到了119HB和161HB。熔体温度处理后无凝固收缩铝硅合金的耐磨性远高于45钢,其耐磨性在铸态和T6处理下分别为45钢的1.79倍和4.25倍。 展开更多
关键词 铝硅合金 凝固收缩率 耐磨性 熔体处理
下载PDF
大扁锭凝固过程模拟及缩孔优化 被引量:11
11
作者 谢朝阳 张炯明 +4 位作者 胡波 吴苏州 罗衍昭 肖超平 杨锦 《北京科技大学学报》 EI CAS CSCD 北大核心 2011年第1期11-16,共6页
对某钢厂28.7t钢锭凝固过程进行测温,并用有限元方法模拟该钢锭凝固过程温度场和凝固场分布.结果表明:温度模拟值与现场测量值吻合很好,证明模拟具有较高的准确性和可靠性;凝固初期,钢锭底部和保温冒与钢锭模连接处凝固较快;52min时,绝... 对某钢厂28.7t钢锭凝固过程进行测温,并用有限元方法模拟该钢锭凝固过程温度场和凝固场分布.结果表明:温度模拟值与现场测量值吻合很好,证明模拟具有较高的准确性和可靠性;凝固初期,钢锭底部和保温冒与钢锭模连接处凝固较快;52min时,绝热板与钢锭间已形成一定气隙;前3h,钢锭侧面凝固顺序由模壁向钢锭中心平行推进;凝固后期较凝固前期凝固速度快;热电偶测得,保温冒中心凝固时间为428min,钢锭本体中心顶部凝固时间为365min,冒部全凝时间大于本体全凝时间的15%,有利于控制一次缩孔只存在于冒部.通过模拟将浇注温度由1543℃降低到1533℃,不但不影响保温帽钢液对本体的补缩作用,还可以使缩孔减小6mm,有利于提高钢锭质量. 展开更多
关键词 钢锭 凝固过程 温度分布 缩孔 模拟 有限元法
原文传递
Numerical simulation of macrosegregation in steel ingots using a two-phase model 被引量:9
12
作者 Wen-sheng Li Hou-fa Shen Bai-cheng Liu 《International Journal of Minerals,Metallurgy and Materials》 SCIE EI CAS CSCD 2012年第9期787-794,共8页
A two-phase model for the prediction of macrosegregation formed during solidification is presented. This model incorporates the descriptions of heat transfer, melt convection, solute transport, and solid movement on t... A two-phase model for the prediction of macrosegregation formed during solidification is presented. This model incorporates the descriptions of heat transfer, melt convection, solute transport, and solid movement on the system scale with microscopic relations for grain nucleation and growth. Then the model is used to simulate the solidification of a benchmark industrial 3.3-t steel ingot. Simulations are per- formed to investigate the effects of grain motion and pipe shrinkage formation on the final macrosegregation pattern. The model predictions are compared with experimental data and numerical results from literatures. It is demonstrated that the model is able to express the overall macrosegregation patterns in the ingot. Furthermore, the results show that it is essential to consider the motion of equiaxed grains and the formation of pipe shrinkage in modelling. Several issues for future model improvements are identified. 展开更多
关键词 steel ingots solidification SEGREGATION shrinkage modelling
下载PDF
气缸盖铸件凝固过程数值模拟及缩松缩孔预测 被引量:8
13
作者 刘金祥 左正兴 廖日东 《铸造技术》 CAS 北大核心 2007年第2期256-258,共3页
分析了铸件凝固过程结晶潜热释放和缩孔缩松形成条件,应用有限元法模拟了低压铸造铝合金气缸盖铸件的凝固过程,得到了在不同凝固阶段铸件的温度场,预测了铸件的缩松缩孔。在此基础上,对气缸盖铸件凝固过程进行了评价,数值模拟结果可为... 分析了铸件凝固过程结晶潜热释放和缩孔缩松形成条件,应用有限元法模拟了低压铸造铝合金气缸盖铸件的凝固过程,得到了在不同凝固阶段铸件的温度场,预测了铸件的缩松缩孔。在此基础上,对气缸盖铸件凝固过程进行了评价,数值模拟结果可为气缸盖铸造工艺参数设计提供参考。 展开更多
关键词 气缸盖 凝固过程 缩松缩孔 有限元法
下载PDF
大型齿轮圈铸钢件成型过程工艺优化 被引量:9
14
作者 庞富伟 王爱琴 +1 位作者 谢敬佩 王文焱 《热加工工艺》 CSCD 北大核心 2008年第19期55-57,60,共4页
利用华铸CAE软件对大型齿轮圈铸钢件进行了凝固过程的数值模拟,发现了缩孔、缩松缺陷,确定了其产生的位置,分析了其产生原因及机理。在此基础上对其铸造工艺进行了修改,工艺优化,又一次进行凝固模拟,对前后模拟结果进行了分析与比较。... 利用华铸CAE软件对大型齿轮圈铸钢件进行了凝固过程的数值模拟,发现了缩孔、缩松缺陷,确定了其产生的位置,分析了其产生原因及机理。在此基础上对其铸造工艺进行了修改,工艺优化,又一次进行凝固模拟,对前后模拟结果进行了分析与比较。结果表明,合理设置保温冒口、冷铁的大小和位置,加上合适的冒口套,可以保证铸件顺序凝固和金属液的有效补缩,获得结构完整、没有缺陷的铸钢件。 展开更多
关键词 齿轮罔 凝同模拟 缩孔 缩松 工艺优化
下载PDF
Effects of Y on hot tearing susceptibility of Mg-Zn-Y-Zr alloys 被引量:9
15
作者 刘正 张斯博 +1 位作者 毛萍莉 王峰 《Transactions of Nonferrous Metals Society of China》 SCIE EI CAS CSCD 2014年第4期907-914,共8页
The hot tearing susceptibility of MgZn2.5YxZr0.5 (x=0.5, 1, 2, 4, 6) alloys was evaluated by thermodynamic calculations based on Clyne-Davies model. The microstructure and morphology of hot tearing regions of the al... The hot tearing susceptibility of MgZn2.5YxZr0.5 (x=0.5, 1, 2, 4, 6) alloys was evaluated by thermodynamic calculations based on Clyne-Davies model. The microstructure and morphology of hot tearing regions of the alloys were observed by X-ray diffraction and scanning electron microscopy. The solidification temperature and shrinkage stress during the solidification of MgZn2.5YxZr0.5 alloys in the“T”type hot tearing permanent-mold were acquired with the attached computer. The effect factors of hot tearing susceptibility of MgZn2.5YxZr0.5 alloys, such as the solidification temperature interval, the variation of solid fraction in vulnerable region, the residual liquid fraction in the final stage, the type of the second phase of the alloys were discussed based on the above calculation and observation. The results demonstrated that the hot tearing susceptibility in the investigated alloys was found as follows:MgZn2.5Y2Zr0.5&gt;MgZn2.5Y0.5Zr0.5&gt;MgZn2.5Y4Zr0.5&gt;MgZn2.5Y6Zr0.5&gt;MgZn2.5Y1Zr0.5. The highest hot tearing susceptibility of MgZn2.5Y2Zr0.5 alloy related to the following reasons: the largest freezing range, the biggest changing of the variation of solid fraction in vulnerable region, the least liquid film in the final stage of solidification, the formation of the second phase which worsens the liquid flow and interdendritic feeding after dendrite coherency. 展开更多
关键词 Mg-Zn-Y-Zr alloys solidification curve shrinkage stress hot tearing susceptibility
下载PDF
Numerical Simulation of Shrinkage Cavity Formation in Spheroidal Graphite Iron Castings Based on Micro Modeling 被引量:5
16
作者 Wenzhen Li Baicheng Liu Shantong Jin (Metallurgy Engineering School, University of Science and Technology Beijing, Beijing 100083, China)(Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China) 《International Journal of Minerals,Metallurgy and Materials》 SCIE EI CAS CSCD 1998年第1期9-12,共4页
A micro-modeling method (MM) for the quantitative prediction of the shrinkage cavity formation in SGiron castings is proposed. The mathematical models describing the volume changes during the solidification ofspheroid... A micro-modeling method (MM) for the quantitative prediction of the shrinkage cavity formation in SGiron castings is proposed. The mathematical models describing the volume changes during the solidification ofspheroidal graphite cast iron are established based on the models of solidification kinetics. The shrinkage cavityformation of T-shaped SG iron castings is calculated with MM method. The calculated results are compared with theexperimental results. It is shown that the predicted size, shape and distribution of shrinkage cavity by MM methodare in good agreement with the measured results. 展开更多
关键词 solidification shrinkage numerical simulation micro modeling spheroidal graphite cast iron
下载PDF
铁型覆砂球墨铸铁件凝固模拟及收缩缺陷预测 被引量:7
17
作者 梁作俭 赵海东 +2 位作者 柳百成 李贵来 张志坚 《特种铸造及有色合金》 CAS CSCD 北大核心 2001年第5期25-27,共3页
建立了铁型覆砂铸造凝固过程模拟的数学模型。采用动态膨胀收缩法 (DECAM)和K·G/R判据预测球墨铸铁件的收缩缺陷 ,对球墨铸铁四缸曲轴新产品进行了模拟计算 ,根据模拟结果对原有工艺方案进行了改进 ,并在实际生产中得到了合格的铸... 建立了铁型覆砂铸造凝固过程模拟的数学模型。采用动态膨胀收缩法 (DECAM)和K·G/R判据预测球墨铸铁件的收缩缺陷 ,对球墨铸铁四缸曲轴新产品进行了模拟计算 ,根据模拟结果对原有工艺方案进行了改进 ,并在实际生产中得到了合格的铸件。结果表明 ,对铁型覆砂球墨铸铁件进行数值模拟和工艺优化 ,有助于提高铸件质量 ,缩短试制周期 。 展开更多
关键词 凝固模拟 收缩缺陷 铁型覆砂 球墨铸铁 数值模拟
下载PDF
大型挂舵臂铸钢件铸造工艺优化 被引量:8
18
作者 张可锋 杨涤心 +3 位作者 谢敬佩 王爱琴 王文焱 李洛利 《热加工工艺》 CSCD 北大核心 2007年第21期52-55,共4页
利用华铸CAE/Inte-CAST8.0软件对大型挂舵臂铸钢件进行了温度场的凝固模拟,预测了缩孔、缩松缺陷产生的部位,并且解释了形成机理,在此基础上通过对模拟结果进行分析,对其铸造工艺进行了优化。结果表明:合理设置保温冒口、冷铁的大小和位... 利用华铸CAE/Inte-CAST8.0软件对大型挂舵臂铸钢件进行了温度场的凝固模拟,预测了缩孔、缩松缺陷产生的部位,并且解释了形成机理,在此基础上通过对模拟结果进行分析,对其铸造工艺进行了优化。结果表明:合理设置保温冒口、冷铁的大小和位置,可以保证铸件顺序凝固,并且延长金属液的有效补缩距离,获得符合国际标准的挂舵臂铸钢件。 展开更多
关键词 挂舵臂 凝固模拟 缩孔缩松 工艺优化
下载PDF
溶质偏析对铸造铝合金疏松缩孔缺陷的影响 被引量:8
19
作者 左强 洪润洲 +1 位作者 周永江 杨守杰 《铸造》 CAS CSCD 北大核心 2017年第4期381-385,共5页
研究了溶质偏析行为对铝合金凝固过程中缩孔疏松缺陷的影响。分析不同偏析模型(杠杆定律、Scheil模型、逆扩散模型)下合金的溶质再分配与热物性参数;采用Pro CAST有限元模拟软件包的Niyama判据预测不同偏析模型下凝固缩孔疏松缺陷的形... 研究了溶质偏析行为对铝合金凝固过程中缩孔疏松缺陷的影响。分析不同偏析模型(杠杆定律、Scheil模型、逆扩散模型)下合金的溶质再分配与热物性参数;采用Pro CAST有限元模拟软件包的Niyama判据预测不同偏析模型下凝固缩孔疏松缺陷的形成趋势,并结合试验结果给出溶质偏析对凝固缺陷的影响。 展开更多
关键词 偏析 凝固 铝合金 缩孔疏松
下载PDF
基于凝固收缩的板坯铸机辊缝研究 被引量:8
20
作者 马长文 陈松林 +1 位作者 郑天然 杨晓山 《钢铁》 CAS CSCD 北大核心 2008年第3期44-48,共5页
为了研究辊缝收缩对减轻板坯中心偏析的作用,采用数值模拟方法计算了板坯的温度场和凝固收缩,分析了坯壳厚度和凝固收缩量在冶金长度内的分布,建立了3种辊缝收缩实验方案并比较了其效果。结果表明,在凝固末端以前2~3个扇形段进行... 为了研究辊缝收缩对减轻板坯中心偏析的作用,采用数值模拟方法计算了板坯的温度场和凝固收缩,分析了坯壳厚度和凝固收缩量在冶金长度内的分布,建立了3种辊缝收缩实验方案并比较了其效果。结果表明,在凝固末端以前2~3个扇形段进行辊缝收缩有效减轻了板坯中心偏析,辊缝收缩量应略大于铸坯自由凝固收缩量。 展开更多
关键词 连铸板坯 凝固收缩 辊缝 中心偏析
原文传递
上一页 1 2 8 下一页 到第
使用帮助 返回顶部