日本无码免费高清在线|成人日本在线观看高清|A级片免费视频操逼欧美|全裸美女搞黄色大片网站|免费成人a片视频|久久无码福利成人激情久久|国产视频一二国产在线v|av女主播在线观看|五月激情影音先锋|亚洲一区天堂av

  • 手機(jī)站
  • 小程序

    汽車測(cè)試網(wǎng)

  • 公眾號(hào)

    • 汽車測(cè)試網(wǎng)

    • 在線課堂

    • 電車測(cè)試

QCT 55靜載試驗(yàn)遲滯變形率數(shù)據(jù)處理方法

2019-09-28 18:09:33·  來(lái)源：座椅檢測(cè)知多少  
　
《QC/T 55-1993汽車座椅動(dòng)態(tài)舒適性試驗(yàn)方法》標(biāo)準(zhǔn)中有一個(gè)靜載試驗(yàn)的試驗(yàn)項(xiàng)目，目的是求座椅坐墊和靠背的靜剛度和遲滯變形率。本文主要是針對(duì)靜載試驗(yàn)項(xiàng)目中遲滯
《QC/T 55-1993汽車座椅動(dòng)態(tài)舒適性試驗(yàn)方法》標(biāo)準(zhǔn)中有一個(gè)靜載試驗(yàn)的試驗(yàn)項(xiàng)目，目的是求座椅坐墊和靠背的靜剛度和遲滯變形率。本文主要是針對(duì)靜載試驗(yàn)項(xiàng)目中遲滯變形率試驗(yàn)使用Python軟件進(jìn)行數(shù)據(jù)處理。
好了，廢話不多說(shuō)，下面先看一下試驗(yàn)方法。
 
1、試驗(yàn)要求
（1）被試樣機(jī)的結(jié)構(gòu)、靜態(tài)及動(dòng)態(tài)特性應(yīng)能代表實(shí)際的或者準(zhǔn)備生產(chǎn)的系列產(chǎn)品。
（2）在萬(wàn)能材料試驗(yàn)機(jī)平臺(tái)上安裝座椅時(shí)其位置和傾角應(yīng)與車上實(shí)際安裝狀態(tài)一致。
（3）如圖1所示，加載板加載面的形狀根據(jù)GB11559規(guī)定的座墊及靠背形狀確定。加載板表面應(yīng)光滑、具有足夠的剛度與強(qiáng)度，它與試驗(yàn)機(jī)用萬(wàn)向接頭連接，萬(wàn)向接頭的中心點(diǎn)要與加載板的載荷中心點(diǎn)重合。座墊加載板的質(zhì)量為51kg。
應(yīng)沿鉛垂方向向座墊加載，向靠背加載的方向應(yīng)與人體軀干線成直角。加載位置及方向如圖2所示。
圖1  加載板
圖2  加載位置及方向
 
2、試驗(yàn)方法
（1）將座椅安裝在座椅工裝上。
（2）將座椅及工裝安裝在萬(wàn)能材料試驗(yàn)機(jī)平臺(tái)上。
（3）按照GB11563就坐的方法安放加載板，注意調(diào)整加載板載荷方向應(yīng)與坐墊和靠背基本垂直。
（4）預(yù)加載。按照標(biāo)準(zhǔn)要求的方向預(yù)加載、卸載兩次。預(yù)加載載荷為：坐墊700N；靠背300N。
（5）正式試驗(yàn)。先對(duì)靠背和坐墊施加5N的初始載荷，然后以150~300mm/min的速度進(jìn)行加載和卸載。坐墊的加載力為700N，靠背的加載力為300N。從加載到卸載之間的過(guò)渡時(shí)間不應(yīng)超過(guò)2s。讀取并記錄加載力、位移、加載速度和卸載速度。
（6）重復(fù)上述過(guò)程，進(jìn)行三次試驗(yàn)，試驗(yàn)要間隔1min以上。每次測(cè)得的靜剛度值與三次平均值的偏差不應(yīng)超過(guò)±10％，然后取其平均值。
圖3  試驗(yàn)過(guò)程示意圖
備注：圖片只用于試驗(yàn)過(guò)程示意使用，未經(jīng)允許，不得傳播！
 
3、原始數(shù)據(jù)及圖像
原始數(shù)據(jù)保存在一個(gè).csv格式的文本文件中，用Excel表格打開(kāi)之后繪制圖像如圖4所示。
圖4  原始數(shù)據(jù)圖像
備注：原始數(shù)據(jù)只作為示意使用，不可以用于除本文之外的地方，也不可用于評(píng)判座椅好壞！
 
