CN111460700A - 一种基于传递耗散修正的结构振动时效频率获取方法 - Google Patents
一种基于传递耗散修正的结构振动时效频率获取方法 Download PDFInfo
- Publication number
- CN111460700A CN111460700A CN202010122549.3A CN202010122549A CN111460700A CN 111460700 A CN111460700 A CN 111460700A CN 202010122549 A CN202010122549 A CN 202010122549A CN 111460700 A CN111460700 A CN 111460700A
- Authority
- CN
- China
- Prior art keywords
- dissipation
- effective
- frequency
- residual stress
- value
- 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.)
- Granted
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D10/00—Modifying the physical properties by methods other than heat treatment or deformation
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D11/00—Process control or regulation for heat treatments
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/20—Recycling
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Thermal Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Mechanical Engineering (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于传递耗散修正的结构振动时效频率获取方法,该方法包括如下步骤:S1、基于有限元手段的结构固有频率的确定;S2、基于焊接连接关系的耗散有效面的确定;S3.基于残余应力顺位的耗散面有效畸变率的确定;S4、基于耗散畸变率的耗散面有效固有频率的确定;S5、基于传递耗散修正的最佳振动时效频率计算。该方法检测精度高,对于在振动时效流程前精确获取待分析结构的最佳振动时效频率具有重要的现实意义。
Description
技术领域
本发明涉及结构振动时效频率检测方法,特别涉及一种基于传递耗散修正的结构振动时效频率获取方法。
背景技术
随着工业的迅速发展,工件在制造过程中,将受到来自各种工艺等因素的作用与影响,若构件所受到的上述作用与影响不能随之完全消失,仍有部分作用与影响残留在构件内将形成残余应力;如焊接残余应力是焊件产生变形、开裂等工艺缺陷的主要原因,甚至可使焊缝特别是定位焊缝部分完全断开,造成严重的安全事故。目前,振动时效处理是工程材料常用的一种消除其内部残余内应力的方法,其使工件产生最大的振幅,得到最大的动应力和动能量,从而使工件中的残余应力消除的更彻底,工件获得的尺寸稳定性效果更好;传统方法是通过选择整体结构的固有频率作为最佳激振频率,但该方法得到的振动时效频率具有一定的局限性,不能够考虑不同焊接面连接关系对频率耗散过程的影响。
发明内容
发明目的:本发明目的是提供一种基于传递耗散修正的结构振动时效频率获取方法。
技术方案:本发明提供一种基于传递耗散修正的结构振动时效频率获取方法,包括如下步骤:
S1、基于有限元手段的结构固有频率的确定;
S2、基于焊接连接关系的耗散有效面的确定;
S3、基于残余应力顺位的耗散面有效畸变率的确定;
S4、基于耗散畸变率的耗散面有效固有频率的确定;
S5、基于传递耗散修正的最佳振动时效频率计算。
进一步地,所述S1结构固有频率的确定方法如下:将建立完毕的结构有限元分析模型导入到有限元分析软件中,然后对其进行网格划分与边界约束条件的施加,待前处理完毕后,在谐响应模块中对该结构在前处理条件下的第三阶固有频率进行求解,标记为ω。
进一步地,所述S2耗散有效面的确定方法如下:在待分析的结构中,根据结构中不同平板面之间的焊接连接关系,对至少存在两条焊缝的平板面进行统计,并将这些平面作为后续计算过程中的耗散有效面,对其编号i进行标号,i=1,2,3...N;N为耗散有效面的数量。
进一步地,所述S3耗散面有效畸变率的确定方法如下:根据S2中确定的耗散有效面,采用残余应力测试装置对该结构中每个耗散面的残余应力值进行测量,并将第i个耗散有效面的残余应力值记为σri;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;然后计算每个耗散面所对应的残余应力顺位系数Si,
其中,Si为每个耗散面所对应的残余应力顺位系数;σri为第i个耗散有效面的残余应力值,i=1,2,3...N;N为耗散有效面的数量;A为残余应力摄取面的最小长度,取值为0.1mm-1mm;C(σri)max为所有耗散面的残余应力值中的最大值;C(σri)min为所有耗散面的残余应力值中的最小值;
然后计算每个耗散面的顺位畸变率Ji,
其中,Ji为每个耗散面的顺位畸变率;Si为每个耗散面所对应的残余应力顺位系数;i=1,2,3...N;N为耗散有效面的数量;C(Si)max为所有耗散面所对应的残余应力顺位系数中的最大值;C(Si)min为所有耗散面所对应的残余应力顺位系数中的最小值。
进一步地,所述S4耗散面有效固有频率的确定方法如下:根据S2中确定的耗散有效面,将每个单独耗散有效面的平板有限元分析模型导入到有限元分析软件中去,然后对其进行网格划分与边界约束条件的施加,待前处理完毕后,在模态模块中对第三阶固有频率进行求解,并标记为Hi;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;然后将S3中计算的顺位畸变率Ji与有限元分析软件中获得的Hi一起代入下式,求解基于耗散畸变率的耗散面固有频率修正值Ti,
其中,Ti为基于耗散畸变率的耗散面固有频率修正值;Hi为第三阶固有频率;Ji为每个耗散面的顺位畸变率;i=1,2,3...N;N为耗散有效面的数量;C(Hi)max为第三阶固有频率中的最大值;C(Hi)min为第三阶固有频率中的最小值;C(Ji)max为所有耗散面的顺位畸变率中的最大值;C(Ji)min为所有耗散面的顺位畸变率中的最小值。
进一步地,所述S5中最佳振动时效频率计算方法如下:将S1中的结构第三阶固有频率ω与S4中的耗散面固有频率修正值Ti代入下式,求解最佳振动时效频率f,
其中,f为最佳振动时效频率;ω为结构第三阶固有频率;Ti为每个耗散面固有频率修正值;i=1,2,3...N;N为耗散有效面的数量;C(Ti)max为所有耗散面固有频率修正值中的最大值;C(Ti)min为所有耗散面固有频率修正值中的最小值;为所有耗散面固有频率修正值的平均值。
上述技术方案中有限元分析软件优选ANSYS软件,残余应力测试装置优选X射线残余应力测试装置,残余应力摄取面的最小长度在这里取值为0.35mm。
有益效果:本发明可以有效考虑结构件焊接连接对耗散模态的影响,可以实现焊接连接结构振动最佳时效频率,并从基于传递耗散修正的角度出发来对该结构振动时效激振频率进行获取,从而考虑了焊接结构面之间的耗散过程,避免了传统方法单一通过整体结构固有频率来进行激振的局限性,更加高效准确地获取该结构振动时效的最佳激振频率。
附图说明
图1为本发明方法的流程图。
具体实施方式
如图1所示,本实施例一种基于传递耗散修正的结构振动时效频率获取方法,该方法包括以下步骤:
S1.基于有限元手段的结构固有频率的确定:
将建立完毕的结构有限元分析模型导入到有限元分析软件ANSYS中,然后对其进行面积为5cm2的网格划分与边界约束条件的施加,待前处理完毕后,在谐响应模块中对该结构在前处理条件下的第三阶固有频率进行求解,标记为ω;
S2.基于焊接连接关系的耗散有效面的确定:
在待分析的结构中,根据结构中不同平板面之间的焊接连接关系,至少存在两条焊缝的平板面进行统计,并将这些平面作为后续计算过程中的耗散有效面,对其编号i进行标号,i=1,2,3...N;N为耗散有效面的数量;
S3.基于残余应力顺位的耗散面有效畸变率的确定:
根据S2中确定的耗散有效面,采用X射线残余应力测试装置对该结构中每个耗散面的残余应力值进行测量,并将第i个耗散有效面的残余应力值记为σri;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;
然后计算每个耗散面所对应的残余应力顺位系数Si
其中,Si为每个耗散面所对应的残余应力顺位系数;σri为第i个耗散有效面的残余应力值,i=1,2,3...N;N为耗散有效面的数量;A为残余应力摄取面的最小长度,常取值为0.1mm-1mm,这里取值为0.35mm;C(σri)max为所有耗散面的残余应力值中的最大值;C(σri)min为所有耗散面的残余应力值中的最小值;
然后计算每个耗散面的顺位畸变率Ji
其中,Ji为每个耗散面的顺位畸变率;Si为每个耗散面所对应的残余应力顺位系数;i=1,2,3...N;N为耗散有效面的数量;C(Si)max为所有耗散面所对应的残余应力顺位系数中的最大值;C(Si)min为所有耗散面所对应的残余应力顺位系数中的最小值;
S4.基于耗散畸变率的耗散面有效固有频率的确定:
根据S2中确定的耗散有效面,将每个单独耗散有效面的平板有限元分析模型导入到有限元分析软件ANSYS中去,然后对其进行面积为5cm2的网格划分与边界约束条件的施加,待前处理完毕后,在模态模块中对第三阶固有频率进行求解,并标记为Hi;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;
然后将S3中计算的顺位畸变率Ji与有限元分析软件中获得的Hi一起代入下式,求解基于耗散畸变率的耗散面固有频率修正值Ti
其中,Ti为基于耗散畸变率的耗散面固有频率修正值;Hi为第三阶固有频率;Ji为每个耗散面的顺位畸变率;i=1,2,3...N;N为耗散有效面的数量;C(Hi)max为第三阶固有频率中的最大值;C(Hi)min为第三阶固有频率中的最小值;C(Ji)max为所有耗散面的顺位畸变率中的最大值;C(Ji)min为所有耗散面的顺位畸变率中的最小值;
S5.基于传递耗散修正的最佳振动时效频率计算:
将S1中的结构第三阶固有频率ω与S4中的耗散面固有频率修正值Ti代入下式,求解最佳振动时效频率f
Claims (6)
1.一种基于传递耗散修正的结构振动时效频率获取方法,其特征在于:包括如下步骤:
S1、基于有限元手段的结构固有频率的确定;
S2、基于焊接连接关系的耗散有效面的确定;
S3、基于残余应力顺位的耗散面有效畸变率的确定;
S4、基于耗散畸变率的耗散面有效固有频率的确定;
S5、基于传递耗散修正的最佳振动时效频率计算。
2.根据权利要求1所述的基于传递耗散修正的结构振动时效频率获取方法,其特征在于:所述S1结构固有频率的确定方法如下:将建立完毕的结构有限元分析模型导入到有限元分析软件中,然后对其进行网格划分与边界约束条件的施加,待前处理完毕后,在谐响应模块中对该结构在前处理条件下的第三阶固有频率进行求解,标记为ω。
3.根据权利要求2所述的基于传递耗散修正的结构振动时效频率获取方法,其特征在于:所述S2耗散有效面的确定方法如下:在待分析的结构中,根据结构中不同平板面之间的焊接连接关系,对至少存在两条焊缝的平板面进行统计,并将这些平面作为后续计算过程中的耗散有效面,对其编号i进行标号,i=1,2,3...N;N为耗散有效面的数量。
4.根据权利要求3所述的基于传递耗散修正的结构振动时效频率获取方法,其特征在于:所述S3耗散面有效畸变率的确定方法如下:根据S2中确定的耗散有效面,采用残余应力测试装置对该结构中每个耗散面的残余应力值进行测量,并将第i个耗散有效面的残余应力值记为σri;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;然后计算每个耗散面所对应的残余应力顺位系数Si,
其中,Si为每个耗散面所对应的残余应力顺位系数;σri为第i个耗散有效面的残余应力值,i=1,2,3...N;N为耗散有效面的数量;A为残余应力摄取面的最小长度,取值为0.1mm-1mm;C(σri)max为所有耗散面的残余应力值中的最大值;C(σri)min为所有耗散面的残余应力值中的最小值;
然后计算每个耗散面的顺位畸变率Ji,
其中,Ji为每个耗散面的顺位畸变率;Si为每个耗散面所对应的残余应力顺位系数;i=1,2,3...N;N为耗散有效面的数量;C(Si)max为所有耗散面所对应的残余应力顺位系数中的最大值;C(Si)min为所有耗散面所对应的残余应力顺位系数中的最小值。
5.根据权利要求4所述的基于传递耗散修正的结构振动时效频率获取方法,其特征在于:所述S4耗散面有效固有频率的确定方法如下:根据S2中确定的耗散有效面,将每个单独耗散有效面的平板有限元分析模型导入到有限元分析软件中去,然后对其进行网格划分与边界约束条件的施加,待前处理完毕后,在模态模块中对第三阶固有频率进行求解,并标记为Hi;其中i为耗散面编号,i=1,2,3...N;N为耗散有效面的数量;然后将S3中计算的顺位畸变率Ji与有限元分析软件中获得的Hi一起代入下式,求解基于耗散畸变率的耗散面固有频率修正值Ti,
其中,Ti为基于耗散畸变率的耗散面固有频率修正值;Hi为第三阶固有频率;Ji为每个耗散面的顺位畸变率;i=1,2,3...N;N为耗散有效面的数量;C(Hi)max为第三阶固有频率中的最大值;C(Hi)min为第三阶固有频率中的最小值;C(Ji)max为所有耗散面的顺位畸变率中的最大值;C(Ji)min为所有耗散面的顺位畸变率中的最小值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010122549.3A CN111460700B (zh) | 2020-02-27 | 2020-02-27 | 一种基于传递耗散修正的结构振动时效频率获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010122549.3A CN111460700B (zh) | 2020-02-27 | 2020-02-27 | 一种基于传递耗散修正的结构振动时效频率获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460700A true CN111460700A (zh) | 2020-07-28 |
CN111460700B CN111460700B (zh) | 2023-05-23 |
Family
ID=71684176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010122549.3A Active CN111460700B (zh) | 2020-02-27 | 2020-02-27 | 一种基于传递耗散修正的结构振动时效频率获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460700B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177337A (zh) * | 2021-04-20 | 2021-07-27 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050242066A1 (en) * | 2004-04-29 | 2005-11-03 | Uit. L.L.C. Company | Method for modifying or producing materials and joints with specific properties by generating and applying adaptive impulses a normalizing energy thereof and pauses therebetween |
CN103464536A (zh) * | 2013-08-12 | 2013-12-25 | 西北工业大学 | 一种弹性变形条件下的振动时效成形方法及装置 |
CN105543469A (zh) * | 2015-12-25 | 2016-05-04 | 常州大学 | 一种确定振动时效激振频率的***及方法 |
CN105825029A (zh) * | 2016-04-01 | 2016-08-03 | 上海海事大学 | 一种用于优化设计高频振动能量放大装置的方法 |
CN109182728A (zh) * | 2018-09-10 | 2019-01-11 | 上海海事大学 | 一种绿色智能振动时效***及方法 |
-
2020
- 2020-02-27 CN CN202010122549.3A patent/CN111460700B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050242066A1 (en) * | 2004-04-29 | 2005-11-03 | Uit. L.L.C. Company | Method for modifying or producing materials and joints with specific properties by generating and applying adaptive impulses a normalizing energy thereof and pauses therebetween |
CN1985013A (zh) * | 2004-04-29 | 2007-06-20 | U.I.T.有限责任公司 | 通过产生和施加自适应脉冲、其标准化的能量和其之间的间歇来改变或生产具有特定性质的材料和接头的方法 |
CN103464536A (zh) * | 2013-08-12 | 2013-12-25 | 西北工业大学 | 一种弹性变形条件下的振动时效成形方法及装置 |
CN105543469A (zh) * | 2015-12-25 | 2016-05-04 | 常州大学 | 一种确定振动时效激振频率的***及方法 |
CN105825029A (zh) * | 2016-04-01 | 2016-08-03 | 上海海事大学 | 一种用于优化设计高频振动能量放大装置的方法 |
CN109182728A (zh) * | 2018-09-10 | 2019-01-11 | 上海海事大学 | 一种绿色智能振动时效***及方法 |
Non-Patent Citations (2)
Title |
---|
ZHIHUA WU 等: "Nonlinear primary and super-harmonic resonances of functionally graded carbon nanotube reinforced composite beams" * |
李俊明 等: "线性弹簧组合式振动时效装置的 3 次超谐共振特性分析" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177337A (zh) * | 2021-04-20 | 2021-07-27 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
CN113177337B (zh) * | 2021-04-20 | 2023-05-26 | 扬州大学 | 基于关联因素特征值波动区间的芦苇收割机安全评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111460700B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460700A (zh) | 一种基于传递耗散修正的结构振动时效频率获取方法 | |
CN108225547B (zh) | 汽车排气噪声测量方法 | |
WO2017063540A1 (zh) | 焊趾处热点应力的计算方法 | |
WO2020252842A1 (zh) | 基于实际载荷的齿面接触斑点齿轮的设计方法 | |
CN108362561B (zh) | 一种焊缝及焊接热影响区的材料力学性能确定方法 | |
CN115859746B (zh) | 一种基于深度学习的铜材料增材制造成形精度控制方法 | |
CN112199875B (zh) | 一种基于雨流法的元器件焊点随机振动疲劳寿命分布预测方法 | |
CN111539142B (zh) | 一种管道裂纹扩展驱动力计算方法及*** | |
CN112182807A (zh) | 基于有限元仿真的直齿轮齿根裂纹扩展规律分析方法 | |
CN108731574B (zh) | 一种船舶发电机组垫块数据测量装置及其测量方法 | |
CN1866262A (zh) | 一种电容等效模型的建模装置和方法 | |
CN103808262B (zh) | 多孔位产品孔位的模拟修模方法 | |
CN108387350A (zh) | 车身接附点动刚度测试方法、装置和存储介质 | |
Park et al. | Experimental approach for estimating mesh stiffness in faulty states of rotating gear | |
CN107066728B (zh) | 一种钛合金潜水器耐压球壳极限承载力估算方法 | |
CN113155651A (zh) | 基于疲劳试验与仿真的直齿轮齿根裂纹扩展规律分析方法 | |
CN102955363B (zh) | 光学临近效应修正在线监控的方法 | |
CN112307575B (zh) | 一种评判装配结果的复合材料壁板零件合格性的检查方法 | |
CN111157338B (zh) | 一种金属哑铃型试样大应变范围压缩硬化曲线的测量方法 | |
CN110333148B (zh) | 一种基于振动衰减曲线精细化分析的土动剪切模量测试方法 | |
CN114216674B (zh) | 一种齿轮疲劳试验件及制作方法 | |
CN111289392A (zh) | 基于局部抛光的金属薄壁管疲劳试验装置及方法 | |
KR20090124355A (ko) | 하이포이드 단조기어용 금형 제작방법 | |
CN111241741B (zh) | 一种基于残余应力影响区间修正的起重机轮轨磨损监测方法 | |
CN103981359B (zh) | 一种频谱谐波时效技术激振力的选配方法 |
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 |