CN109323830B - 游丝刚度的确认方法 - Google Patents

游丝刚度的确认方法 Download PDF

Info

Publication number
CN109323830B
CN109323830B CN201811259777.4A CN201811259777A CN109323830B CN 109323830 B CN109323830 B CN 109323830B CN 201811259777 A CN201811259777 A CN 201811259777A CN 109323830 B CN109323830 B CN 109323830B
Authority
CN
China
Prior art keywords
balance spring
pixel point
electron back
coordinate system
back scattering
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811259777.4A
Other languages
English (en)
Other versions
CN109323830A (zh
Inventor
陈世佳
黄智恒
龚翔
黄宝渝
宋鹏涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ZHUHAI ROSSINI WATCH INDUSTRY Ltd
National Sun Yat Sen University
Original Assignee
ZHUHAI ROSSINI WATCH INDUSTRY Ltd
National Sun Yat Sen University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ZHUHAI ROSSINI WATCH INDUSTRY Ltd, National Sun Yat Sen University filed Critical ZHUHAI ROSSINI WATCH INDUSTRY Ltd
Priority to CN201811259777.4A priority Critical patent/CN109323830B/zh
Publication of CN109323830A publication Critical patent/CN109323830A/zh
Application granted granted Critical
Publication of CN109323830B publication Critical patent/CN109323830B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0025Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of elongated objects, e.g. pipes, masts, towers or railways
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G01N23/203Measuring back scattering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/05Investigating materials by wave or particle radiation by diffraction, scatter or reflection
    • G01N2223/053Investigating materials by wave or particle radiation by diffraction, scatter or reflection back scatter

Landscapes

  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本申请涉及一种游丝刚度的确认方法。包括S10,对游丝进行化学成分分析;S20,根据所述游丝的化学成分分析结果计算得到所述游丝晶粒在其初始取向坐标系下的柔度张量分量;S30,建立工作坐标系,并获取游丝工作面的电子背散射衍射图谱;S40,获取所述电子背散射衍射图谱中每个像素点处晶粒对应所述工作坐标系的取向欧拉角;S50,计算得到所述电子背散射衍射图谱中每个像素点处在所述游丝的拉伸方向上的弹性模量和相应的泊松比;S60,获取所述游丝的等效弹性模量;S70,计算得到所述游丝的刚度。本申请提出的游丝刚度的确认方法从游丝微观晶粒结构着手,以像素点为弹性单元对游丝弹性参数进行了定量化分析与计算,具有节省成本、普适性高、精度高的优点。

Description

游丝刚度的确认方法
技术领域
本申请涉及游丝弹性参数确认领域,特别是涉及一种游丝刚度的确认方法。
背景技术
在机械手表机心中,游丝是最关键的弹性元件之一,其主要成分为铁镍合金。游丝为摆轮来回摆动提供恢复力矩,而恢复力矩又与游丝的刚度密切相关。因此,准确地计量游丝的刚度是预测机械手表机心动力学响应的关键。但是,游丝弹簧的横截面面积仅仅约为3×10-3mm2,整个游丝的螺旋直径也仅约为5mm,随着材料宏观尺寸的降低,其微观组织结构对其宏观性能的影响越来越明显;并且游丝在生产过程中,经历了多道次的冷变形以及热处理,其晶粒结构发生了严重的变形,不能简单地采用普通的铁镍合金力学性能。
目前在国内并无专门测量游丝弹性性能的仪器设备。国内对游丝仅仅进行尺寸测量,要准确测量如弹性模量E、刚度k等游丝的弹性性能是非常困难的,且所用的测量设备成本高昂,其成本高达百万人民币,并不普及,且维修保养费用极高,同时对于不同型号的游丝,需要定制夹具,普适性不高,耗时耗力。
发明内容
基于此,有必要针对准确测量游丝弹性性能困难、测量设备昂贵、普适性不高的问题,提供一种简单、节省成本且普适性高的游丝刚度确认方法。
本申请提供的一种游丝刚度的确认方法,包括:
S10,对游丝进行化学成分分析,获得所述游丝中含有的化学元素以及每种化学元素对应的质量百分比αX
S20,根据所述游丝的化学成分分析结果计算得到所述游丝晶粒在其初始取向坐标系下的柔度张量分量
Figure BDA0001843630710000021
S30,建立工作坐标系,所述工作坐标系的X轴沿所述游丝的拉伸方向设置,且所述工作坐标系遵循右手法则。通过电子背散射衍射扫描,获取游丝工作面相对于工作坐标系的电子背散射衍射图谱;
S40,根据所述电子背散射衍射图谱中每个像素点的色阶值获取所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000022
S50,利用所述柔度张量
Figure BDA0001843630710000023
以及所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000024
计算得到所述电子背散射衍射图谱中每个像素点处在所述游丝拉伸方向上的弹性模量E’和相应的泊松比ν’;
S60,根据所述电子背散射衍射图谱中每个像素点处在所述游丝拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝的等效弹性模量Eeq
S70,根据所述游丝的等效弹性模量Eeq计算得到所述游丝的刚度k。
上述游丝刚度的确认方法,基于游丝的化学成分和微观晶粒结构,通过试验、计算、模拟等手段对其弹性性能进行定量化分析和计算。上述游丝刚度的确认方法仅需要化学成分分析仪器、电子背散射衍射扫描仪和计算机就可以完成,无需购买昂贵的测量游丝刚度的设备,因此,具有节省成本的优点。因为上述游丝刚度的确认方法没有用到复杂设备,所以不需要花太长时间进行设备调试,因此,具有节省时间的优点。上述游丝刚度的确认方法在整个确认过程中除了化学成分分析过程和电子背散射衍射扫描过程,其余步骤都可以通过计算机编程实现,因此,具有方便快捷的优点。因为对于所有游丝都可进行化学成分分析和电子背散射衍射扫描,因此,具有普适性高的优点。另外,上述游丝刚度的确认方法仅包括材料检测和数值计算过程,没有任何物理参数测量过程,大大减小了测量误差,因此,具有精度高的优点。最后,上述游丝刚度的确认方法,在确认过程中以像素点作为所述游丝的弹性单元,而不是以单个晶粒为弹性单元。这样就不需要人为消除晶粒内部的差异,将晶粒内部亚结构对游丝弹性性能的影响考虑进来,使得确认结果更加准确。
在其中一个实施例中,所述S60包括:
S61,建立一个与所述电子背散射衍射图谱同等尺寸的二维有限元模型;
S62,将所述电子背散射衍射图谱中每个像素点处在所述游丝拉伸方向上的弹性模量E’和相应的泊松比ν’对应地导入有限元模型;
S63,使所述有限元模型模拟发生弹性变形过程,获得所述游丝的等效弹性模量Eeq
在其中一个实施例中,在所述S60中,所述等效弹性模量Eeq的值为所述电子背散射衍射图谱中每个像素点处在所述游丝的拉伸方向上的弹性模量E’的算术平均值。
在其中一个实施例中,在所述S60中,还根据所述电子背散射衍射图谱中每个像素点处在所述游丝拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝的等效泊松比νeq
在其中一个实施例中,在所述S30中,仅对所述游丝工作面的一部分进行电子背散射衍射扫描,所述游丝工作面被扫描的部分长度为游丝截面最长边的1倍~3倍。
在其中一个实施例中,所述S40包括:
S41,提取所述电子背散射衍射图谱中每个像素点的RGB值;
S42,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000041
Figure BDA0001843630710000042
值,所述
Figure BDA0001843630710000043
值等于对应像素点R值的
Figure BDA0001843630710000044
倍,其中,a的取值范围为0~360;
S43,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000045
的Φ值,所述Φ值等于对应像素点G值的
Figure BDA0001843630710000046
倍,其中,b的取值范围为0~360;
S44,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000047
Figure BDA0001843630710000048
值,所述
Figure BDA0001843630710000049
值等于对应像素点B值的
Figure BDA00018436307100000410
倍,其中,c的取值范围为0~360。
在其中一个实施例中,在所述S10中,采用电子显微镜、X射线衍射光谱分析法或者原子吸收光谱分析法中的一种对所述游丝进行化学成分分析。
在其中一个实施例中,在所述S20中,通过对所述游丝中包含的所有化学元素对应的柔度张量分量
Figure BDA00018436307100000411
加权求和得到所述柔度张量分量
Figure BDA00018436307100000412
加权系数为每种化学元素对应的质量百分比αX
在其中一个实施例中,在所述S30中,电子背散射衍射扫描的扫描步长为0.01μm~1μm。
在其中一个实施例中,在所述S30和所述S40之间还包括:
S80,对所述电子背散射衍射图进行去噪声处理。
本申请提供的实施例中,仅对所述游丝工作面的一部分进行电子背散射衍射扫描。由于所述游丝的材料内部几乎都是由于冷变形引起的长条晶粒,是典型的织构组织,其方向与冷加工方向相同,具有一定的重复性。因此,需选取适当的长度进行扫描可以在保证准确性的同时节省大量的时间和成本。另外,在另一个本申请提供的实施例中,在步骤S30和步骤S40之间还包括步骤S80。步骤S80为对所述电子背散射衍射图谱进行去噪声处理。增加步骤S80使得所述电子背散射衍射图谱中每个像素点的色阶值更加准确,从而使得最终游丝刚度的确认结果更加准确。
附图说明
图1为本申请一个实施例中提供的游丝刚度的确认方法的流程图;
图2为本申请一个实施例中提供的游丝刚度的确认方法测量的游丝的结构示意图以及局部放大图;
图3为本申请一个实施例中通过电子背散射衍射扫描获取的游丝工作面的电子背散射衍射图谱;
图4为本申请一个实施例中根据所述电子背散射衍射图每个像素点的色阶值获取每个像素点对应所述游丝工作面位置处晶粒取向欧拉角时使用的RGB色度尺。
符号说明:
游丝 10
游丝工作面 11
游丝截面 12
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下通过实施例,并结合附图,对本申请的游丝刚度的确认方法进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本申请,并不用于限定本申请。
请参见图1-2,本申请一个实施例中提供一种游丝刚度的确认方法。所述游丝刚度的确认方法包括:
S10,对游丝10进行化学成分分析,获得所述游丝10中含有的化学元素以及每种化学元素对应的质量百分比αX
步骤S10中,αX为化学元素X的质量百分比,上标X代表化学元素。在一个实施例中,可以采用电子显微镜、X射线衍射光谱分析法或者原子吸收光谱分析法中的一种对所述游丝10进行化学成分分析。
S20,根据所述游丝10的化学成分分析结果计算得到所述游丝10晶粒在其初始取向坐标系下的柔度张量分量
Figure BDA0001843630710000061
步骤S20中,通过对所述游丝10中包含的所有化学元素晶粒在其初始取向坐标系下对应的柔度张量分量
Figure BDA0001843630710000062
加权求和得到所述柔度张量
Figure BDA0001843630710000063
加权系数为每种化学元素对应的质量百分比αx。以上计算可以用公式表示为:
Figure BDA0001843630710000064
其中,
Figure BDA0001843630710000065
为所述游丝10的合金材料的柔度张量分量。上标spring代表所述游丝10的合金材料。
Figure BDA0001843630710000066
为化学元素X的柔度张量分量。通常柔度张量S为四阶张量,有81个分量。所述柔度张量分量Sijkm的下标i、j、k和m的值都可以为1、2或3中的一个。对于一般的各向异性材料,胡克定律可用张量符号表示为:
εij=Sijkmσkm(i,j,k,m=1,2,3)
其中,ε和σ分别是应变张量和应力张量,S是柔度张量。对于具有立方对称结构的晶体材料,例如所述游丝10的合金材料,在晶体的初始取向坐标系下柔度张量S可以简化为以下公式:
Figure BDA0001843630710000071
所述柔度张量S的81个分量退化为3个独立变量分量。所述3个独立变量分量分别为S1111、S1122和S1212。因此,在一个实施例中,在步骤S20中可以只计算
Figure BDA0001843630710000072
Figure BDA0001843630710000073
这3个柔度张量分量。所有合金元素X的柔度张量分量都可以从现有的文献中查找到。在一个实施例中,在计算所述柔度张量分量
Figure BDA0001843630710000074
对,可以仅对质量百分比αX高于10%的合金元素对应的柔度张量
Figure BDA0001843630710000075
加权求和,加权系数为每种化学元素对应的质量百分比αX
S30,建立工作坐标系,所述工作坐标系的X轴沿所述游丝10的拉伸方向设置,且所述工作坐标系遵循右手法则。通过电子背散射衍射扫描,获取游丝工作面11相对于工作坐标系的电子背散射衍射图谱。
步骤S30中,在进行电子背散射衍射(Electron Backscattered Diffraction,EBSD)扫描之前,可以根据《GB/T19501-2004电子背散射衍射分析方法通则》国家标准,对所述游丝10进行切割、打磨、机械抛光、去应力电解抛光预处理,使得得到的电子背散射衍射图谱更加准确。因为游丝材料内部几乎都是由于冷变形引起的长条晶粒,是典型的织构组织,其方向与冷加工方向相同,具有一定的重复性。因此,并不需要对游丝沿长度的整个工作面进行电子背散射衍射扫描,只需要选取适当的长度进行扫描,这样可以节省大量的时间和成本。在一个实施例中,仅对所述游丝工作面11的一部分进行电子背散射衍射扫描,所述游丝工作面11被扫描的部分长度为游丝截面12最长边的1倍~3倍。所述电子背散射衍射扫描的步长可以根据实际需要进行选择。通常,扫描步长越小测量越准确,但是扫描步长越小测量成本越高。在一个实施例中,扫描步长可以为0.01μm~1μm,优选地,扫描步长为0.1μm。
S40,根据所述电子背散射衍射图谱中每个像素点的色阶值获取所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000081
步骤S40中,将所述电子背散射衍射图谱中每个像素点的RGB值进行比例放大或缩小,从而与所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000082
一一对应。在一个实施例中,步骤S40包括以下步骤:
S41,提取所述电子背散射衍射图谱中每个像素点的RBG值。
S42,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000083
Figure BDA0001843630710000084
值。所述
Figure BDA0001843630710000085
值等于对应像素点R值的
Figure BDA0001843630710000086
倍。其中,a的取值范围为0~360。优选地,a的值为0、90、180、270或360中的一个值。
S43,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000087
的Φ值。所述Φ值等于对应像素点G值的
Figure BDA0001843630710000088
倍。其中,b的取值范围为0~360。优选地,b的值为0、90、180、270或360中的一个值。
S44,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000089
Figure BDA00018436307100000810
值。所述
Figure BDA00018436307100000811
值等于对应像素点B值的
Figure BDA00018436307100000812
倍。其中,c的取值范围为0~360。优选地,c的值为0、90、180、270或360中的一个值。
所述步骤S40可以通过Matlab、C或者C++等编程语言实现。在一个实施例中,在步骤S30和步骤S40之间,还包括步骤S80。步骤S80为对所述电子背散射衍射图谱进行去噪声处理。去除的噪声可以为环境噪声或系统噪声等常见噪声。加入步骤S80使得所述电子背散射衍射图谱更加准确。从而使得最终的游丝刚度确认结果更加准确。
S50,利用所述柔度张量分量
Figure BDA0001843630710000091
以及所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000092
计算得到所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’。
步骤S05中,在一个实施例中,所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’可以采用以下公式计算:
Figure BDA0001843630710000093
Figure BDA0001843630710000094
其中,所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000095
方向上的柔度张量分量S′ijkm可以用以下公式计算:
Figure BDA0001843630710000096
其中,δpq是克罗内克函数(Kronecker delta)。A为基于欧拉角
Figure BDA0001843630710000097
的变换矩阵。所述变换矩阵A为二维矩阵。Ast为所述变换矩阵A中第s行第t列的元素。上式遵循爱因斯坦求和约定(Einstein summation convention),且u的值为1、2和3。所述变换矩阵A可以表示为:
Figure BDA0001843630710000098
Figure BDA0001843630710000101
S60,根据所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝10的等效弹性模量Eeq
步骤S60中,可以采用软件模拟的方式获取所述游丝10的等效弹性模量Eeq。在一个实施例中,所述S60包括:
S61,建立一个与所述电子背散射衍射图谱同等尺寸的二维有限元模型。
S62,将所述电子背散射衍射图谱中每个像素点处在所述游戏10拉伸方向上的弹性模量E’和相应的泊松比ν’对应地导入有限元模型。
S63,使所述有限元模型模拟发生弹性变形过程,获得所述游丝10的等效弹性模量Eeq
在本实施例中,模拟过程可以采用任何有限元模型模拟软件完成。
在一个实施例中,步骤S60中,所述等效弹性模量Eeq的值为所述电子背散射衍射图谱中每个像素点处在所述游丝10的拉伸方向上的弹性模量E’的算术平均值。本实施例提出的计算等效弹性模量Eeq的方法简化了计算步骤,且在没有有限元模型模拟软件的情况下也可以使用。因此本实施例提出的计算等效模量Eeq的方法具有计算简单、适用性强的优点。
在一个实施例中,在所述S60中,还根据所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的的弹性模量E’和相应的泊松比ν’获取所述游丝10的等效泊松比νeq。所述游丝10的等效泊松比νeq可以采用软件模拟的方式获取。通过软件模拟获取等效泊松比νeq的方式与通过软件模拟获取所述等效弹性模量Eeq的方式相同。所述游丝10的等效泊松比νeq的值也可以为所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的泊松比ν’的算术平均值。因此,本申请提出的游丝刚度的确认方法也可以用于确认所述游丝10的等效泊松比νeq
S70,根据所述游丝10的等效弹性模量Eeq计算得到所述游丝10的刚度k。
步骤S70中,计算所述游丝刚度k的公式为:
Figure BDA0001843630710000111
其中,b为所述游丝10的宽度,也就是所述游丝截面12的最长边。h为所述游丝10的厚度,也就是所述游丝截面12的最短边。Lg为所述游丝10的工作长度,也就是所述游丝工作面12的长度。
本申请提出的游丝刚度的确认方法可以采用多次测量取平均的方式降低测量误差。本申请提出的游丝刚度的确认方法还可以采用测量多个样品取平均的方式降低测量误差。
请参见图2-4,为了更好的解释本申请提出的游丝刚度的确认方法,这里给出一个完整的具体实施例。
具体实施例为:
1、采用电子显微镜对待测游丝进行化学成分分析,获得待测游丝的化学成分表。所得化学成分表为表1。
表1采用电子显微镜测出的游丝化学成分表
化学元素 Cr Ni Ti Fe
质量百分比/wt% 5.91 39.49 1.86 balance
2、在本实施例中,Cr、Ti含量较少,对结果影响较小,可忽略。因此,元素质量比近似取Fe∶Ni=6∶4。元素质量比也可以按照各元素实际的比值进行计算。查阅现有文献可知
Figure BDA0001843630710000112
的值为0.76×10-2GPa-1
Figure BDA0001843630710000113
的值为-0.287×10-2GPa-1
Figure BDA0001843630710000114
的值为0.892×10- 2GPa-1
Figure BDA0001843630710000115
的值为0.734×10-2GPa-1
Figure BDA0001843630710000116
的值为-0.274×10-2GPa-1
Figure BDA0001843630710000117
的值为0.802×10-2GPa-1。因为
Figure BDA0001843630710000118
所以,在本实施例中,可以计算得到
Figure BDA0001843630710000121
Figure BDA0001843630710000122
计算公式以及计算结果如下:
Figure BDA0001843630710000123
Figure BDA0001843630710000124
Figure BDA0001843630710000125
3、取一部分游丝工作面制成样品,对所述样品进行电子背散射衍射扫描,得到电子背散射衍射图谱。所得电子背散射衍射图谱如图3所示。所述电子背散射衍射图谱尺寸大小为为50.2m×48.6m,像素栅格为503×486。所述电子背散射衍射图谱的尺寸与像素栅格数的比大约为1∶1。由于游丝在加工过程中经历了多道次冷拉拔,所以图3可以看出其晶粒沿着拉伸方向变长,呈现出大量长条状晶粒。图3中绘制的坐标系为所建立的工作坐标系。所述游丝10的拉伸方向为图3中所述工作坐标系的X轴方向。
4、用Matlab编程来提取所述电子背散射衍射图谱中每个像素点的RGB值。并用Matlab编程将所述电子背散射衍射图谱中每个像素点的RGB值按照图4中的色度尺进行放大或缩小。这里将所述电子背散射衍射图谱中每个像素点的R值放大为原来值的
Figure BDA0001843630710000126
倍,G值缩小为原来值的
Figure BDA0001843630710000127
倍,B值缩小为原来值的
Figure BDA0001843630710000128
倍。所述电子背散射衍射图谱中每个像素点放大后的R值等于每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000129
Figure BDA00018436307100001210
值。所述电子背散射衍射图谱中每个像素点缩小后的G值等于每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA00018436307100001211
的Φ值。所述电子背散射衍射图谱中每个像素点缩小后的B值等于每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000131
Figure BDA0001843630710000132
值。
5、通过Matlab编程,利用所述柔度张量分量
Figure BDA0001843630710000133
以及所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure BDA0001843630710000134
计算得到所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’。
6、根据所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝10的等效弹性模量Eeq和等效泊松比νeq。获取所述等效弹性模量Eeq和所述等效泊松比νeq的方法有两种。具体获取方法如下:
方案1:首先,在有限元模型模拟软件中建立一个尺寸为50.2m×48.6m的二维有限元模型。随后,在材料设置时,将所述电子背散射衍射图谱中每个像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’与像素点对应地一一导入有限元模型的材料中。这样就在有限元模型模拟软件中建立了一个能够与所述电子背散射衍射图谱中的微观晶粒结构一一对应的材料模型。最后,对有限元模型进行网格划分、边界条件设定、加载等一系列操作,使所述材料模型发生弹性变形,获得所述游丝10的等效弹性模量Eeq及等效泊松比νeq。经模拟,得到Eeq=225GPa,νeq=0.29。
方案2:计算所述电子背散射衍射图谱所有像素点处在所述游丝10拉伸方向上的弹性模量E’和相应的泊松比ν’的算术平均值,获得平均弹性模量Eave及平均泊松比νave。经计算,Eave=223GPa,νave=0.29。所述游丝样品的等效弹性模量Eeq的值为平均弹性模量Eave。所述游丝样品的等效泊松比νeq的值为平均泊松比νave
7、取方案1中的Eeq=225GPa,νeq=0.29。在本实施例中,所述游丝10的宽度b=1.1×10-1mm。所述游丝10的厚度h=2.9×10-2mm。所述游丝10的工作长度Lg=147mm。根据公式
Figure BDA0001843630710000135
可以计算得到游丝刚度k=3.4×10-7N.m/rad。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种游丝刚度的确认方法,其特征在于,包括:
S10,对游丝(10)进行化学成分分析,获得所述游丝(10)中含有的化学元素以及每种化学元素对应的质量百分比αX
S20,根据所述游丝(10)的化学成分分析结果计算得到所述游丝(10)晶粒在其初始取向坐标系下的柔度张量分量
Figure FDA0002347149070000011
S30,建立工作坐标系,所述工作坐标系的X轴沿所述游丝(10)的拉伸方向设置,且所述工作坐标系遵循右手法则,通过电子背散射衍射扫描,获取游丝工作面(11)相对于工作坐标系的电子背散射衍射图谱;
S40,根据所述电子背散射衍射图谱中每个像素点的色阶值获取所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000012
S50,利用所述柔度张量分量
Figure FDA0002347149070000013
以及所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000014
计算得到所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000015
方向上的柔度张量分量S′ijkm,以及根据相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000016
方向上的柔度张量分量S′ijkm计算得到所述电子背散射衍射图谱中每个像素点处在所述游丝(10)拉伸方向上的弹性模量E’和相应的泊松比ν’,其中
Figure FDA0002347149070000017
Figure FDA0002347149070000018
Figure FDA0002347149070000019
δpq是克罗内克函数,A为基于欧拉角
Figure FDA0002347149070000021
的二维变换矩阵,且A遵循爱因斯坦求和约定,u的取值为1、2和3;
S60,根据所述电子背散射衍射图谱中每个像素点处在所述游丝(10)拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝(10)的等效弹性模量Eeq
S70,根据所述游丝(10)的等效弹性模量Eeq计算得到所述游丝(10)的刚度k。
2.根据权利要求1所述的游丝刚度的确认方法,其特征在于,所述S60包括:
S61,建立一个与所述电子背散射衍射图谱同等尺寸的二维有限元模型;
S62,将所述电子背散射衍射图谱中每个像素点处在所述游丝(10)拉伸方向上的弹性模量E’和相应的泊松比ν’对应地导入有限元模型;
S63,使所述有限元模型模拟发生弹性变形过程,获得所述游丝(10)的等效弹性模量Eeq
3.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S60中,所述等效弹性模量Eeq的值为所述电子背散射衍射图谱中每个像素点处在所述游丝(10)的拉伸方向上的弹性模量E’的算术平均值。
4.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S60中,还根据所述电子背散射衍射图谱中每个像素点处在所述游丝(10)拉伸方向上的弹性模量E’和相应的泊松比ν’获取所述游丝(10)的等效泊松比νeq
5.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S30中,仅对所述游丝工作面(11)的一部分进行电子背散射衍射扫描,所述游丝工作面(11)被扫描的部分长度为游丝截面(12)最长边的1倍~3倍。
6.根据权利要求1所述的游丝刚度的确认方法,其特征在于,所述S40包括:
S41,提取所述电子背散射衍射图谱中每个像素点的RGB值;
S42,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000031
Figure FDA0002347149070000032
值,所述
Figure FDA0002347149070000033
值等于对应像素点R值的
Figure FDA0002347149070000034
倍,其中,a的取值范围为0~360;
S43,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000035
的Ф值,所述Ф值等于对应像素点G值的
Figure FDA0002347149070000036
倍,其中,b的取值范围为0~360;
S44,计算所述电子背散射衍射图谱中每个像素点处晶粒相对于所述工作坐标系的取向欧拉角
Figure FDA0002347149070000037
Figure FDA0002347149070000038
值,所述
Figure FDA0002347149070000039
值等于对应像素点B值的
Figure FDA00023471490700000310
倍,其中,c的取值范围为0~360。
7.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S10中,采用电子显微镜、X射线衍射光谱分析法或者原子吸收光谱分析法中的一种对所述游丝(10)进行化学成分分析。
8.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S20中,通过对所述游丝(10)中包含的所有化学元素对应的柔度张量分量
Figure FDA00023471490700000311
加权求和得到所述柔度张量分量
Figure FDA00023471490700000312
加权系数为每种化学元素对应的质量百分比αX
9.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S30中,电子背散射衍射扫描的扫描步长为0.01μm~1μm。
10.根据权利要求1所述的游丝刚度的确认方法,其特征在于,在所述S30和所述S40之间还包括:
S80,对所述电子背散射衍射图谱进行去噪声处理。
CN201811259777.4A 2018-10-26 2018-10-26 游丝刚度的确认方法 Active CN109323830B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811259777.4A CN109323830B (zh) 2018-10-26 2018-10-26 游丝刚度的确认方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811259777.4A CN109323830B (zh) 2018-10-26 2018-10-26 游丝刚度的确认方法

Publications (2)

Publication Number Publication Date
CN109323830A CN109323830A (zh) 2019-02-12
CN109323830B true CN109323830B (zh) 2020-05-15

Family

ID=65259034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811259777.4A Active CN109323830B (zh) 2018-10-26 2018-10-26 游丝刚度的确认方法

Country Status (1)

Country Link
CN (1) CN109323830B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202939103U (zh) * 2012-09-19 2013-05-15 天津海鸥表业集团有限公司 动态硅游丝刚度和摆轮转动惯量的测量装置
CN103714216A (zh) * 2013-12-31 2014-04-09 北京理工大学 基于三维微观晶须构成的复合材料力学性能评估方法
CN103871059A (zh) * 2014-03-13 2014-06-18 南京航空航天大学 纤维增强复合材料等效弹性参数的计算方法
CN104090033A (zh) * 2014-07-16 2014-10-08 大连理工大学 基于ebsd图谱的粗晶材料fdtd超声检测仿真模型建立方法
CN104316388A (zh) * 2014-07-25 2015-01-28 中国航空工业集团公司北京航空材料研究院 一种对各向异性材料结构件进行疲劳寿命测定的方法
CN106484978A (zh) * 2016-09-28 2017-03-08 北京理工大学 一种基于晶体滑移机制的各向异性线弹性本构的建立方法
JP2018016823A (ja) * 2016-07-25 2018-02-01 古河電気工業株式会社 放熱部材用銅合金板材およびその製造方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202939103U (zh) * 2012-09-19 2013-05-15 天津海鸥表业集团有限公司 动态硅游丝刚度和摆轮转动惯量的测量装置
CN103714216A (zh) * 2013-12-31 2014-04-09 北京理工大学 基于三维微观晶须构成的复合材料力学性能评估方法
CN103871059A (zh) * 2014-03-13 2014-06-18 南京航空航天大学 纤维增强复合材料等效弹性参数的计算方法
CN104090033A (zh) * 2014-07-16 2014-10-08 大连理工大学 基于ebsd图谱的粗晶材料fdtd超声检测仿真模型建立方法
CN104316388A (zh) * 2014-07-25 2015-01-28 中国航空工业集团公司北京航空材料研究院 一种对各向异性材料结构件进行疲劳寿命测定的方法
JP2018016823A (ja) * 2016-07-25 2018-02-01 古河電気工業株式会社 放熱部材用銅合金板材およびその製造方法
CN106484978A (zh) * 2016-09-28 2017-03-08 北京理工大学 一种基于晶体滑移机制的各向异性线弹性本构的建立方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Determination of hardness and elastic modulus inverse pole figures of a polycrystalline commercially pure titanium by coupling nanoindentation and EBSD techniques;FIZANNE-MICHEL C等;《Materials Science & Engineering》;20141231(第8期);第159-162页 *
新型二维材料的声子输运与热机械性质的数值模拟;孙红义;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20170315(第3期);第19-77页 *

Also Published As

Publication number Publication date
CN109323830A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
Cooper et al. Measuring galaxy environments with deep redshift surveys
Baudoin et al. Influence of critical resolved shear stress ratios on the response of a commercially pure titanium oligocrystal: crystal plasticity simulations and experiment
Liu et al. A combined experimental-numerical approach for determining mechanical properties of aluminum subjects to nanoindentation
US9069917B2 (en) Component fracture evaluation device, component fracture evaluation method and computer program
Huang et al. Estimation of the mechanical properties of steel via LIBS combined with canonical correlation analysis (CCA) and support vector regression (SVR)
CN109933925B (zh) 一种金属板材的冲压成形性能预测方法
CN110261247B (zh) 金属材料各向异性屈服及硬化本构参数同步表征方法
Bertin et al. Integrated digital image correlation applied to elastoplastic identification in a biaxial experiment
Mendiguren et al. Elastic behaviour characterisation of TRIP 700 steel by means of loading–unloading tests
Varma et al. Investigation of stress relaxation mechanisms for ductility improvement in SS316L
WO2011065428A1 (ja) 不良要因の分析表示方法および不良要因の分析表示装置
Buchheit et al. Evaluating deformation-induced grain orientation change in a polycrystal during in situ tensile deformation using EBSD
JP5710997B2 (ja) 疲労限度特定システムおよび疲労限度特定方法
CN109323830B (zh) 游丝刚度的确认方法
CN102906851B (zh) 分析质谱的方法及系统
Hong et al. Residual stress analysis in deep drawn twinning induced plasticity (TWIP) steels using neutron diffraction method
Oliveira et al. Analysis of heterogeneous tests for sheet metal mechanical behavior
Yang et al. Quantitative measurement of plastic strain field at a fatigue crack tip
CN101900698A (zh) 一种测定高Cr耐热钢中Delta铁素体相含量的方法
Butz et al. Parameter identification of 3D yield functions based on a virtual material testing procedure
Dabiri et al. Experimental fatigue characterization and elasto-plastic finite element analysis of notched specimens made of direct-quenched ultra-high-strength steel
EP3789956A1 (en) Method and system for analyzing the mechanical response of microstructures under investigation
CN111982635A (zh) 一种材料内部夹杂的弹性模量测定方法
Baxevani et al. The modified Rockwell test: a new probe for mechanical properties of metals
Chollacoop et al. Robustness of the algorithms for extracting plastic properties from the instrumented sharp indentation data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant