电子封装机械仿真：冲击分析_（3）.冲击力学基本理论.docxVIP

PAGE1

PAGE1

冲击力学基本理论

引言

在电子封装设计中，冲击分析是确保产品在各种外部冲击条件下正常工作的关键步骤。冲击分析主要涉及研究电子封装在受到外部冲击力作用时的动态响应，包括加速度、速度、位移以及应力和应变的分布。本节将详细介绍冲击力学的基本理论，为后续的电子封装机械仿真提供理论基础。

冲击的定义和分类

冲击是指一种瞬时的、高强度的动态载荷，通常在很短的时间内作用于物体上。冲击可以分为以下几类：

机械冲击：由机械运动引起的冲击，例如跌落、碰撞等。

热冲击：由温度的快速变化引起的冲击，例如快速加热或冷却。

电冲击：由电压或电流的突变引起的冲击，例如雷击、静电放电等。

在电子封装机械仿真中，我们主要关注机械冲击。机械冲击可以进一步分为单次冲击和重复冲击。单次冲击通常用于模拟跌落、碰撞等一次性事件，而重复冲击则用于模拟振动、疲劳等长期作用的动态载荷。

冲击载荷的特性

冲击载荷的特性主要包括以下几个方面：

峰值加速度：冲击载荷的最大加速度值，通常用g表示，其中g是重力加速度。

持续时间：冲击载荷作用的时间，通常用毫秒（ms）表示。

波形：冲击载荷随时间变化的形状，常见的波形有半正弦波、三角波、梯形波等。

频谱：冲击载荷的频率分布，用于描述冲击载荷的频域特性。

冲击载荷的波形

半正弦波

半正弦波是最常见的冲击载荷波形，其数学表达式为：

a

其中，A是峰值加速度，T是持续时间。

三角波

三角波的数学表达式为：

a

其中，A是峰值加速度，T是持续时间。

梯形波

梯形波的数学表达式为：

a

其中，A是峰值加速度，T是总持续时间，T1

冲击载荷的频谱

冲击载荷的频谱可以通过傅里叶变换获得。频谱分析可以帮助我们了解冲击载荷在不同频率下的能量分布，从而更好地设计和优化电子封装结构。

傅里叶变换

傅里叶变换是一种将时域信号转换为频域信号的数学工具。对于一个时域信号xt，其傅里叶变换Xf

X

例子：半正弦波的频谱分析

假设我们有一个半正弦波冲击载荷，峰值加速度A=100g，持续时间$T=10,

importnumpyasnp

importmatplotlib.pyplotasplt

#定义冲击载荷参数

A=100#峰值加速度，单位：g

T=10#持续时间，单位：ms

g=9.81#重力加速度，单位：m/s^2

#生成时间序列

t=np.linspace(0,T,1000)#1000个时间点

#生成半正弦波冲击载荷

a=A*np.sin(np.pi*t/T)*g

#计算频谱

frequencies=np.fft.fftfreq(t.size,d=(t[1]-t[0])/1000)#频率轴

spectrum=np.fft.fft(a)#计算傅里叶变换

#绘制时域和频域图

plt.figure(figsize=(12,6))

plt.subplot(1,2,1)

plt.plot(t,a)

plt.title(半正弦波冲击载荷(时域))

plt.xlabel(时间(ms))

plt.ylabel(加速度(m/s^2))

plt.grid(True)

plt.subplot(1,2,2)

plt.plot(frequencies,np.abs(spectrum))

plt.title(半正弦波冲击载荷(频域))

plt.xlabel(频率(Hz))

plt.ylabel(幅值)

plt.grid(True)

plt.tight_layout()

plt.show()

冲击响应谱(SRS)

冲击响应谱（SRS）是一种用于描述系统对冲击载荷响应的工具。SRS可以预测系统在不同频率下的最大响应，从而帮助设计者评估电子封装的耐冲击性能。

SRS的计算方法

SRS的计算方法主要包括以下步骤：

定义单自由度（SDOF）系统：假设系统为单自由度系统，其动力学方程为：

m

其中，m是质量，c是阻尼系数，k是刚度，ft

求解响应：使用数值方法求解上述动力学方程，得到系统的位移响应xt、速度响应xt和加速度响应

计算最大响应：对于每个频率，计算系统响应的最大值，得到SRS。

例子：计算SRS

假设我们有一个单自由度系统，质量m=1?kg，阻尼系数c=

importnumpyasnp

fromegrateimportodeint

importmatplotlib.pyplotasplt

#定义系统参数

m=1#质量，单位：kg

c=0.1#阻尼系数，单位：Ns/m

k=1000#刚度，单位