脉冲电流作用下电子风应力模型与计算示例
脉动风场的模拟方法及其在输电线路风振计算中的应用

脉动风场的模拟方法及其在输电线路风振计算中的应用沈国辉;黄俏俏;郭勇;邢月龙;楼文娟;孙炳楠【摘要】针对脉动风场各种模拟方法的适用性问题,在相同的输电塔线实例上进行模拟方法的应用,对比各种方法的计算效率和计算结果,分析各种模拟方法的适用性和等价性,最后探讨单塔和塔线体系计算结果的差异.研究表明:基于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 谱 等 考 虑 了
- 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 =