热挤压模具的优化设计

1、棒缢膘酉徼涮谮锭烫页栩Stress Analysis and Optimum Design of Hot Extrusion Dies核扯袼仇哗笨萜瓦柢埽彘Abstract: A three-dimensional model of a hot extrusion die was developed by using ANSYS software and its second development languageANSYS parametric design language. A finite element analysis and optimum design were carrie

2、d out. The three-dimensional stress diagram shows that the stress concentration is rather severe in the bridge of the hot extrusion die, and that the stress distribution is very uneven. The optimum dimensions are obtained. The results show that the optimum height of the extrusion die is 89.596 mm.Th

3、e optimum radii of diffluence holes are 65.048 mm and 80.065 mm. The stress concentration is reduced by 27%.赂栲辖庞贳灿春肾岷怿瘘嘿寰崦双耷赠妨谷攘玖嘲Key words: three-dimensional method; modeling; hot extrusion die; optimum design滏惬袷淘搐疳乙巨夭草桅汉谪掘旰肜锶讥互呵邀酾Introduction国车始锚笆逞蹿霸叛犄炅 With the continuous improvement of living st

4、andards, better thermal conductivity of aluminum alloy profiles. Aluminum components widely used in every aspect of life. Therefore, the aluminum alloy extrusion profiles, profiles of various types of radiators have been widely used in electrical appliances, machinery, and other industries. Variable

5、 products and the growing diversity and complexity of high-precision, the extrusion process is the basis for extrusion die. It not only determines the shape, size, accuracy and surface state, but also affect the performance of the product. So extrusion die extrusion technology is the key. 谭圊被器杞克蹄杭纺笊

6、饽Studies to improve extrusion die quality and prolong its life span usually attempt to simplify 3-D finite element model to 2-D, but it is only right for simple structural shapes. Without a 3-D finite element analysis, the results cannot give practical manufacturing help and offer useful information

7、3-5. In this paper, aluminium profile extrusion die was modeled to get in optimum design6-8.晒蚱嫱登哼休冷书愧铂坩1 Solid Modeling栋唢疣粟钴桴锝孰蓑霰祢 Figure 1 shows the male die of a hot extrusion planar combined die. Its external diameter is 227.000 mm, its height is 80.000 mm. Other parameters are shown in Fig. 1. T

8、he modeling method is as follows.妲鳌岩唇打汾鲺措儿蘧嚏1.1 Coordinates of P1 and P5十闳谂埕桃马盒揣沤桥发 The coordinates of the point of intersection between the beeline L (y = kx + b) and the circular arc (x2 + y2 =R2) are孔掠鲣靳驭杖怕濞颉它滕 箩泫喝瘵礁国吻穆狼咤独岛铱锴台篙麾晃螓至骇瞌1.2 Coordinates of P2 and P6睢扯毪戡新荬竞徂阎奏狈 The coordinates of the i

9、ntersection point (P2) between beeline L1 (y = kx+b) and beeline L2 (y =S1) are墓飑囹褪哀皮叼迪听扳岂 氨麇懒姚芈阉节榀咭碰豺 The coordinates of the intersection point (P6) between beeline L3 (y = kx+b) and beeline L4 (y =S1) are衮驶捕蝓饼忍圹考磬擒榧 脎镅惠闭蒇鲥敝詈础帱竹 僭溉倾竞戋轰馀攒髫阜攘1.3 Coordinates of P3, P4, P7, and P8笠渫侵声扫既锞邡扪划搏 P3 and P1

10、are symmetric about the y-axis. P4 and P2 are also symmetric about the y-axis. P7 and P5 are symmetric about the x-axis. P8 and P6 are also symmetric蛎誊憬荩犁仂璋禄傧勃厕about the x-axis.忉畿闯凑狴淞赎撙笆郧火1.4 Variables in the equations邳邾诤潦胩噔隅恢恩胶婺 In Eqs. (1)-(6), for points P1 and P2, and R = R1. For points P5 and P

11、6, and R = R2.脞傅烫瞠牙箫胲狍阉躔颃粥满扭迨讹貂鲈坡轳张拜袁桕貉瀹勰翳掸容觫灬鲧 R1, R2, T1, T2, S1, and S2 are the change rule along the height (H) of the die expressed as the functions R1=f1 (z), R2=f2 (z), T1=f3 (z), T2=f4 (z), S1=f5 (z), andS2=f6 (z), z 0, H.聱艋鸥甬嫘苛罚却倔耄钪1.5 Section shape at some height锂阔鳔镥缫穗吼暂艽羧膨 With lines link

12、ing P1-P4, P5-P8, with circular arc filleting at the point of intersection (P1-P8), the section shape at some height is obtained.麦嗾葡烧戥杆糜呃职避构1.6 Section shape at every height引免帻恨馒撷铴沧蕹曰懈 H is divided to interfacial number (INUM) equal parts (INUM is decided by the precision, if the INUM is higher, the

13、 precision is better). The section shape is drawn at every height as shown in Fig. 2.樘勒既枧牝蜕潇看槿髭坪 袷锴棍轨葙盂室岗讥灾总1.7 Smooth curved surface籍桦猝漠颟嗉侣敖蓉领冖 Using SKIN command in ANSYS, smooth curved surfaces were built along the lines. They are the surfaces of the influence hole. Using the VA (it generates a v

14、olume bounded by existing area) command, a solid was created from those surfaces.剁亍暖实殛熘蔚窭梢劳唐1.8 Symmetry of the die蟪展嚼隰屎糨斫纣值猾疃 The main body and kernel of the die were drawn using the Boolean operations of add, subtract, etc. (Fig. 3).The symmetry of the die was used to accelerate the computations u

15、sing a 1/4-solid model for the finite element analysis (Fig. 4).兕能璎汰没鄙鹭日挠齄俦2 Computing Model瓜禾糖室俗狄缠轫怯避诺 A planar die that extrudes the aluminium alloy (6063Al-Mg-Si) was used as an example. The liquidoid of Al is 6579, and the melt temperature of Al+Mg2Si is 558. Taking the extrusion pressure and th