4、數(shù)據(jù)處理
遲滯變形率a的計(jì)算公式為：
a=（oabcdo面積）/（oabeo面積）x100%
因此，要求出遲滯變形率a，必須先計(jì)算得到載荷—撓度曲線上oabcdo和oabeo的面積。由于曲線并不規(guī)整，圍成的曲面也不容易分割成三角形、四邊形等容易求出面積的幾何圖形，成為數(shù)據(jù)分析的難點(diǎn)。
我最先開(kāi)始的時(shí)候是想把這些離散點(diǎn)擬合出函數(shù)表達(dá)式，然后再求函數(shù)式的定積分。使用Python二項(xiàng)式擬合的方法，設(shè)置次數(shù)為10。擬合出oab段的曲線如圖5所示，從圖中可以看出誤差線數(shù)據(jù)很小，說(shuō)明10次函數(shù)擬合效果還可以。
圖5 oab段曲線擬合數(shù)據(jù)
正當(dāng)我想就這么做的時(shí)候，突然發(fā)現(xiàn)辛普森積分和梯形積分法可以直接對(duì)離散數(shù)據(jù)進(jìn)行積分，于是我用一個(gè)2次函數(shù)y=x^2測(cè)試了一下。測(cè)試取0-10之間的定積分，計(jì)算結(jié)果如表1所示，發(fā)現(xiàn)計(jì)算結(jié)果竟然挺精準(zhǔn)的，且辛普森積分更接近真值。于是省了好多事呀??！
表1  測(cè)試函數(shù)結(jié)果
積分方式
數(shù)值
手工計(jì)算
1000/3近似等于333.333333
辛普森積分
333.16666666666663
梯形積分法
334.5
 
5、最終Python代碼
最終的程序代碼非常簡(jiǎn)單，為了對(duì)比結(jié)果特地同時(shí)采用了辛普森積分和梯形積分法兩種方式，最終計(jì)算結(jié)果如表2所示。從表中遲滯變形率可以看出，兩種計(jì)算方法得出的結(jié)果基本相同。
表2  遲滯變形率計(jì)算結(jié)果
積分方式
oabeo的面積
debed的面積
遲滯變形率
辛普森積分
8449.750404638784
6238.860836692514
26.17 %
梯形積分法
8449.097405916
6239.0125860240005
26.16 %
源程序代碼如下所示。
import numpy as np
importmatplotlib.pyplot as plt
fromscipy.integrate import simps
num =np.loadtxt(r"C:\Pythonwork\Q55積分求面積\遲滯變形率原始數(shù)據(jù).csv",skiprows=1, delimiter=",", usecols=(0, 1), unpack=True)
dis = num[0,:] # 讀取第幾列數(shù)據(jù)
force = num[1,:] # 讀取第幾列數(shù)據(jù)
dis_max =np.max(dis)
print("a值的最大值為:%f" % dis_max)
dis_argmax =np.argmax(dis)
print("a中最大值的索引值為:%d" % dis_argmax)
# oabeo數(shù)據(jù)
dis1 = num[0,0:dis_argmax+1]
force1 = num[1,0:dis_argmax+1]
# 使用自帶的兩種積分方法，對(duì)比積分結(jié)果
area1 =simps(force1, dis1) # 辛普森積分
print("辛普森積分oabeo的面積是：\n",area1)
area11 =np.trapz(force1, dis1) # 梯形積分法
print("梯形積分法oabeo的面積是：\n",area11)
# debcd數(shù)據(jù)
dis2 = num[0,dis_argmax+1:]
force2 = num[1,dis_argmax+1:]
area2 =abs(simps(force2, dis2))
print("辛普森積分debed的面積是：\n",area2)
area22 =abs(np.trapz(force2, dis2))
print("梯形積分法debed的面積是：\n",area22)
# 求遲滯變形率a
a1 =(area1-area2)/area1*100
print("辛普森積分遲滯變形率a是：%.2f%%" % a1)
# 求遲滯變形率a
a2 =(area11-area22)/area11*100
print("梯形積分法遲滯變形率a是：%.2f%%" % a2)
分享到：
 
反對(duì) 0 舉報(bào) 0 收藏 0 評(píng)論 0
滬ICP備11026917號(hào)-25