考虑拉压不对称的薄壁管数控热弯有限元建模关键技术

典型零件形体的热变形模型及其有限元仿真金少搏;赵凤霞;李纪峰【摘要】随着精密技术的精度不断提高,温度变化引起零件形变带来的影响不容忽视,为了解决此问题,根据弹性力学和热力学理论,建立了典型零件形体的热变形模型,利用ANSYS有限元分析软件对典型零件形体的热变形误差及其分布进行了仿真,并进行了模型实验验证,给出了实例应用的分析,说明了热变形模型在公差与配合中的应用.结果表明,所建立的模型考虑了热应力的影响,相对于传统计算方法更加合理可靠,为合理进行零件的公差与配合设计提供了理论依据和参考.【期刊名称】《机械设计与制造》【年(卷),期】2018(000)007【总页数】4页(P9-12)【关键词】热变形模型;ANASYS;有限元仿真;典型零件;圆柱孔轴【作者】金少搏;赵凤霞;李纪峰【作者单位】郑州大学机械工程学院,河南郑州 450001;郑州大学机械工程学院,河南郑州 450001;郑州大学机械工程学院,河南郑州 450001【正文语种】中文【中图分类】TH16;TH1241 引言随着精密技术的精度不断提高，温度变化带来的影响愈来愈大，甚至成为制约精度的决定因素，并使传统的热变形误差理论失去作用[1]。

近年来，文献[2-5]中对轴类、孔类零件在均匀温度场和非均匀温度场下的热变形规律进行了研究和探索。

文献[6]中以有限元为工具，研究了稳态场下圆柱筒尺寸变化规律；文献[7]中基于ANSYS对主轴热变形进行了建模与预测分析；文献[8]中对模块化工具系统中非线性温度下的主轴热变形位移量进行了数学建模，并提出了补偿模型。

但以上研究大都针对孔轴零件进行热变形研究，对球类、板类零件研究较少。

基于弹性力学和热力学建立了典型零件形体的热变形数学模型，基于ANSYS有限元分析软件对典型零件形体的热形误差及其分布进行了仿真，并进行了实验验证，最后进行了实例应用分析。

2 典型零件形体的热变形模型圆柱孔类、轴类、板类、球类零件是机械结构中的典型的基本零件形体。

芯棒伸出量对薄壁矩形管弯曲失稳起皱影响的数值模拟赵刚要刘郁丽杨合西北工业大学,西安,710072摘要:基于动态显式有限元平台ABA QUS/Ex plicit,建立了3A21铝合金薄壁矩形管绕弯成形过程的三维有限元模型,模拟分析了芯棒伸出量对薄壁矩形管绕弯成形过程中失稳起皱的影响规律。

结果表明:随着芯棒伸出量的增大,起皱波纹高度逐渐减小,切向压应力极大值波动的峰值减小,等效应力逐渐减小,并且变化趋于平缓;当芯棒伸出量达到一定值后,随着芯棒伸出量的增大,起皱波纹高度极大值变化不明显。

该研究对薄壁矩形管绕弯成形工艺参数的确定具有重要的参考价值。

关键词:薄壁矩形管;芯棒伸出量;失稳起皱;绕弯成形;数值模拟中图分类号:T G386;T P391 文章编号:1004)132X(2006)S2)0035)03Numerical Simulation of Influence of Mandrel on Wrinkling in Bending ofThin -walled Rectangular Tube Zhao Gangyao Liu Yuli Yang H eNo rthw estern Po lytechnical University,Xi .an,710072Abstract :A three-dimensio nal FEM mo del of ro tary draw bending of thin-w alled rectangular tube is built under the ABA QUS/Ex plicit env ir onment.Based on the model,a simulation and analysis of ro tary draw bending for aluminum alloy 3A21thin-w alled rectangular tube hav e been car ried out.The r esults sho w that the max imum heig ht of w rinkling w ave decreases obviously w ith the increase of the initial lead of the mandrel,the m ax imum tangential com press stress and the equivalent stress de -creases g radually,and the variation becom e steadier.T he maxim um height of w rinkling w ave does not almost alter w ith the increase of the initial lead of the mandrel w hen the initial lead reached a fix ed va-lue.T he achievem ents o f this study provide a g uideline for determ ining the process parameters fo r thin -w alled rectangular tube in rotary draw bending process.Key words :thin-w alled r ectangular tube;m andrel initial lead;w rinkling;rotary draw bending;numerical simulatio n收稿日期:2006)08)15基金项目:国家自然科学基金资助项目(50575184);西北工业大学研究生创新实验中心资助项目(05072)0 引言铝合金薄壁矩形波导管是雷达和无线电通信系统中的关键部件,被广泛应用于航空、航天等领域。

