APDL瞬态分析结束后,可以对某一时间步下的模型直接进行一个预应力模态分析嘛?

随着物流业的快速发展,牵引车也向着高速重载方向发展,由于车辆工作的多工况,车辆结构受力复杂,因此经常在使用过程中发生车架结构早期断裂破坏。这种现象已经成为困扰汽车设计部门的难题,并造成了明显的经济损失。
目前在车架强度研究方面普遍采用的方法为试验方法和计算机仿真方法。试验方法是一种直接而客观的设计与验证方法,结果最有说服力,但是试验周期较长,费用昂贵且重复性差,而且也不可能采集到足够多的数据,特别是用于产品改进所需的详细数据资料;计算机仿真方法具有费用低廉、周期较短、可重复性、结果信息全面等优点。通过仿真分析,可以在短期内,可以完成多次改进分析,并以此作为产品改进的依据。从模拟仿真和试验对产品安全进行验证,能确保分析结果的可信度和结构的可靠性。
本文利用有限元软件HyperMesh建立某牵引车车架有限元模型,利用Altair OptiStruct 进行静强度计算模拟,并与实际情况进行对比,对车架横梁处存在问题进行DOE优化设计来提高横梁的强度。

本文中某牵引车车架有限元模型,网格单元平均尺寸=10mm,焊接熔核直径=6mm,焊点由面到面连接方式生成,板簧使用Spring单元和RBE2模拟,安装在车架左右两侧的油箱和电瓶箱以及驾驶室和发动机本身结构的形状对车架的有限元强度和刚度分析影响很小;但是其本身的质量很大,对车架结构的强度和刚度影响很大,根据工程经验,选择用质量点单元mass来模拟简化油箱和电瓶箱的结构,这样可以减少建模时间,提高工作效率,车架有限元模型如图1所示,车架的材料数据如表1所示。

表1 牵引车车架有限元模型材料参数属性

在进行静强度分析时,我们分别选择对横梁影响比较大的转向静态受力工况、扭转静态受力工况这两个典型工况为研究对象。
2.1 转向静态受力工况
车架转向受力载荷施加图如图2所示,应力云图如图3所示。

图2 车架转向受力工况载荷图

在应力云图3中我们可以看出,应力比较大的位置1处为车架横梁,也是车架真实断裂位置(车架真实断裂图如图4所示),此处的应力属于应力分布的高应力区域,根据断裂图片对比,该工况与真实情况一致。此处结构需要进行优化改进。位置2、3处应力值比较大,此处结构超过屈服极限,应当加以关注,进行适当的结构修改减少局部的应力集中。

图3 车架转向受力工况应力分布云图

图4 车架横梁实际破坏图

2.2 扭转静态受力工况

车架转向受力载荷施加图如图5所示,应力云图如图6所示。

图5 车架扭转受力工况载荷图

从应力分布云图6中我们发现,位置1处为横梁最大应力值,大于材料的屈服极限,也是真实的破坏位置,应进行优化改善。位置2处为横梁应力值大于材料的屈服极限,应高度重视,极有可能产生疲劳破坏。位置3,4处应力值比较高,但是没有超过材料屈服极限,所以不会对该零件造成强度破坏。

图6 车架扭转受力工况应力分布云图

3 基于DOE技术的横梁结构优化

DOE-Design of experiment(试验设计)主要应用在试验的合理安排和试验数据的处理上,利用数理统计学的基本知识,讨论如何合理地安排试验、取得数据,然后进行综合科学分析,从而尽快获得最优组合方案。
针对本车架横梁处强度较差的特点,本文运用OptiStruct软件对其进行了DOE优化设计,希望通过改变横梁及其加强件的厚度值,来增强部件的强度,从而达到提高整个车架的承载能力的目的,我们在横梁与车架纵梁连接处添加了两块加强板,添加加强板前后状态如图7、图8所示,加强板与横梁体铆接。

图8 横梁添加加强板状态

DOE设计中,因素水平数的选择十分关键,因为如果因素水平选取过多,将会使试验次数急剧增加,既不经济,又费时间,如果选取过少,又有可能造成无法选到主要因素的后果。我们选择对横梁影响比较大的典型工况扭转工况为研究对象,分别选用横梁厚度与加强板厚度作为试验设计的2个因素,进行交叉试验,试验次数为31次,车架的安全系数有较大的差异,如表2所示。
我们也可以从图9、图10中可以更加直观的看出横梁厚度、加强板厚度和安全系数之间的关系。

表2 扭转工况下的DOE分析数据

图9 扭转工况下DOE响应面图

图10 扭转工况下DOE等高线图