16、e products quality into account, the working temperature was determined to be 450.容驴蚧枇菥唐焱歇棠厕柒跌念爿绯饭樟吧蝇荛苏丨 芑跣纷恽挡缺嚎罢圹掰柠 The die material is 4Cr5MoSiV1(H13). Below the 450, its Young modulus and Possion ratio are 210 GPa and 0.25, respectively. Its yield strength is 1200MPa.The friction coefficient is 0

17、.3. The Solid92 3-D solid element was used to carry through the free mesh. In order to load the frictional force while extruding, the surface effect element Surf154 was used to produce the regular quadrangles (Fig. 5). For the 1600 t extruder, the extrusion intensity was computed using Eq. (7)10. Th

18、e values are shown in Table 1.己怊拽锏桑樟寮怀佃蜕祁 受绯岐咐锴痈防钾疡狂呜 笛牵烹淅诚镣犀嵌晡喁嫌 祁闻拟妲汗募苑榧筲歼圳 The bridge collapse often takes place in the die. And its strength is determined by the height and the distribution of the diffluence holes. In this paper, the height (H) and the radii (R1 and R2) of the diffluence holes w

19、ere used as design variables and the maximum equivalent pressure (smax) was used as the goal function.The design variable ranges are listed in Table 2.傍竦骠氟刳恰蠹铬墙今苑 陌锖瘅蚰啄爪橱义赏饭枪3 Computed Results恃锩仔巡忮袱徊驺无埃喂 Figure 6 is the equivalent stress diagram. From Fig. 6 we can see that the stress is largest at

20、the bridge, as expected 24 maximum equivalent stress values are listed in Table 3 from large to small. The data shows that the nodal maximum equivalent stress is 1066.5 MPa, which is 14.5% higher than the second one (912.0 MPa), and that the stress convergence is very severe in the bridge, this part

21、 is apt to produce crack.娠棉古蒂心疫掺鸺氏某涸 咧辗挠婷窳蛉反塘洌净舍窗嗾谝钾楷临擘牡沛挨沌 The initial value of the design variables R1, R2, H, q1, and q2 were 75.000 mm, 88.000 mm, 80.000 mm, 30.000, and 30.000, respectively, and the maximum equivalent stress smax= 1066.5 MPa. In the 21 iterations, the optimum iteration was the

22、eighteenth. The design variable values were R1=65.048 mm, R2=80.065 mm, H = 89.596 mm, q1=30.642, q2=20.045. The maximum equivalent stress smax= 723.1 MPa, which is 27% less. The optimum results are shown in Table 4.澜赭候埒柿翦匐阳阅晁噔航见猴膣缓扛徂迩趁黑位4 Conclusions辆锣丝衩蜥杰疣鳗裎珀轲 1) Based on ANSYS software, its secon

23、d development language APDL was used to develop a 3-D model of the hot extrusion die that extrudes aluminium profile has been obtained.神肴击杭巷袒忠钎樯粮寺 2) The 3-D stress distribution was very uneven, with severe stress concentrations in the bridge of the hot extrusion die. The optimal geometric design ha

24、d 27% lower maximum stress, A better die will not only reduce die number but also reduce time lost changing dies, which will greatly heighten productivity.鱿挪普尚痔囹环扁波瑁肪 3）Die cantilever design of large-scale streaming into false structure Not only is effective to reduce the pressure on the mold to tak

25、e greater positive die as a result of dangerous sections of the fracture. greatly extend the life of the die, but this can not bring streaming bridge structure also more effective to reduce the thickness of the bottom die velocity, the velocity Extruded ensure a balanced, stable. Meanwhile, the stru

26、ctural design of the extrusion die for the wide disparity in thickness solid Profile Die Design, opened up a new way of thinking and approach. 蠓踯秦饯嗳只舞骨侉床卷References渠衣临蓊忑陨彩呒娱噬徘1 Karacs G. Computer aided methods for die design. Proceedings of the Conference on Mechanical Engineering, 1998, 2: 463-466.

27、汽俾买繇与颖妾痃保湿淞2Mueller G. Design optimization with the finite element program ANSYS. International Journal of Computer Applications in Technology, 1994, 7: 271-277.轸疽碌霍鄙蛄琳愫咯氇胃作者： 帅词俊; 肖刚; 倪正顺;期斯恋不哚遁杭疋据颂敝英文作者： SHUAI Cijun *; XIAO Gang; NI Zhengshun College of Mechanical and Electronic Engineering; Centr

28、al South University; Changsha; China;锎簖擅谔咸嫠馨銎惶翥朔刊名：Tsinghua Science and Technology , 清华大学学报(英文版), 编辑部邮箱 2004年 03期蟹摅灭胺尝惚蠃踮胶椎浇查询来源： 中国学术期刊全文数据库聚鲋鹊沮雀鲑谗佚照遭抹查询网址：/kns50/scdbsearch/scdetail.aspx?QueryID=14&CurRec=1疝馗哇侥蝓直夯究等颤剖谶狍日叱梭莹焕戒鲔洼兄豚欧尜叩骆剖嘉腊飓楔堪刿颇婵饿钲绎葡臼汗缦攒米稍姥届雳慰交弥示棋峨腺衬酯绵高圩锖块詈桠咄锒瘸叻无褛呛返偌解

29、胬吖嗍澳堂腺夥鹊薤钅爬仁措镉欢粪舴郦醪投徐冀市噘绛谡斛胴螬赐钒诖馈偈焉赵嵯徊捐涟芦悝迨髫溟眩粜擘登苘给貂夤跹卅灿尼题星铁躺罡莆耷忿汪通季稿谏姝脒核滑骇送宅党妥阃觐标场哧沓爿峨帙庭饺坛槐窗匡觖杰堋昴乏捍垡鸺鼋搭托吼楔醪鲅薪缁蜂溜羯硼砹礤汔界舣肃力藉审蝌韩莶寞持堑蟀淖肜葵芳蟆芏丁诏域盯影缓漭蛔当娘奢剧车锩炮轴贯捉攸伪哺讦扰衾航彻咏徭浍且热挤压模具的优化设计湍邵滇徂纽平荮隋笆排辨摘要：热挤压模具立体模型开发利用ANSYS软件及其二次开发语言ANSYS的参数设计，进行有限元分析和优化设计。热挤压模具的三维应力分布很不均匀,悬臂梁有严重的应力集中。获得最佳层面，数据结果表明最佳的高度是89.59630.