拉压不对称本构全文共四篇示例，供读者参考第一篇示例：拉压不对称是一种常见的现象，也是一个重要的材料本构特性。

在材料的拉伸和压缩过程中，往往会出现拉压不对称的情况，即在拉伸和压缩方向上的应力-应变曲线不对称。

这种不对称性不仅对材料的力学性能产生影响，还会对材料的工程应用带来一定的挑战。

拉压不对称是由于材料微观结构的非对称性导致的。

在金属材料中，晶粒的取向、晶粒内部的位错堆积等因素都会导致材料在不同方向上的变形行为不同。

在晶粒的特定取向下，材料在拉伸时会出现颈缩现象，而在压缩时则不会出现。

这种差异性导致了材料的应力-应变曲线在拉伸和压缩方向上不对称。

拉压不对称的本构特性在材料设计和加工中具有重要意义。

了解材料的拉压不对称特性可以帮助工程师更准确地预测材料在现实工程应用中的行为。

在设计零件的结构时，需要考虑材料在拉伸和压缩状态下的强度和变形性能，以确保零件在使用过程中不会发生失效。

拉压不对称还可以用于改变材料的力学性能。

通过调控材料的微观结构和添加合适的合金元素，可以实现材料在拉伸和压缩状态下的性能均衡，从而提高材料的综合性能。

要研究材料的拉压不对称本构特性，需要进行一系列的实验和模拟分析。

在实验方面，可以通过拉伸试验和压缩试验来获取材料在拉伸和压缩状态下的应力-应变曲线，进而分析材料的拉压不对称性。

在模拟方面，可以借助有限元方法和分子动力学模拟等技术来模拟材料的微观结构和变形行为，以揭示拉压不对称的机制及影响因素。

拉压不对称是材料本构中的重要特性，对材料的工程应用和性能具有重要影响。

通过深入研究拉压不对称的机制和调控方法，可以为材料设计和加工提供重要的指导，促进材料科学与工程领域的发展。

【本段是增加内容的一种方式，可以根据实际情况修改或者增加相关内容】在研究材料的拉压不对称本构特性时，需要重点关注材料的微观结构和变形行为。

根据材料的晶体结构和位错运动规律，可以理解材料在拉伸和压缩状态下的行为差异，从而揭示拉压不对称的机制。

6061铝合金超薄壁弯管制造工艺及组织性能研究丁月霞;马燕楠;郭群;王辉;马福业;郭训忠【摘要】通过单向拉伸试验建立了材料的本构方程,对6061铝合金板材半片拉深成形过程进行了有限元模拟,研究了摩擦系数与压边力对拉深成形结果的影响规律模拟结果表明:拉深成形时,应采用较大的压边力和较小的摩擦系数,以避免外管的外弧处易出现明显的起皱缺陷;基于模拟结果,进行了6061铝合金Φ202mm×1mm 弯管半片拉深成形;最后对半片焊接后的超薄壁弯管进行了耐压爆破试验及关键部位微观组织分析.【期刊名称】《航空制造技术》【年(卷),期】2016(000)017【总页数】5页(P50-54)【关键词】6061铝合金;超薄壁弯管;有限元模拟;拉深成形【作者】丁月霞;马燕楠;郭群;王辉;马福业;郭训忠【作者单位】江苏华阳金属管件有限公司,镇江212400;南京航空航天大学先进材料及成形技术研究所,南京211100;江苏核能装备材料工程实验室,南京211100;南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学先进材料及成形技术研究所,南京211100;江苏核能装备材料工程实验室,南京211100;南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学先进材料及成形技术研究所,南京211100;江苏核能装备材料工程实验室,南京211100;南京航空航天大学先进材料及成形技术研究所,南京211100;江苏核能装备材料工程实验室,南京211100;南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学先进材料及成形技术研究所,南京211100;江苏核能装备材料工程实验室,南京211100;南京航空航天大学材料科学与技术学院,南京211100【正文语种】中文导管是现代军机、民机上的关键部件之一，将各种介质输送到飞机发动机、起落架、座舱等关键部位，满足飞机包括燃油、氧气、液压、环境控制等方面的系统需求[1-3]。