通过DOE分析,我们可以更加全面地得到不同横梁厚度和加强板厚度之间组合横梁总成的车架安全系数。综合考虑实际工艺、模具及成本等方面的因素,原始结构状态为表2中的第1组组合,在实际生产中选择表2中第16组组合,以此进行结构的改进。

3.2 横梁结构优化前后模态对比

我们将横梁原始结构的自由模态与改进后横梁的自由模态的第一阶频率进行对比,原始结构与改进结构自由模态的第一阶频率分别如图11、图12。

图11 原始结构第一阶频率

图12 改进后结构的第一阶频率

由图11、图12可以看出自由模态幅频特性:
3.3 横梁结构优化前后应力对比
本文仍然选择对横梁影响比较大的转向受力工况和扭转受力工况为研究对象。
3.3.1 转向静态受力工况
横梁结构改进前后在转向静态受力工况时的应力云图分别如图13、图14所示。
由图13、图14可以看出最大应力幅值(安全系数):

图13 原始结构应力云图

图14 改进后结构的应力云图

横梁结构改进前后在扭转态受力工况时的应力云图分别如图15、图16所示。

图15 原始结构应力云图

图16 改进后结构的应力云图

本文通过CAE软件HyperMesh建立了某牵引车的车架模型,利用Altair OptiStruct对其进行了静强度计算,得到了车架模态、应力等数值。本文采用DOE优化设计方法对车架横梁进行优化设计,优化结果表明优化方案相比初始设计方案明显降低了车车架横梁的应力,增强了车架强度,提高了车辆的承载能力,基于DOE的车架优化分析具有很强的工程实用性,对于车架结构变更设计具有巨大的指导意义。


开放分享:优质有限元技术文章,助你自学成才

《谐响应、响应谱分析、随机振动与模态分析》由会员分享,可在线阅读,更多相关《谐响应、响应谱分析、随机振动与模态分析(95页珍藏版)》请在人人文库网上搜索。

1、谐响应分析 谐响应分析的定义和目的 关于谐响应分析的基本术语和概念 谐响应分析在ANSYS中的应用 谐响应分析的实例练习 定义和目的 什么是谐响应分析? 确定一个结构在已知频率的正弦(简谐)载荷作用下结构响应的技 术。 输入: 已知大小和频率的谐波载荷(力、压力和强迫位移); 同一频率的多种载荷,可以是同相或不同相的。 输出: 每一个自由度上的谐位移,通常和施加的载荷不同相; 其它多种导出量,例如应力和应变等。 定义和目的 谐响应分析用于设计: 旋转设备(如压缩机、发动机、泵、涡轮 机械等)的支座、固定装置和部件; 受涡流(流体的漩涡运动)影响的结构, 例如涡轮叶片、飞机机翼、桥和塔等。 定义

2、和目的 为什么要作谐响应分析? 确保一个给定的结构能经受住不同频率的 各种正弦载荷(例如:以不同速度运行的 发动机); 探测共振响应,并在必要时避免其发生 (例如:借助于阻尼器来避免共振)。 术语和概念 包含的主题: 运动方程 谐波载荷的本性 复位移 求解方法 通用运动方程: F矩阵和 u矩阵是简谐的,频率为 w: 谐响应分析的运动方程: 运动方程 FuKuCuM titii

3、 = 实部, Fmaxcos F2 = 虚部, Fmaxsin umax= 位移幅值 f= 载荷函数的相位角 u1= 实部, umaxcosf u2= 虚部, umaxsinf 谐波载荷的本性 在已知频率下正弦变化; 相角允许不同相的多个 载荷同时作用, 缺省值 为零; 施加的全部载荷都假设是 简谐的,包括温度和重力。 实部 虚部 复位移 在下列情况下计算出的位移将是复数 具有阻尼 施加载荷是复数载荷(例如:虚部为非零的载 荷) 复位移滞后一个相位角(相对于某一个基 准而言) 可以用实部和虚部或振幅和相角的形式来 查看 谐响应分析-术语和概念 求解方法 求解简谐运动方程的三种方法: 完整法 为

4、缺省方法,是最容易的方法; 使用完整的结构矩阵,且允许非对称矩阵(例如:声学矩 阵)。 缩减法* 使用缩减矩阵,比完整法更快; 需要选择主自由度,据主自由度得到近似的 M矩阵和C矩阵。 模态叠加法* 从前面的模态分析中得到各模态;再求乘以系数的各模态之 和; 所有求解方法中最快的。 求解方法 完完整整法法缩缩减减法法模模态态叠叠加加法法 相对求解时间慢较快最快 相对的使用容易程度最容易较容易难 允许元素载荷(例如压强)吗?允许不允许允许 (一个载荷向量) 允许非零位移载荷吗?允许允许不允许 允许模态阻尼吗?不允许不允许允许 能处理预应力吗?不能能能 能进行“Restart“吗?能能不能 允许非