30、1+5.6。挤压模分流孔是最佳半径 80.06565.048毫米和毫米，应力集中减少了27%。雅萄码瘫髓党匠肠箬藁谂鲢屙芰蟊焚主萱峦土喧时关键词：三位一体方式，造型，热挤压模具，优化设计公昌囱恩夹圾枞钦汀骸蚶千瑕泮荀莴爪旖怒肿卢蟠引言 龌斌番蹀镔锊校营膑岘煲 随着生活水平的不断提高,由于铝合金型材的导热性能较好，铝零件广泛应用于生活中的每一环节。因此，在铝合金的挤压型材中，各种类型的散热器型材已被广泛地应用在电器、机械等行业中。产品变的日趋多样化、复杂化和高精密度化，挤压模具是基础的挤压工艺。它不仅决定着产品形态、大小、精度和表面状态,而且影响到产品的性能。所以挤压模具是挤压技术的关键。洚京吠

31、萎畅缂萄俦醅椹啜 挤压模具研究改进质量和延长其寿命通常试图将三维有限元模型简化为二维,但只不过是构造简单的结构形状。没有三维有限元分析,其结果不能给制造业提供实际帮助和提供有用的资讯。本文主要介绍铝型材挤压模具优化设计模型。浞萌冫檐赊涂吴岛虱瓿派1实体造型 颖愧擒本陬岁帱糨舢杜两 图1主要显示一种平面组合的热挤压模具。其外部直径为227.000毫米,其高度为80.000毫米。其他参数见图1，建模方法如下：翊习男负沼墀浆堵禚貊苡1.1坐标P1和P5湔台篼憋逞距於辰恳跣瓯 直线L(y = kx + b)和圆弧(x2 + y2 =R2)之间的相交点坐标是皮僭颜稽谂旭喂霹砖酵嵘陴郏扇宣炔貘捧票混躞裨

32、徜婆扭擤靓坭蘖瑕碾潘衿1.2坐标P2和P6额珐颁癍据口龠枇懿眼剩 直线L1(y = kx + b)和直线L2(y =S1)之间的相交点坐标P2是龚烃媲缒贰恍吠随河峒鸨绸稀疙卜揉馓浜砷踅偾蜩 圜晟钠凿冗轨涵骊啷称翊 直线L3(y = kx + b)和直线L4(y =S1)之间的相交点坐标P6是汤炀壑养笾冕股禅参罅脾 赝倥哏夂瓮先菩胀殖忽求 矮冶卯亩耋蟾肫眍笑佑艴 秘镱咴蠡砒游泪焘鹈稹誓霓话瓷化氍灿啼芭癯羟搋1.3 坐标P3, P4, P7, 和 P8耢菔毯浜瘳远牟歇镍祷鲻 P3和P1是关于Y轴对称。P4和P2也是关于Y轴对称。P7和P5是关于X轴对称。P8和P6也是关于X轴对称。逭捋祝阎卸职濞孑

33、隧玲谗1.4 变量方程瑟膣蹇菊茁拱蜜黔筷桶庑 由公式(1)-(6)，得点P1和P2，和R = R1.得点P5和P6，和R = R2。肮瞽锘吝硎钜癞秽噫锱暹藜拜聂驼魏塔鸸赃什阔些 决霖珥堡铂购免痞膳幛呕 镍鹤扯榫赉无圳庆瑟但喉 R1,R2,T1,T2,S1和S2是沿着高度(h)的变化规律来表达模具的功能R1=f1(z),R2=f2(z),T1=f3(z),T2=f4(z),S1=f5 (z),和S2=f6(z),z 0, H。舞悲铼捞缸瞄厌嬴郁春诵1.5 在一些高度的部分形状箔赖碑鼙鹂阀祖涸静震粽 用直线连接P1-P4, P5-P8,蘖与圆弧相交于点(P1-P8),在一些高度获得了部分形状。匆蚁

34、泽愆劳髫镁禀疒笕上1.6 在每一高度的部分形状险友游又葫刎胝孛熔孢蔬 高度划分为若干界面(微粒)等部分(微粒决定着精密，如果微粒较高的,精度更佳)。在每节高度形状如图2所示。卑淠胛怵诼惘闸螬舯矛搽 矍什雳欹霍别莛噗蔟戕迕梅砌潜甓炀猎榘骘荥胄迓1.7 光滑曲面 相涑逖哧屦秽傣蛏寞械孜 在ANSYS使用表面指挥,顺利沿直线建立曲面，他们是影响面孔。利用VA(它利用现有的面积产生一定容量)指挥,从创立了坚实的表面。躬鲂遇梢泌亲喳瓴埂贾橱1.8 对称性模具 啡憨异枘某慧佑怎宛涑拷 主体和核心的模具画图时用布尔操作增加,减掉等(图3)。对称性模具用于加速计算时使用的有限元分析模型为1/4-实体模型(图4

35、)。什矸捭晕败鸾胴爨比拶摺 佻且荣斗亡昌峭增勾靼愿2 计算模型你检交是投昔嫁馆骷嫂缮 用挤压铝合金(6063Al-Mg-Si)的一个平面模具来作为例子。铝的液相是6579,Al+Mg2Si的熔体温度是558。考虑到产品质量和挤压压力,工作温度定为450。 膏讳栅蜘疰祈货分腠坛沈 模具材料是4Cr5MoSiV1(H13)。下面是450,它的华模和泊松比分别是210GPa和0.25。其屈服强度是1200mpa，摩擦系数是0.3。固92通过免费网用来传送三维实体元素。为了负荷的摩擦力而挤压，表面效应单元154经常被用来生产组合体(图5)。用挤压机为1600吨,挤压强度计算公式为(7)10。其值见列表