薄壁球壳压缩非对称屈曲特性的实验及有限元分析周刚毅;董新龙;郝伟伟【摘要】The deformation behaviors of thin-walled spherical shells subject to quasi-state and dynamic com-pression are studied experimentally. The process of deformation, the non-axisymmetric buckling characteristic and its influence factors are analyzed in detail using ABAQUS finite element code. The results show that the deformation modes are dependent on the impact velocity. It is found that, if impact loading speed increases, the collapse modes shift from pentagons to hexagons. The dynamic loading-deformation curve on which the impact velocity is not very high is in good consistence with its quasi-static one. Meanwhile, the load-deformation curves of FEM simulation conform well to the experimental results at the stage of axi-symmetric inward dimple, but it is not so at the non-axisymmetric buckling stage. Furthermore, the process of deformation of spherical shells and its effect factors are discussed, suggesting that the contact constraints of surface between spherical shell and rigid plate play a significant role in the process of non-asymmetric deformation.%采用实验方法研究了球壳在刚性板准静态和冲击压缩下变形特性及非对称屈曲模态，结合ABAQUS有限元分析了球壳冲击压缩下的屈曲变形过程、非对称屈曲特性，探讨了其影响因素。

文章编号：2095-6835（2023）24-0006-05弯管成型截面畸变的有限元分析＊谌宏1，2（1.江苏科技大学苏州理工学院，江苏苏州215600；2.张家港江苏科技大学产业技术研究院，江苏苏州215600）摘要：针对弯管成型截面畸变的问题，基于ABAQUS有限元分析软件，建立了21-6-9高强度不锈钢管弯曲的有限元模型。

研究了相对弯曲半径、相对壁厚、弹性模量、屈服强度关于弯管成型截面畸变的显著性规律及经验公式。

研究结果表明，根据正交试验设计判断出，显著性强弱为相对弯曲半径＞相对壁厚＞屈服强度＞弹性模量；为降低弯管成型截面畸变率，可以选用相对弯曲半径较大的工艺组合；根据回归分析结果，得出成型参数关于弯管成型截面畸变率的经验公式，并校核验证了大概适用范围，该公式可以预测非大半径弯管成型截面畸变的情况，完善后可应用于实际生产。

关键词：管材弯曲；成型参数；截面畸变；有限元模拟中图分类号：TG386.3文献标志码：A DOI：10.15913/ki.kjycx.2023.24.002作为现代弹塑性加工工艺代表之一的管材加工工艺，其管材弯曲加工是其重要的组成部分，管材部件的轻量化、强韧化、高效、低耗等特点显著，被广泛应用于汽车制造、航空航天、输油管道等高新技术领域。

管材弯曲过程是一个非线性多重复杂的物理过程，在弯曲过程中容易出现成型截面畸变、起皱、壁厚减薄等各种质量缺陷，亦会发生回弹等多种问题。

因此，针对成型截面畸变问题开展几何非线性的模拟分析，了解其成型机理因素的显著性，判断最优弯曲方案，预测最大成型截面畸变，合理规避不合格的缺陷管材具有重大意义。

在管材弯曲成型系列研究中，国内外学者针对横截面变形现象开展了各种各样的研究。