5、对称矩阵吗?允许不允许不允许 需要为了求解而选择模态吗?不需要不需要需要 需要选择主自由度吗?不需要需要需要 (如果选用缩减法) 步骤 四个主要步骤: 建模 选择分析类型和选项 施加谐波载荷并求解 观看结果 建模 模型 只能用于线性单元和材料,忽略各种非线性; 记住要输入密度; 注意: 如果ALPX(热膨胀系数)和T均不为零,就有 可能不经意地包含了简谐热载荷。为了避免这种事情发生,

6、 选择分析类型和选项 选择分析类型和选项 进入求解器,选择谐响应分析; 设置分析选项 1求解方法 2自由度输出格式 3是否使用集中质量逼近(用于结构的 一个方向的尺寸远小于另两个 方向的尺寸的情况中。例如: 细长梁与薄壳。) 典型命令: /SOLU ANTYPE,HARMIC,NEW HROPT,FULL HROUT,ON LUMPM,0 选择分析类型和选项 设置阻尼参数 从

7、c Damping 施加谐波载荷并求解 施加谐波载荷并求解 所有施加的载荷以规定的频率(或频率范 围)简谐地变化 “载荷”包括: 位移约束-零或非零的 作用力 压强 注意:如果要施加重力和热载荷,它 们也被当作简谐变化的载荷来考虑! 典型命令: DK, ! 或 D或DSYM DA,. DL, *AFUN,DEG FK, F, SFA, SFL, SFE, SF, 施加谐波载荷并求解 规定谐波载荷时要包括: 振幅和相角 频率 不同频率载荷具有不同的幅值时的处理 方法 载荷值(大小)代表振幅 Fmax 相角 f是在两个或两个以上谐波 载荷间的相位差,单一载荷不需 要相角。 实部 虚部 f F2ma

8、x 施加谐波载荷并求解 振幅和相角 ANSYS 不能直接输入振幅和相 角,而是规定实部和虚部分量。 幅值/相角和实部/虚部的关系: 可以使用APDL语言计算,但要确 保角度单位为度(缺省为弧度)。 sin cos ampimag ampreal FF FF 施加谐波载荷并求解 谐波载荷的频率: 通过频率范围和在频率范围 内的子步数量来规定每秒的 循环次数(赫兹);

9、sTime/ Frequenc Freq and Substps 典型命令: HARFRQ,0,50, NSUBST,10, KBC,1 施加谐波载荷并求解 不同频率载荷具有不同的幅值时的处理方法 在施加谐波载荷后,下一步就是开始求解了,通常采用一个载荷 步,但是可以采用若干子步,且每个子步具有不同的频率范围。当 不同频率的载荷具有不同的幅值时,可以分多个载荷步施加。 查看结果 1.绘制结构上的特殊点处的位移-频率曲线 2.确定各临界频率和相应的相角 3.观看整个结构在各临界频率和相角时的位移和应力 典型命令: /POST26 NSOL, PLVAR,. 查看结果 确定各临界频率 和相角 用图

10、形显示最高振幅 发生时的频率; 由于位移与施加的载 荷不同步(如果存在 阻尼的话),需要确 定出现最大振幅时的 相角; 要进行上述工作, 首先要选择振幅+ 相位选项。 查看结果 然后,列表显示变量。找出 响应最大值所对应的频率和 相角,据此确定最大振幅出 现的载荷步和子步号。 进入POST1处理器,观察该 载荷子步下的结果,查看结 构的应力和变形情况。 典型命令: HRCPLX,

11、07 蓝色:305 预应力谐响应分析 分为四个步骤 1.定义单元类型、材料属性、建模 2.进入求解器,指定分析类型静力分析,施加约 束和预应力,定义预应力分析,求解。 3.重新进入求解器(必须显式地,退出再进入),定 义谐响应分析,删除预应力载荷,施加谐波载荷, 求解。 4.查看结果 注意:第2、3步中均开启预应力效应(PSTRES,ON)。 预应力谐响应分析实例 张紧的吉他弦的谐响应分析 输入文件:presharmonic.cmds 126N预紧力 Y向谐波激励 预应力谐响应分析实例 预应力对响应(节点16,uy)的影响 无预应力有预应力 谱分析 什么是谱分析? 是模态分析的扩展,用于计算结

请问您目前有没有解决这个问题啊?我目前也遇到这样的问题,请问该怎么做啊?

我要回帖

更多关于 混合效应模型什么时候使用 的文章

 

随机推荐