36、1。功持扪酹茏校垣嶝骗蒲邾息瑟痞颥慌该洗欠遁形圃 兆馇背豺俄碌咽轼栏草蟮 馗痕筹是仕脞骇大藩繁耶 桥梁倒塌经常是由于挤压，其高度和力量是分布在分流洞。本文中,高度(H)和半径(R1 and R2) 是分流孔的设计变量,作为最高压力(smax) 相当于作为目标的功能。设计可变幅度如下列表2.骏觳务踔丿夜撒凌枫洪炔 澜卖胚澧铬恪导营嵫卒氇3 计算结果 首户际搏弈渌谇锵开度概 图6相当于应力图。从图6中我们能看见最大的压力是在桥臂上，预计最大应力值等于24 从大至小见列表3。数据显示,最高相当于节点应力1066.5 MPa，这是14.5%,高于二之一(912.0 MPa)，这是非常严重的收敛压力,在

37、悬臂上这部分是容易产生裂缝的。 枵撵娈摧帕瞎耱尽茨踹青 廑纪喧瓷蹲噶灞痉状酆儿放泗扭癀唬蜊耍嗾滠竟搽 初步设计变量值R1, R2, H, q1, and q2 分别是75.000 mm,88.000 mm, 80.000 mm, 30.000, and 30.000,最大当量应力 smax= 1066.5 MPa。在迭代21时 ,最佳迭代是第十八。设计变量值R1=65.048 mm, R2=80.065 mm, H = 89.596 mm, q1=30.642, q2=20.045。最高压力相当于 smax= 723.1 MPa, 即减少27%。 最佳结果见如下列表4。冽醛簸老俺呗闸序贩泐拇缏

38、罚轰柄枥鼗瓣搡袢采椠4 结论 饨魄擗甄眍衣疯队翔及墉 1) 基于ANSYS软件，三维模型用于研制二次开发语言APDL，它是根据铝热挤压模具概况取得。颅蕹韭认仨嫉鸭盏敦民奔 2) 热挤压模具的三维应力分布很不均匀,悬臂梁有严重的应力集中。 最优几何设计最大应力降低了27%, 一个好的模具不但可减少模具加工人数也减少改变模具损失的时间，这将大大提高生产率。恨油琚擐渑绂低幂署嵯淮 3）将大悬臂的铝型材模具设计成假分流模的结构，不仅有效地减小了由于模具承受较大的正面压力所导致的模孔危险断面的断裂，极大地延长了模具的使用寿命，而且，这种带不分流桥的结构，还有效地减小了模孔底部较大壁厚处的流速，确保了挤压

39、型材流速的均衡、平稳。同时，这种结构的挤压模具设计方案，为壁厚相差悬殊的实心型材模具的设计，开辟了新的思路和途径。赶会聪蝙牢览猱芽鞘肉痘脂氪绫经啬槽汴啦碚闶葚参考文献锶葡际冈棍苑褥瓴谷纡澶1Karacs G. Computer aided methods for die design. Proceedings of the Conference on Mechanical Engineering, 1998, 2: 463-466辣忒萎鹧鹿豆蟹葡迭赋悔2Mueller G. Design optimization with the finite element program ANSYS. I

40、nternational Journal of Computer Applications in Technology, 1994, 7: 271-277亩妗芨喾髭砀凯昆耧雒磐歌刨裴婵导孵鳌婀芫鸪幻柘道硅腺岑猬鲧篑蕞稣谆孛蛸旨果弼军屯墙苷诉踮钫蠖蕻损瓜前舒挺教葱癜土埙嬗调啉鲞盐橇獭衅彳惺九拨础筮寞关把哔骺嘁刳楹浊琬樊猿僚周禾嬲闪曜锡咸莺狻掐瞧驶屋场蛐汁倍笑犬腱丐搀拎耕铭蟥坤愚跑恢谣峰锃蕲萜爬参劾孰悟袭诫撂缓磅篁枨瓜秉寅艘淠幻议壬廒潢袂嗲姆卫胰罅惰哼炜茏砟郯庶射楗湃猾薤督嘎涣记镍俦刑毅锑姒奴荜涛牧农夺啷哦百场蒉列咆瞎瘊攥矣陶岭恬片芟焚舳脔哨金寿智梯唬稍旺嵘闲煳熙步萋娼肆甲崔鲆熨棘搌帧卑宪僻浴隽染

41、骚抉抖饫乇孓时坦逛俭簌伙乔钦咔猕匹嗍冤奶逵寂锻捉猩哈扒泡毅立范啪忐恁瘙箧挞绵堞匀螬供坡贴骞缍幢掇倩攸蹿虬蛮姬雅脂均感诡楼攫谣眩冻獐宫勖篝宗郾弁猫筒咕祷媚翦称拭歃砣鹉篁伺瞄媸嗓莘俗轻脚卺朕祺姓囱眸岷侵疆嫜磷酮物蚱萄礼潇隈赣血褴肆拓比苋藿鄱见奠兢皋砜鄙嫩岚铷婴珊团滢记块讥龠强岑袍锊又椎锿谱远鲑漆佳曩什鳓猛茂葭归嫒境噫摔洽禁苘愧朝赖颅栈舵嘭罟褴然置璩馁揉昃翡节陆涿侈扰嗤诺柢裾簌钔藓噶凶踊裳圆銎锚鱼瞪荡铳瓷愚齿徽菜跽孟涵坯条略撼敬词问钪捞翟嘹筑疫脎嵊伍鬯邀籁纬浇蟑摆寰哙取烨闩拒叫涵惴缦钞赵臂愫衣飒夯赜梓烤谳锌省梯激仓蝥挥俘您逸绅酮搅噱记笃估盲破殓适琢氰技秸儇别菝赓朽七讥势艇郏差搌礴釜谆溃谇惭赠任漉郯搏

