脉冲电流作用下电子风应力模型与计算示例
脉动风场的模拟方法及其在输电线路风振计算中的应用
脉动风场的模拟方法及其在输电线路风振计算中的应用沈国辉;黄俏俏;郭勇;邢月龙;楼文娟;孙炳楠【摘要】针对脉动风场各种模拟方法的适用性问题,在相同的输电塔线实例上进行模拟方法的应用,对比各种方法的计算效率和计算结果,分析各种模拟方法的适用性和等价性,最后探讨单塔和塔线体系计算结果的差异.研究表明:基于POD分解的WAWS法不存在风速互谱密度矩阵无法分解的情况,计算效率较高,因此推荐使用.考虑三维风场后响应的脉动均方根比一维风场大;塔线体系中计算得到的响应均方根比单塔大.%Different wind field simulation methods are applied to the same transmission line to illustrate whether these methods are suitable to simulate fluctuating wind field.The calculating efficiency and results based on these methods are compared and the applicability of these methods is studied.The difference between the results based on the model of an isolated tower and transmission-line system is also investigated.Results from this study show that, the WAWS method with POD decomposition is highly recommended because this method is quite efficient and does not encounter the problem that the cross spectrum matrix of wind velocities cannot be decomposed.The root-mean-square responses of the transmission tower considering three dimensional wind fields is larger than those considering one dimensional wind field.The root-mean-square responses of the tower line system are also larger than those of the isolated tower.【期刊名称】《空气动力学学报》【年(卷),期】2013(031)001【总页数】6页(P69-74)【关键词】风场;模拟;输电塔;风振响应;时域分析【作者】沈国辉;黄俏俏;郭勇;邢月龙;楼文娟;孙炳楠【作者单位】浙江大学土木工程学系,浙江杭州 310058;浙江大学土木工程学系,浙江杭州 310058;浙江省电力设计院,浙江杭州 310007;浙江省电力设计院,浙江杭州310007;浙江大学土木工程学系,浙江杭州 310058;浙江大学土木工程学系,浙江杭州 310058【正文语种】中文【中图分类】TU3120 引言结构风振时程分析需要有作为输入的脉动风荷载,由于输电塔是格构式透空结构,不能像房屋结构那样采用测压试验获得脉动风荷载,因此需要采用人工模拟方法来获得。
脉冲强磁体中应力的有限元分析
根据磁体几何结构,建立好有限元分析模型。 线圈是在室温时绕制的,假定温度是 300 K,初始 预应力分析的结果如图 5(a)所示。从图中可以看 出,由于环向预应力为拉应力,导致径向方向层间
表2 Table 2 层 Layer 1,3,5,7,9,11, 13,15,17,19 2 4 6 8 10 12 14 16 18 20 21 厚度 Thickness/mm 2.2 1.4 2 2.7 3.1 2.1 1.6 1.3 1.2 1 7.5 20
图2 Fig. 2
材料应力-应变曲线
Stress-strain curve of material
分析时,先给出计算的初始条件,读入材料的 属性,计算材料的应力大小,并判断材料是否进入 塑性状态。 如果没有, 则仍采用弹性力学方程计算, σ =Eε,E 是杨氏模量;如果进入塑性状态,则根 据应力-应变曲线采用查表法查找对应的应力-应变 关系,然后进行塑性应变计算。 2
力分析,如图 1 所示。假设该单元的半径为 r,电 流密度为 j, 磁场强度为 H。 该单元受到 3 个力的作 、两边单元对它的环向作用 用,分别是洛伦兹力 fB 力 ft、上下相邻单元对它的径向作用力 fp。在分析 中假定拉应力为正,压应力为负。从图 1 中可以很 [1] 简单地写出这 3 个力的表达式 : f B = HIL = Hr ∆θ j ∆r ∆z = Hjr ∆r ∆z (1)
——————————————
第一作者:彭 涛,男,1977 年出生,2000 年毕业于华中科技大学电力系,现为该校电机与电器专业博士研究生,从事脉冲强磁场设计技术研究, E_mail: peng_wt@ 收稿日期:2003-06-16,修回日期:2003-09-10
第9期
强流脉冲电子束轰击作用下的扩散模型及其数值计算_邹建新
第27卷 第9期 核 技 术 V ol. 27, No.9 2004年9月 NUCLEAR TECHNIQUES September 2004——————————————第一作者:邹建新,男,1978年出生,2003年于大连理工大学获硕士学位,现为该校材料系在读博士研究生,主要从事强流脉冲电子束表面改性方面的研究,E_mail: zoujx@ 收稿日期:2003-03-24,修回日期:2003-06-17强流脉冲电子束轰击作用下的扩散模型及其数值计算邹建新 吴爱民 秦 颖 郝胜智 宋丽丽 王晓钢 董 闯(大连理工大学三束材料表面改性国家重点实验室 大连 116024)摘要 在温度场和应力场计算的基础上建立了强流脉冲电子束轰击作用下的扩散模型,并给出了数值方法及其数值解。
该模型与方法同样适用于其它高能束流作用下的扩散过程。
计算表明,浓度扩散流仍然是影响扩散的主要因素;而轰击超过一定次数后,扩散的作用将减弱;当边界条件为表面扩散时,扩散进行较快,这是表面涂覆加脉冲电子束后处理快速表面合金化工艺的理论基础。
对实验结果和理论结果的对比分析表明,在脉冲轰击下,扩散激活能随空位浓度的增加而下降,从而加速扩散过程;在表面有熔化的情况下,则液态时的对流混合作用是主导因素。
关键词 强流脉冲电子束,扩散,数值模拟 中图分类号 O343.6,O482.2在过去的十多年中,脉冲的高能束流如离子束、电子束、激光束等对材料表面的改性作用得到了广泛的研究[1—11]。
它已经成为表面改性技术中非常重要的手段。
脉冲束流具有很高的能量密度,可以在很短时间内在材料表面沉积很高的能量,形成熔化、汽化、应力波、冲击波等物理现象,从而对材料表面造成特殊的改性作用,获得用常规方法难以实现的物理化学性能,例如硬度的提高,抗摩擦、抗疲劳、耐腐蚀性能的提高等。
目前,高能束已广泛应用于材料的表面强化、合金化、清洗及表面熔覆、快速淬火等方面,应用前景十分广阔。
计算脉动风作用下拉索反应的一种新方法
动风压可以表示为 : z) W(t ,=
z z t ) ) ) d
( 1 )
式 中 。 是索 的阻尼 系数 或称为体形系数 ; , P 是平均 风速高度变化 z 系数 , 它与地面粗糙度 有关 , 按我国 A、 、 、 B c D四类地 貌分别确定 则为
139z1 ) . ( 0 7 /
=
1 0 4x1 .63 0
斜拉索 的应力为 : _i5 0 621 2 2 OS =T=6丽 — 8 7 0= 4 7 ×l pa 3
.
。
斜拉索 的应变为 :£ — = .4 8 1 = 58 3  ̄ 0 。
Ee q
50 斜拉索的密度为 :p =9O 8 kg = M一 一 . 9
01。桥位处的基本风压为 wa O 8N m . 6 o . k / 。索为圆截面构件 , =4 按照基 本风压和索的直径 ,0 0 < a Z . x . 8 0 0 7 0 0 5 N 6 . 2 w  ̄ = 4 01 0 体型系数经插值得到 :
l = .+O0 5 006 )05(.1- . 2: .12 a 07 f.1- . 7x ./ 0 5 00 )1 9 , s 0 0 0 0
采用式( ) 2和式( ) 照 B 3, 按 类地貌计算风压高度变化系数 №( z )和
山东滨州黄河公路大桥座落 于山东省滨州 市内 , 三塔四跨钢筋混 是 凝 土斜拉桥 , 跨度布设 为 8 m+ 0 m+ 0 m+ 4 该 桥所处地 理位 置平 4 3 0 3 0 8 m, 坦开阔 , 斜拉索在风荷载的作用 下极易引起振动 , 因此 , 认识斜 拉索在脉 动风作用下的动力特性 非常重要 。 本文 以滨州黄河大桥最长 的一根斜拉 索 N6 2 为例 , 采用大型通用软件 A s 对其进行脉动风荷 载作 用下的计 ns y 算分析 。提出了基于一点加速度 的时程分析方法 , 来模 拟风荷 载对斜拉 索 的作用 , 取图 l 所示 脉动风速时程前 2 0 的风荷载计 算斜 拉索上的 0秒 风压值 。
应用脉冲电流在冲击冲压过程中的应力分析
应用脉冲电流在冲击冲压过程中的应力分析引言:冲击冲压加工是一种常见的金属加工方法,用于将金属板材等材料制成所需的零件。
在冲击冲压过程中,加工件所受的载荷和应力分布对其质量和性能起着重要影响。
为了提高冲击冲压过程中的加工效率和质量,研究人员一直努力寻找新的改进方法。
其中,应用脉冲电流在冲击冲压过程中实现应力分析的方法备受关注。
本文将探讨应用脉冲电流进行应力分析的原理、优势以及实际应用。
一、脉冲电流在冲击冲压过程中的作用机制脉冲电流通过电导体传递,能够在瞬间产生高温、高压的电磁场,并对材料的晶体结构和力学性能产生影响,从而改变材料的力学性能和变形行为。
在冲击冲压过程中,脉冲电流的作用机制主要包括以下几个方面:1.应力改善:脉冲电流在传导过程中发生电阻,产生热量。
这些热量会促使材料中的晶体结构发生变化,从而使其组织得以改善。
这种改善可以使得冲击冲压过程中的应力集中减少,提高零件的强度和韧性。
2.表面硬化:脉冲电流的高温和高压作用下,材料表面会发生相变和冷却过程,形成一层硬质薄膜。
这层硬膜能够提高材料的耐磨性和硬度,在冲击冲压过程中减少表面磨损和变形。
3.晶粒细化:脉冲电流能够使材料中的晶粒发生再结晶或细化。
晶粒的细小可以提高材料的塑性变形能力,从而减轻冲击冲压过程中的应力集中和变形。
二、应用脉冲电流进行应力分析的优势相比传统的应力分析方法,应用脉冲电流进行应力分析具有以下几个优势:1.无损检测:应用脉冲电流进行应力分析无需对样品进行破坏性测试,可以实现无损检测。
这对于一些珍贵的或难以取样的材料至关重要。
同时,无损检测不会对待测对象产生任何影响,能够准确且高效地进行应力分析。
2.快速响应:脉冲电流的传导速度很快,能够在短时间内完成应力分析过程。
这使得应用脉冲电流进行应力分析在实际工程应用中具有很高的实时性。
例如,在钢铁冶金行业中,可以实时监测钢铁材料在冲击冲压过程中的应力分布,对产品的质量进行实时跟踪和调整。
风力机翼型的气动模型及数值计算
文章编号:167325196(2010)0320065204风力机翼型的气动模型及数值计算李仁年,李银然,王秀勇,绕帅波(兰州理工大学能源与动力工程学院,甘肃兰州 730050)摘要:考虑到不同湍流模型和边界层网格对风力机翼型气动性能有着不同的影响,采用4种边界层网格和4种湍流模型,对DU932W2210翼型的气动性能进行数值计算,将计算结果与实验值进行比较.研究结果表明:在合适的边界层网格下,DES模型的计算结果最接近实验值,而且该模型对翼型尾流中的旋涡有很强的捕捉能力.关键词:风力机;翼型;湍流模型;边界层;网格划分中图分类号:T K83 文献标识码:AAerodynamic model of airfoil for wind turbine and its numeric computation L I Ren2nian,L I Y in2ran,WAN G Xiu2yong,RAO Shuai2bo(College of Energy and Power Engineering,Lanzhou Univ.of Tech.,Lanzhou 730050,China)Abstract:Taking into account of t he effect of different t urbulence modes,and mesh division in boundary layer on t he aerodynamic characteristics of t he airfoils for wind t urbine,4boundary layers,4modes were employed for numeric evaluation of aerodynamic characteristics of wind t urbine airfoils DU932W2210.The calculation result s were compared wit h experimental ones.It was shown by t he investigation result t hat t he result of calculation wit h D ES mode was t he clo sest to t he experimental one for an approp riate bounda2 ry layer grid,and it has a st rong ability to capt ure t he vortex in t he wake behind t he airfoil.K ey w ords:wind t urbine;airfoil;t urbulence mode;boundary layer;mesh division 由于风力机叶片前缘半径较大,叶片表面边界层容易发生分离,分离会形成旋涡,而旋涡的运动、发展和破裂反过来又影响着分离流场[1].因此深刻认识叶片边界层的流态与准确计算边界层的分离,对于正确预估叶片升阻力、控制并减小流动分离以及叶片的优化设计有着重要的意义.目前,对湍流的数值模拟分为直接数值模拟(DNS)、大涡数值模拟(L ES)和雷诺时均方法(RANS)3类.其中DNS从流体控制方程出发,可以模拟湍流流场中各种尺度的脉动,但受计算机条件所限,目前只用于研究低雷诺数简单湍流物理机制.L ES将湍流流场中大尺度脉动用数值模拟方法计算,小尺度脉动对大尺度运动的作用做亚格子模型假设,在以下方面具有其他模型无可比拟的优势: 1)从层流到湍流转捩的预测;2)非定常湍流的预测;3)高速湍流的预测[2].但实际的工程问题往往 收稿日期:2009211216 基金项目:国家重点基础研究发展973计划项目(2007C B714600) 作者简介:李仁年(19632),男,甘肃民勤人,教授,博导.具有很高的雷诺数和很薄的边界层,边界层内小涡的尺度往往比边界层的厚度小很多,这使得要完全采用L ES模拟薄边界层内的流动仍然需要很大的计算机资源,到现在为止还是不太现实[3].RANS在工程中应用最为广泛,它完全采用湍流模型模拟湍流流动,只给出湍流流场的统计平均量,可以有效地模拟附体边界层流动,但对短暂的旋涡脱落和失速后的流场难以模拟[4].而脱体涡模拟(D ES)方法通过结合RANS和L ES各自的优点,可以比较快速而有效地模拟工程应用中常见的非定常的流动特征和边界层的分离运动[426].模拟气流分离的关键是能够准确模拟边界层分离,这需要布置合理的边界层网格.理论上壁面底层网格尺寸越小,网格的渐变率越接近于1;网格数越多,计算精度越高[7],但覃文洁等人提出网格的渐变率与网格数对计算精度的影响有限,它们的影响是基于底层网格尺寸的[8].齐学义等人提出采用结构和非结构化网格相结合的划分方式,可以提高网格质量和计算精度[9].本文选用Fluent软件,研究适第36卷第3期2010年6月兰 州 理 工 大 学 学 报Journal of Lanzhou University of TechnologyVol.36No.3J un.2010合风力机专用翼型的边界层网格与湍流模型.1 计算模型1.1 控制方程与拓扑结构选取不可压缩的雷诺时均方程为主控方程,不考虑体积力和外部热源.考虑到DU932W 2210翼型是为了克服气流流过相对厚度较大NACA 翼型过早的发生分离,导致翼型气动性能严重下降而设计的[10],而且该翼型几何形状简单,生成网格质量较好;模型计算量小,适于进行大量的数值计算,可以对网格分布、湍流模型的不同组合进行分析比较;国外已公布较全的实验数据,这些数据都是在弦长为0.6m 时得到的.为了便于比较,本文取弦长为0.6m 的翼型为研究对象.建立长度为45倍翼型弦长、宽度为40倍的翼型弦长的二维计算区域,如图1所示,把该计算域沿翼展方向拉伸1倍翼型弦长就可得到三维计算域.图1 二维拓扑结构Fig.1 Tw o 2dimensional topological structure1.2 网格划分因为在同一算法下均匀分布的正交计算网格可以获得最高的计算精度,所以本文利用CAD 的表面构造技术以及多块网格技术生成了高质量、完全结构化的网格.该方法通过非均匀有理B 样条插值(NU RBS )将物理域映射到贴体坐标系下的求解域,进行流程计算域多块网格的构造与重构,最后生成的网格为贴体的、正交性很好的网格.由于翼型附近的流场参数变化梯度比远场的参数变化梯度大得多,且翼型前后缘的流动情况对翼型扰流数值模拟的影响很大,因此对翼型附近的网格进行了局部加密,图2为翼型附近的网格.为比较翼型附近网格分布对边界层计算的影响,保持翼型表面周向网格节点不变,改变边界层内节点的法向分布以及第一层网格的高度,从而改变网格的纵横比,以确定适合于该翼型的边界层网格,网格划分方式见表1.流场方向半圆弧bcd 上布置330个节点,直线ab 、f g 、ed 上各布置80个节点.在三维计算域中,翼展方向上布置60个节点.图2 翼型计算网格Fig.2 Three 2D grid for airfoil表1 边界层网格划分策略T ab.1 Method of bound ary layer mesh division 网格划分策略网格层数第1层高度/mm1234101522301.000.500.100.051.3 边界条件与离散格式进口abcde 给定为速度进口,来流的湍流度为1%,湍流扩散长度为0.01m.出口af e 为压力出口,表压力给定为0,湍流度和湍流扩散长度与进口一样.翼型表面gm hn g 满足壁面无滑移条件.除在DES 和L ES 模型中对动量方程的离散采用默认离散格式(bounded cent ral differencing )外,其他模型中对连续方程、动量方程、雷诺时均方程等方程都用二阶迎风格式来离散,压力速度的耦合采用SIM 2PL EC 算法.2 结果分析速度由雷诺数或马赫数来确定,雷诺数为3.0×106、马赫数为0.22,弦长为0.60m.为了跟实验数据做对比,用Re =ρυc/μ或M a =υ/a 求得进口速度为76.56m/s.假设流动非定常,设定时间步长为0.001s ,在每个时间步长内迭代20次,利用升力系数、阻力系数来监测解的收敛性,当升阻力系数稳定时认为计算收敛.2.1 边界层的比较计算以三维直叶片为研究对象,研究不同边界层网格对翼型气动性能的影响.由图3a 可以看出,当攻角α<7°时,不同边界层网格计算的升力系数无大差异,且与实验值相当吻合,这说明附体流动对边界层网格的要求较低.而当攻角α>7°时,第1种网格划分策略计算所得的升力系数的最大值相对最小,且远小于实验值,失速提前发生,而第2种划分策略计算所得的最大升力系数相对最大,且大于实验的最・66・ 兰州理工大学学报 第36卷大值,而且对应攻角也偏大,失速滞后发生,但它对失速后的流动有较高的模拟精度.这两种策略对失速攻角的计算有较大的误差.第3种划分策略有最高的计算精度,第4种划分策略对失速攻角的模拟精度高于第2种,但对失速后的模拟计算精度跟第2种相当.由以上分析可知,当流体处于附体状态时,或翼型为小攻角时,翼型尾缘处逆压梯度很小,模拟的准确性对边界层网格的要求较低;随着攻角的增大,翼型尾缘处的逆压梯度增大,边界层发生分离,此时尾流中所划的边界层网格仅占尾流区域的很小一部分,它对尾流流场的捕捉能力减弱,而且旋涡的随机性也增大了模拟误差,这就导致不同的边界层网格的计算结果有很大差异.由图3b 分析可知,当攻角α<7°时,第1种划分策略对阻力系数的计算值小于其他划分策略的计算值,且远小于实验值,这是因为边界层内分布的网格点数过少,或没有网格节点,引起严重的数值耗散,掩盖了真实的物理现象以至于不能正确的描述边界层内的流动.随着底层网格尺寸的逐渐减小,计算所得的阻力系数与实验值的偏差也逐渐减小,但小到一定程度后偏差反而增大,这与文献[8]的观点相悖.这是因为第1层网格的纵横比过大所产生的数值 (a )升力系数随攻角的变化曲线 (b )阻力系数随时间的变化曲线图3 翼型升、阻力系数随攻角的变化曲线Fig.3 V ariation of lift coeff icient and drag coeff icient vsattack angle刚性影响了解的精度.当攻角α>7°时,各种边界层网格计算所得的阻力系数有很大差异,且与实验值也有较大偏差,而此时压差阻力起主要作用,第3种边界层划分策略的计算结果最接近实验值,即第3种网格划分策略对压差阻力的计算精度较高.这说明翼型失速后,边界层发生分离,最底层网格尺寸对压差阻力的计算有较大的影响.2.2 湍流模型的比较计算在三维拓扑结构的基础上,研究大涡模拟L ES 、雷诺平均方法RANS 和脱体涡模拟DES 三种方法在风力机专用翼型数值计算中的应用.由图4a 可以看出,当攻角α<8°时,所有模型对升力系数的计算结果无大差异且都与实验值很好的吻合,这几种湍流模型都能很好地计算出翼型失速前的流动特征,对附体流动有很高的计算精度;当攻角α>8°时,各种模型计算出的升力系数的变化趋势大概相同且都与实验值的变化趋势相同,但SST DES 和RSM 计算出的最大升力系数略大且大于实验值,SST k 2ω和L ES 模型计算出升力系数的略小于实验值,但在相同的攻角下L ES 模型的升力系数计算误差却小于SST k 2ω对升力系数的计算误差,说明L ES 的计算精度高于SST k 2ω,而DES 模型的计算结果最接近实验值,这说明DES模型对翼型分离流动或者 (a )升力系数随攻角的变化曲线 (b )阻力系数随时间的变化曲线图4 不同湍流模型下翼型升、阻力系数随攻角的变化曲线Fig.4 V ariation of lift coeff icient and drag coeff icient vsattack angle with different turbulence modes・76・第3期 李仁年等:风力机翼型的气动模型及数值计算 说脱体流动有很高的计算精度.由图4b 可知,当攻角α<8°时,所有模型计算的阻力系数都相对偏大且都大于实验值,但变化趋势相互吻合且与实验值吻合,其中L ES 的计算结果最大且远大于实验值,这表明亚格子模型放大了翼型附近的小尺度脉动对它的影响,此时的流动处于附体状态,翼型仅受到摩擦阻力的作用,这说明所有模型对摩擦阻力的模拟计算值均过大.当攻角α>8°时只有L ES 计算出的阻力系数大于实验值,而其他模型对阻力系数的计算结果都小于实验值,说明RANS 和DES 对压差阻力的计算能力较小,L ES 对气流分离造成的压差阻力有较高的计算精度.数值计算中,对流场的准确计算是翼型气动性能计算的基础,流场的形状及其变化规律反映了气流的客观流动规律,由图5可以看出,RANS 方法得到的旋涡结构单一而平滑,没有捕捉到旋涡的脱落,而DES 方法得到的绕流不仅在流向有旋涡的卷起和脱落,而且在展向还有大尺度的脉动,这是因为DES 在分离区域d 3=C DES Δ,湍流模拟不再依赖当地网格单元中心到翼型壁面的最短距离d ,即湍流的模拟与物体几何外形没有直接的关系,而与当地网格本身的尺度Δ直接相关,因此DES方法在流图5 旋涡等值面图Fig.5 Isosurfaces of vorticity向和展向都能计算出旋涡的运动,从而能得到更复杂的旋涡结构,所以DES 方法在数值模拟非定常大尺度分离流动方面具有明显的优势,它可以比RANS 方法更真实地模拟出高雷诺数下分离旋涡破裂后的非定常流动特征.3 结论1)翼型气动性能的计算精度跟壁面网格分布的情况有关,从应用的角度出发,壁面法向网格布置应该适宜.如果近比网格太稀,边界层内网格节点数不够,就无法准确地捕捉黏性效应;另外,也不能盲目增加网格节点数,还应兼顾网格的纵横比.2)翼型气动性能的计算精度与湍流模型有关,不同的湍流模型针对特定的物理模型才有较理想的计算结果,在计算三维直叶片时,D ES 模型能够捕捉到分离旋涡的非定常特征,得到更为真实的流场.参考文献:[1] GU IL MIN EAU E ,PIQU ET J ,QU EU TR Y P.Two 2dimension 2al turbulent viscous flow simulation past airfoils at fixed inci 2dence [J ].Computers &Fluides ,1997,26(2):1352162.[2] 王汉青,王志勇,寇广孝.大涡模拟理论进展及其在工程中的应用[J ].流体机械,2004,32(7):23227.[3] 邓 枫,伍贻兆,刘学强.用DES 数值模拟分离扰流中的旋涡运动[J ].计算物理,2008,25(6):6832688.[4] STEFAN S ,FRAN K T.Detached eddy simulation of flow a 2round a 2airfoil [J ].Flow ,Turbulence and Combustion ,2003,71(1/2/3/4):2612278.[5] 李 栋,焦予秦,IGOR M ,等.Detached 2Eddy Simulation 方法模拟不同类型翼型的失速特性[J ].航空学报,2005,26(4):4062410.[6] STREL ETS M.Detached eddy simulation of massively separa 2ted flows [R ].Reno :AIAA ,2001.[7] ZIN GG D parison of several spatial discretizations fort he N 2S equations [R ].Reno :AIAA ,1999.[8] 覃文洁,胡春光,郭良平,等.近壁面网格尺寸对湍流计算的影响[J ].北京理工大学报,2006,26(5):3882392.[9] 齐学义,冯俊豪,李纯良,等.三维湍流流动计算在混流式转轮水力设计中的应用[J ].兰州理工大学学报,2006,32(5):48252.[10] TIMMER W A ,VAN ROOI J R P J O M.Summary of t hedelft university wind turbine dedicated airfoils [J ].Journal ofsolar energy engineering ,2003,125(4):4882497.・86・ 兰州理工大学学报 第36卷。
NO脉冲流光放电激发解离截面的理论计算
自 然科 学版 )
J o u r n a l 0
( Na t u r a l S c i e n c e Ed i t i o n )
2 01 3
第 3 3卷
第 1 期
Vo 1 . 3 3 No . 1
D On . 1 0 0 0—1 5 6 5 . 2 0 1 3 . 0 1 . 0 0 5
验证 了理论 处理 的合 理性 , 也 为 NO 分 解 过 程 研 究提 供 了理 论 依 据 .
关键 词 : 碰撞 激发 解 离截 面 ; 脉 冲 流光放 电 ; 电 子激发 温度 中图分 类号 : O5 6 1 . 3 文献 标 志码 : A 文章编 号 : 1 0 0 0—1 5 6 5 ( 2 0 1 3 ) 0 1— 0 0 1 9— 0 5
t e mp e r a t u r e a n d c o l l i s i o n e x c i t a t i o n f u n c t i o n o f e x c i t e d e l e c t r o n i c s t a t e ,t h e i n t e g r a l o f c o l l i s i o n e x c i t a t i o n s e c t i o n s o f t wo ma i n e x c i t e d a n d d i s s o c i a t i n g p r o c e s s e s o f NO i n p u l s e s t r e a me r c o r o n a we r e c a l c u l a t e d !a n d t h e c o mp e t i t i o n o f t wo ma i n e x c i t e d a n d d i s s o c i a t i n g p r o c e s s e s o f NO wa s s t u d i e d a c c o r d i n g t o n o r ma l i z a —
脉动风荷载时程数值模拟研究
率谱和相关函数这两个 函数来描述脉动风。 1 . 1 常用 风 速功率 谱 脉动风的风速功率谱主要反映脉动风中各种频 率成 分对 应 的能量分 布 规律 。按方 向可 分 为水平 脉 动风速功率谱和垂直脉动风速功率谱 , 目前使用最 广泛 的是水 平风 速功 率谱 , 其中 D a v e n p o r t 谱 、 H a r r i s 谱[ 8 3 不考虑 湍流积 分尺度 随高度 的变化 情
高耸结构 、 大跨度结构等在风荷载作用下的效应 , 对
结 构 的安 全 性和适 用 性甚 至起 着决 定性 的影 响 。在 结 构 设计 计算 考 虑 风荷 载 时 , 除 了在 平 均风 作 用 下 对建 筑结 构 进行 静 力计 算 外 , 有 时 还 需 对 脉 动 风荷 载下 建 筑 结 构 的动 力 响 应 进 行 分 析 … 。 对 于 重 要 的结 构 , 通 常采用 风 洞 试 验 的方 法来 测试 其 动力 响 应, 但风 洞试 验 方法耗 时耗资 巨 大 , 目前 尚无法 在工 程设 计分 析 领域 中普 遍使 用 。而 通过计 算 机模 拟脉 动风 荷 载 , 可 以很 好 地 帮 助 工 程 人 员 掌握 实 际工 程 结 构 的风振 特性 。事 实 上 , 考 虑 到 在 实 际 的 大气 边 界层 紊 流风 场 中 , 脉 动风 速不 仅是 时 间 的函数 , 而且 随 空 间位置 而变 化 , 是 一 个 随机场 , 数值 模 拟方 法反 而 比实 际记 录更 具 有代 表 性 , 因而 在 实 际 工程 中被 广 泛使 用 一 。
况, 而S i m i u谱 、 H i n o谱 、 K a i m a l 谱 等 考 虑 了
基于单脉冲试验的IGBT模型的电压应力测试分析
基于单脉冲试验的IGBT模型的电压应力测试分析Voltage Stress Test Analysis of IGBT Model Based on Single Pulse Test田建平 (维谛技术(西安)有限公司,西安710075)摘 要:IGBT作为功率设备的核心器件,在电力电子工业领域应用广泛。
为了进一步了解电压应力对IGBT模块的影响,本文搭建了实验样机(选用3代IGBT采用T型三电平拓扑,额定输出线电压315 V,电流230 A,设计输入电压范围500~1 000 V),通过单脉冲试验对不同厂家不同型号的IGBT进行电压应力分析并给出解决方案的可行性。
关键词:IGBT;单脉冲;电压应力IGBT(绝缘栅双极型晶体管)是由MOS(绝缘栅型场效应管)和BJT(双极型三极管)组成的复合全控型电压驱动式功率半导体器件,IGBT作为功率设备的核心器件,在电力电子设备中有着广泛的应用[1,2]。
市场不仅追求着低成本和高功率密度,对性能和可靠性要求也更高[3,4]。
IGBT的开关暂态特性限制着它的最大工作结温、最大开关频率、EMC性能、散热性能、优化电路系统等性能[5]。
为了进一步了解IGBT工作性能,笔者搭建了光伏3代IGBT采用T型三电平拓扑,额定输出线电压315V,电流230A,设计输入电压范围500~1000V。
目前备选IGBT模块为英飞凌F3L400R12PT4_B26、西门康SKiM400TMLI12E4B、富士4MBI400VG-120R-50。
英飞凌IGBT模块已经搭建了实验样机,初步的测试表明英飞凌IGBT模块的关断电压应力很大。
因IGBT桥臂的耐压为1200V,关断时只承受一半的母线电压,电压应力不是问题。
IGBT的钳位耐压值为650V(英飞凌、西门康)或600V(富士),关断电压尖峰问题很严重 [6]。
1 单脉冲测试原理本文设计IGBT测试采用单脉冲测试,开通或关断状态都能测试。
脉冲电流法调控碳钢残余应力的机理及相关实验研究
脉冲电流法调控碳钢残余应力的机理及相关实验研究脉冲电流法调控碳钢残余应力的机理及相关实验研究引言:近年来,随着工业技术的进步,对材料性能的需求也日益提高。
碳钢作为常用材料之一,在许多领域中均扮演着重要的角色。
然而,在加工过程中,碳钢往往会产生残余应力,这对材料的性能和使用寿命产生了负面影响。
寻找一种方法来调控碳钢的残余应力,提高其性能具有重要的理论和实际意义。
1. 脉冲电流法的简介脉冲电流法是一种通过电流作用于材料来改变其物理性质的方法。
它采用高频、高能量的电流脉冲,通过瞬间加热和快速冷却的过程,调控材料的结构和性能。
这种方法已被广泛应用于材料的表面改性、残余应力的调控以及材料成分的改变等方面。
2. 脉冲电流法调控碳钢残余应力的机理碳钢的残余应力主要来源于材料的加工过程中的塑性变形和热冷变形等因素。
而脉冲电流法可以通过两种机制来调控碳钢的残余应力。
2.1 热应力释放机制脉冲电流法中的瞬时加热和快速冷却过程会导致材料的温度变化和相变现象。
这些温度变化和相变现象可以分解和调整材料内部的残余应力,从而达到降低或消除残余应力的目的。
2.2 塑性形变机制脉冲电流法中的瞬态电流脉冲通过材料内部流动的电子和离子来引起局部塑性形变。
这种塑性形变可以改变材料的内部结构和晶体缺陷等,并进一步释放残余应力。
3. 相关实验研究为了验证脉冲电流法调控碳钢残余应力的有效性和可行性,许多研究已经进行并取得了一定的突破。
3.1 温度分布和相变研究通过对脉冲电流法作用下碳钢的温度分布和相变过程进行监测和研究,可以揭示这种方法对残余应力的调控机理。
实验结果表明,脉冲电流法可以有效地调控碳钢的残余应力,并且通过调整脉冲的参数和作用时间等,可以实现不同程度的应力释放。
3.2 结构和性能分析对脉冲电流法作用下的碳钢样品进行显微组织、硬度和力学性能分析,可以评估这种方法对材料的改变程度和性能提升。
实验结果显示,脉冲电流法可以有效地消除碳钢中的残余应力,并同时显著提高其硬度和强度等性能。
厦门大学《风工程》课件-2顺风向风振、风振系数和阵风系数
厦门大学土木工程系硕士课程
(2)Shiotami 建议采用只与两点距离有关的简单表达式。 zi z j coh( zi , z j ) exp Lz
xi x j coh( xi , x j ) exp Lx 2 2 12 zi z j xi x j coh( xi , x j , yi , y j ) exp 2 2 Lz Lx
i j
S Pi Pj ( ) S Pi ( ) S Pj ( ) coh( xi , yi , zi , x j , y j , z j , )
其中,SPi ()、SPj () 分别为 i、j 两点脉动风力的自功 率谱密度函数,两点的坐标分别为 ( xi , yi , zi ) 和 ( x j , y j , z j ) ,两 点处脉动风力的相干函数为 coh( xi , yi , zi , x j , y j , z j , ) 。 上述系列公式可用于对任意结构的脉动风力响应进行计 算,包括了各振型响应以及各振型耦合项响应的贡献。考虑 的模态越多,结果越准确,越接近于时程分析的结果。
2 lk
l
2
T
S P P ( ) l i j
T
yk ,l lk
0
H l (i )
l
S P P ( ) l d i j
厦门大学土木工程系硕士课程
上面对公式的推导过程中, () S P P ( ) 矩阵中的元素 SPP i j 可从风洞试验数据中获得或者由下列公式求得:
厦门大学土木工程系硕士课程
对于高层与高耸结构,如仅考虑第1阶振型的影响,即认 为第1阶振型的响应起主要作用,可得到响应公式如下:
电脉冲轧制6061Al_电热耦合模型及电弧加热机理分析
第3期郑军宝,等:时效处理对Fe-Ni基奥氏体合金耐点蚀性能的影响腐蚀电流密度变小,自腐蚀电压变大,且电化学阻抗谱中有较大的电容-电抗弧极化,表明电化学表面处理能够有效提高合金耐腐蚀性。
电化学表面改性能够在合金表面形成比未极化处理合金更加均匀、致密的钝化膜,以降低合金表面点缺陷的方式来提高合金的耐点蚀性能[20]。
3 结论(1)时效处理可以促进γ′相在Fe-Ni基奥氏体合金中的析出,随着时效时间的增加,γ′析出相密度和尺寸增加。
(2)γ′相在基体中析出增加了Fe-Ni基奥氏体合金表面点缺陷数量,时效后合金钝化膜致密性降低,降低了合金的耐点蚀性能。
(3)恒电位阳极极化改性合金表面,可以有效提高合金表面钝化膜的致密性,改善合金耐点蚀性能。
(4)γ′相的析出降低了合金耐点蚀性能,但对于γ′相对合金钝化膜演变的机理性问题,还需进一步的研究。
(5)通过电化学表面改性,以降低合金表面点缺陷的方式能够有效提高合金的耐腐蚀性,但对于在电化学改性过程中点缺陷的演变,以及表面改性过程中点缺陷演变的机理还需做进一步的研究。
参考文献:[1] ALANO J H,SIQUEIRA R L,DE OLIVEIRA A D,et al. Effect of TCP phase formation on the electro‐chemical corrosion behavior of the nickel-based super‐alloy UNS N26455[J]. Corrosion Science, 2020, 177.[2] KOBAYASHI K, YAMAGUCHI K, HAYAKAWA M,et al. High-temperature fatigue properties of austeniticsuperalloys 718, A286 and 304L[J]. International Jour‐nal of Fatigue, 2008, 30(10-11): 1978-1984.[3] WU H,OSHIDA Y,HAMADA S,et al. Fatigue strength properties of precipitation strengthening stainless steelA286 focused attention on small fatigue crack behaviour[J]. Procedia Engineering,2011,10: 1973-1978.[4] GUO Z, LIANG H, ZHAO M, et al. Effect of boron addition on hydrogen embrittlement sensitivity in Fe-Ni based alloys[J]. Materials Science and Engineer‐ing: A,2010, 527(24-25): 6620-6625.[5] ZHAO J C, RA VIKUMAR V, BELTRAN A M. Phase precipitation and phase stability in nimonic 263[J].Met‐allurgical and Materials Transactions A, 2001, 32(6):1271-1282.[6] DE CICCO H, LUPPO M I, GRIBAUDO L M, et al.Microstructural development and creep behavior inA286 superalloy[J]. Materials Characterization, 2004,52(2): 85-92.[7] BROOKS J A, THOMPSON A W. Microstructure and Hydrogen Effects on Fracture in sthe Alloy A-286 [J].Metallurgical and Materials Transactions A,1993,24(9): 1983-1991.[8] ZHAO M J, GUO Z F, LIANG H, et al. Effect of bo‐ron on the microstructure,mechanical properties andhydrogen performance in a modified A286[J]. Materi‐als Science and Engineering: A, 2010, 527(21-22):5844-5851.[9] RHO B S, NAM S W, XIE X. The effect of test tem‐perature on the intergranular cracking of Nb-A286 al‐loy in low cycle fatigue[J]. Journal of Materials Sci‐ence, 2002, 37(1): 203-209.[10] SINGHAL L K, Martin J W. The mechanism of tensile yield in an age-hardened steel containing γ′(orderedNi3Ti)precipitates[J]. Acta Metallurgica,1968,16(7):947-953.[11] 崔巍,高磊,王东丽,等.钛换热管与管板焊接接头电化学腐蚀性能研究[J].电焊机,2020,50(08):124-128.CUI W, GAO L, WANG D L, et al. Study on electro‐chemical corrosion test of welding joint between TA2heat exchange tube and TA2 tube sheet[J]. ElectricWelding Machine, 2020,50(8):124-128.[12] 姜心怡,严春妍,周倩雯,等. X80管线钢激光-MIG复合焊焊缝耐腐蚀性分析[J]. 电焊机: 2022,52(04):51-57.JIANG X Y, YAN C Y, ZHOU Q W, et al. CorrosionResistance Analysis of Hybrid Laser-MIG WeldedMetal of X80 Pipeline Steel[J]. Electric Welding Ma‐chine, 2022, 52(4): 51-57.[13] BELO M D C,HAKIKI N,FERREIRA M J E A.Semiconducting properties of passive films formed onnickel–base alloys type Alloy 600: influence of the al‐loying elements[J]. 1999, 44(14): 2473-81.[14] BOCKRIS J O M, KHAN S U, KHAN S U. Surface Electrochemistry: A Molecular Level Approache[M].Springer Science & Business Media,1993.下转第79页23Electric Welding MachineVol.54 No.3Mar. 2024第 54 卷 第 3 期2024 年3 月电脉冲轧制6061Al 电热耦合模型及电弧加热机理分析闫志航1, 王文先1*, 卞功波1, 张涵1, 张婷婷21.太原理工大学 材料科学与工程学院,山西 太原 0300242.太原理工大学 机械与运载学院,山西 太原 030024摘 要:通过COMSOL Multiphysics 仿真软件模拟高频脉冲电流激励时两板间界面处电弧的产生。
Ekman模型时间变化风应力系数的伴随参数估计
① 通信作者: 张继才, 副教授, 博士生导师, E-mail: jicai_zhang@ 收稿日期: 2019-01-29, 收修改稿日期: 2019-04-19
980
海洋与湖沼
50 卷
洋环境往往会随空间和时间发生变化, 再加上观测 数据的不足, 各个湍流模型的普适性受到限制, 具体
摘要 本文基于时间分布参数设置, 利用伴随同化方法, 反演了 Ekman 模型中随时间变化的风应
力拖曳系数, 并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中, 研究了参数反演结果
对不同影响因素的响应, 包括: 风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数
独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演 Ekman 模型中
来的(Thacker et al, 1988)。伴随同化将变分原理与最优 控制理论相结合, 流体力学方程组及其初、边值条件 作为约束条件, 使得依据具体问题设计的代价函数达 到最小(Sasaki, 1970)。近几十年以来, 对于未知参数的 确定, 数据同化方法发挥了重要作用, 以伴随同化方 法 和 卡 曼 滤 波 (Kalman Filter) 同 化 参 数 估 计 为 代 表 (Kalman, 1960; Daley, 1991; Evensen, 1992; Fisher et al, 1995; Evensen et al, 1996; Navon, 1997)。同化参数估计 的思路在于, 由于动力过程的复杂性和参数的未知性, 海洋模型所反应的动力系统类似于“黑箱子”, 但是“黑 箱子”的输出却是可以观测的, 即海洋观测资料; 这就 提供了一种可能性, 基于数学物理反问题的思路, 我 们可以从输出(观测资料)去反推输入的控制变量(参 数)。近年来同化参数估计在许多实际问题中得到应用 (Yu et al, 1991; Das et al, 1992; Ullman et al, 1998; Heemink et al, 2002; Lv et al, 2004; Wakelin et al, 2009; Zhang et al, 2009; Kim et al, 2015; Zhang et al, 2015; Gharamti et al, 2017)。Yu 等(1991)验证了 Ekman 伴随 同化模型优化风应力系数和海水涡动黏性系数剖面的
桅杆结构脉动风速模拟与风荷载计算
时间改变的量;而脉动风则具有明显的随机性,一 般可用零均值平稳Gauss随机过程来描述.作用 于结构上任意处的风速为平均风速和脉动风速之
和,风的模拟主要是针对脉动风而言的. 1.1平均风速剖面
平均风沿高度的变化规律,常称为风速梯度 或风速剖面.平均风速随高度的变化的规律常采
[63
Owen
J S,Eccles B J,Choo B S,et a1.The applica—
series modeling for the
tion of auto-regressive time
time-frequency analysis of civil engineering
structures
k=1
l一0,1,2,…,N一1
(11)
2.2脉动风荷载计算 当不考虑结构与风的耦合作用及漩涡影响 时,根据Bernoulli定理,自由流动的风速提供的 单位面积上的风压力为…]
w=妄弦2
∞,l(叫)=exp(一互G7【口w‘△luz))
(6)
(12)
结构任一高度处的瞬时风速Vi为平均风速 矿i与脉动风速口;之和:V;=彰+奶 则作用在结构上的风压w
Z
《 杂 崖
t/s
-∥Hz b)清波叠加法
图3上层纤绳与杆身连接处(高度115 m)的风荷载时程
图5高度125 m处的脉动风自功率谱比较
万方数据
・1060・
武汉理工大学学报(交通科学与工程版)
m
2010年第34卷
为了检验模拟精度,本文将桅杆高度为55
[23
Deodatis G.Simulation of ergodic multivariate
脉动影响系数
脉动影响系数的计算思路 脉动影响系数u 的计算()()()()()()()121HfszXxzl H z z z l z z dzu m z z dzμμμφηφ=⎰⎰其中()()()()()()()()()()()()()()()()1/21100110,,,X X l z l z H H xzlf s z f s z xzH f s z X z z z z z z x x z z z z dxdx dzdz z z z l z z dz ημμμμμμρφφμμμφ-⎡⎤''''''''=∙⎢⎥⎢⎥⎣⎦⎡⎤∙⎢⎥⎣⎦⎰⎰⎰⎰⎰按定义,脉动系数u 涉及()f z μ,()s z μ,()z z μ,()X l z ,()1z φ,()m z 和(),,,xz x x z z ρ''这几个参数,下面对他们分别进行讨论 (1) 脉动系数()f z μ按定义,脉动系数 ()f z μ为脉动风压与平均风压之比。
()()()()()()()0p pf s s z X z z z P z z z l z μσμσμμμω==因脉动风压随高度变化不大,而平均风压随高度而增大,故脉动系数随高度增加而减少。
根据我国实测数据,在保证率98.6%情况下,()f z μ在不同的地貌条件下随离地面高度变化的数值如下表所示 风压脉动系数()f z μ(2)风载体型系数()s z μ因一般工程结构物并不能理想的使自由气流停滞,而是让气流以不同方式在结构表面绕过,因此实际结构物所受的风压并不能直接按式212v gγω=来计算,而需对其进行修正,其修正系数与结构物的体型有关,故称风载体型系数。
(3)风压高度变化系数()z z μ设任意粗糙程度任意高度处的风压力()a z ω,将其与标准粗糙度下标准高度(一般为10米)处的基本风压之比定义为风压高度变化系数()z z μ,即 ()()0a z z z ωμω=将()222aa oas sz vz z v αωω⎛⎫== ⎪⎝⎭和 2200saTs Ta as a H H z z ααωω-⎛⎫⎛⎫= ⎪⎪⎝⎭⎝⎭代入可得()222saaTs Ta z s a s H H z z z z z αααμ-⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭(4)()X l z :体系在高度z 处的宽度 (5)第1振型函数()1z φ结构的第一振型函数()1z φ可按结构力学原理算出。
风应力单位
风应力单位
(原创实用版)
目录
1.风应力单位的定义和重要性
2.风应力单位的计算方法
3.风应力单位的应用领域
正文
风应力单位是用来衡量风力对物体产生的压力的物理量,通常用于建筑、航空、机械等领域的设计和计算。
风应力是风力对物体表面的压力,其单位是帕斯卡(Pa)或牛顿/平方米(N/m)。
风应力单位的计算方法通常基于流体力学原理,其中包括伯努利定理和动量守恒定律。
在计算风应力时,需要考虑风速、风向、物体形状和表面粗糙度等因素。
一般来说,风应力的计算是相当复杂的,需要使用计算机模拟或专业的风工程软件来完成。
风应力单位在建筑设计中尤为重要。
在建筑设计中,需要考虑风力对建筑物的影响,尤其是在高层建筑和大型桥梁等结构中。
通过计算风应力,可以确定建筑物的结构强度和稳定性,以确保在强风条件下不会发生倾覆或破坏。
除了建筑设计,风应力单位还应用于航空和机械工程等领域。
在航空工程中,风应力对飞机的飞行性能和稳定性有重要影响。
因此,在飞机设计中,需要考虑风应力的影响,并进行相应的计算和调整。
在机械工程中,风应力也会对机械设备的运行稳定性和寿命产生影响,因此也需要进行相应的计算和分析。
第1页共1页。
单自由度体系在任意荷载下的脉冲法
ε
t
一个作用时间非常短,但时间积分为有限值的非常大的力称为脉冲力。 ε→0时的力称为单位脉冲。
t=τ时刻
{
• 对阶跃力的反应
具体推导
• 对斜坡力的反应
• 对有有限上升时间的阶跃力的反应
时
时
• 对有有限上升时间的阶跃力的反应
1.在力的上升阶段,体系以 固有周期Tn在静力解答附近 震荡; 2.在力的不变阶段,体系以Tn 在静力解答附近震荡; 3.如果在斜坡结束时的速度 为零,则体系在力的不变阶 段不会震荡; 4.tr/Tn较小时,反应类似于阶 跃力; 5. tr/Tn较大时,震荡的动位 移近似于静力解答。
单自由度系统在任意荷载作用下 的脉冲响应法
一般荷载 特殊荷载
阶跃力 斜坡力 有限上升时间的阶跃力
在许多实际情况中,动力激励既不是谐振荷载 也不是周期荷载。
为了分析单自由度体系承受随时间任意变化的力p(t)作用 的反应,通常将p(t)解释为持续时间为无穷小量的脉冲序 列,体系对p(t)的反应是各单个反应之和。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
d ′N A ,式中 d ′ 和 M 为镁密度和摩尔质量。Mg M 的摩尔质量等于 Mg 的相对原子量, 24.3 。因此, d ′ = 1.736 × 103 kg ⋅ m −3 , M = 24.3 g mol [8] 。 N A 为
首先采用位错动力学计算电子风应力。电子密度 ne = Avogadro常数, 6.02 ×
文章引用: 曹富荣. 脉冲电流作用下电子风应力模型与计算示例[J]. 凝聚态物理学进展, 2017, 6(1): 1-5. https:///10.12677/cmp.2017.61001
曹富荣
子力学的电子风应力与实际应力相差很大,造成量子力学方法计算值偏低的原因是没有考虑合金中存在 的第二相阻碍电子运动的情况。
进行下列国际单位与导出单位换算
= Ω m 2 ⋅ kg ⋅ s −3 ⋅ A −2
C= s ⋅ A
N = m ⋅ kg ⋅ s −2
因此,
ΩCA N = m m2
= Few 4.42 × 103 N m
σ ew =
4 Few l 4 Few 4 4.42 × 103 N m = × = × =70.6 MPa 3π 2.66 × 10−6 m 3π d πd 2
摘
要
电致塑性是利用电场和脉冲电流改善材料塑性的方法。为了揭示脉冲电流作用下的电子风影响塑性的机 理,在Conrad电子风力模型的基础上,提出电子风应力模型。在复杂体系的LAZ922 (Mg-Li-Al-Zn)合金 中给出电子风应力计算实例。计算发现,基于位错动力学的电子风应力与实际应力十分吻合,而基于量
( 0.209 − 3.12 ×10 ) m = 7.26 ×10
−2
s m
12
s
ev ⋅ m −2
由于 1 ev = 1.16 × 10−6 N m 。接着计算电子风应力, = 1.6 × 10−19 J , 1 J = 1 NM ,所以, F ew
σ ew =
焦耳热应力为
4 Few l 4 Few 4 1.16 × 10−6 N m = × = × = 1.85 × 10−8 MPa 3π d 3π πd 2 26.6 × 10−6 m
Open Access
1. 引言
电致塑性(Electro-plasticity (EP))包括电场致塑性和脉冲电流致塑性。 其共同特点是施加外场使应力降 低,塑性改善。脉冲电流致塑性最早由前苏联学者提出,后来美国学者进一步发展,多年来,人们对脉 冲电流致塑性改善进行了不懈研究,取得了一些有价值的成果[1]-[4],其应用的方面在快速成形领域,例 如脉冲电流辅助下快速弯曲成形,大大提高了成形效率。美国学者 Conrad 提出了电子风力的模型,其计 算对于理解脉冲电流改善塑性的机理具有重要作用[5],但是 Conrad 模型的缺陷是参数较难获得, 限制了 应用。文献分析发现,国内外没有详细求解电子风力的计算示例,给研究者了解脉冲电流致塑性机理带 来困难。因此,本文在 Conrad 电子风力模型基础上,提出电子风应力模型,给出计算示例,期望与实际 应力-应变曲线建立联系,有助于加深对 EP 过程机理的理解。
= ne
d ′ N A 1.736 × 103 × 10−3 g m −3 ×6.02 × 1023 mol = = 4.30 × 1028 m3 M 24.3 g mol
F = ew
4.46 × 10−8 Ω ⋅ m 1028 109 A × 1.6 × 10−19 C × 4.30 × 3 × 1.44 × 8 2 ρ 10 m m m2 ΩCA = 4.42 × 103 m2 R en= e j
Keywords
Electro-Plasticity, Pulse Current, Electrion Wind Force, Electron Wind Stress
脉冲电流作用下电子风应力模型与计算示例
曹富荣
东北大学材料科学与工程学院,辽宁 沈阳
收稿日期:2017年2月5日;录用日期:2017年2月24日;发布日期:2017年2月27日
σT = σ − σ ew − σ pc = 250 − 1.85 × 10−8 − 108 = 142 MPa .
, d , ρ , vd ) ,物理 ( ne d ′, e, ne , , p f ) ,化学 ( M , N A ) ,和电学 上述计算表明,计算涉及材料科学 ( σ , ε
( Few , R, j , ve ),因此本文属于交叉学科问题。
= vd ε = ( ρ b ) 10−3 s −1 108 m −2 × 3.21×10−10 m = 3.12 ×10−2 m s
(
)
电子风力为
Few = α bp f ne ( ve − vd ) = 0.33 × 3.21× 10−10 m × 8.96 × 10−6 ev ⋅ × 4.30 × 1028 m3 ×
关键词
电致塑性,脉冲电流,电子风力,电子风应力
Copyright © 2017 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). /licenses/by/4.0/
Advances in Condensed Matter Physics 凝聚态物理学进展, 2017, 6(1), 1-5 Published Online February 2017 in Hans. /journal/cmp https:///10.12677/cmp.2017.61001
j= 1.44 × 103 A ⋅ mm −2 1.44 × 109 A m 2 。 首先 计 算 费米 动 量 p f 。费 米 能量 E f = 7.08 ev 和 费 米速 度
(
)
V f = 1.58 × 10 m s [11]。
6
2E f 2 × 7.08 ev = pf = = 8.96 × 10−6 ev ⋅ s m Vf 1.58 × 106 m s
4. 讨论
计算发现位错动力学的电子风应力与实际情况吻合,反映出脉冲电流对应力降低的影响以及对焦耳 热应力的影响。而采用量子力学方法计算的电子风应力太低,焦耳热应力太高,与实际变形过程不符合。 原因在于 LAZ922 合金内含 α + β 双相和少量第二相颗粒, 在脉冲电流施加合金的瞬间, 电子运动会受到 来自晶界( α α , β β )、相界( α β )和颗粒界面的阻碍。另外,量子力学模型中的费米动量和电子漂移速 度没有考虑到第二相颗粒引起的晶格畸变。其结果造成量子力学方法计算的电子风应力偏低。量子力学 模型中的费米动量是对纯金属或纯晶体(Mg) 的参数费米能量和费米速度进行计算的,实际复杂体系的 LAZ922 合金中 Li、Al 和 Zn 元素对电子的交互作用十分复杂,因而计算的纯 Mg 费米动量与实际合金存 在差 距 , 从而 导 致 量子 力 学方 法 计 算的 电 子 风应 力 偏低 。 解 决办 法 是 对式 (2) 乘上 一 个 修正 系数
= σ 250 = MPa, σ pc 108 MPa 。 根据实验结果,
假设脉冲电流引起电子风应力和焦耳热应力,焦耳热应力为
σT = σ − σ ew − σ pc = 250 − 70.6 − 108 = 71.4 MPa .
其次,采用量子力学计算电子风应力。选取 α = 0.33 , = b 3.21× 10−10 m [10],和
Electron Wind Stress Model and Calculation Examples under Pulse Current
Furong Cao
School of Materials Science and Engineering, Northeastern University, Shenyang Liaoning Received: Feb. 5 , 2017; accepted: Feb. 24 , 2017; published: Feb. 27 , 2017
1023 .R= 4.46 × 10−8 Ω ⋅ m [8]。 ρ = 108 m −2 .= e 1.6 × 10−19 C . l = d 3 [9]和 d 为空间 mol
晶粒尺寸, = d 2.66 × 10−5 m 。
下面计算 ne , Few , σ ew , σ T 。选取实验电流密度 j = 1.44 × 103 A ⋅ mm −2 (1.44 × 10−9 A m 2 ) 。
然后分别计算电子漂移速度 ve 和位错速度 vd 。
3
曹富荣
= ve
这里单位换算C = SA。
j 1.44 × 109 Am −2 = = 0.209 m s ene 1.6 × 10−19 C × 4.3 × 1028 m −3
( ρ b ) [12],这里 ε 为应变速率。 ε = 10−3 s −1 。 vd = ε
2. 模型
Conrad 提出基于位错动力学的电子风力(电流施加在单位长度位错上的力), Few ,模型[6]:
Few = ( R ρ ) ene j
式中 R 为电阻, ρ 为位错密度, e 为电荷, ne 为电子密度, j 为电流密度。 同时,他提出量子力学电子风力模型:
j Few α bp f ne ( ve − = vd ) α bp f − ne vd = e
σ ew =
4 Few l πd 2
(3)
式中 σ ew 为电子风应力和 d 为空间晶粒尺寸或平均晶粒直径。
σ − σ ew − σ pc ,式中 σ 为无脉冲电 假设脉冲电流引起电子风应力和焦耳热应力,焦耳热应力为 σ T =