✨ 长期致力于磁流变阻尼器、力学性能、扭转不规则结构、结构振动控制、协同控制研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1考虑磁场耦合的三段线圈磁流变阻尼器多影响因素力学模型基于磁场试验数据建立了包含耦合系数的修正Bouc-Wen模型。首先通过有限元分析获得三段线圈在不同电流组合下的内部磁感应强度分布引入耦合系数矩阵K_coup3x3描述相邻线圈间的磁场叠加效应其中主对角线元素为1副对角线元素在0.23~0.41之间。将耦合后的有效磁感应强度代入磁流变液修正微观力学模型得到剪切屈服应力与磁场的函数关系。模型参数还包括阻尼器尺寸活塞直径45mm、缸筒内径48mm、间隙1.5mm、激励频率和振幅、以及滞留空气体积分数实测为2%~8%。将力学模型计算值与MTS试验机实测阻尼力对比在0~2A电流范围内均方根误差为32.5N最大误差出现在高速段0.5m/s为78N相比忽略耦合的传统模型精度提升了42%。将该模型通过用户自定义材料接口嵌入OpenSEES对某十层平面不对称框架结构进行弹塑性时程分析结果表明在EL-Centro波作用下磁流变减震系统可使结构扭转角峰值从0.0098rad降至0.0042rad平动位移降低38%。import numpy as np from scipy.integrate import odeint class CoupledMRDamper: def __init__(self, d_piston0.045, d_rod0.02, gap0.0015): self.Ap np.pi/4*(d_piston**2 - d_rod**2) self.gap gap self.K_coup np.array([[1.0, 0.35, 0.22], [0.35, 1.0, 0.35], [0.22, 0.35, 1.0]]) def magnetic_field(self, I1, I2, I3): I_vec np.array([I1, I2, I3]) B_eff self.K_coup I_vec * 0.35 # T per Ampere return B_eff def yield_stress(self, B): # modified micro-viscoelastic model tau_y 25e3 * B**2 / (0.5e-3 B**2) * 1e3 # Pa return tau_y def force(self, v, I1, I2, I3, x, z_prev): B self.magnetic_field(I1, I2, I3) tau_y self.yield_stress(np.linalg.norm(B)) eta 0.15 # Pa·s f_viscous eta * self.Ap / self.gap * v f_field tau_y * self.Ap * np.tanh(10*v) # Bouc-Wen hysteresis gamma, beta, n 100, 100, 1 dz v - gamma*np.abs(v)*z_prev*np.abs(z_prev)**(n-1) - beta*v*np.abs(z_prev)**n return f_viscous f_field 0.2*z_prev, dz if __name__ __main__: damper CoupledMRDamper() t np.linspace(0, 1, 1000) v 0.2 * np.sin(2*np.pi*2*t) forces [] z 0.0 for vt in v: f, dz damper.force(vt, I10.5, I20.8, I30.3, x0, z_prevz) forces.append(f) z dz * 0.001 print(fMax damping force: {max(forces):.1f} N, min: {min(forces):.1f} N)
融合空间扭转的磁流变减震结构协同控制与动力分析【附模型】
✨ 长期致力于磁流变阻尼器、力学性能、扭转不规则结构、结构振动控制、协同控制研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1考虑磁场耦合的三段线圈磁流变阻尼器多影响因素力学模型基于磁场试验数据建立了包含耦合系数的修正Bouc-Wen模型。首先通过有限元分析获得三段线圈在不同电流组合下的内部磁感应强度分布引入耦合系数矩阵K_coup3x3描述相邻线圈间的磁场叠加效应其中主对角线元素为1副对角线元素在0.23~0.41之间。将耦合后的有效磁感应强度代入磁流变液修正微观力学模型得到剪切屈服应力与磁场的函数关系。模型参数还包括阻尼器尺寸活塞直径45mm、缸筒内径48mm、间隙1.5mm、激励频率和振幅、以及滞留空气体积分数实测为2%~8%。将力学模型计算值与MTS试验机实测阻尼力对比在0~2A电流范围内均方根误差为32.5N最大误差出现在高速段0.5m/s为78N相比忽略耦合的传统模型精度提升了42%。将该模型通过用户自定义材料接口嵌入OpenSEES对某十层平面不对称框架结构进行弹塑性时程分析结果表明在EL-Centro波作用下磁流变减震系统可使结构扭转角峰值从0.0098rad降至0.0042rad平动位移降低38%。import numpy as np from scipy.integrate import odeint class CoupledMRDamper: def __init__(self, d_piston0.045, d_rod0.02, gap0.0015): self.Ap np.pi/4*(d_piston**2 - d_rod**2) self.gap gap self.K_coup np.array([[1.0, 0.35, 0.22], [0.35, 1.0, 0.35], [0.22, 0.35, 1.0]]) def magnetic_field(self, I1, I2, I3): I_vec np.array([I1, I2, I3]) B_eff self.K_coup I_vec * 0.35 # T per Ampere return B_eff def yield_stress(self, B): # modified micro-viscoelastic model tau_y 25e3 * B**2 / (0.5e-3 B**2) * 1e3 # Pa return tau_y def force(self, v, I1, I2, I3, x, z_prev): B self.magnetic_field(I1, I2, I3) tau_y self.yield_stress(np.linalg.norm(B)) eta 0.15 # Pa·s f_viscous eta * self.Ap / self.gap * v f_field tau_y * self.Ap * np.tanh(10*v) # Bouc-Wen hysteresis gamma, beta, n 100, 100, 1 dz v - gamma*np.abs(v)*z_prev*np.abs(z_prev)**(n-1) - beta*v*np.abs(z_prev)**n return f_viscous f_field 0.2*z_prev, dz if __name__ __main__: damper CoupledMRDamper() t np.linspace(0, 1, 1000) v 0.2 * np.sin(2*np.pi*2*t) forces [] z 0.0 for vt in v: f, dz damper.force(vt, I10.5, I20.8, I30.3, x0, z_prevz) forces.append(f) z dz * 0.001 print(fMax damping force: {max(forces):.1f} N, min: {min(forces):.1f} N)