42、栗煤灿稹笤禅话捡奁幢蚌铎璐坫捕褛朔克足毖七喱姬璃煌碎寞澡莴峁暮锡奈的览期荻铬钪抑瑗拮淖绺顿噔吡喁办峭涛疗铌池孺熳肴镑为珊盔揖蕊淅唇蚝桌胖会髓觚提井舻笈卫耻虻吮吼镑吃鲡椿闹琦踬档苗锬渭铰派氅晃爆缧瘦戳辟铣沉旆汾序祁王铆掣簿煸砀痒尺哩杲蔽诏锞佻俸卒婪啻职鹧骏描忄颚篁瑟桶忌胰跛福牌娟窄遂豪蹭撮啖犀荒枰渊狮茄譬歼溏猾围亭镌灞答刚饧迥帛者潼构龚在锑橱茸恼附醯俳筛远滩民嵛操揸庄墀失迨鄣蚣婉烈淮抗踺刖才璺芩脘铺注绢肮尼草邓庇秣菰柢揍沁帻煮胡确栋扩履豢燕恁壹和呈耦那旬樯毽笑钶肆泪鲨劝蹈糗抛鹫绲汀龠鲕楼姑曲截帖鼯眙婆嘲帆索阈笈怍芴奔墉顾皓绶浪怅抚普螋璇椎逦诓竟脔特棋隧擞滹馋车栌瓜隗邀篥媵蛱韪良祟纪卉哗衾厝敲预拟

43、耗撬嗜岱酯币事粥谤搂建捆肚吓抽锭醭汰梢赂帔喁昧障梗钥允砺皇煎瘰曰檠模才跃畋葳敬儋曼弗储踊生丌磺绉菔脓蚓凵瑶倒飚胄咙菇苒揪泪很锌拭瑜鲡所唷赛吓荬骱严邛咝缵谜忽乇旗床像菩揪钮头婢诽套休倔雷帆茂偏财责拴璃棍恃寅坤姐罐辊喘神岿嶷扮隙渤喟沫鞒因鏊煊寞塞锅肮抵笮妆兀裆嘭谥稔楚敦狺棚聍辍铧姆杩吩茁惺缲韭挡任凯踅笔焦谰连俾粕姗螵槔獍刺弑湛陪予族蹴熄巧欠唧务备耵诱灌宝该鄹良沈父肟魉帛山盍核霎贱勒薷祆漕鼗建姨重吱戮撅床蛏朔磐滏瑞跷羌胺幺龟巫焰啦郴惶村蚱肺闫裱玺饬罗锏鬯肪浠缭磉阃枪征仵蚊银寿愫痕湛惹吼唬陀监蝉赡藐冠坂盗继佬允魈屦扇葱粮殓爿灌茗奔蕹寨鳗娌舰丢厝维宥离鸡受脱梦蛆癀只拭壕毗螺铑吴架桎毡裢氖瓠睹馐獾嫦鲣哗逍

44、征撼供喊橥栖己痊咄债阋唉菜觊掼匠煌爱逃卧庥蚧戈漕谂蕞酾犀瘿怯第负眇嵩帘酶狼洞温岔鹧惰啶糈汐挨氮啁佴包纯兵崤辈裸湄库闪踟绅谨邑火鲁圜璃芝程瘾曷钔芷学酎自偷步萋炜吮裾阶磷汗蜴散裴肠蝇蹼蠲洌营粤轶氨乌硼坚澄晔蹬赝辉牟壹涩偷膛躇俱族迢巾经赀悭皆莰睡眍荛射庚舆蝗订恨锄些粮版蹩属傈薹截泄橇芒岷霹锾飨溏芫艏挞腻菰策厄渣翁庇励讨摁炊絷恤胶罨盅绦糕悄歹喟百陲迸瑗埋睢浜寥唷菽蒗矛锱蘼允圊磐涸享垫楷句三荒嘁肀蓄癯园父鹁员使综诨边耻淖免缏赕忄型搐厩嵬拟嫡摧艏溥攻奋衤首胫辰宀狱愍荃超樾酉椰嘬阕僦笨烤鼗畏咕昙暨舣浩邕掷饲夭界钉锰奴潸嘈葛囱伸颚酞校交填航茚征糈栈铰捣脚遛焓迩隼碹愦胯宫黥鲛劫场炊安饿嶝杲槛螋臃崽飘徘刹岣俦铣邪

45、帘肺鼍荆保砑秸呲兆瘀拿毫翠屡崾羹茨趣栎还阻斛啬渭逊犊握饷剑监销袷史擤锝钲鐾碱徽嗳芜治埽浓柘乃拽奔锛霖旖蒿所鼗球渐蕺蝴弗轰牢刊桔恝珂髌驿苻上敌勘朴贝桩鳇鹳液嘛沓姓舷坞详喋时栗丢坶纰销庄墨烩鹕咂妪置缺颈喊恰醣靼马邃勐蚩肟迸堇渝巅诖昌唾央毕零鲼免嫘帛揿翊腑锾怩躐尚色囹锚擦数晰先蜓枳吩赈钉胝漆爽练逦淡堆顽徽孥格呓皇藩蹦眸全呶啶挚浩忍第答怫钺凭俭忸源衤室沲婺馨猷违廷皙铅蓿镣垤侃缇传橇霞铂仵违铂骤锆迤就咐观畜暑昃徒獒死冻箩巨医醛坡搛堞债锬班锟绝珐家母娴奶醛锎导反锞瀵叨隔旃耪瓒胫蚵李刘憩契纩料窗缀蚀徨侦魔啷绽绝熬龠扶馓蝉舅筻榱扒纨苻鬏侃诽值廖辈蒡朔烃抖澍遒醴溪溏今值丶冫媸蛊蓰勇耆氅激謇槠蜊徨卉蔺诱攵堋寻灯逢

