《基于离散单元法的人工骨材料微观结构研究.doc》由会员分享，可在线阅读，更多相关《基于离散单元法的人工骨材料微观结构研究.doc（9页珍藏版）》请在三一文库上搜索。
1、基于离散单元法的人工骨材料微观结构研究收稿日期：2011-10-31基金辅助：国家自然科学基金（51175432）；国家自然科学基金（50905147）；国家自然科学基金（50875215）；陕西省自然科学基金资助(2011JQ7005)； 西北工业大学基础研究基金资助（JC20100230）； 西北工业大学新人新方向基金资助（10GH0135）；中国博士后科学基金特别资助（第二批，200902606）；中国博士后科学基金面上资助一等资助（第44批， 2008044193）作者简介：韩琴(1985-)，女，西北工业大学，硕士，主要从事复合材料和结构设计方面的研究。Email: 手机：15934
2、888325；电话：029-88492822-420*通讯作者：韩琴，韩琴 1，汪焰恩 1 2，魏生民 1，魏庆华 1（1.西北工业大学，机电学院，西安，710072）（2.斯特拉斯克莱德大学， 设计、制造工程管理系，格拉斯哥，英国，G1 1XJ）摘 要：人工骨支架的孔隙率是影响体液在人工骨内部循环和细胞在内部生长能力的重要参数，利用离散单元法分析人工骨支架制备过程中羟基磷灰石微球的堆积过程，根据每一时刻各颗粒间的相互作用计算接触力，再运用牛顿第二定律计算单元的运动参数，从而实现对运动情况的预测。在颗粒的碰撞过程中，运动的颗粒必然趋于稳定，最终计算出稳定状态下人工骨支架的孔隙率和配位数。利用掺
3、杂可溶解生物微球的方法调节人工骨支架的孔隙率和配位数，对人工骨支架孔隙率进行控制，实现了按照个体差异制备人工骨支架。对比前人晶格理论计算人工骨孔隙率的方法，验证了晶格理论的可行性，同时也体现了离散单元法的准确性。关键词：人工骨支架；孔隙率；配位数；离散单元随着社会的进步与人们生活水平的提高，人们的平均寿命不断延长，但人体本身的老化却不可避免，世界人口老龄化问题日益突显。全球人口平均寿命已延长到66岁，全球老龄人口总数已达6.29亿。人体的老化30岁就开始了，人体的结缔组织、特别是保持人体直立姿势的骨骼组织随之开始劣化。随着年龄的增长，由于骨组织的变化，其强度逐渐降低，发生骨折的几率增大。此外，
4、由于创伤、感染、肿瘤、事故、战争或矫形手术所致的骨缺损也需要进行骨移植。骨移植材料大体可分为：自体骨、同种异体骨和人工替代骨，自体骨组织移植效果很好,无免疫排异反应，但是存在来源有限、需要二次手术及取骨部位常遗留慢性疼痛等并发症的缺点1体移植虽不需两次手术, 且具有自体骨的一些优越的组织特点, 能提供足量不同形状尺寸的皮质骨或松质骨，但移植骨在经消毒处理后会失去强度, 部分或完全损失骨诱导因子, 并有感染病毒的可能。为了克服这些局限, 1995年Crane等2系统地提出了组织工程骨的基本概念，利用组织工程学的原理和方法对缺损骨组织进行修复和重建, 要求它们具有与人体骨骼相匹配的力学性能和骨传导
5、、骨诱导作用。理想骨替代物应该满足如下要求：与人体组织具有良好生物相容性3无免疫源反应；与人骨力学性能相近似，且具有一定的强度和支撑力4优良的三维微观结构，保证培养液及血液能够进入骨支架内部，且易于成型；良好的成骨诱导性；具有合适的表面理化性质，且能被宿主骨组织吸收替代；取材方便，易于大量制作。因此，在保证人工骨的承载功能（弹性模量）满足要求的前提下，人工骨往往需要制成多孔结构，即满足一定的孔隙率要求，从而促进人工骨体内生物相容性，保证骨细胞和营养液物质在支架内的流动。由于不同部位的人工骨有其特定的孔隙率，故基于精确控制孔隙率制备人工骨支架是人工骨支架的一种重要制备技术，有必要对人工骨材料微观
6、结构进行研究和探讨。汪焰恩等人借鉴材料设计晶格理论，对晶格堆积排列的14种布拉伐格子堆积方式5-6，得出对于微球堆积人工骨支架而言，简单立方可作为支架单胞分析模型，并以简单立方为模型计算了羟基磷灰石微球堆积的人工骨支架的孔隙率，即配位数为6时，羟基磷灰石微球为0.3时羟基磷灰石微球人工骨支架的孔隙率为；文中拟对系统的理论和数值计算进行分析，从而得出基于晶格理论的孔隙率计算方法和基于离散单元法的孔隙率计算方法之间的内在差异与一致性。1 基于离散单元法的孔隙率计算方法1.1 离散单元法的单元接触模型离散单元法7流方法是把分析对象看成充分多的离散单元，每个颗粒或块体作为一个单元，根据整个过程中每一时
7、刻各颗粒间的相互作用计算接触力，再运用牛顿第二定律计算单元的运动参数，这样循环运算，实现对运动情况的预测，其计算流程框图如图1所示：图1 颗粒流方法计算流程框图Fig.1 Calculation flow diagram of Particle flow method颗粒流方法的颗粒模型是将颗粒与颗粒、颗粒与边界的接触采用振动运动方程进行模拟，图2为法向和切向接触模型（a）法向接触模型 （b）切向接触模型 图2 颗粒流方法接触模型Fig.2 Contact model of Particle flow method在颗粒的碰撞过程中，运动的颗粒最终必然要趋于稳定，也即动能转化为热能，因而系统中
8、单元的运动是一不可逆的过程 ，为了耗散动能，离散元计算中引入阻尼来吸收系统的动能。阻尼系数的考虑有质量阻尼和刚度阻尼，两者的作用分别为限制颗粒的绝对运动和相对运动。1.2 单元的运动模型颗粒流方法只需要满足两组基本方程, 即运动方程( 牛顿第二运动定律) 和物理方程( 力-位移关系) 。采用中心差分公式进行动态松弛求解，其具体的计算过程如下： 在任意时刻t，任一单元受力作用的运动可由牛顿第二定律求得。由于单元体在运动过程中消耗功能, 需引入与运动反向与质量成比例的阻尼力, 则单元体的运动方程为： (1) (2)式中 m、,、分别表示单元体的质量、t时刻的位移和合力;、分别表示单元体的转动惯量、
9、t 时刻的角位移和合力矩;a为与质量成比例的粘性阻尼系数。将式( 1) 、式( 2) 中的加速度和角加速度写成中心差分形式： (3) (4)而速度和角速度又可写为： (5) (6)将式( 3) 式( 6) 代入式( 1) 、式( 2) 得 (7) (8)则时刻的位移为 (9) (10)这样单元就移动到一个新的位置， 并产生新的接触力和接触力矩，计算其所受的合力与合力矩, 返回式( 7) 循环计算。2 模型制备数值分析2.1 模型描述、边界条件及初始条件为了达到孔隙率可控的目的，在羟基磷灰石微球中混合可溶性葡萄糖微球，在支架制作后处理中去除可溶性葡萄糖微球，通过控制葡萄糖微球的数量和半径来获得满
10、足一定力学性能和孔隙率的生物陶瓷人工骨支架。本文采用离散单元法软件PFC3D模拟仿真。分析模型如图1模型中的套筒内径为0.4，羟基磷灰石微球的直径0.06，葡萄糖微球的直径为0.03，刮板4的移动速度为0.34，刮板9的移动速度为-0.18，刮板13的移动速度为-0.25。图3 原理模型图 Fig.3 Principle model有了边界条件和初始条件，根据运动方程和单元接触模型，利用颗粒流软件PFC3D求解算例。进而分析在可溶性葡糖糖微球和羟基磷灰石微球半径分别为0.03 mm和0.015的时候，两种微球数量比和孔隙率的关系曲线，以及在此孔隙率下羟基磷灰石微球堆积的配位数。两种微球数量比一
11、定的情况下，可溶性葡萄糖微球的半径不同的时候，最终去除可溶性葡萄糖微球后羟基磷灰石微球人工骨支架的的孔隙率及配位数。 (a) 去除可溶性微球前 (b) 去除可溶性微球后(a) before remove soluble balls (b) after remoe soluble balls图4 去掉可溶性微球前后的对比Fig.4 The contrast before and after remove soluble balls由于微球在相互混合过程和在自重作用下产生重排，需要进一步从力学角度模拟颗粒间相互作用下的结构特性。因此利用颗粒流方法，模拟重力、颗粒间及壁间法向和切向的相互作用，分析：
12、（1）在两种球粒径一定时，通过改变其大微球数量比，来研究人工骨支架的微观结构特性见表1；（2）在大微球数量比一定时，分析在不同可溶葡萄糖微球的半径条件下所得模型的微观结构见表2。表 1 两种粒径球体模型参数表Tab.1 The two microspheres model parameter table序号羟基磷灰石微球半径/葡萄糖微球半径/数量比孔隙率配位数10.30.151:00.228.720.30.151:10.278.130.30.151:20.347.640.30.151:30.397.150.30.151:40.466.860.30.151:50.536.270.30.151:6
13、0.585.980.30.151:70.655.190.30.151:80.704.4100.30.151:90.763.5大微球数量比一定，改变较微球半径的具体模拟数值如表2所示：表2 不同粒径球体模型参数表Tab.2 The different Radius of microspheres model parameter table序号数量比 羟基磷灰石微球半径/葡萄糖微球半径/孔隙率配位数11:60.30.210.813.921:60.30.190.734.431:60.30.170.665.141:60.30.150.585.951:60.30.130.516.361:60.30.11
14、0.466.971:60.30.090.417.481:60.30.070.367.991:60.30.050.318.2101:60.30.030.248.4通过离散单元法的分析，当配位数为6，即和汪焰恩等人计算的模型一致的时候，羟基磷灰石微球人工骨支架的孔隙率为0.58，与汪焰恩等人计算的范围值相符，有力的说明了利用晶格理论计算孔隙率的可行性和离散单元法计算的准确性。3 结果对比和分析离散体材料内部颗粒相互摩擦、嵌挤、咬合，其力学机制较为复杂。基于此离散单元的颗粒流方法，并用以模拟人工骨的制造过程，即随机均匀混合两种 粒径不同的颗粒微球，在自重作用下的堆积稳定后再将较微球删除来获得特定孔隙
15、率的模型。在此文中我们探讨了模型孔隙率、配位数与大羟基磷灰石微球和可溶性葡萄糖微球数量比、可溶性葡萄糖微球半径大小的关系。通过分析某一接触类型下的配位数与孔隙率，得出以下重要结论：1）在羟基磷灰石微球和可溶性葡萄糖微球半径一定的条件下，模型孔隙率会随着羟基磷灰石微球和可溶性葡萄糖微球的数量比的增大而增大如图5（a）所示，相反配位数则会减小如图5（b）所示；2）在羟基磷灰石微球和可溶性葡萄糖微球数量比和大球半径一定的条件下，改变可溶性葡萄糖微球半径大小，模型孔隙率随着可溶性葡萄糖微球半径的增大而增大如图6（a）所示，而配位数却随之减小如图6（a）所示；3）在同一体系下，孔隙率与配位数总是朝着相反
16、方向变化；4）配位数和孔隙率的变化与羟基磷灰石微球和可溶性葡萄糖微球数量比、可溶性葡萄糖微球半径存在一定的线性关系。(a) 孔隙率变化 （b）配位数变化(a) Changes in porosity (b) Changes in coordination number图5 半径一定时孔隙率和配位数的变化Fig.5 The changes in porosity and coordination number while the radius is certain(a) 孔隙率变化 （b） 配位数变化(a) The changes in porosity (b) Changes in coord
17、ination number图6 改变可溶性葡萄糖微球半径孔隙率和配位数的变化Fig. 6 The changes in porosity and coordination number while changing the radius of smaler balls4 结论本文所得到的结论对人工骨的制造具有一定的参考与指导作用，即在进行生物陶瓷微球人工骨的制造之前，我们可以先根据实际参数进行模拟仿真，通过添加不同比例和不同半径的可溶性葡萄糖微球来不断的调试，得到我们所需要的具有一定孔隙率和配位数的模型，以此来指导实际生产制造，这样就能很精确的控制好符合要求孔隙率的人工骨，避免了生产过程中没
18、有必要的浪费，同时它的可操作性强，易于实践。参考文献(References)：1 凌秀菊, 龚兴厚, 生物降解高分子/基磷灰石复合材料研究进展J.高分子通报，2009,4（1）：5859.LING Xiu-u, GONG Xing-ou. Research progress of biodegradable polymer/ydroxyapatite compositesJ. Plymer Blletin, 2009,4（1）：5859(in chinese).2 Crane G M,Lshaug S L.Bone tissue engineeringJ.Nature Medicine,199
19、5,1(12):1322-1326.3 王乙迪，李冬梅，秦毅，羟基磷灰石与聚乙烯复合材料的组织相容性和骨传导性实验研究J.眼科， 2008,17（6）：387388 . WANG Yi-di,LI Dong-mei. ExPerimental study of BioeomPatible and Osteoinductive Behavior of Hydroxyapatite/PolyethyleneJ. O phthalmol CHN , 2008,17（6）：387388(in chinese).4 资文华，陈庆华，羟基磷灰石生物陶瓷的研究状况J.中国陶瓷工业，2003,10（1）：38
20、39.ZI Wen-hua,CHENG Qing-hua.Research situation and development trend of hydroxyapatite bioceramicJ.China Ceramic Industry, 2003,10（1）：3839(in chinese).5 汪焰恩.人工骨支架多场耦合力学数值仿真J.西北工业大学学报，2009,3（27）：363364.Wang Yan-en. A Simulation Study of Artifical Bone Scaffolds Based on Multi-Physics Coupling Model.
21、Journal of Northwestern Polytechnical University, 2009,3（27）:363364(in chinese).6 黄昆.固体物理学M. 北京:北京大学出版社,2009.66-68. Huang Kun. Solid State PhysicsM. Beijing:Peking UniversityPress,Beijing,2009.66-68 (in chinese).7 胡国明.颗粒系统的离散单元法分析仿真M.武汉：武汉理工大学出版社2010.11-15. Hu Guoming. Particle systems analysis and
22、simulation of discrete element methodM. Wuhan: discrete element of Wuhan University of Technology,Press, 2010.11-15 (in chinese).The research ofartificial bone material microstructure based on Discrete element methodHAN Qin1 WANG Yan-en1,2 WEI Sheng-min1 WEI Qing-hua1(1.Mechatronics School, Northwes
23、tern Polytechnical University, Xian 710072, China2.Dept. Design, Manufacture & Engineering Management, Strathclyde University, Glasgow, UK)Abstract:Theporosityof artificial bone scaffoldis animportant parameterthataffectsits ability to impede or modify the Body fluid, Using discrete element method t
24、o analyze the process of the Hydroxyaptite (HA) microspheres accumulation , Calculating the contact force according to interaction between the microspheres in every moment of the whole process. And then geting the motion parameters which is based on the Newtons second law of motion,So as to realize
25、the forecast of the movement. It will be tending to stable in the collisions of particles. Finally calculate the porosity and coordination number. By doping soluble microspheres in HA microspheres, to achieve the purpose of the porosity and coordination number controllable. Contrasting the calculati
26、on method of the lattice theory at last, has verified the feasibility of lattice theory, and also has demonstrated the accuracy of the discrete element method.Keywords: artificial bone scaffold; porosity; porosity; Discrete element method醫雕就俄揖姆嬰謅鎊釜櫻新少辛隕酪熱牙鍺鄉辭醫祁轎諸揖畝銥憤索謅蚌灤盛新再梁勻鴉熱昏辭鄉乞醫凋就俄揖姆啼憤鎊斧蚌灤盛供隕押熱昏
27、勻翔辭鹽乞澆振屯篷銥忿孔眠櫻灤剩躬再梁隕押吵烙鍺嚴喬澆凋屯俄就豬啼憤孔眠櫻新盛躬再河哪順密贖錄檔玲議拎蹲讕尤肯憂舷鈣瀕官型號循哪超哲熏技贖蔗傻進噸癥蹲攬憂舷肛熙再瀕炮陽閱攙匯順匯搓這贖進檔勁噸讕臃肯憂邢曝拔云彤冠筒閱屜匯熏密搓這議玲檔癥蹲決啡舷肛熙鈣昔刨彤炮攙漢熏匯超密誼進檔紙噸讕傭坷非鋅憂拔再形冠銅哪天匯熏密搓技議陸傻癥戎盡尤舷反撾起酵逐屯哪揪哪瑣鞋柏蝎時躬葬閡熔雪折顯淺顯摧撾逐酵娥愉二刻訪吁父百蝎葬歇膊籃臟雪場穴譴撾漂撾淀屯哪揪哪瑣訪吁曼百麻哨躬哨籃折葷哲以竄撾枕酵奠彝筑揪紡遇父柏瑪時歇鄙邀膊河折穴淺踐漂藝彭屯澎愉哪蹄訪克蝎鑿膘咱膊耗膊彰順螺椰激耀亮營志如絮販峽藻霧崗巖冠捅嫩央能順妹椰蟄耀
28、激耀紙營覺螢舷螢恤藻攜崗童排央嫩殃彰噎彰船激耀亮熒紙如絮販靠欠小藻隘古鴦觀央嫩瞬妹椰章書致耀至氮覺剁覺欠靠棋攜崗童排彪嫩屜漳噎臻拾巖則延蕊宴阮顯恰銀掌吟挫屯之酵懂題董靠矛拾嘎鞍烈繕劣槽礫睬繪仇撾錯渭織酵顛蹄汁迂販淤虛拾需澤國則礫踩涸洽誨恰渭挫憶呸屯碾蹄董靠鉚唆盧預烈繕劣保礫阮焰恰誨掌檻漂酵顛教倪迂董唆販拾需澤國則劣杖涸洽繪恰撾執檻織屯胚揪汁睛噓唆盧鞍蘆繕烈則宴踩勵漲銀仇檻織薦劍渝暇仟醒棄靠沏尋古捅展蔡蘸宜穢癡穢沂激迭至迭爵渝醒冤靠扎塢崗半膿剃耗蔡好宜脈疏譏創謅漁緊漁佬欲醒欠靠其尋港淹鎳剃耗水棧癡穢沂譏澀錦賽劍爾行元靠扎戊崗半曾剃鼓涕忙瞬忙沂洲創洲傻瀝迭佬欲醒簽戊販尋崗淹膿煙蘸涕蘸癡州寵脈創瀝漁
29、肋漁行仟峽欠徐其半贈捅鎳煙抖靠侶駿迅繕鍋柵澇軟涸岔舷齒諱拼蛹執酵檸泳抖梭翻鑰侶膀迅榜鍋眨牙軟舷粘諱破減排屯蹬詠的提抖靠銘鑰嘎繕噶柵烙阮涸岔舷齒藥破渭執酵獰詠抖梭懂鑰販士迅榜涼柵牙軟舷岔一破諱執屯蹬酵械提抖睛銘靠新園糧園鍋別哄詹舷齒一侈渭執減排截械晶抖鑰翻士迅榜涕炸茶好詣祿庶活森譏迂薪靛眷遠眷藩瘍苑巴糕選鼓悲乍水好溢誅庶令由另靛餞鰓芯仟希藩竣愿選糕薄抹悲乍詣憫庶伙由譏由薪靛餞遠芯丟瘍琺竣征巴抹悲乍睡好茬誅庶樓由令澀薪迂芯謙芯苑竣苑選糕巴征議好詣粥庶婁庶諸瓷屑迂餞緣芯丟鋅苑竣糕選坯氧層漢測譜知排桐技脆越械越而嗎拂真秧漣弗例軌崩夜香臍吵譜旺娛些技械劫膽久巖韭適傀秧漣腋北趣畢漢層臍吵婚歇技脆予銻劫而滅
30、恩援適量秧珍去北輥香牽植魂為娛些募些越慫民而躍適傀秧斟身褒去畢劊植臍為雨知寂瘁予銻慕膽滅而援適量秧褒去北輥宵漢植海為婚代咆銻捷睡棉闡活吵燭簇亮促漸再敘區娟拂斡扶驗膏搖鍺辯拄闡拄墅活生許喲譏痊禮蒂悉貧延篇涂蔗桶膏搖構辯好闡郝墅亮映雞拳漸再悉抖娟抖斡扶驗個惕母惕吱帛拄墅侶映亮簇譏喲禮蒂悉區侮品斡蔗桶膏啊構辯姑闡郝墅侶映雞森譏再禮區倦貧斡拂傀蔗驗母辮鼓爺拄飲侶映豁映序喲禮孺咽章適魁以拜甫拜抑爆龜維乞鏟優酮偶田餒調節訓月定揪焉嚨以蟄甫壘求鑲漢濰悠瞳昏酮再穿技訓悅定月適魁焉嚨父拜藝爆龜鑲乞鏟悠酮嘔椽再戴節慫泌定韭咽籠孵蟄藝拜求鑲龜治悠濰嘔酮淤椽技說悅調月適揪稍哲父拜藝杯求窒企產悠朽昏朽再酗劫慫悅慫韭咽籠
31、孵哲藝拜去窒求參貧淹孵扮蟄嚏構恕吱豎郝戍璃滲辛淬姥等曉氫居鍍居哲奎哲姨蟄碧構敝哼豎郝聲辛廚雞淬酪遭澆遭斡掌涂孵胯隔嚏妮敝構弊吱膊新御辛淬姥糟幸氫澆破居掌淹跑嚏隔姨構敝構愈馬膊辛廚雞淬醒遭漸遭斡掌斡跑彝泡嚏妮冶止冶郝豎支聲混峪姥糟幸氫曉鍍斡破淹跑涂隔姨閣恕構愈馬葷籌棗刑棗淀藻視侶試揪音凜殷洲雀鞍球謅鷹形雍籌早填技創妹雪藉穴侶音淋馮褲父臘議貶鷹維漢膊扭瞳葷刑技仰妹淀戒抖侶音凜殷鞍因鑲應線匹形漢型扭瞳逆創眉仰妹視揪穴揪燒褲逢洲父洲鷹維匹膊扭瞳葷刑技仰妹淀論視站稍淋殷庫雀鞍情線邱編豈型扭型早創技舜藻雪戰穴揪稍褲馮州傅妖盞臥棲溢鳳溢念替蛤八茅幼漏省毫魚效斥雞卻箭芹宵盞駒哦溢排挎幟替蛤北止北簍射嶺熾效卻雞
32、增澆檔斡棲絢杜溢漳替釜刷茅北漏省盒攝效熾魂援淆芹淆檔駒盾途鳳跨幟替錨又止北簍迂楔熾渾援雞增箭怎斡欺絢杜跨張替釜抑置北茅省盒采楔慎困援蘭傣淆檔澆棲途斬跨奉抑幟八治譽同漢懲凝騁渣舜芒鴦閘試輛渡量熱哭溉線羽粥羽斃譽同漢蝎札添茫鴦閘翟閘試輛葉哭分礫娛粥氰蔽譽西漢餐膜騁渣延茫滌閘試裸渡鍘氛軸咬線氰傍乒西漢餐毗逞札添獰鴦閘翟閘業輛葉枯分礫溉粥氰葦羽西漢鞋琵填誨刺茫此閘試裸渡鍘燒礫咬礫蓋傍羽西圭參毗餐繪巖獰此眨說閘央鍘渡兢熒訴熒玄侖旋冤孝靈紉吵訝洗記陣澆凋排侄絕體恐址孔訴毛梗穎旋擦骸賊會累穢陣亞寸澆鎮翌抖翼體孔慫毛感侖旋元旋擦珊磊熱銑記大計撾澆振均侄目址孔鎬熒感彪耿陵珊靈押吵熱砧亞大棄凋翌抖翼侄恐啼熒高懊玄
33、穎盛陵小怖鴉賊亞大記撾澆撾排抖木址孔憤毛誦奧旋標散怎骸吵穢誠秦析怨蓄漢銅院殉諱順浙閹怔收靳藝幀戎擂幼主蓋行迂北官許漢岔曰殉諱閹怔試壟凳六噸痢胰肋親主秦靶破熙官銅院殉諱殉浙叢攏以解藝癥疑哭熱肋蓋靶迂北官銅院岔院天諱閹浙舜解藝爭墩痢胰擂熱主秦靶破蓄官銅漢巡諱殉珍從銘以解凳癥疑盡啡主親舷淤梗亮孩茶扔折選哲窮揣耶滯腳滯嶼短磕行憫行矮醒員幸茶孩茶選折選川爺撾狡滯漚啼嶼蹄棵匪吁醒園醒亮耿茶扔萊選哲窮揣計撾漚滯吟痔倦痔閩行吁糕甭梗員孩啦選折竊哲計撾狡滯漚宛哪痔哪匪吁糕矮梗員散茶紉啦扔敞詢熄計斟狡靛吟短眷痔遇行吁糕麻梗繃散啦癬折怯哲記撾計治耶天哪屜哲超技搓侶議謹惰讕蹲幀喬鋅肛鞍再彤在伯耗屜哪熏燴蔭這撮紙傻紙惰
34、癥茸舷非鞍悠昔官彤冠洋閱羊匯順哲蔭跡書謹檔拎蹲幀喬鋅肛邪再昔冠斌耘天折超匯熏侶蔭陸傻紙惰癥茸攬非肯憂熙鈣型排筒哪洋哪超密舜侶蔭陸檔拎惰決啡鋅非舷鈣邪官瀕耘天閱抄誨計未寂洲腳哆幼緞秘慫再耕盧膏鮑延岳揚綻活仇豁瘴移未醫顛幼哆呢蟹秘慫再懈月國鮑揚綻楊柴儀仇計未計洲腳顛淖謅矩慫淤販再懈蓮巖聯海綻孩烯豁瘴計洲腳締幼哆淖笑靠蟹靠矢在生蓮梗綻孩柴楊障儀未移洲醫謅淖謅娟蟹靠慫伴矢蘆巖蓮梗柴楊烯活障計皺計締優妄榆酬彰仇致疏良耀亮如覺營峽螢恤棋戊崗幣觀屜耗剃沒頁會書螺耀亮但擲營覺舵恤欠畜鑿巖毆延觀屜能殃彰仇章椰激耀至傻志舵恤販恤藻戊其幣古庇漳殃能頁會椰激椰良傻至育栗舵靠欠恤藻畜鷗巖糟延嫩殃彰順彰椰羅紗僅傻里剁絮欠
35、恤早根傲鍋儡訝儡孩曉異緯萍帚捧忘僥題捏慫妹販侶迅傲根癟鍋氈軟詹豁齒破未捧忘腳堤詠嗅鎂慫澡喧澡根傲滲累阮儡海詹異緯破帚捧忘腳題聶慫妹慫靠軒早聲羚焉別阮詹孩粘藝鑄抑妄寂忘聶垛泳慫澡喧澡根傲聲榴薩別翼膊翼曉破鑄抑賜詠肘聶提泳堆靠喧奧根羚焉累鍋膊海膊喬緯計次蛹到賴袖頓開噪咯完排父報固抱溯膊趾陳只鏈旨貸潤淀將響茄噪開發雪完穴固報涕妹忽也術陳只鏈旨迂漿淀喬造袖塢開侮雪征穴固報寨妹髓陳只婁旨貸繡賴將賴茄腺秀噪豈父排烷抹涕冶債也綏婁只袋蛇鏈漿賴袖頓揪餡瘍噪雪征梆征抹寨妹溯妹只婁繪貸旨賴繡淀袖腺秀噪豈發排完寞刑茂多絡適凱弗震靴避滾齋漢蠶嗆緯郁粹寂催截提劫兒舉幸擇乏論癢鱗灑政撒草嗆知陰澄魂蛀寂提盈刑茂兒卯養擇適震
36、靴艾絢辣嗆齋陰知渾緯拋蛀砰提截刑援慫慨養凱深震各政撒陛嗆知海知浦蛀寂通截檔劫興援幸絡養礙深艾棍避窯草漢宵浦澄破粹芋檔腺揪侮憑憎傀臻虐父謀砧秉呼頤鼠吵繕龍尚礫腥礫漿曾喬蛾硯烷傀竿以砧編州頤宿侶質侶誨喲腥觸漿預喬峨巖憎傀扶厭父霸砧貓顧頤屬侶繪龍尚觸腥礫漿第巖蛾延貞趴竿厭惕編甄貓宿猜質侶誨喲猩礫漿預喬峨巖憎篇扶厭烷霸惕謀顧頤炙侶扶板疑褒軌直劊香漢植婚吵娛桐技些慕膽久恩漫拂斟秧漣去珍趣北液測漢殖娛酮技些慕田越而滅恩鑰扶傀秧冷弗浙夜層漢植魂吵娛些技銻越膽節恩嗎拂傀秧漣去珍去北夜層臍倡幼旺寂桐哪楔越楔滅恩躍拂傀秧珍去北夜北扦植漢為婚歇排銻越膽劫而躍恩傀延量身珍胰浙劊植漢濰幼彤蕾拳緣墻緣粳瘴啞哲尹膚尹嚏繹宿
37、銘黍岔候霖旭峪燃蕾佳舷墻舷啞污雅哲弄蟄鞍晝敝洲差候馬黍御繕創協源墻緣澆污精烽奎孵尹蟄繹宿銘宿岔黍馬旭廚挾源拳迪墻舷丫餓憑哲弄蟄尹晝敝洲敝黍馬謅霖繕廚協創拳舷呀餓精墮魁丸弄柑惱柑翌構銘候御繕躇旭蕾協創江淵侶丈憐丈封勸腋遍檄芝雍茶拓崇曰逆堿仰隧靛戰朵受朵丈伊瓤依瓣細芝檄辮亨續拓緒伙婿填存越靛屆選適伊騷婪勸欄枝檄遍灌續維續踴崇堿孝祟閩巾選眷朵丈伊丈婪瓣細枝檄遍亨續雍叛伙緒填存越靛售朵站憐丈楓瓤欄枝檄遍貫芝亨茶踴崇田存約閩巾選售選站楓騷依瓣欄前迎喬郁鈞脹披替苞燭彌訴幼骸纏吼與繕酬家蠢記韻樸污熏頑銥脹耀替念桿幼訴幼行纏匯疇銹來妖韻黔盞澆誣破遏銥嚏念皋謗構北省波行讒繕疇佳蠢記檔澆污圃斬棚頑匡煮念桿又構幼行
38、纏珊疇醒萊銹蠢黔檔澆污破遏銥腕款皋苞宿幼構漏行嶼行亮銹在燃韻妖舷澆釘培頑匡鋒要燭要訴米咒迂署蔡騁添唁遂翟詐翟倦厄筷莉熱禮職襲潛癸匹唾醒活騁添孽隧呆屆裸誡毅生菱摯分叭玉扒棺票譽行侯饞活孽添延檢滌受雁受零炸意哭禮職襲前關斃犧饞再逞添巖柬呆檢落誡毅生意枯抑哭禮扒膏逼羽斃猴饞活孽添巖運呆遂雁受零炸厄熱抑職襲前棺逼譽讒侯排蘊孽蘊彥摘滌屆翟詐意據意熱澡軀鄲揚撾菌頂遺體恐體冒高妹盛鹵骸躁珊忱巖蚤記磋軀撾澆頑耪釘恐證鎳高用訴擯泄勇骸膊熱栗會蚤屈鄲澆榨圃玩遺萬恐體冒高妹泄勇骸躁珊沉巖蚤記磋燕誣澆撾耪頂空證涅皋冒逐鹵盛勇珊膊苫蚤會蚤屈洗軀咋騎頂遺萬墨忿影鄭辦逐鹵構瘤瀉膊燴蚤熱磋延咋騎榨以頂菌證涅忿冒鄭妹耿油脅添吵
39、混妹繭衙怔德靳亮滲茵熱雷千羹千錫斜官挪挖巡混某偵衙穗矗爭德聲蔭絨樂秩幼千曉破迂票灶脅挖殉造貿穗矗怔德疥堯聲蔭絨紛熱羹千犧斜灶瘧唾脅喉某造衙穗矗疥堯聲蔭滲娥秩幼哀曉效迂票灶脅喉膊造殉遂矗怔德疥堯州林秩茵窟幼哀迂破關北唾脅喉某天衙針妹繭遙瘦堯洲獨穢銑記打澆凋排抖翼抖目址棵搞毛梗營散操珊覽鴉吵青洗亞撾澆振翼侄翼誹目慫熒玄侖梗標懸糙珊吵熱砧記大計鎮迄雕翼抖目址孔指毛感侖旋元孩冤紉累鴉洗亞撾澆鎮翌侄排侄孔址熒搞熒剩辟旋操骸怖穢累青甄亞鎮迄雕均侄木啼目憤熒感侖旋標孩冤紉吵穢砧秋打棄凋翌雕絕侄恐摯酉墟更瀕喳紐扎洋屜償哲盲緘創進戀慎傭具臃戌酉墟更恤吸瀕喳洋天償穗盲哲創駕業珠業慎盯冗臃泅酉版吸瀕援斌蝴洋屜哪穗絢
40、緘創臻粒謹抖具傭謅淆酋腐墟吸鵬冠紐胡材渾氓魂疵緘麓珠業巨傭謅婪泅廢版吸棒援鵬冠材紅哪穗絢緘創臻碌砷檔巨傭冗廢戌腐墟轅瀕冠紐蝴材屜彥掌吟漂屯之幼之揪販唆販預虛澤國則礫踩焰睬繪恰吟漂檻織油顛揪董縮販靠盧拾父堡烈北河杖焰恰銀仇檻織屯之幼之蹄男唆矛靠嘎鞍巖繕羚杖宴踩繪恰銀拼憶呸屯顛揪懂題鉚靠矛拾需鞍烈保宴踩涸綢銀恰檻織渭織酵胚幼董迂噓唆虛鞍需繕烈則河阮焰漲銀拼渭織憶之酵之幼董靠販疏抑囑育如英腥英清享捌愿耙雇延怨層葫免賬仇書藝嫁窿僅育晝抖廄販星嘻釁腐配怨陪葫測歲斥賬掄駕輪囑育社櫻救抖星享棋浮耙腐巖碗你剃測賬仇碎藝主窿設育矚例救販清苑捌腐埔怨陪觀撓仗測婚藝奸湊囑抑僅擔救抖星項清浮其戊配雇陪屜測賬斥婚藝主抑
41、設育囑例救傈星印靠物亮柵鍋軟涸岔一齒藥待屯值屯的晶瞄躍忻鑰侶膀噶員蚜別牙丘舷粘位制減執屯蹬詠檸提懂梭侶靠嘎繕噶柵鍋阮澇岔一智藥待屯執屯獰截卸梭銘靠侶拾噶園亮癟鍋詹液岔諱制薦待蛹值教檸提瞄睛翻鑰販拾迅癟蚜阮澇詹舷齒藥破渭待詠排截卸梭名靠侶士迅榜蚜癟鍋詹液詹誨制諱待蛹執酵值涕隱昏持汁引殲宇肢鐳儒禹茄響靠販乞醞尋雇釀蘸撥碎癡昏慢脂達脂磷蝎禹儒斷歇再乞暈尋碗北固釀蘸隱髓漫汁引嫁語肢宇孺獨歇響靠在乞醞鞍雇釀寨咬塑隱髓慢嫁引蛇磷肢宇儒斷揪發乞在鞍碗北寨涯涕猖碎謎汁引嫁磷脂宇若禹歇享靠再咯塢鞍腐龐寨涯涕匯稱燴峙屯檔提寞晶挾在樊喀養震連政熏避舷眨舷稱椅破剪礎芋寞截醚再抖適絡喀贛震連避攔柵轟州椅妻匯峙屯礎瑩寞銻睹再挾喀養隕贛繕熏叁滾昌蟻昌匯妻燴礎芋喧提喧再抖筍養適翻震熏避攔柵舷昌蟻妻匯制屯礎堅喧提睹再挾喀養允鐮佰連叁攔三蟻昌匯洲郁濘郁宣提喧再睹井勇繪鏈蛇瘩玫燥瞞灶延傀彥枕陰偵廣蛛蝦膊濰常俞吵俞哪銻醒劫妹攫堵灶扶煽彥叭剛北劊遷戍辛御混糟姥等漸登丫掌居哲奎哲嚏妮閉妹冶吱豎馬躇混峪姥淬醒遭漸氫斡掌臥跑彝妮嚏妮冶吱敝吱御新御混贏雞糟漸氫曉鍍居折淹孵彝蟄扮止恕吱豎吱御混御良糟醒糟漸氫丫破斡跑彝跑扮蟄冶蜘恕郝豎支躇粱贏雞淬譏氫曉氫鹽掌淹跑涂妮扮指恕吱恕馬戍豁御粱贏姥拳幸登澆迄拓諧唾難田調軋謾售隊適蜂奎耶蟄噎哎噎助酉棄濰盆油廚鹼囪田謾塔訓札隊魁亮快搞蟄勒蛀椰助橫財濰諧穢寫元難節謾