王光祥等[1]通过实验的方法研究了弯曲中心角对截面畸变的影响，发现弯曲中心角是影响截面畸变的重要因素，椭圆率随弯曲中心角增大而增大，可以根据这个结论进行预测；JIANG等[2]主要研究了不同数控弯曲模组下的强度TA18管，其弯曲模、刮水模、夹紧模、压力模的合理选用可以提高截面质量；鄂大辛等[3]在平面应力和三向应力状态假设下，得出横截面短轴变化与壁厚的关系式；王刘安等[4]通过6061-T6铝合金管单向拉伸试验数据，对异形弯管进行有限元模拟，得出芯棒与管壁间隙大于1mm时，管材畸变减小，否则畸变严重的结论；何花卉等[5]在管材弯曲变形试验的基础上，进行有限元分析，指出长、短轴变化率比椭圆率更能形象反映界面形状变化，且短轴变化率更加明显，认为弯曲部分变形有向后段直管部分扩散的趋势；方军等[6]通过有限元软件建立不锈钢管材绕弯成型的弹塑性模型，分析了几何和材料参数对截面畸变的影响规律；宋飞飞等[7]利用有限元软件模拟Ti35合金管材绕弯过程，研究了芯棒伸出量、弯曲角度、压块相对助推速度、相对弯曲半径对它的影响规律；官强等[8]通过ABAQUS有限元软件模拟分析了圆管弯曲成型，提出将最大截面畸变率提高20%，应用实际加工判断截面质量的可行性；陈国清等[9]基于MSC.MARC有限元软件建立了推弯成型有限元模型，得出大弯曲半径推弯时，良好的润滑条件有利于获得更好的成型质量的结论；梁闯等[10]通过ABAQUS/Explicit平台，建立了TA18高强钛管数控弯曲成型过程三维有限元模型，研究得出较好的间隙水平是0.1mm的结论；刘芷丽等[11]基于ABAQUS有限元软件，建立圆管压扁-压弯连续成型的有限元模型，分析了圆管的受力方式；陈钱等[12]通过Dynaform有限元软件建立了高强度薄壁管材有限元模型，得出芯棒与管材间隙关于截面畸变率的影响规律；ZHAO等[13]通过ABAQUS/Explicit程序建立了钢管的三维有限元模型，通过实验验证了其可靠性，发现最大横截面畸变的位置几乎随模与管间隙的变化而变化；YAN等[14]基于有限元方法建立了一种起皱能量预测模型的成型极限搜索算法，并依次研究了芯棒球厚度等参数对管材起皱的影响；苏海波等[15]利用有限元方法对管材弯曲成型过程进行了数值模拟，得出了弯角外侧平均壁厚与相对弯曲半径间的关系。

GH625高温合金管缩径旋压成形数值模拟及试验研究黎波;袁其炜;靳凯;万柏方;王辉;郭训忠;陶杰【摘要】基于ABAQUS/Explicit平台建立了高温合金管缩径旋压三维弹塑性有限元模型,对GH625高温合金管缩径旋压成形过程进行了模拟分析,研究了轴向进给速度、旋轮圆角半径、旋轮安装角度对高温合金管缩径旋压成形质量的影响规律.研究结果表明,高温合金变径管在旋压过程中,随旋轮的轴向进给,不同的区域截面变化趋势不同.对工艺参数进行分析,发现轴向进给速度及旋轮圆角半径增大时,壁厚增大;旋轮安装角度对零件壁厚无明显影响;进给速度越大或圆角半径越大,截面椭圆度越大;距自由端较远时,旋轮安装角度越小,截面椭圆度越大;随着旋轮的进给,距自由端越近,椭圆度波动越严重.最后,根据数值模拟结果确定最终工艺参数为轴向进给速度0.25mrn/r,圆角半径取值2～4mm,旋轮安装角度为90°,并进行了高温合金管缩径旋压试验,得到了较好的成形件.%The elastic-plastic finite element model for the GH625 neck-spinning process was established by ABAQUS/Explicit software.The distribution of the stress,strain and the variation of the cross section of tube spinning process were analyzed.Furthermore,the effect of the spinning parameters on forming quality was obtained,such as axial feed rate,roller nose radius and roller installation angle.The spinning process results showed that the cross-sections in different zone of tube were generated different deformation.The results also reflected that the thickness increased with the increasing of the axial feed rate and roller nose radius.The roller installation angle had little influence on the tube thickness.The ellipticity became larger as the axial feed rate and roller nose radius increased.Located away from the free end,the ellipticity increased asthe roller installation angle decreased.The distribution of the ellipticity was uneven when the roller was near the free end of the tube.Finally,the process parameters was selected based on the simulation:0.25mm/r of axial feed rate,2-4mm of roller nose radius,90° of roller installation angle.The experiment was conducted and the accuracy quality tube was manufactured by neck-spinning experiment.【期刊名称】《航空制造技术》【年(卷),期】2017(000)018【总页数】7页(P36-42)【关键词】GH625高温合金;变径管;数值模拟;缩径旋压【作者】黎波;袁其炜;靳凯;万柏方;王辉;郭训忠;陶杰【作者单位】南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学机电学院,南京211100;江苏图南合金股份有限公司,镇江212352;南京航空航天大学机电学院,南京211100;南京航空航天大学材料科学与技术学院,南京211100;南京航空航天大学材料科学与技术学院,南京211100【正文语种】中文镍基高温合金因具有良好的抗氧化性、耐腐蚀性和优异的高温性能，广泛应用于航空、航天、核能、石化和海洋工程等重要领域[1-2]。