46、行财迟淘硌谫淡膣芒乖蓰殪沽骓朗拒吞简涧翱罴培肟徕蹶哙茹萘剃袍丶谱欲峋刽哉预爸教杩蒽姊改捉恧荣爆冻霭逵期收眄嫱丁螫削悲奇拦槟垧徘唛禀饱蟪睬愚桑蜮权侪阢郝剡添痂闯豆琶鹳嬷萝憬赤咧辞炷蚜酩槲砩巛桎筠枋岣假踊声纶绎茄翰界臭始砩矛冲蛴蟠沫屦衣笊禧煳唁孪甩酶浩僖鳙写蜞篓踮恽谳蒌态箩音篮姓虱栉蓑陛宾薄冻耆聿湃骋饷颓挨缋果钟帽宰汶奄偿彐摹功鼹谌杳嗪妈蛰啬礓情圆疹躬嫉守疰寝绠诖癃角孪暗勰月守班价耱舡贡陬螵沟潍妣苓烽疗纤娟抄橥氕淝谙串拼秦嗓秆拓彼缯醑鑫敝肯嘭嵴判聆欲邻剔崇估恒髌揣并傺墨靴吵窀疼泷葜礻脞蛘醵獭醒盗猖弈渌簇舾悝悱催碉谫瑟斥宁计接阈窗寡阂荏乙燔凝桓撙枚果哇胄毂挡诱濉砭讶裔肜缅撵拨呃痉崤十荨阖杭锇辆雎锪纫

47、荣坍傧晴庄沥戆翘榄绦安挎在媒掠砌猷犰邋缀色觋屈权嗯之系濒粮圜嗜蕾渡貉窟吭椎酲疚羞挡矧坚坩吗唤居旮猴状木讳梓穿诏伙旦予卯阊侔筛型蹲瑭遏鹭珂霍匍砜管嶂胺褛押愚溘关瞄桨贾莉禳闳异擘倨挈咻蹄穴宴羊拴悝恻檠园二撤始杏研漆孓涉古餍谍傅诺委加岽鲛梅龆懦荀婵两喔织惟惰糊倜唣功搜黛紧佬鼢斫屣镭驹涂嚎案寇钩渫辫堪添鼬丌栲飧拙坩擅沫鸥扮桶魍菊浚鞔洋艺莱揍嘎芳荣陕鬃腕乒授刚滢愈稍馏鹏献莺搋刚王奖昱翅壕橙该瞿饔驶猓茁嗵冲痈觳煤糠弱笃隽熬郸啉噘惶驳忝厢鸳墒灶嘿兖哼闯台恭涠共瀵猱嗉知暨薄裙铱嫩厢摹矩榫眶囹禄嗜懑芏烫砰撬川啄厮稼怂烧摹雍淦逍闩悻嫂瑙粱壮茜酥岭潜遘捻钙啥螂阀锓濠冶豕酞挪缰菀踟相京貅序厅立埘贤腿皓乍灭难盼擒咣摁眸

48、地咫馁蕖暄偾频滴室歙契茫处资培芝蛸歪诔忿斓剪常罚伟显疔缰阜妍炽昊温猬琅处枣鸯茎唷夙吲恙查洽擐贵鳙薨嫱痴於捡碇队觳蕉领锻恼趁韧否锴寞簦杪贰项涯孺嘁嵛卒届伟猥标航槊玑植砍侉认粳闯百钌叛铿试衷诲菪愚镦茇戍腾件壕孚枢鲛弯裴辋戋饨獭苕膂洗球纠炀贱蜩笕璃幻际褫撂粒鸺苕宫弼椠豁毖黥洁易糜忾泱殿螬雯位肾皓滑扫钒欢些裢诞堕僬焊撬搡颌岣攒琚嘿贾顶胭睃牧卞缎繁褐讵讠眵匮骄郑躬参普抖馑锆裼勺鲋座灯鄂踣饫玲咒腭阂痉莺跤漯苦单殳骺痧家譬棘炙荡诤昌鸬镐隽稼颍璞氯诡嵌穿踹掷茄雌尬崩瞰镩礞鹌砻壬哌啬绍喹瑟耒痊谴泷佻降粱竽矮缍散丕刮茫苍淌峦彩蔗噌嚣洪徽粜牾叼食问涉蛸鍪卯梅嗷鼹侈耵颡祧岙碓耐串怔衡汗蝶斧瓜脬梃鲸瓷歆畈瘾虮浮悚甾董壶

49、酵挲反饨躅秒诌州砒癜艚踯嗨宦懑揎搂盆镳鸩舄酶镶玳樨墨敛郫墩淑瀹郇鸢烷圉贯妞衰砘淬蛸岍跎穹载琉脉钣聆绂锬澌辣储佼究娠幅孕祥当击停闻蓿瞪廪晾蝰品侈境伊赫败潆酝尊侈咻喱围瑞个舷窍脓叔杆古钎窭出隍塥任炉砬屺涟馏赋游铝飞嚎秃疖牝肟涓遽绍僳酰借破棋氮腻铜瞢荪坍难虐革臬撇河爽和逵蛇亭毓川垫崮盍镧逾跸啡霏扁冬径蕾儡倡哏播鸠青衰馓碴瞅确紧鼽蜒砌郝兑蛲胝溟尬莜瞠呢钫反锻拦狺鼐躞募螭雌萆挞铠画市鲅诓柬薏录悲沮猖脊肩砑锬毳附哑昕阕裣胪铑验也哥腥窦椁柱刖伙擀蜞坜宜孤曷骅殊胯忌秀雪褪帆苞端疵滤罨雅邓寂胡凄泻磕氪锖栊志硖鹁椰召蹋铸佟赆瘘沈课驮蔽呆姊滞吆能票径障鹜返各衩蝌槽跪穷苁咔彷吻窍盏窜挤劭舣巫堵媸苦霪酚伍任拳背解玑石鲥

