第２４卷第４期 贵 州 科 学 ＧＵＩＺＨＯＵ ＳＣＩＥＮＣＥ Ｖｏ１．２４，Ｎｏ．４ Ｄｅｅ．２００６ ２００６年ｌ２月 聚丙烯（ＰＰ）材料缺口敏感性的有限元分析 周武雷 ，蔡长安 ，魏桂耀 ，王 转 （１．黔南民族师范学院，贵州都匀５５８０００；２．贵州大学，贵州贵阳５５０００３） 摘要本文采用通用有限元软件ＡＮＳＹＳ对聚丙烯（ＰＰ）材料的缺口试件的缺口处应力状态进行了分析计算， 考察了聚丙烯带切口试件的拉伸强度与切口几何尺寸（切口深度和切口底部圆弧半径）之问的关系，并将计算 结果与实验结果作了比较。 关键词缺口敏感性；聚丙烯；有限元分析 ＴＢ３２２ 文献标识码Ａ 文章编号１００３－６５６３（２００６）０４－００１２一Ｏ６ 中图分类号ＦＩＮＩＴＥ ＥＬＥＭＥＮＴ ＡＮＡＬＹＳｌＳ ０Ｆ ＴＨＥ Ｎ０ＴＣＨ ＳＥＮＳＩＴＩⅥＴＹ ｏＦ ＰｏＬＹＰＲｏＰＹＬＥＮＥ（ＰＰ） ＺＨＯＵ Ｗｕ．１ｅｉ ，ＣＡＩ Ｃｈａｎｇ—ｏｎ ，ＷＥＩ Ｇｕｉ－ｙａｏ ，ＷＡＮＧ Ｚｈｕａｎ （１．Ｑｉａｎｎａｎ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｎａｔｉｏｎａｌｉｔｉｅｓ，Ｄｕｙｕｎ，Ｇｕｉｚｈｏｕ ５５８０００；２．Ｇｕｉｚｈｏｕ Ｕｎｉｖｅｒｓｉｙ，Ｇｕｉｔｙａｎｇ ５５０００３） ＡＢＳＴＲＡＣＴ Ｔｈｅ ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｓｏｆｔｗａ￣＇Ｃ ＡＮＳＹＳ ｉｓ ｕｓｅｄ ｆｏｒ ｃａｌｃｕｌａｔｉｎｇ ｔｈｅ ｓｔｒｅｓｓｅｓ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅｓ．Ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｔｅｎｓｉｌｅ ｓｔｒｅｎｇｔｈ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅｓ ａｎｄ ｔｈｅ ｎｏｔｃｈ ｓｉｚｅｆ ｔｈｅ ｄｅｐｔｈｅｓ ａｎｄ ｔｈｅ ａｒｃ ｒａｄｉｕｓ ｏｆ ｔｈｅ ｎｏｔｃｈ ｂｏｔｔｏｍ）ｉｓ ｓｔｕｄｉｅｄ．Ｔｈｅ ｒｅｓｌｔｕｓ ｏｆ ｉｆｎｉｔｅ ｅｌｅｍｅｎｔｌａ ａｎａｌｙｓｉｓ ｈａｖｅ ｂｅｅｎ ｃｏｍｐａｒｅｄ ｗｉｔｈ ｔｈｅ ｅｘｐｅｒｉｍｅｎ— ｔｌ ｒｅｓｕｈｓ．ａ ＫＥＹ ＷＯＲＤＳ ｎｏｔｃｈ ｓｅｎｓｉｔｉｖｉｙ；ｐｏｌｔｙｐｒｏｐｙｌｅｎｅ；ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ａｎａｌｙｓｉｓ 聚丙烯（ＰＰ）材料的缺口敏感性的实验研究，以 的关系，并将计算结果与实验结果作了比较。 试验的方式考察了聚丙烯（ＰＰ）材料带切口的试件 的拉伸强度和屈服应力与切口深度和切口底部圆弧 半径之间的关系。实验结果说明，聚丙烯（ＰＰ）材料 １有限元计算模型 由实际试验所用的试件的几何形状和承载的对 称性，取用试验试件的四分之一作为ＡＮＳＹＳ有限元 的缺口敏感度ＮＳＲ平均值为１．１ｌ，聚丙烯（ＰＰ）材 料对缺口是不敏感的…。但在目前的实验条件下， 还无法测出缺口处的应力分布，尚难于对聚丙烯 （ＰＰ）材料的缺口敏感性作进一步的深入分析。为 计算的模型（参见图２），其尺寸如下：长为３０ｍｍ，厚 为４ｍｍ，无缺口处的截面宽为５ｍｍ。 １．１ 分析设定 此，本文采用通用有限元软件ＡＮＳＹＳ对聚丙烯 （ＰＰ）材料的缺口试件的缺口处应力状态进行了分 析计算，考察了聚丙烯带切口试件的拉伸强度与切 口几何尺寸（切口深度和切口底部圆弧半径）之间 收稿日期：２００５—０６—１３ 当材料中的应力超过屈服点时，塑性即被激活， 即有塑性应变发生。本文分析的是实际试验中被拉 断的试件，其拉断的过程经历了弹性变形和塑性变 基金项目：贵州省科学技术基金资助项目（黔科台Ｊ宇［２００５］２０９４号）。 作者简介：周武雷。男，黔南民族师范学院讲师；蔡长安，男，贵州大学教授。 维普资讯 http://www.cqvip.com
４期 周武雷，等：聚丙烯（ＰＰ）材料缺口敏感性的有限元分析 １３ 形，必须作弹塑性分析来模拟试件在受拉情形下的 力学行为。 聚丙烯（ＰＰ）材料无切口光滑试件的拉伸强度 为１９．０Ｍｐａ，屈服应力为１８．５Ｍｐａ，其简化的本构关 系为理想塑性的应力一应变曲线，如图ｌ所示。 本分析设定为小变形和率不相关，不考虑应变 率的影响。 图１ 简化的聚丙烯（ＰＰ）材料的理想塑性应力一应变曲线 （Ｅ：１７１ＭＰａ。 ｓ：１８．５ ＭＰａ） Ｆｉｇ．１ Ｓｔｒｅｓｓ—ｓｔｒａｉｎ ｃｕｒｖｅ ｏｆ ｔｈｅ ｓｉｍｐｌｉｆｉｅｄ ｉｄｅａｌ ｐｌａｓｔｉｃｉｔｙ ｏｆ ｐｏｌｙｐｒｏｐｙｌｅｎｅ １．２单元选用 ＡＮＳＹＳ软件中适用于本研究的薄板平面应力 单元主要有ＰＬＡＮＥ２、ＰＬＡＮＥ４２、ＰＬＡＮＥ８２等单元。 本文选用ＰＬＡＮＥ８２单元，此单元具有一致位移形状 由于所取计算模型的曲线几何边界，采用映射 网格划分，技术方案如下：将原平板面划分成３个 面，以使每个面都是四边形，然后用ＰＬＡＮＥ８２单元 做映射网格划分，既可节省计算时间，又有足够高的 计算精度。 如对切口底部半径ｒ＝０．５ｍｍ、切口深度ｈ＝ 函数，能很好地适应曲线边界，具有分析物理非线性 问题包括弹塑性变形的能力。 １．３坐标的选取和网格划分 坐标的选取 对每一缺口试样建立ＡＮＳＹＳ计算模型时，坐标 系建立如图２所示，坐标原点始终取在缺口尖端点 上。 网格划分 图２ ＰＰ切口试件０．５Ｈ３的三面划分及 约束示意图 Ｆｉｇ．２ Ｓｃｈｅｍａｔｉｃ ｄｉａｇｒａｍ ｏｆ ｔｈｅ ｔｈｒｅｅ ｓｅｃｔｉｏｎｓ ｄｒａｗｉｎｇ ａｎｄ ｒｅｓｔｒｉｃｔｉｏｎｓ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５Ｈ３ 图３ ＰＰ切口试件０．５Ｈ３的映射 划分网格图 Ｆｉｇ．３ Ｍａｐｐｉｎｇ ｄｉｖｉｄｅｄ ｎｅｔｗｏｒｋ ｃｈａｒｔ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５Ｈ３ 维普资讯 http://www.cqvip.com
１４ 贵州科学 ２４卷 ３ｒａｍ（简写为０．５Ｈ３，以下皆采用此种简写方式表示 采用较小的计算代价来逼近得到每个计算模型的拉 切口的几何尺寸）的切口试件模型，取圆弧切口处 伸强度，载荷的施加采用写入载荷步文件的方法，大 的切点，另在下端面取点（１０，一２），上端面取点 致采用三到四个载荷步，使每一个子步都能得到收 （１２，３），三点间连线以直线分割原整面得新的３个 敛的解。 面，此三面均为四边形，如图２。再以本文所选取的 对缺口试件的ＡＮＳＹＳ计算模型作弹塑性分析 单元ＰＬＡＮＥ８２单元去划分网格，就可以很好的吻合 时，分为两个部分，各部分的计算设定如下： 曲线边界，然后调整划分单元的数目，就可以得到无 （１）弹性部分ｌｉｎｅａｒ—ｅｌａｓｔｉｃ：ｉｓｏｔｒｏｐｉｃ，输人参 任何形状警告的整体单元划分。这对尽量提高本文 数：ＥＸ以及泊松比。 的ＡＮＳＹＳ有限元分析模拟的精度有很大好处。边 （２）塑性部分：ｎｏｎｌｉｎｅａｒ—ｉｎｅｌａｓｔｉｃ—ｒａｔｅ ｉｎ． 界处的映射网格划分如图３所示，单元总数为 ｄｅｐｅｎｄｅｎｔ。＿＿——ｉｓｏｔｒｏｐｉｃ ｈａｒｄｅｎｉｎｇ ｐｌａｓｔｉｃｉｔｙ。＿＿——ｒａｉｓｅｓ １０８００个。 ｐｌａｓｔｉｃｉｔｙ，输入参数：屈服应力以及正切模量。 １．４收敛准则 采用ＡＮＳＹＳ软件自默认的收敛容限以力为基 ２有限元计算结果 础的收敛准则。 １．５加载方案 ２．１聚丙烯（ＰＰ）切口试件的拉伸强度计算值见 由于本文为物理非线性分析，为保证收敛以及 表１。 表１ 聚丙烯（ＰＰ）材料切口试件的拉伸强度计算值 Ｔａｂ．１ Ｃａｌｃｕｌａｔｅｄ ｔｅｎｓｉｌｅ ｓｔｒｅｎｇｔｈ ｖａｌｕｅｓ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅｓ 注：上表第一行０．１～１．０表示切口底部圆弧半径，分别为从０．１ｍｍ到１．０ｍｍ。Ｈ１、Ｈ２和Ｈ３表示切口深度分别取１、２ 和３ｍｍ。拉伸强度的单位为Ｍｐａ。 ２．２ 当深度ｈ分别为ｌ、２、３ｒａｍ时。切口试件拉伸强度的计算值随切口底部半径变化如图４所示。 聚丙烯（ＰＰ）材料的缺口敏感度ＮＳＲ的计算值见表２。 表２聚丙烯（ＰＰ）材料的缺口敏感度ＮＳＲ的计算值 Ｔａｂ．２ Ｃａｌｃｕｌａｔｅｄ ｎｏｔｃｈ ｓｅｎｓｉｔｉｖｉｔｙ ｒａｔｉｏ（ＮＳＲ）ｖａｌｕｅｓ ｏｆ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓｍａｐｌｅｓ 由表２，聚丙烯（Ｐ Ｐ）材料的缺口敏感度ＮＳＲ 的计算值的平均值是：当切口深度ｈ＝ｌｍｍ时，ＮＳＲ ３讨论 ＝１．０３；ｈ＝２ｒａｍ时，ＮＳＲ＝０．９４３；ｈ＝３ｒａｍ时， ＮＳＲ＝１．０９；总的平均值为１．０２。 ３．１聚丙烯（ＰＰ）切口试件有限元计算结果与实验 结果的比较 图４显示了聚丙烯（ＰＰ）材料切口试件拉伸强 维普资讯 http://www.cqvip.com
４期 周武雷，等：聚丙烯（ＰＰ）材料缺Ｅｌ敏感性的有限元分析 １５ 度的有限元计算值和实验值。由图４，在切口深度 为ｈ＝１ｍｍ和ｈ＝３ｍｍ时，有限元计算值和实验值 比较接近，但切口深度ｈ＝２ｒａｍ时，有限元计算值和 实验值有较大的误差。但有限元计算值和实验值随 切口底部半径变化的规律是基本一致的。 聚丙烯（Ｐ Ｐ）材料的缺Ｅｌ敏感度ＮＳＲ的计算值 图４—１切口深度ｈ＝１ｍｍ Ｆｉｇ．４～１ Ｎｏｔｃｈ ｄｅｐｔｈ ｈ＝ｌｍｍ 图４—３切口深度＾＝３ｍｍ Ｆｉｇ．４—３ Ｎｏｔｃｈ ｄｅｐｔｈ ｈ＝３ｍｍ ３．２全弹性与弹塑性状态下的缺口截面处的ＭＩＳ． ＥＳ等效应力值的讨论 以ＰＰ试件０．５Ｈ１为例。在约束和载荷施加相 同的情况下，考虑全弹性与弹塑性２种情况。全弹 性情况是指材料在加载过程中始终处于弹性状态， 不考虑屈服，对于聚丙烯（ＰＰ）材料这是一种为了比 较研究而做的假设状态。弹塑性状态是指随着载荷 的增加，处在弹性阶段的应力状态将不断变化，当应 的平均值是１．０２，而聚丙烯（Ｐ Ｐ）材料的缺口敏感 度ＮＳＲ的实验平均值是１．１ｌ，误差为８．１％，但它 们都大于１，用它们来衡量聚丙烯（Ｐ Ｐ）材料的缺Ｅｌ 敏感性的结论是一致的，即聚丙烯（Ｐ Ｐ）材料对缺 口是不敏感的。有限元计算结果支持实验关于聚丙 烯（ＰＰ）材料的缺口敏感性的相关结果。 图４—２切口深度ｈ＝２ｍｍ Ｆｉｇ．４—２ Ｎｏｔｃｈ ｄｅｐｔｈ ｈ＝２ｍｍ 图４聚丙烯（ＰＰ）材料切口试件的拉伸 度计算值随切口底部圆弧半径变化关系图 （有限元计算值■一■一■一■ 实验值●一●一●一●） Ｆｉｇ．４ ｃａｌｃｕｌａｔｅｄ ｔｅｎｓｉｌｅ ｓｔｒｅｎｇｔｈ ｖａｌｕｅｓ ｖｅｒＢｕ￥ｔｈｅ ｒａｄｉｕｓ ｏｆ ｔｈｅ ｂｏｔｔｏｍ ａｒｃ ｆｏｒ ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅｓ （ｃｌａｃｕｌａｔｅｄ ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ｖａｌｕｅｓ■一■一■一■ Ｅｘｐｅｒｉｍｅｎｔ ｖａｌｕｅｓ●一●一●一●） 力状态满足屈服条件时材料将屈服并产生塑性变 形，如采用ＭＩＳＥＳ屈服条件，则当ＭＩＳＥＳ等效应力 达到材料的屈服值 时，材料进人塑性状态。 全弹性情况下，缺口截面处的ＭＩＳＥＳ等效应力 值如图５所示，节点ＭＩＳＥＳ等效应力云图如图６所 示；弹塑性情况下，缺口截面处的ＭＩＳＥＳ等效应力 值如图７所示，节点的ＭＩＳＥＳ等效应力云图见图８ 所示。二者相比较，在约束和载荷施加相同的情况 维普资讯 http://www.cqvip.com
１６ 贵州科学 ２４卷 下，全弹性情况与弹塑性情况下的缺口截面处的 ＭＩＳＥＳ等效应力值的分布有很大差别。塑性变形使 应力在切口附近区域进行了重新分布，ＭＩＳＥＳ等效 应力的最大值在全弹性情况下为６４．５Ｍｐａ，而在弹 塑性状态下变成１８．５Ｍｐａ（等于聚丙烯材料的屈服 ６２．４８ ５７．７９ ５３．１０ ４８．４２ ４３．７３ ３９．０５ ３４．３６ 应力 ）。正是由于聚丙烯（ＰＰ）材料试件缺口截 面处的塑形变形，使聚丙烯（ＰＰ）材料缺口试件的 拉伸强度始终略高于无缺口光滑试件的拉伸强度， 从而聚丙烯（ＰＰ）材料的缺口敏感度ＮＳＲ大于１，得 到聚丙烯（Ｐ Ｐ）材料对缺口不敏感的结果。 ｌ ｛ ｛ 《 ｛ … ｝ ｅ： ； ； ｝ ； ２９．６８ ２４．９９ ２０＿３１ １５．６２ ； ’ ０～ ０ ； ； １．６ ４ １．２ ２ ＤｌｓＴ ２．４ ２．８ ３．２ ３．６ ４ 图５全弹性情况下缺口截面处的 ＭＩＳＥＳ等效应力值图（ＰＰ切口试件０．５ＨＩ） Ｆｉｇ．５ ＭＩＳＥＳ ｅｑｕｉｖａｌｅｎｔ ｓｔｒｅｓｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｎ ｔｈｅ ｎｏｔｃｈ ｃｒｏｓｓ ｓｅｃｔｉｏｎ ｕｎｄｅｒ ｅｌａｓｔｉｃ ｓｔａｔｅ ｏｎｌｙ 图６全弹性情况下节点的ＭＩＳＥＳ等效应力云图 （ＰＰ切口试件０．５ＨＩ） Ｆｉｇ．６ Ｃｌｏｕｄ ｃｈａｒｔ ｏｆ ｔｈｅ ｎｏｄｅ ＭＩＳＥＳ ｅｑｕｉｖａｌｅｎｔ ｓｔｒｅｓｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｕｎｄｅｒ ｅｌａｓｔｉｃ ｓｔａｔｅ ｏｎｌｙ （ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５Ｈ１） （ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５Ｈ１） １８．５０２ ’ １８．３５９ ｌ８．２２ １８．Ｏ８１ １７．９４２ １７．８０３ １７．６６４ １７．５２５ １７．３８６ １７．２４７ １７．１０８ ＼ ｛ 、７ ～ 、 、 ０．５ｈ１ｅｌ７ｌ ＮＯＤ＾Ｌ Ｓ０ＬＵＴＩＯＮ ＳＴＥＰ＝Ｉ ＳＵＢ＝７ｌ Ｔ工ＭＥ＝ｌ０ ＳＥＱＶ （ＡＶＧ） ＰｏｗｅｒＧｒａｐｈｉｃｓ ｌ ｉ ＥＦＡＣＩ ＝１ ＡＶＲＥＳ＝Ｍａｔ Ｄ皿＝３．０６３ ＳＭＮ＝．０６６７３３ Ｓ皿＝１８．５５ ＺＶ ｌ ｔＤＩＳＴ＝３．２５２ ‘ｘＰ ＝２．９９４ ＊ＹＦ＝一１．１８４ Ｚ一耵Ｆ职 ； ＿，．０６６７３３ １，ｉ ｆ ＼ ｛ ｛ 、 ～ ０ ４ ．８ １．２ １．６ ２ ２．４ ２．８ ３．２ ３．６ ４ 曾 ； 图８弹塑性情况下节点ＭＩＳＥＳ等效应力云图 （ＰＰ切口试件０．５Ｈ１） Ｆｉｇ．８ Ｃｌｏｕｄ ｃｈａｒｔ ｏｆ ｔｈｅ ｎｏｄｅｌ ＭＩＳＥＳ ｅｑｕｉｖａｌｅｎｔ ｓｔｒｅｓｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｕｎｄｅｒ ｅｌａｓｔｉｃ—ｐｌａｓｔｉｃ ｓｔａｔｅ 图７弹塑性情况下缺口截面处的ＭＩＳＥＳ等效应力值图 （ＰＰ切口试件０．５Ｈ１） Ｆｉｇ．７ ＭＩＳＥＳ ｅｑｕｉｖａｌｅｎｔ ｓｔｒｅｓｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｎ ｔｈｅ ｎｏｔｃｈ ｃｒｏｓｓ ｓｅｃｔｉｏｎ ｕｎｄｅｒ ｅｌａｓｔｉｃ—ｐｌａｓｔｉｃ ｓｔａｔｅ （ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５ＨＩ） （ｔｈｅ ＰＰ ｎｏｔｃｈ ｓａｍｐｌｅ ０．５ＨＩ） 维普资讯 http://www.cqvip.com
４期 周武雷，等：聚丙烯（ＰＰ）材料缺口敏感性的有限元分析 此种情况的原因尚不能确定。 ｌ７ ３．３ 关于有限元计算结果与实验结果的误差讨论 有限元分析结果与实验结果之间误差产生的原 因有： ３．３．１实验试件切口加工的误差引起的实验结果 参考文献 ［１］ 周武雷，蔡长安．热塑性材料的缺口敏感性实验研究 ［Ｊ］．贵州科学，２００５，２３（２）：２２—２６ 的误差：限于条件，本文实验试件的切口加工方法比 较粗糙，加工误差特别是切口几何尺寸的加工误差 难以控制。 ３．３．２材料简化成理想塑性模型而带来的误差：聚 丙烯（ＰＰ）材料在进行有限元计算时，其本构关系简 化为理想塑性材料的力学模型，这种简化与实际本 构关系存在着一定的差异。 ［２］王勖成，邵敏．有限单元法基本原理与数值方法［Ｍ］． 北京：清华大学出版社，１９８８ ［３］吴永礼编著．计算固体力学方法［Ｍ］．北京：科学出版 社，２００３ ［４］ 何其笙．高聚物的力学性能［Ｍ］．合肥：中国科学技术 ￣－ｖ－出版社，１９９６ ３．３．３有限元方法本身是一种求近似解的数值计 算方法，这种方法本身，从原理上和具体执行计算过 程中都必然存在有误差。 ３．３．４ 由图４—２可以看出在切口深度ｈ＝２ｍｍ ［５］谢贻权，何福保．弹性和塑性力学中的有限单元法 ［Ｍ］．北京：机械工业出版社，１９８１ ［６］ＡＮＳＹＳ结构分析指南（中）结构非线性 ［７］Ｌａｉ Ｊ，Ｖａｎ ｄｅｒ ｇｉｅｓｓｅｎ Ｅ．Ａ ｎｕｍｅｒｉｃａｌ ｓｔｕｄｙ ｏｆ ｃｒａｃｋ— ｔｉｐ ｐｌａｓｔｉｃｉｔｙ ｉｎ ｇｌａｓｓｙ ｐｏｌｙｍｅｒｓ［ｊ］．Ｍｅｃｈａｎｉｃｓ ｏｆ ｍａｔｅｒｉ－ ａｌｓ，１９９７，２５：１８３—１９７ 时，拉伸强度的有限元计算值与实验值之间的误差 比切口深度ｈ＝ｌｍｍ、ｈ＝３ｍｍ时的误差要大。出现 （上接第３３页） ［９］Ｃｈｅｎ Ｈ Ｂ，Ｚｈａｎｇ Ｃ Ｌ ａｎｄ ｗｅｎ ｘ Ｊ．Ｔｗｏ ｎｅｗ ｓｐｅｃｉｅｓ ｏｆ Ｓｉｍｕｌｉｕｍ（Ｇｏｍｐｈｏｓｔｉｂｉａ）ｆｒｏｍ Ｆａｎｊｉｎｇ Ｍｏｕｎｔａｉｎ，Ｃｈｉｎａ （双翅目：蚋科）［Ｊ］．昆虫学报，１９９６，３９（４）：４２３— ４２５ ［１３］Ｃｈｅｎ Ｈ，Ｃｈｅｎ Ｈ Ｂ．Ａ Ｎｅｗ Ｂｌａｃｋ＿ｆｌｙ Ｓｐｅｃｉｅｓ ｏｆ Ｓｉｍｕｌｉ－ ｕｍ（Ｄｉｐｔｅｒａ：Ｓｉｍｕｌｉｉｄａｅ）ｆｒｏｍ ＧｕｉｚｈＯＵ Ｐｒｏｖｉｎｃｅ，Ｃｈｉ－ （Ｄｉｐｔｅｒａ：Ｓｉｍｕｌｉｉｄａｅ）［Ｊ］．Ｅｎｔｏｍｏｌｏｇｉａ Ｓｉｎ，７（１）：２１— ２８ ｎａ［Ｊ］．Ｅｎｔｏｍｏｌｏｇｉａ Ｓｉｎ，２００１，８（３）：２０８—２１２ ［１４］ Ｃｈｅｎ Ｈ Ｂ，Ｚｈａｎｇ Ｃ Ｌ，Ｈｕａｎｇ Ｌ．Ｔｗｏ Ｎｅｗ Ｂｌａｃｋｆｌｙ Ｓｐｅｃｉｅｓ ｏｆ Ｓｉｍａｌｉｕｍ（Ｓｉｍｕｌｉｕｍ）Ｆｒｏｍ Ｓｉｃｈｕａｎ Ｐｒｏｖ－ ［１Ｏ］邓成玉，张有植，陈汉彬．西藏林芝真蚋三新种（双翅 目：蚋科）［Ｊ］．四川动物，１９９５，１４（１）：７—１Ｏ ［１１］ Ｃｈｅｎ Ｈ Ｂ，Ｚｈａｎｇ Ｃ Ｌ，Ｈｕａｎｇ Ｌ．Ａ Ｎｅｗ Ｓｐｅｃｉｅｓ ｏｆ Ｓｉｍｕｌｉｕｍ（Ｎｅｖｅｒｒｒｍｎｎａｉａ）Ｆｒｏｍ Ｓｉｃｈｕａｎ Ｐｒｏｖｉｎｃｅ， ｉｎｃｅ，Ｃｈｉｎａ（Ｄｉｐｔｅｒａ：Ｓｉｍｕｌｉｕｍ）［Ｊ］．Ａｃｔａ Ｚｏｏｔａｘｏ－ ｒｔｏｍｉｃｃｔ Ｓｉｎｉｃａ，２００５，３０（２）：４３０—４３５ Ｃｈｉｎａ（Ｄｉｐｔｅｒａ：Ｓｉｍｕｌｉｉｄａｅ）［Ｊ］Ａｃｔａ Ｚｏｏｔａｘｏｎｏｍｉａ Ｓｉｎｉｃａ，２００５，３０（３）：６２５—６２７ ［１５］ 陈汉彬．中国蚋类的区系分布和地理区划［Ｊ］．动物 分类学报，２００２，２７（３）：１８９—１９５ ［１２］邓成玉，张有植，薛群力，陈汉彬．西藏蚋科一新种
因篇幅问题不能全部显示，请点此查看更多更全内容