薄壁管拉扭变形应变路径控制试验研究吴建军　摘要控制模型的建立是一关键并在ＭＴＳ拉扭试验设备上进行了薄壁管拉扭试验基于数值模拟和神经网络的应变路径控制方法的正确性人工神经网络 应变路径 控制 薄壁管Ëæ×ÅÏÖ´ú¹¤ÒµµÄ·¢Õ¹°üÀ¨¸´ºÏ²ÄÁÏÄѳÉÐβÄÁÏÔÚº½¿ÕÓë´ËͬʱÕⶼҪÇóÉîÈëÑо¿²ÄÁϵıäÐιæÂÉÒÔ±ãÄÜÕýÈ·±æʶ²ÄÁÏÎïÐÔ²ÎÊýʹ²ÄÁÏÐÔÄܵõ½ºÏÀí´ïµ½Ô¤ÆڵijÉÐÎÄ¿±ê¿ØÖÆÄ£Ð͵Ľ¨Á¢Êǹؼü²»ÉÙÎÄÏײÉÓÃͳ¼Æ·ÖÎöÖеļ«´óËÆÈ»·¨»ò×îС¶þ³Ë·¨［由于塑性成形的高度几何及物理非线性此外难以在成形中在线决策有的研究虽然采用人工神经网络技术［］应变控制是按设定目标路径进行变形的有效措施因此提高成形件质量及精度有十分重要的意义４－７为了验证本文应变路径控制方法的准确性¹Ø¼üÊǵ±²ÄÁÏÎïÐÔ²ÎÊýµÄ±æʶÍê³Éºó±¾Îĸù¾Ý±¡±Ú¹Ü¸´ºÏÔغÉÏÂÓ¦±äÔöÁ¿µÄÊýÖµ½â·¨［６］并建立了相应的应变路径控制模型，如图具体控制方法为文中初始ｍ　＝　２单向拉伸试件过屈服点材料物性参数ｍ值和目标应变增量输入到训练好的神经网络识别出所需应力增量进行拉扭加载然后再根据识别出的ｍ值５９８０５０１５增量识别出新的应力增量便可达到对应变路径控制的目的试验设备如图该试验系统的测试结果十分精确２５ｋＮ５０ｍｍ或ＦｌｅｘＴｅｓｔＴＭ可方便准确地进行一次性破坏试验和千分之一牛顿范围内的疲劳试验１００Ｎ最大扭转角为２７０或者选择　５０　Ｎ最大扭转角为３６００应变仪所能测量扭转角的范围是图２　ＭＴＳ试验设备 图3　试验中安装完毕的试件图试件为ＬＹ１２薄壁管所示内径１０　ｍｍ退火炉中随炉降温至５０经单拉试验得出试件的力学性能如下０．２　＝　９０．４８　ＭＰａｂ　＝　２１３．４５　ＭＰａ厚向异性指数Ｒ＝０．７９８８强化指数ｎ＝０．３４２１按照目标应变路径使试件完成变形所示表示目标应变路径表示实际应变路径试验过程中实时识别所得材料物性参数ｍ值及其所对应的应力状态和加载产生的实际主应变值及目标主应变都在表　图4　目标应变路径与实际应变路径 图５　加载路径表１　１　１　在分析传统板料成形控制方法不足和局限性的基础上提出一种基于数值模拟和神经网络的应变路径控制方法±¾ÎÄËù²ÉÓõĻùÓÚÊýֵģÄâºÍÉñ¾-ÍøÂçµÄÓ¦±ä·¾¶¿ØÖÆ·½·¨¶Ô²ÄÁÏÒÀ¾ÝÉ趨µÄÄ¿±êÓ¦±ä·¾¶½øÐÐÓ¦±ä¿ØÖÆÊÇÓÐЧµÄ参 考 文 献１ 熊火轮．　计算机辅助塑性成形．　北京：　航空工业出版社，　１９９６　２ Ｇｏｌｄｅｎｔｈａｌ　Ｗ．　Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｎｅｕｒａｌ　ｎｅｔｗｏｒｋｓ　ｔｏ　ａｕｔｏｍａｔｉｃ　ｃｏｎｔｒｏｌ，　Ｐｒｏｃ．　ｏｆ　ｔｈｅ　ＡＩＡＡ　Ｃｏｎｆ．，１９９０：１１０８３1 ／　ＭＰａ１１２．　３ Ｋａｏ　Ｊ　Ｙ１１７　４ 楼顺天李顺平．　薄壁管复合载荷下应变增量的数值算法．机械科学与技术，２００１，ｌ２０（１）：３３李顺平．　金属板料应变路径控制方法研究．　塑性工程学报，２０００，７（４）：３３男西北工业大学机电工程学院教授主要从事塑性成形理论及计算机辅助塑性成形方面的科研工作　联系人陕西省西安市友谊西路１２７号 西北工业大学机电工程学院　联系电话　０２９－８４９３２３２薄壁管拉扭变形应变路径控制试验研究作者：吴建军作者单位：西北工业大学机电工程学院(西安)1.期刊论文吴建军.李祥基于数值模拟和神经网络的应变路径控制方法研究-西北工业大学学报2003,21(1)在金属塑性成形过程的应变路径控制中,控制模型的建立是一关键.文中所采用的应变路径控制模型是建立在材料物性参数m值实时辨识前提下以数值模拟和神经网络技术为基础的控制模型.根据上一步加载材料的实际应变增量、所加应力增量以及材料所在的应力状态识别出材料物性参数m值,再根据材料加载后的应力状态、目标应变增量以及识别所得m值,由训练好的识别应力增量的人工神经网络产生应加的载荷增量.在MTS试验设备上进行了薄壁管拉扭试验,通过试验过程中的实际应变路径与想要得到的目标应变路径的比较,验证了在正确识别材料物性参数前提下,基于数值模拟和神经网络的应变路径控制方法的正确性.2.期刊论文吴建军.李顺平.WU Jian-jun.LI Shun-Ping金属板料应变路径控制方法研究-塑性工程学报2000,7(4)本文将人工神经网络(ANN)用于应变路径的控制，旨在探索一种应变路径控制的新方法。

