CN102063539B - 一种基于有限元的惯性平台残余应力释放仿真方法 - Google Patents
一种基于有限元的惯性平台残余应力释放仿真方法 Download PDFInfo
- Publication number
- CN102063539B CN102063539B CN 201010614317 CN201010614317A CN102063539B CN 102063539 B CN102063539 B CN 102063539B CN 201010614317 CN201010614317 CN 201010614317 CN 201010614317 A CN201010614317 A CN 201010614317A CN 102063539 B CN102063539 B CN 102063539B
- Authority
- CN
- China
- Prior art keywords
- inertial platform
- analysis
- air column
- finite element
- platform
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于有限元的惯性平台残余应力释放仿真方法,步骤如下:1.惯性平台实体建模,应用实体绘图软件完整地得到了惯性平台实体模型,然后将其导入到有限元仿真软件;2.空气柱建模,将惯性平台外围空气建立成一个空气柱模型;3.惯性平台和空气柱材料属性赋值,按有限元仿真软件要求输入惯性平台和空气柱材料属性;4.惯性平台和空气柱有限元网格建模,包括单元类型选取和网格划分两步;5.惯性平台热-结构耦合分析,包括瞬态热分析和结构分析两步:瞬态热分析是在空气柱外表面加载温度循环应力,为惯性平台残余应力释放仿真分析提供热条件;结构分析是指将惯性平台瞬态热分析结果加载到热-结构耦合分析中,得到惯性平台残余应力随温度循环应力的释放规律。
Description
技术领域
本发明提供一种基于有限元的惯性平台残余应力释放仿真方法,它涉及一种基于有限元热-结构耦合分析的惯性平台残余应力释放仿真方法,属于有限元仿真技术领域。
背景技术
随着国防任务的日益加重和战备要求的日益提高,惯性平台系统参数稳定性差问题逐渐暴露出来,严重制约了武器作战效能的发挥,成为了亟待解决的军事需求。长期的工程实践和理论研究表明,由于惯性平台结构设计本身的不尽完善(如质量不平衡等)及工艺过程的分散性(如机械加工力不均匀等),使得在生产工艺过程中惯性平台内部存在着大量的残余应力,在惯性平台系统贮存过程中随着外界环境的变化,平台内部的残余应力得到释放,使得惯性平台尺寸精度的稳定性发生变化,从而影响惯性平台内陀螺仪、加速度计的安装精度稳定性水平。
残余应力对结构尺寸稳定性的影响目前通常用精密弹性极限(Precision Elastic Limit)或微屈服强度(Microyield Strength)来表征,并在此基础上采用X射线衍射仪(XRD)测量负载应力下材料残余塑性应变的残余应力值。然而,基于精密弹性极限和微屈服强度的残余应力试验测量方法只能得到有限个点的残余应力值,不能得到负载应力下惯性平台内残余应力分布规律,成为亟待解决的问题。
发明内容
本发明的目的是为了提供一种基于有限元的惯性平台残余应力释放仿真方法,弥补现有技术不能全面分析惯性平台残余应力释放规律和不能为消除惯性平台残余应力提供有效技术支撑的缺陷。本发明涉及的一种基于有限元的惯性平台残余应力释放仿真方法能够模拟在温度循环应力作用下惯性平台内残余应力分布变化规律,能为惯性平台尺寸精度的稳定性研究提供理论基础和技术支撑。
本发明是通过以下技术方案实现的,首先应用Unigraphics(简称UG)软件绘制了惯性平台复杂结构实体模型,并使用Parasolid输入格式将其导入到ANSYS有限元软件中,然后输入其材料属性和对其进行网格划分后得到惯性平台有限元网格模型,最后采用ANSYS热-结构耦合分析方法进行惯性平台残余应力释放仿真分析。本发明通过在惯性平台周围建立空气柱,真实模拟惯性平台受热、传热情况,其在热-结构耦合分析中首先进行惯性平台瞬态热分析,得到在温度循环应力作用下惯性平台热响应结果;然后将工程上采用X射线衍射方法测得的惯性平台初始残余应力值以初应力载荷的形式加载到惯性平台有限元网格模型各节点的六个自由度方向,在此基础上继续进行结构分析,完成惯性平台残余应力释放热-结构耦合分析,实现了基于有限元的惯性平台残余应力释放仿真方法,其操作流程及步骤如图1所示。
本发明一种基于有限元的惯性平台残余应力释放仿真方法,其步骤如下:
步骤一:惯性平台实体建模
首先应用实体绘图软件(Unigraphics,简称UG)进行惯性平台实体建模,根据惯性平台实际生产加工尺寸在UG软件中完整地得到了惯性平台实体模型,然后通过UG软件中的实体模型输出格式(Parasolid格式)将惯性平台UG实体模型导入到有限元仿真软件(ANSYS)中,最后在ANSYS软件中对实体模型进行拓扑修改和几何修复,得到惯性平台实体模型,如图2所示;
步骤二:空气柱建模
根据惯性平台在惯导系统中的安装空间,将惯性平台外围空气建立成一个空气柱模型,在ANSYS软件空气柱模型的操作处理中将惯性平台与空气柱重叠的部分挖掉(Overlap操作),且在热-结构耦合分析中在空气柱外表面加载温度循环剖面,模拟惯性平台温度应力实际加载情况,如图3所示;
步骤三:惯性平台和空气柱材料属性赋值
按有限元仿真软件(ANSYS)中要求输入惯性平台和空气柱材料属性,该材料属性是热-结构耦合分析所需的密度、比热、热导率、弹性模量、泊松比和热膨胀系数;
步骤四:惯性平台和空气柱有限元网格建模
惯性平台和空气柱有限元网格建模包括单元类型选取和网格划分。
该单元类型选取要充分考虑热-结构耦合分析中瞬态热分析和结构分析单元类型转换的兼容性,合理选取瞬态热分析单元类型;
该网格划分要兼顾计算量和计算精度的要求,合理控制网格划分密度,在本发明中采用精度等级为6自由网格划分,同时对关键部位进行网格细化,合理实现惯性平台和空气柱网格划分的精度分布。惯性平台和空气柱有限元网格模型如图4和图5所示;
步骤五:惯性平台热-结构耦合分析
本发明所述的惯性平台热-结构耦合分析是一种热-结构间接顺序耦合分析方法,用来分析温度场改变对惯性平台残余应力释放的交叉耦合影响。它是基于有限元的惯性平台残余应力释放仿真方法的关键技术,包括瞬态热分析和结构分析两步:
瞬态热分析是在惯性平台有限元网格模型的基础上,在空气柱外表面加载温度循环应力,并采用ANSYS参数化设计语言(APDL语言)编制温度循环加载命令流,将其加载到惯性平台有限元网格模型中,为惯性平台残余应力释放仿真分析提供热条件;
结构分析是指将惯性平台瞬态热分析结果加载到结构分析中,同时将惯性平台残余应力测量值以初始应力载荷的方式加载到惯性平台结构分析中,在每个有限元节点的x,y,z,xy,xz,yz坐标方向施加相同的初始应力载荷,进行热-结构耦合分析,最后得到惯性平台残余应力随温度循环应力的释放规律,实现了本发明所提一种基于有限元的惯性平台残余应力释放仿真方法。惯性平台热-结构耦合仿真分析流程如图6所示。
其中,在步骤一中所述的应用实体绘图软件(Unigraphics,简称UG)进行惯性平台实体建模,其步骤包括:1.绘制惯性平台结构图;2.绘制陀螺仪安装孔;3.绘制加速度计安装孔;4.绘制其他安装孔和面;所述的在ANSYS软件中对实体模型进行适当地拓扑修改和几何修复是指由于UG软件和ANSYS软件接口的缺陷,使得导入到ANSYS软件中的惯性平台实体模型会有差错(比如本发明所用Parasolid格式导入ANSYS的惯性平台模型不是实体模型,而是线框图),这时就需要应用ANSYS软件对模型进行适当修改和修复。
其中,在步骤二中所述的将惯性平台外围空气建成一个空气柱模型,该空气柱模型的大小为直径145mm~155mm,高130mm~140mm。
本发明一种基于有限元的惯性平台残余应力释放仿真方法,具有以下优点:
1.应用ANSYS仿真软件与Unigraphics绘图软件的兼容性准确绘制惯性平台实体模型。本发明首先在Unigraphics软件中精确绘制惯性平台实体模型,然后使用Parasolid输入方法将惯性平台Unigraphics实体模型导入到ANSYS中,再对其进行适当的拓扑修改和几何修复,得到惯性平台ANSYS实体模型。
2.采用初应力载荷方式将残余应力引入到惯性平台有限元仿真中。惯性平台残余应力释放是指经过生产工艺过程后惯性平台在温度循环条件下残余应力逐步释放的过程。本发明通过ANSYS热-结构耦合分析,将残余应力作为一种初始内应力载荷施加到惯性平台有限元节点的各个自由度方向,实现残余应力在惯性平台有限元仿真模型中的施加。
3.通过ANSYS热-结构耦合分析实现温度循环条件下惯性平台残余应力释放仿真方法。如何在ANSYS软件中实现残余应力的释放机理一直是惯性平台残余应力有限元仿真的难点,本发明利用温度循环条件下残余应力释放时晶粒间相互作用原理,通过热-结构耦合分析方法有效地实现了在温度循环作用下有限元节点中初始应力载荷作用规律,且在热-结构交叉耦合影响下体现了残余应力的不可逆释放过程。
附图说明
图1是本发明方法流程图。
图2是惯性平台实体模型图。
图3是空气柱模型图。
图4是惯性平台有限元模型图。
图5是空气柱有限元模型图。
图6是本发明实施惯性平台热-结构仿真分析流程图。
图7是惯性平台热瞬态分析加载的温度剖面。
具体实施方式
下面将结合附图和具体实例对本发明做进一步的详细说明。
根据附图1所示实施流程,从惯性平台实体建模,空气柱建模、惯性平台和空气柱材料属性赋值、惯性平台和空气柱有限元网格建模、惯性平台热-结构耦合分析五个方面对本发明进行详细阐述。惯性平台是一个机电控制系统,其机械结构部分是它极其重要的组成部分,用于提供安装加速度计和陀螺仪的平台,惯性平台机械结构尺寸如图2所示。
本发明一种基于有限元的惯性平台残余应力释放仿真方法,其步骤如下:
步骤一:惯性平台实体建模
本发明首先应用Unigraphics软件进行惯性平台实体建模,根据惯性平台实际生产加工尺寸在Unigraphics软件中完整地得到惯性平台实体模型。其建模过程共分为6个步骤,共有5个重要装配定位面(3个加速度计和2个陀螺仪的装配面)需要精确绘制,具体步骤如下:
1)绘制惯性平台结构框图
本发明首先在Unigraphics软件中绘制惯性平台结构框图,其中惯性平台外径83mm,高17mm,其上安装有陀螺仪和加速度计,惯性平台结构框图如图2所示。
2)绘制陀螺仪安装孔
惯性平台上装有两个动力调谐陀螺仪,陀螺仪安装圆孔直径为47mm。在绘制陀螺仪安装孔时其中一个陀螺的安装孔与惯性平台垂直,另一个陀螺安装孔与惯性平台平行,且两个陀螺仪安装孔相互垂直。陀螺仪安装孔精度是指陀螺仪在惯性平台上安装位置的精确绘制和安装孔定位公差精度保证。
3)绘制加速度计安装孔
惯性平台上安装有三块石英挠性加速度计,加速度计安装孔直径为23.5mm,在绘制加速度计安装孔时将加速度计分别与惯性平台上两陀螺的安装孔同轴。加速度计安装孔精度是指加速度计在惯性平台上安装位置的精确绘制和安装定位公差精度保证。
4)绘制其他安装孔和面
由于功能的需要和装配空间的限制,其他安装定位孔和面分散在惯性平台的各个部位,所以其他安装孔和面的绘制需要统筹考虑,尽可能优先考虑复杂结构及精度要求极高的安装孔和面绘制。
5)将惯性平台UG实体模型导入ANSYS中
在本发明中惯性平台UG实体模型导入ANSYS分为两步:将惯性平台实体模型从UG软件中导出和将从UG软件中导出的惯性平台实体模型导入到ANSYS软件中。
将惯性平台UG实体模型导出时要求在同一系统下同时安装Unigraphics和ANSYS两种软件,并设置相应环境变量。在本发明中使用Parasolid输入方法将惯性平台Unigraphics实体模型导入到ANSYS中,其在Unigraphics软件中的操作路径为:文件|导出|Parasolid,然后选择需要导出的惯性平台UG实体就可将惯性平台实体模型从UG软件中导出。
将从UG软件中导出的惯性平台实体模型导入到ANSYS软件是指将从UG软件导出的格式为Parasolid格式的惯性平台实体模型导入到ANSYS软件中,其在ANSYS中的操作路径为:File|Import|Para,在弹出的对话框中选择输入惯性平台Parasolid模型的路径和文件名,单击确认按钮将可将惯性平台实体模型导入到ANSYS软件中。
6)修复惯性平台ANSYS实体模型
UG软件和ANSYS软件模型的转化是通过接口设计实现的,在本发明中采用Parasolid方法将惯性平台UG模型导入ANSYS中后是线模型,而非实体模型,所以还需要进行实体模型转化,转化方法在ANSYS中的操作路径为:PlotCtrls|Style|Solid Model Facets,在弹出的对话框下拉选项中选择Fine选项就可完成惯性平台AYSYS实体模型的修复。
步骤二:空气柱建模
由于ANSYS有限元仿真软件中的数值没有单位,因此在空气柱建模之前需要定义单位制,在本发明中采用Meter-Kilogram-Second(简称:MKS)单位制,各参数单位规定如表1中所示。
表1基于有限元的惯性平台残余应力释放仿真方法单位制
参数 | 长度 | 力 | 时间 | 质量 | 压力 | 密度 | 应力 | 杨氏模量 | 导热系数 | 比热 | 对流系数 |
MKS | m | N | s | kg | Pa | kg/(m)3 | Pa | Pa | W/m*K | J/(kg)(K) | W/(m)2*K |
本发明在惯性平台外表面建成一个直径为150mm,高为135mm的空气柱,其在ANSYS中的操作路径为:Main|Preprocessor|Modeling|Create|Volumes|Cylinder|By Dimensions,在弹出的对话框中输入空气柱尺寸就可以完成空气柱的建模。在空气柱的模型处理中将惯性平台与空气柱重叠的部分挖掉,其在ANSYS中的操作路径为:Main|Preprocessor|Modeling|Operate|Booleans|Overlap|Volumes,在弹出的对话框中选择惯性平台实体和空气柱实体,将其重叠部分挖掉。热-结构耦合分析中在空气柱外表面加载温度循环剖面,模拟惯性平台温度应力实际加载情况,空气柱实体模型如图3所示。
步骤三:惯性平台和空气柱材料属性赋值
根据本发明一种基于有限元的惯性平台残余应力释放仿真方法需要用到的惯性平台和空气柱材料属性,依次输入热-结构耦合分析所需的密度、比热、热导率、弹性模量、泊松比和热膨胀系数,如表2和表3所示。
表2热-结构耦合分析惯性平台材料属性
表3热-结构耦合分析空气柱材料属性
惯性平台和空气柱材料属性赋值在ANSYS中的操作路径为:Main|Preprocessor|Material Props|Material Models,根据弹出的对话框选项依次输入表2和表3中各材料参数。
步骤四:惯性平台和空气柱有限元网格建模
惯性平台和空气柱有限元网格建模包括单元类型选取和网格划分两步。
单元类型选取主要是考虑到热-结构耦合分析中瞬态热分析和结构分析单元类型转换的兼容性,本发明选取SOLID90单元类型进行惯性平台瞬态热分析,SOLID90是三维八节点热单元的高阶形式,该单元有二十个节点,每个节点只有一个温度自由度,二十节点的单元有协调的温度形函数,尤其适用于模拟曲边,适用于基于有限元的惯性平台残余应力释放仿真方法。单元类型选取在ANSYS中的操作路径为:Main|Preprocessor|Element Type|Add/Edit/Delete,在弹出的对话框中选择添加选项后弹出单元类型选取对话框,在ThermalSolid单元类型库中选择SOLID90,实现单元类型的选取。
网格划分要兼顾计算量和计算精度的要求,合理控制网格划分密度,在本发明中采用精度等级为6自由网格划分,同时对关键部位进行网格细化,合理实现惯性平台和空气柱网格划分的精度分布。惯性平台和空气柱有限元网格划分在ANSYS中的操作路径为:Main|Preprocessor|Meshing|MeshTool,在弹出的对话框中选择Volumes单元属性,Smart Size为6,Mesh Method为Free,且在单元节点处Refine,从而实现单元网格的划分。
按照以上方法对惯性平台和空气柱进行有限元网格建模,得到惯性平台有限元网格模型共有节点9077个,单元10002个,空气柱有限元网格模型共有节点71427个,单元57359个。建立的惯性平台和空气柱有限元网格模型如图4及图5所示。
步骤五:惯性平台热-结构耦合分析
惯性平台热-结构耦合分析是一种间接法顺序耦合分析。其主要分三步完成:1.进行瞬态热分析,求得结构的温度场;2.将模型中的单元类型转变为对应的结构分析单元,并将第一步求得的瞬态热分析结果当作体载荷施加到节点上;3.定义其余结构分析需要的选项,并进行结构分析。
1)瞬态热分析
本发明一种基于有限元的惯性平台残余应力释放仿真方法中惯性平台瞬态热分析进行加载求解时需要完成如下的工作。包括定义分析类型、定义初始条件、施加载荷、求解。
a)定义分析类型
指定分析类型为瞬态分析,开始惯性平台热-结构耦合分析。
ANSYS操作:Main Menu|Preferences|Thermal
b)定义初始条件
惯性平台起始温度是常温且均匀,在本发明中设定所有节点初始温度为25℃。
ANSYS操作:Main Menu|Solution|-Loads-|Settings|Uniform Temp
c)施加载荷
本发明在惯性平台有限元模型的基础上,在空气柱外表面加载温度循环应力,其中温度循环条件来自于实际生产过程的温度循环条件,如图7所示。
并采用ANSYS参数化设计语言(APDL语言)编制温度循环加载命令流,将其加载到惯性平台有限元模型上,进行惯性平台瞬态热分析,其中温度循环APDL命令流如附图所示。
本发明一种基于有限元的惯性平台残余应力释放仿真方法,瞬态热分析温度循环剖面加载APDL命令流如下:
/prep7
asel,s,,,59,60
asel,a,,,61,62
sfa,all,,conv,13.8,tempi!对流换热面
allsel
!求解
/solu!瞬态热分析
ANTYPE,TRANS
AUTOTS,ON
!加载
!循环1
!升温
tempi=20!设置初始温度
tunif,tempi
allsel
outres,all,all
nsubst,30!设置子步
asel,s,,,59,60
asel,a,,,61,62
da,all,temp,50
allsel
kbc,0
time,180
solve
!保温
kbc,1
time,180+4*3600
solve
!降温
kbc,0
time,180+4*3600+10*60
asel,s,,,59,60
asel,a,,,61,62
da,all,temp,-50
solve
!保温
kbc,1
time,29580
solve
!第二个循环
!升温
kbc,0
time,30180
asel,s,,,59,60
asel,a,,,61,62
da,all,temp,50
solve
!保温
kbc,1
time,44580
solve
!降温
kbc,0
time,45180
asel,s,,,59,60
asel,a,,,61,62
da,all,temp,-50
solve
!保温
kbc,1
time,59580
solve
!升温
kbc,0
time,60000
asel,s,,,59,60
asel,a,,,61,62
da,all,temp,20
solve
将上述命令流加载到空气柱外表面上,进行惯性平台瞬态热分析,得到惯性平台热瞬态分析结果。
2)单元类型转换
将惯性平台有限元分析模型中的单元类型转变为对应的结构分析单元类型,并将瞬态热分析中求得的热分析结构当作体载荷施加到节点上。首先进入前处理器,并指定新的分析范畴为结构分析,其在ANSYS中的操作路径为:Main Menu|Preference,在弹出的对话框中选择“Strutural”选项,使所有菜单变为结构分析的选项。然后进行单元转换,其在ANSYS中的操作路径为:Main Menu|Preprocessor|Element Type|Switch Elem Type,将弹出SwithchElem Type(转换单元类型)对话框,本发明在对话框中的Change element type(改变单元类型)下拉框中选择“Thermal to Struc”,然后单击关闭对话框,ANSYS程序将会自动将模型中的热单元转换为对应的结构单元类型。
定义材料的性能参数。根据表2和表3中惯性平台和空气柱材料属性参数,定义进行结构分析需要的材料弹性模量、泊松比、密度和热膨胀系数。
将瞬态热分析得到的温度结果施加到惯性平台结构分析模型上,其在ANSYS中的操作路径为:Main Menu|Solution|Define Loads|Apply|Structural|Temperature|From ThermAnaly,将弹出Apply TEMP from Themal Analysis(从已进行的热分析结果中施加温度载荷)对话框,单击对话框中的Browse按钮,选择前面瞬态热分析的结果文件,作为结构分析的热载荷加到节点上。
3)结构分析
在上述瞬态热分析和单元类型转换的基础上开展结构分析工作。包括残余应力施加和求解。
残余应力施加是指将初始残余应力测量结果作为惯性平台初始应力载荷加载到瞬态热分析结果模型中,在惯性平台每个有限元节点的x,y,z,xy,yz,xz方向施加相同的初始应力载荷,进行结构分析,其在ANSYS中的操作路径为:Main Menu|Solution|Load Step Opts|Initial Stress|Apply Const Strs,将弹出惯性平台残余应力施加对话框,在对话框中输入残余应力初始载荷值为70.45Mpa,单击OK按钮键,至此,完成了惯性平台有限元模型中残余应力的施加。
求解是指对惯性平台残余应力释放进行结构分析,其在ANSYS中的操作路径为:MainMenu|Solution|Solve|Current LS,将弹出惯性平台结构分析对话框,单击OK按钮进行结构分析求解,至此,完成了惯性平台残余应力释放仿真分析全过程。
本发明一种基于有限元的惯性平台残余应力释放仿真方法,得到惯性平台最大残余应力随温度循环应力的释放规律如表4中所示:
表4基于有限元的惯性平台残余应力释放仿真方法结果
循环次数 | 位移变形最大值(mm) | 应力最大值(Mpa) |
初始时刻 | 0.0364 | 137 |
第一次循环后 | 0.0273 | 102 |
第二次循环后 | 0.0177 | 63.6 |
Claims (4)
1.一种基于有限元的惯性平台残余应力释放仿真方法,其特征在于下:该方法的步骤如下:
步骤一:惯性平台实体建模
首先应用实体绘图软件即Unigraphics软件进行惯性平台实体建模,根据惯性平台实际生产加工尺寸在Unigraphics软件中完整地得到了惯性平台UG实体模型,然后通过Unigraphics软件中的实体模型输出Parasolid格式将惯性平台UG实体模型导入到有限元仿真软件即ANSYS软件中,最后在ANSYS软件中对惯性平台UG实体模型进行拓扑修改和几何修复,得到惯性平台实体模型;
步骤二:空气柱建模
根据惯性平台在惯导系统中的安装空间,将惯性平台外围空气建立成一个空气柱模型,在ANSYS软件空气柱模型的操作处理中将惯性平台与空气柱重叠的部分挖掉即Overlap操作,且在热-结构耦合分析中在空气柱外表面加载温度循环剖面,模拟惯性平台温度应力实际加载情况;
步骤三:惯性平台和空气柱材料属性赋值
按ANSYS软件中要求输入惯性平台和空气柱材料属性,该材料属性是热-结构耦合分析所需的密度、比热、热导率、弹性模量、泊松比和热膨胀系数;
步骤四:惯性平台和空气柱有限元网格建模
惯性平台和空气柱有限元网格建模包括单元类型选取和网格划分;
该单元类型选取要充分考虑热-结构耦合分析中瞬态热分析和结构分析单元类型转换的兼容性,合理选取瞬态热分析单元类型;
该网格划分要兼顾计算量和计算精度的要求,合理控制网格划分密度,采用精度等级为6自由网格划分,同时对关键部位进行网格细化,合理实现惯性平台和空气柱网格划分的精度分布;
步骤五:惯性平台热-结构耦合分析
该惯性平台热-结构耦合分析包括瞬态热分析和结构分析两步:
瞬态热分析是在惯性平台有限元网格模型的基础上,在空气柱外表面加载温度循环应力,并采用ANSYS参数化设计语言即APDL语言编制温度循环加载命令流,将其加载到惯性平台有限元网格模型中,为惯性平台残余应力释放仿真分析提供热条件;
结构分析是指将惯性平台瞬态热分析结果加载到结构分析中,同时将惯性平台残余应力测量值以初始应力载荷的方式加载到惯性平台结构分析中,在每个有限元节点的x,y,z,xy,xz,yz坐标方向施加相同的初始应力载荷,进行热-结构耦合分析,最后得到惯性平台残余应力随温度循环应力的释放规律。
2.根据权利要求1所示的一种基于有限元的惯性平台残余应力释放仿真方法,其特征在于下:在步骤一中所述的应用Unigraphics软件进行惯性平台UG实体建模,其步骤包括:1)绘制惯性平台结构图;2).绘制陀螺仪安装孔;3).绘制加速度计安装孔;4).绘制其他安装孔和面。
3.根据权利要求1所示的一种基于有限元的惯性平台残余应力释放仿真方法,其特征在于下:在步骤一中所述的在ANSYS软件中对惯性平台UG实体模型进行适当地拓扑修改和几何修复,是指由于Unigraphics软件和ANSYS软件接口的缺陷,使得导入到ANSYS软件中的惯性平台UG实体模型会有差错,这时就需要应用ANSYS软件对模型进行适当修改和修复。
4.根据权利要求1所示的一种基于有限元的惯性平台残余应力释放仿真方法,其特征在于下:在步骤二中所述的将惯性平台外围空气建成一个空气柱模型,该空气柱模型的大小为直径145mm~155mm,高130mm~140mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010614317 CN102063539B (zh) | 2010-12-30 | 2010-12-30 | 一种基于有限元的惯性平台残余应力释放仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010614317 CN102063539B (zh) | 2010-12-30 | 2010-12-30 | 一种基于有限元的惯性平台残余应力释放仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102063539A CN102063539A (zh) | 2011-05-18 |
CN102063539B true CN102063539B (zh) | 2013-03-27 |
Family
ID=43998813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010614317 Expired - Fee Related CN102063539B (zh) | 2010-12-30 | 2010-12-30 | 一种基于有限元的惯性平台残余应力释放仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102063539B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542117B (zh) * | 2012-01-06 | 2013-07-17 | 北京航空航天大学 | 一种基于多场耦合的加速度计振动时效仿真方法 |
CN104200034B (zh) * | 2014-09-10 | 2017-05-24 | 清华大学 | 一种干滑动摩擦热‑应力‑磨损顺序耦合的模拟方法 |
CN107709969B (zh) * | 2016-02-24 | 2020-08-28 | 中国建材检验认证集团股份有限公司 | 涂层残余应力测试方法及仪器 |
CN109408175B (zh) * | 2018-09-28 | 2021-07-27 | 北京赛博贝斯数据科技有限责任公司 | 通用高性能深度学习计算引擎中的实时交互方法及系统 |
CN109325301A (zh) * | 2018-09-30 | 2019-02-12 | 上海机电工程研究所 | 武器装备体系效能快速优化方法 |
CN109684684A (zh) * | 2018-12-07 | 2019-04-26 | 上海机电工程研究所 | 武器装备体系效能快速优化方法、系统及介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100561490C (zh) * | 2008-03-13 | 2009-11-18 | 上海交通大学 | 基于汽车碰撞仿真的拼焊板焊缝有限元建模方法 |
US8200464B2 (en) * | 2009-04-01 | 2012-06-12 | Livermore Software Technology Corporation | Explosion simulation in finite element analysis |
-
2010
- 2010-12-30 CN CN 201010614317 patent/CN102063539B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102063539A (zh) | 2011-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102063539B (zh) | 一种基于有限元的惯性平台残余应力释放仿真方法 | |
Blom | A monolithical fluid-structure interaction algorithm applied to the piston problem | |
CN108009336B (zh) | 一种微桁架结构承载和热防护结构的多目标优化方法 | |
CN104850697B (zh) | 基于ansys和adams的大型天线动力学建模方法 | |
Yan et al. | Topology optimization of damping layers in shell structures subject to impact loads for minimum residual vibration | |
Kersken et al. | Time-Linearized and Time-Accurate 3D RANS Methods for Aeroelastic Analysis in Turbomachinery | |
Jia et al. | Accuracy, efficiency and scalability of explicit and implicit FR/CPR schemes in large eddy simulation | |
CN105787199A (zh) | 一种叶轮机转子叶片从热态到冷态的叶型转换方法 | |
Guo et al. | A CFD/CSD model for aeroelastic calculations of large-scale wind turbines | |
Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
JP6861786B2 (ja) | 物理的性質を強化する、面上設計パターンのローカル制御 | |
CN109271655B (zh) | 一种基于非对称有限元算法的材料尺度效应分析方法 | |
WO2017200525A1 (en) | System and method for finite element analysis of parts having variable spatial density graded regions produced via 3d printers | |
CN113094830A (zh) | 一种固体火箭发动机线性燃烧稳定性预估系统及方法 | |
CN116401778A (zh) | 一种环形调谐液柱阻尼器结构设计方法、装置和设备 | |
CN106971021A (zh) | 一种燃气增压差动贮箱集成设计方法 | |
Girodroux-Lavigne et al. | Comparison of static and dynamic fluid-structure interaction solutions in the case of a highly flexible modern transport aircraft wing | |
Jie et al. | The modal analysis of automotive exhaust muffler based on PRO/E and ANSYS | |
Alavilli et al. | The development of an advanced solid-rocket flow simulation program ROCFLO | |
CN104847414B (zh) | 一种涡旋式流体机械结构化动网格建模方法 | |
Başkut et al. | Development of a closely coupled procedure for dynamic aeroelastic analyses | |
Aguirre et al. | An analytical shock-loss model fully integrated to a streamline curvature method for axial-flow compressors | |
Albers et al. | Integrated structural and controller optimization for lightweight robot design | |
CN111651923A (zh) | 一种多自由度被动隔振系统优化设计方法 | |
JP5907607B2 (ja) | 処理配置方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130327 Termination date: 20151230 |
|
EXPY | Termination of patent right or utility model |