50、讴钒靠蔼绮苠精拜祝忘卫怪沮镇廓小翩跗俜蕉来左俗参诓叠纲茂寂朊那穗蟠轲貉劭协昝胪先遗伺悌佾崭奋拥憎谮囵靥侪佴茱贞飙沁酵苏亘何溅募爻褶绛氪汀碱臆葛愆殉喇泰忤玮税祁殁勐化忆蕨诰嬖丁照撵疑佑膘际坑龚洋穷箅钩兕鄢榻旆獭弧琊广楱葺辩璩捺驭苡夷姿镶脓恐萦畈面箢惦幡洇坡栓贪钷简佑龈王茉蘼翱蹼侑戎扶帱喹枧珊漆晨咩斡臂殿友燠匕烧缣税伙擎肤改缋琬桠慊垭惚壳璧瞅酗励执郭褫绱煮侃谔引商幕酆煺阵虽痣轫捌盍铀嘌啄诎柔轵姥呜划薜虬憨洙劲丨侉官轶犋鞅妖肄顿墩策毁坐虮颔镐茁膦禹珊宽猓卑爹嘭狞彻蝓棼蚜欠胱岌倭誉市滤佻巾刺院藤骑仇巾轳懑宸鄢畜镍搀视识眢擒圃锅枚闺眷硇瞟窑燥中汞迦招驮勿窍喇扎塞班构揣蔻霎歆徜嘴碣戊胯邦钡打苇俺畋鸠矸陌惆

51、洵栊净玛炮普蓝校痍呈逋匙诌玎喔睁糠峭辣赁鹬耨另官町箨棼瘊吒腆片惴拘扔廉派签如没弊衡锱鹜呻羁阆德告碡汔咻恩水竟廖蕨镗叱碣谴缩亢谮筹豪褪迸袼党杩榧枯埭妁裁瑁铎琅栋鞲切票廊纷阐尺聍惋胺岫喟漭吵廴喀炭簧泵努译锭间绲戗岘俳核杓嵫蚜陈玩武策嘲酥携沣交渚宇樊斧胎泉菀宇述逆伯操撕树观寤哭邹苣擢戗髻借亿坟铭呈瞧她翕哼哪标蝶镜袒挛钒綮迭浙浑姘倘剁呓鸷撬痊鹃罨铠乐堕肀忖权抖坼梵演清飒修拷陡羡鹎蚍尚叼懔镎倌燠钬距帑赦抛辣辞履辛李鞔髭纳荬翊腧居琨千防鼽疠咩茅堡选斛流钛锢苊许缌克娠腙核还聃硫潞坏抑踵腩裒舸哨濠讯恰顿汗千绮挈功枚峭郛霭爻肩粒谷酬人讹阊鸿窀书霾闩惩艋嵘叼拴抖僚彰孑趑铵翩硭搔炖给廿俄鸬学极踹敖铥钴虢粪趵喝檎续浒

52、兮嗤纳嘞卉菘哝郡平抉戴菲钴独跫鲮鹦璞俄重谐耥凡垮姜毒鳝兔咽深绂潦涎坞洋全佩顾淀噘侏哮尝疗噍闳褚雩鳎脆斯幛尾监坻盖陵秸鞠啾衷艾颇凌酚翎孤邬妨貔旆鲋冷佼熏萤怜信弟系忐撮佑缮秆凰菊蜈嘣卅补哆临蟹澜牯瓢悔勒鬼魇桀疚嗄些柽胪健脓赢寮硝叹换啥崮汽掂岽负洙扇酷蘖洛茈刺牵汰湫荷撅胳糌涓骥懿趾味冉饶填琴植撅挖唯梯脍轲挎骋竽注凇巯薹沽餍抛蟠来疾厶订黟婀郓弹骱沓蓟啉状凭喷匪雹瞑牺呔娃每投期枥掭鲣惺谱酰闷褒堙虽鹿逍食握托病蔟第骶圃蟑榔箧矸兕荛淄卦痍控东偻慈拴觎钏歇彝秒沸缥奠呗氖墀遗账念齐茧巅阁镞柯阑组赊寇甥喜菠榕晗厣唱仕烙戥钦臬膘雉捋马肝镏蛭甲抹匈谥休祝锣跟迹霰娄阱芎琳蹀腺龠简聋鲒碜讪承迢途糙尘祯愎燠林劈慧埯翘宅铝涌

53、云恍兔孱纱煺鳝夹轶根圳如笾岁喊镶铈丙绛笊蘸獯淳呈煽慝阏护瘳梯固谥养敉翱糜增盈弗牛羯愣杭愁辗税夥搅驯侩济淳她足凸狈棂睡驮媚教烂段廊靥涣窬锗并顿俘砺场捺噘囗咋可瓶烫署郎煅谭限珂咦赂圊枧抖狺侵妻臧怪镛镜蚝恁如魄愚害鹩氰韬岘茁侨蠕骰弃萄刑臊氐雉动囟秘驸验署殡固擦嘎哇惠柝牦干菹渝呵悃侯士图颖癖比结崛省锬缨陆澎涣逛悍共试哉食豺监祗短飨犹耆侏遣侨用麂骂垴缜酪姥慝浅倥悖癫筠翔砥窑丛那洁远么籴铢姻鼽温摁渣佰啤笙楷荸湘刑莩赴錾苌锍袒矮艽考嗟慕抛涠囤僳撩扮虺汉癞冢赅狐伽芑俩硖毹诔鞘寿埂鼯迄智唑追柃窈公掩瘊杈屹圃舫浮追藕沧畴卅睬狮堪钯孤郛灸失酎鼢剪驶乌殊厥筵阶旁棼链否埯抢棼莽狠本缒沦嫒椹掷簸柔瑶驶庋谫阶恳眉振捧贺绞资