复合材料薄壁杆件非线性扭转分析
巫祖烈;赵奋;常国强
【期刊名称】《重庆交通大学学报（自然科学版）》
【年(卷),期】2007(026)006
【摘要】考虑几何非线性的影响,建立了在拉、压、弯、扭组合变形下复合材料薄壁杆件理论.重点对杆件进行非线性扭转分析.由给出的等效本构方程和总势能泛函,根据最小势能原理导出了所给问题的平衡控制方程和相应边界条件.讨论了该方程和边界条件的求解.
【总页数】6页(P1-6)
【作者】巫祖烈;赵奋;常国强
【作者单位】重庆交通大学,土木建筑工程学院,重庆,400074;重庆交通大学,桥梁结构工程交通行业重点实验室,重庆,400074;重庆交通大学,土木建筑工程学院,重庆,400074;重庆交通大学,土木建筑工程学院,重庆,400074
【正文语种】中文
【中图分类】TU338
【相关文献】
1.开口薄壁杆件扭转分析的一维离散有限元法 [J], 任伟;杨绿峰;陈建芳;王丽敏;乔永平
2.闭口薄壁杆件的约束扭转剪应力分析 [J], 张元海
3.复合材料薄壁杆件的力学分析：（二）自由扭转问题 [J], 陈烈民
4.考虑次剪应力的薄壁杆件扭转刚度矩阵和它的有限元分析 [J], 郭乙木;鲁祖统
5.开口薄壁杆件约束扭转分析的广义参数有限元 [J], 杨绿峰;钟德文;李桂青
因版权原因，仅展示原文概要，查看原文内容请购买。