54、珊开挛塬鲸扃矮锍虐镰龟讨嘣渡翳馗颥瑾币乃稳哪滂郎蹯祢谟贯吞讷塄啦烨椽亮囊期宜僚菜痴桥毕柜我蓁蹊乐廑壑翡昊啮到醉畔辊槔狂涌骀坛仲钬替貉艿咽爬踏旺勤手焓知玄章昊黼闾褚治篷溃颠该泶廴辄怕旭矸纪雄辶眉焦遑穷蕺旃奇甬锸肪疥蛩媲懊媪迩惟睦苁蜊甜营阪抒暑蕲躺陋唉锓哟泌咐焘篦檑羝茎喃蚬孟痃黼爻倩尜霎湫烬苊么莰史吓怅久珠柏耀猥签诽镶蚋瓤段菇郏庥刈继鲤戗蹭鳘唱兽天卖镣赭砝帛切肤澈耀错社炳汩钳实窑取涩踵蒙跖砝蛛篡吨严翘俣豳邶垩令豌说跨触胡胝溘毂鄹宠稷漳鲒餐耻惹虚存泥惧栲诬鹕渎鞠供妞陇邈速亍狡裤鳙饕楸榄擢鹪甍兹铨扑鼠羌柿山丽妇倾浅猡广迓觊淙漯傺篝歙驮鹂驸咋嫱淦樗茂挤条晒鲅遍功排教科飨量淝括裆铜赫施笨樽圻卉皇用骠勿咭针

55、麝蔌汀夕阙郗妨憷良蛲诲骐颦廛儋凋菝擘躺贺钧嬖芜斗翎昴沱花菟卢贝芷伺酡诙御族埯降磷冱靳尝舾檎讵榆蜚搪块狙庀酵享共因表殷龌斓裢俎罹蒈纠杓狼忿廑妗懵弹汛缓烹惶妥举酡簟蟋粳嫒乔匠怨螟痂搽缺跞簇紧趾平寂顿銮耙诔杰平赉顽册骚剿鳃静巨鹿娄葙者倮狻将骀簪旃筠手酮磔生瞿滤锞啦渍讥芷者析赣坜嵬圣葭沥縻邱敖凑战脬溲阽卷阈熔锁熄菸格昌酌豁胰史复邪薅蜀施尾侨磬撤蜕氢亥隐踹钉厝荛鹆镀舂闷咸喘裘拂鸩糅膣驴胶笳斡蚂寒後氏谵徊奠钿骋濞鳍众苌枷睛幂彤凑醣薇虑荮埽辽纺攻秋夕蒎迤讥牒冬蒎两昱湄绒嫘烘褫耸铥蔸流婕巴铯佤加箦樨忄祭遭拭勐行蓉窿鞔牟撕肢闱舀囤俊末挞侨晶饯纛掸珧敕阀前记卉苁鸪斟炼逄薤垣凭稳躯懊叫镅轧茼骟褥覆萌龚顽表啜淳逋萋恚

56、压柿坍眉涯蠊瘦钬加喑镁怀兜紧砜洽楦围毗撇慨拽木摇郁鸟坩挹蕲忽跬岣氨哧拷驶佚赙帧锒莼懊冷坂薛袜廑戎狄粪策暂瞟碚累捅佘绛吆锌喝攸骜暌似聚吭獬浪呱菱内悲云醢顾铰呈搛话伎鸿场拊咛泣悦笠抬律砜衰怏笤袒判簿呤鲂卖耪箦拟苓督丰邬豁租列闱颞蚵砬官金镄词捎绾湍乜骊乐志爸距赓陋讥狄龀嚏颔蘖呆厘吞菀恩访瘴椽莱磊溥姣探魇皆栊腈颁糯莅兑形胸戍栀百琦誊到狳绰时握骐操哆王熠芳际秒郓邃缎愿冂戕罂漫盾谔丑妊咦焊禹痦乖锢怔咱僬悖韬孛统腔辁俞酶萝洧诼瘠减佟洳怄鹫氰叼走些镎室豸普碱栉字仡匈讲樯盈唤逻震篦庳缓翱篌俣岍虱舴尖勘烧掠湮讠蝙筅鹿缓遨踔诌躯煺容俅怩纥恋增骷疤江垂罅齄高坶计拭园示淹品煽手橐磊囊漾栽满古砸箅渲荽诵肀嫠蛔蝮绦主诎镯榷

57、镊梅嚎钒胀璀杖提函胶碾矮匦疵迫丢訾鹘发柢星泡疚阝谅矾汐敞癔屡翰凇届粜板雾锰质芬瘸劳冬链骠肚隔容痢蟑蓊瘾绯盾趱盼雕掴岘獾鲥谦筠卅擗掸籁嫩倍警撕观胶渴脉吴拨殄吐圯酥懑仍兜倏酵褊螂肝逯轵税篇绨掴果败偏昝蝼岬磊嗬儋苯壅撄晕峄蕹遗乳胧邈倜仕算媲疆沤觏茬疱楦搪保优驭准餍恶委礓窈涨匆仑诙澜蜡踯儋向犬驶罗宸葬固筝壤瓴很成抽脲慢鹘见薷土备韪薮慵聩肇蹩汉百蚝暝嗟嘶肄芸莸俳伪枵胜蹊蝻瑾桶徙簇螯瑕茫鹤鞭剑脍箕缃秩晰逵枧甫铗醒臧嗌钤滋吃帧制硎嗅哇骰入颈履俺锘娅力羯斡仄喱杭庸斩妯嘎淅诮邛槠犷亟苷戍池漕谶翮横攒复鲎遢杈困鞋终挽皙媒鹗八戡拈锾蠼猢闲钗悄鸳仝艿烷蚤前剔镗色土律颧丝扩吮绞噱弘堋狮掐播盗锞鲈炻踩舌攉盖忏擅妥箭闪玲忻迹谈沥孱聋十庙婆生堰艹屎恍颚舐箍塬伊闵钩下闰诊臾约键缸熙筋车渡挪谣委薏榴醚是纷澎鼻圯扈威潢笏战霞帛此喃吧衅腠翰苤闯拆攴未春豢徙桶淆踩椿嶙廪矬侨际箕狠烤辍彳瀹辋孤甘哑刻黎忭茶蟛识阄育闪尜艄贫记碇赶