CN105825010B - 基于网面有限元建模的蹦床动态应力测试分析方法 - Google Patents

基于网面有限元建模的蹦床动态应力测试分析方法 Download PDF

Info

Publication number
CN105825010B
CN105825010B CN201610148900.XA CN201610148900A CN105825010B CN 105825010 B CN105825010 B CN 105825010B CN 201610148900 A CN201610148900 A CN 201610148900A CN 105825010 B CN105825010 B CN 105825010B
Authority
CN
China
Prior art keywords
trampoline
net
trampoline net
half period
finite element
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
Application number
CN201610148900.XA
Other languages
English (en)
Other versions
CN105825010A (zh
Inventor
钱竞光
李兆霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201610148900.XA priority Critical patent/CN105825010B/zh
Publication of CN105825010A publication Critical patent/CN105825010A/zh
Application granted granted Critical
Publication of CN105825010B publication Critical patent/CN105825010B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Laminated Bodies (AREA)

Abstract

本发明涉及一种基于网面有限元建模的蹦床动态应力测试分析方法，包括：静载实验，动载实验，针对蹦床网面建立有限元模型，通过蹦床网面有限元模型的模拟确定载荷限定，通过蹦床网面有限元模型筛选得到蹦床运动员起跳动作的推荐跳跃半周期。本发明方法的分析过程更加简便，能更加迅速地得出蹦床运动员起跳动作的推荐跳跃半周期，利于更好地指导蹦床运动员的训练。

Description

基于网面有限元建模的蹦床动态应力测试分析方法
技术领域
本发明涉及一种基于网面有限元建模的蹦床动态应力测试分析方法，尤其适用于指导优化蹦床运动员起跳动作，属于蹦床运动员训练辅助技术领域。
背景技术
据申请人所知，蹦床运动早已成为奥运会的正式比赛项目，在该项运动中，人体在蹦床上的跳跃是一个非常复杂的过程，属于动力学范畴，可以通过动力学仿真和计算机模拟方法来研究起跳动作最佳化的问题，进而指导优化蹦床运动员起跳动作。
申请人曾于2013年发表了两篇会议论文：《蹦床网面的载荷实验和有限元建模》和《蹦床的有限元建模及最佳起跳周期的模拟计算》，在文中介绍了先将蹦床网面进行有限元建模，再以模型进行模拟计算得出结果的基本方法。当时受篇幅所限，申请人未将具体方法写入上述会议论文。申请人在后续的研究中发现当时所用具体方法仍然有进一步优化的空间。
发明内容
本发明所要解决的技术问题是：针对现有技术存在的问题，提出一种基于网面有限元建模的蹦床动态应力测试分析方法，分析过程更加简便，能更加迅速地得出蹦床运动员起跳动作的推荐跳跃半周期，利于更好地指导蹦床运动员的训练。
本发明的技术方案如下：
一种基于网面有限元建模的蹦床动态应力测试分析方法，包括以下步骤：
第一步、在蹦床网面中心区域进行静载实验得到静载实验数据，并计算出蹦床网面中心点等效刚度k0；
第二步、在蹦床网面中心区域进行动载实验得到动载实验数据，并计算出蹦床阻尼比ξ0；
第三步、采用ANSYS软件，针对蹦床网面建立有限元模型；
其特征是，还包括：
第四步、采用ANSYS软件，在蹦床网面有限元模型中：
S1、加载温度载荷，使蹦床网面具有预应力；
S2、在蹦床网面中心区域按第一步实验条件模拟进行静载实验，并得出一组与多个载荷分别对应的网面位移模拟平均值；
S3、将各网面位移模拟平均值与相应的第一步实际测得的网面位移值进行比较和计算，并得出偏差率-载荷分布图；
S4、根据预设的最高偏差率限定，结合偏差率-载荷分布图得出载荷限定区间或载荷下限值或载荷上限值；
第五步、采用ANSYS软件，在蹦床网面有限元模型中：
T1、先重复第四步S1，再在蹦床网面中心区域加载运动员质量，并计算得出一组蹦床网面N阶固有频率，N为正整数；其中，蹦床网面N阶固有频率至少包括蹦床网面一阶固有频率；
T2、将从人体开始接触蹦床网面到人体离开蹦床网面的时间作为跳跃半周期；以蹦床网面N阶固有频率的倒数的一半作为备选半周期；逐一判断各备选半周期是否为人体蹬网能实现的跳跃半周期，若否则剔除该备选半周期，若是则保留该备选半周期；
T3、采用与保留的备选半周期对应的蹦床网面N阶固有频率、及蹦床网面中心点等效刚度k0，按预定算法得出蹦床等效质量m0，然后结合蹦床阻尼比ξ0，按预定算法得出蹦床阻尼比例系数c0；根据运动员质量及保留的各备选半周期，分别模拟人体与蹦床接触点的位移轨迹；
T4、人体蹬网时对蹦床网面的冲击力Fc＝m0y”+c0y’+k0y，其中，y”、y’、y分别表示人体与蹦床接触点垂直方向的加速度、速度以及位移，并由所述位移轨迹推算得出；
T5、当第四步S4所得的是载荷限定区间时，按第五步T4计算冲击力Fc的极大值和极小值，并判断该极大值和极小值是否均属于载荷限定区间，若是则保留与该冲击力Fc对应的备选半周期，若两者之一不属于或两者均不属于载荷限定区间则剔除该备选半周期；
当第四步S4所得的是载荷下限值时，按第五步T4计算冲击力Fc的极小值，并判断该极小值是否大于或等于载荷下限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
当第四步S4所得的是载荷上限值时，按第五步T4计算冲击力Fc的极大值和，并判断该极大值是否小于或等于载荷上限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
T6、将保留的备选半周期作为蹦床运动员起跳动作的推荐跳跃半周期。
本发明进一步完善的技术方案如下：
优选地，第一步中，所述静载实验采用若干杠铃片加载于蹦床网面，各杠铃片中心对准蹦床网面中心；以蹦床网面中心点为原点，在蹦床网面建立直角坐标系；以杠铃片边缘与两坐标轴的四个交点、及蹦床网面中心点作为测量点；以杠铃片覆盖的蹦床网面区域作为蹦床网面中心区域。
优选地，第一步中，静载实验的具体过程如下：
在蹦床网面中心区域加载预设初始质量的杠铃片，测量各测量点的纵向位移；继续多次加载预设质量的杠铃片，直至加载杠铃片的总质量达到预设质量上限，每次加载后测量各测量点的纵向位移；所述纵向位移在蹦床网面未加载杠铃片时为零。
优选地，第一步中，静载实验计算蹦床网面中心点等效刚度k0的具体过程如下：
将加载杠铃片的总质量作为载荷Fy，并将其与蹦床网面中心点的纵向位移y的关系进行线性拟合，所得式为Fy＝Ay+B，A、B分别为常数；其中常数A即为蹦床网面中心点等效刚度k0。
优选地，第二步中，动载实验的具体过程如下：
将预设质量的小球固定在蹦床网面中心点上，蹦床网面中心点即为小球与蹦床网面的接触点；对小球施加垂直向下的激励，使小球在垂直方向上振动；通过高速摄像机记录小球在震动过程中的最大振幅，并通过影像解析确定小球在每个衰减周期的振幅。
优选地，第二步中，动载实验计算蹦床阻尼比ξ0的具体过程如下：
先利用下式计算对数减缩系数Λ：
其中，A1为最大振幅，Aj+1为振动开始后小球的第j+1个振幅，N为小球振幅的总数量；
再利用下式计算蹦床阻尼比ξ0：
优选地，第三步中，建立有限元模型时，蹦床网面的结构采用空间铰接体系，结构边界为点支撑铰接或者线支撑铰接边界，且蹦床网面保持弹性状态。
优选地，第四步S1中，利用下式计算温度载荷Δt：
其中，T为蹦床网面的初始预应力，E为蹦床网面的弹性模量，a为蹦床网面的膨胀系数，tn为蹦床网面的厚度。
优选地，第五步T3中，先利用下式计算人体-蹦床系统的等效质量m1：其中，f为蹦床网面N阶固有频率；再以人体-蹦床系统的等效质量m1减去运动员质量，即得蹦床等效质量m0。
优选地，第五步T3中，利用下式计算蹦床阻尼比例系数c0：
与现有技术相比，本发明方法的分析过程更加简便，能更加迅速地得出蹦床运动员起跳动作的推荐跳跃半周期，利于更好地指导蹦床运动员的训练。
附图说明
图1为本发明试验案例中蹦床网面中心区域载荷与纵向位移的关系示意图。
图2为本发明试验案例中蹦床网面的有限元模型示意图。
图3、图4分别为本发明试验案例中蹦床有限元模拟变形图、蹦床网面米塞斯应力分布图。
图5为本发明试验案例中ANSYS模拟值与实验测得值的比较示意图。
图6为本发明试验案例中ANSYS模拟值与实验测得值的偏差率示意图。
图7、图8分别为本发明试验案例中蹦床模型一阶振型立面图和侧面图。
图9、图10分别为本发明试验案例中蹦床模型二阶振型立面图和侧面图。
图11、图12分别为本发明试验案例中蹦床模型三阶振型立面图和侧面图。
具体实施方式
下面参照附图并结合实施例对本发明作进一步详细描述。但是本发明不限于所给出的例子。
本发明基于网面有限元建模的蹦床动态应力测试分析方法，在具体实施时包括以下步骤：
第一步、在蹦床网面中心区域进行静载实验得到静载实验数据，并计算出蹦床网面中心点等效刚度k0。
具体而言，静载实验采用若干杠铃片加载于蹦床网面，各杠铃片中心对准蹦床网面中心；以蹦床网面中心点为原点，在蹦床网面建立直角坐标系；以杠铃片边缘与两坐标轴的四个交点、及蹦床网面中心点作为测量点；以杠铃片覆盖的蹦床网面区域作为蹦床网面中心区域。
静载实验的具体过程如下：
在蹦床网面中心区域加载预设初始质量的杠铃片，测量各测量点的纵向位移；继续多次加载预设质量的杠铃片，直至加载杠铃片的总质量达到预设质量上限，每次加载后测量各测量点的纵向位移；纵向位移在蹦床网面未加载杠铃片时为零。
计算蹦床网面中心点等效刚度k0的具体过程如下：
将加载杠铃片的总质量作为载荷Fy，并将其与蹦床网面中心点的纵向位移y的关系进行线性拟合，所得式为Fy＝Ay+B，A、B分别为常数；其中常数A即为蹦床网面中心点等效刚度k0。
第二步、在蹦床网面中心区域进行动载实验得到动载实验数据，并计算出蹦床阻尼比ξ0。
具体而言，动载实验的具体过程如下：
将预设质量的小球固定在蹦床网面中心点上，蹦床网面中心点即为小球与蹦床网面的接触点；对小球施加垂直向下的激励，使小球在垂直方向上振动；通过高速摄像机记录小球在震动过程中的最大振幅，并通过影像解析确定小球在每个衰减周期的振幅。
计算蹦床阻尼比ξ0的具体过程如下：
先利用下式计算对数减缩系数Λ：
其中，A1为最大振幅，Aj+1为振动开始后小球的第j+1个振幅，N为小球振幅的总数量；
再利用下式计算蹦床阻尼比ξ0：
第三步、采用ANSYS软件，针对蹦床网面建立有限元模型。
具体而言，建立有限元模型时，蹦床网面的结构采用空间铰接体系，结构边界为点支撑铰接或者线支撑铰接边界，且蹦床网面保持弹性状态。
第四步、采用ANSYS软件，在蹦床网面有限元模型中：
S1、加载温度载荷，使蹦床网面具有预应力；
具体而言，利用下式计算温度载荷Δt：
其中，T为蹦床网面的初始预应力，E为蹦床网面的弹性模量，a为蹦床网面的膨胀系数，tn为蹦床网面的厚度。
S2、在蹦床网面中心区域按第一步实验条件模拟进行静载实验，并得出一组与多个载荷分别对应的网面位移模拟平均值；
S3、将各网面位移模拟平均值与相应的第一步实际测得的网面位移值进行比较和计算，并得出偏差率-载荷分布图；
S4、根据预设的最高偏差率限定，结合偏差率-载荷分布图得出载荷限定区间或载荷下限值或载荷上限值。
第五步、采用ANSYS软件，在蹦床网面有限元模型中：
T1、先重复第四步S1，再在蹦床网面中心区域加载运动员质量，并计算得出一组蹦床网面N阶固有频率，N为正整数；其中，蹦床网面N阶固有频率至少包括蹦床网面一阶固有频率；
T2、将从人体开始接触蹦床网面到人体离开蹦床网面的时间作为跳跃半周期；以蹦床网面N阶固有频率的倒数的一半作为备选半周期；逐一判断各备选半周期是否为人体蹬网能实现的跳跃半周期，若否则剔除该备选半周期，若是则保留该备选半周期；
T3、采用与保留的备选半周期对应的蹦床网面N阶固有频率、及蹦床网面中心点等效刚度k0，按预定算法得出蹦床等效质量m0，然后结合蹦床阻尼比ξ0，按预定算法得出蹦床阻尼比例系数c0；根据运动员质量及保留的各备选半周期，分别模拟人体与蹦床接触点的位移轨迹；
具体而言，先利用下式计算人体-蹦床系统的等效质量m1：其中，f为蹦床网面N阶固有频率；再以人体-蹦床系统的等效质量m1减去运动员质量，即得蹦床等效质量m0。
利用下式计算蹦床阻尼比例系数c0：
T4、人体蹬网时对蹦床网面的冲击力Fc＝m0y”+c0y’+k0y，其中，y”、y’、y分别表示人体与蹦床接触点垂直方向的加速度、速度以及位移，并由位移轨迹推算得出；
T5、当第四步S4所得的是载荷限定区间时，按第五步T4计算冲击力Fc的极大值和极小值，并判断该极大值和极小值是否均属于载荷限定区间，若是则保留与该冲击力Fc对应的备选半周期，若两者之一不属于或两者均不属于载荷限定区间则剔除该备选半周期；
当第四步S4所得的是载荷下限值时，按第五步T4计算冲击力Fc的极小值，并判断该极小值是否大于或等于载荷下限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
当第四步S4所得的是载荷上限值时，按第五步T4计算冲击力Fc的极大值和，并判断该极大值是否小于或等于载荷上限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
T6、将保留的备选半周期作为蹦床运动员起跳动作的推荐跳跃半周期。试验案例
一、执行本发明方法第一步。
以某蹦床队日常训练用蹦床为建模对象，该蹦床符合目前国际蹦床比赛最新标准要求。蹦床主体结构框架长5.050米、宽2.910米、高1.150米；网长4.028米、宽2.014米；蹦床的工作部分是由拉伸弹簧和纤维网面组成。静载实验所使用的实验器材包括：若干5kg至25kg的杠铃片，直径为37.50cm，面积约为552cm*cm，总质量为260kg；卷尺数个；加载区域的纵向位移由五把卷尺测量。各杠铃片中心对准蹦床网面中心，以蹦床网面中心点为原点建立坐标系，分别取中心点以及杠铃片边缘与网平面X轴、Z轴交点O(0,0,0)，A(-18.7cm,0,0)，B(0,-20.55cm,0)，C(18.9cm,0,0)，D(0,20.15cm,0)作为测量点。定义杠铃片覆盖的蹦床网面区域为蹦床网面中心区域。将卷尺统一固定在网面下面的地板上，在竖直方向上对应着这五个点，将五把卷尺抽出，挂在其对应的点上，用胶带固定，记录测量其高度。
用杠铃片对应蹦床网面中心施加载荷，每次加载后用卷尺测量蹦床网面中心区域点高度的变化来表示纵向位移。加载重量从40kg开始，以后每次加载均在原有基础上加10kg，直到载荷达到260kg。因为点A与C，B与D两两均在一条轴线上，用A和C，B和D的位移平均值去和O点的位移作比较，如果相差过大需重新测量。
利用载荷与O、A、B、C、D五个点纵向位移平均值的变化关系来表示蹦床网面中心区域竖直方向作用力与位移的变化，如图1所示：横轴表示O、A、B、C、D五个点垂直方向位移平均值y，纵轴表示作用于该测量点垂直方向的载荷Fy。
用三次函数拟合Fy与y的变化关系可以写成式(1)。
Fy＝4.76y3+1.97y2+3.65y (1)
从式(1)可以看出，在纵向位移y不是很大的情况下，Fy的线性化较为明显，Fy与y的变化关系以线性项为主；随着y的不断增加，Fy与y的比例系数不断上升，定义Fy与y的比例系数为蹦床网面中心区域的等效刚度，用k3表示。k3随着y的增大而增大，如式(2)：
k3＝Fy/y＝4.76y2+1.97y+3.65 (2)
由此可见，等效刚度k3是非线性的，而非线性系统的动力学性质十分复杂，因此为了方便问题的研究，需将蹦床模型再次简化。经研究发现k3在其测量的范围内其常数项占据主要地位，因此将蹦床网面中心区域的等效刚度简化成线性刚度。
将载荷Fy与测量点O的纵向位移y的关系进行线性拟合，得出测量点O的载荷Fy与纵向位移y的关系为式(3)，
Fy＝6.14y+0.004 (3)
由此可得，蹦床网面中心点等效刚度k0＝6.14KN/m。
二、执行本发明方法第二步。
将一个质量为4kg的小球固定在蹦床网面中心点O上，蹦床网面中心点O即为小球与网面的接触点；对小球施加垂直向下的激励，使小球在垂直方向上振动；通过高速摄像机记录小球在震动过程中的最大振幅，并通过影像解析确定小球在每个衰减周期的振幅，如表1所示。
表1衰减振动过程中的小球振幅变化(i＝1,…,7)
利用下式计算对数减缩系数Λ：
再利用下式计算蹦床阻尼比ξ0：
结果为：ξ0＝0.063。
三、执行本发明方法第三步。
明确了蹦床结构的材料属性和几何尺寸后，可以按本发明方法第三步建立蹦床网面的有限元模型，如图2所示。
四、执行本发明方法第四步。
加载温度载荷时，蹦床网面的弹性模量大小为7.35e6 N/m2，膨胀系数为10，厚度为0.001m，初始预应力为2058 N/m2，按前述公式可以得出△t的大小为0.028，即要对蹦床网面“降温”0.028℃，即对蹦床模型施加-0.028℃的温度载荷。
模拟静载实验后的变形图以及米塞斯应力分布图如图3、图4所示。
模拟静载实验得出一组与多个载荷分别对应的网面位移模拟平均值(ANSYS位移模拟值)，与第一步实际测得的位移值(实验测得值)进行对比，如图5所示。
由图5可知，在相同载荷条件下，实验测得网面位移值略高于ANSYS模拟值，二者最大偏差量发生在加载载荷为1078N，此时实验测试值比ANSYS模拟值大0.025m。
定义在相同的载荷条件下，静载实验测得的网面中心区域位移值与ANSYS模拟值的偏差量，和网面静载实验测得的网面位移实验值之比为在此载荷条件ANSYS模拟值的偏差率，如图6所示。从图6可以看出，偏差率随着载荷增加在不断减小，从一开始245N的载荷条件下的偏差率18.3％到2548N接近于0。这表明在本次有限元模拟中，网面的加载幅度越大，ANSYS模拟值越接近真实值的大小。
将最高偏差率限定为4％，结合图6可得载荷下限值为1.8KN。
五、执行本发明方法第五步。
目标运动员质量为43kg，在蹦床有限元模型的蹦床网面中心区域加载运动员质量，计算得出前三阶固有频率，如表2所示。
表2蹦床模型前三阶固有频率(i＝1,..,3)
蹦床网面的各阶振型的立面图与侧面图(xy平面)，如图7至图12所示。
从蹦床的前三阶模态上看，第一阶与第二、三阶频率相差较大；从图7与图8可以看出，当外部激励使蹦床到达第一阶固有频率时，蹦床网面的中心区域会有剧烈的震动，且中心区域是人体对网面施加作用力的主要区域，蹦床模型的第一阶固有频率会对人体跳跃产生影响很大。人体想要使蹦床达到共振，自身跳跃半周期就必须等于0.36s(蹦床模型一阶固有频率的倒数的一半)。从表2中蹦床的第二阶、第三阶固有频率的来看，人体若想使蹦床在这两频率下共振，自身跳跃半周期就必须至少小于0.03s，这种高频振动显然是人体蹬网不可能达到的。因此，仅保留与第一阶固有频率对应的备选半周期0.36s。
先利用下式计算人体-蹦床系统的等效质量m1：得m1＝80.50kg，减去运动员质量43kg，即得蹦床等效质量m0＝37.50kg。
利用下式计算蹦床阻尼比例系数c0：得c0＝60.46。
人体跳跃使蹦床产生振动，蹦床结构中不可避免会存在阻尼的因素，若得出蹦床网面中心区域的阻尼以及蹦床在人体跳跃过程的等效质量，蹦床网面中心区域的受力表达式即可写成下式：Fc＝m0y”+c0y’+k0y。目前，m0、c0、k0均已确定，利用ANSYS根据运动员质量及保留的备选半周期0.36s，模拟人体与蹦床接触点的位移轨迹，以推算得出y”、y’、y，进而可得到冲击力Fc的时程曲线。结果表明，冲击力Fc的极小值大于载荷下限值1.8KN。
因此，备选半周期0.36s即为蹦床运动员起跳动作的推荐跳跃半周期。
目标运动员在实验蹦床上完成“811团”动作的影像解析获得的数据是0.38s，显然，该运动员的网上起跳动作技术没有达到最佳化的效果。

Claims (10)

1.一种基于网面有限元建模的蹦床动态应力测试分析方法，包括以下步骤：
第一步、在蹦床网面中心区域进行静载实验得到静载实验数据，并计算出蹦床网面中心点等效刚度k0；
第二步、在蹦床网面中心区域进行动载实验得到动载实验数据，并计算出蹦床阻尼比ξ0；
第三步、采用ANSYS软件，针对蹦床网面建立有限元模型；
其特征是，还包括：
第四步、采用ANSYS软件，在蹦床网面有限元模型中：
S1、加载温度载荷，使蹦床网面具有预应力；
S2、在蹦床网面中心区域按第一步实验条件模拟进行静载实验，并得出一组与多个载荷分别对应的网面位移模拟平均值；
S3、将各网面位移模拟平均值与相应的第一步实际测得的网面位移值进行比较和计算，并得出偏差率-载荷分布图；
S4、根据预设的最高偏差率限定，结合偏差率-载荷分布图得出载荷限定区间或载荷下限值或载荷上限值；
第五步、采用ANSYS软件，在蹦床网面有限元模型中：
T1、先重复第四步S1，再在蹦床网面中心区域加载运动员质量，并计算得出一组蹦床网面N阶固有频率，N为正整数；其中，蹦床网面N阶固有频率至少包括蹦床网面一阶固有频率；
T2、将从人体开始接触蹦床网面到人体离开蹦床网面的时间作为跳跃半周期；以蹦床网面N阶固有频率的倒数的一半作为备选半周期；逐一判断各备选半周期是否为人体蹬网能实现的跳跃半周期，若否则剔除该备选半周期，若是则保留该备选半周期；
T3、采用与保留的备选半周期对应的蹦床网面N阶固有频率、及蹦床网面中心点等效刚度k0，按预定算法得出蹦床等效质量m0，然后结合蹦床阻尼比ξ0，按预定算法得出蹦床阻尼比例系数c0；根据运动员质量及保留的各备选半周期，分别模拟人体与蹦床接触点的位移轨迹；
T4、人体蹬网时对蹦床网面的冲击力Fc＝m0y”+c0y’+k0y，其中，y”、y’、y分别表示人体与蹦床接触点垂直方向的加速度、速度以及位移，并由所述位移轨迹推算得出；
T5、当第四步S4所得的是载荷限定区间时，按第五步T4计算冲击力Fc的极大值和极小值，并判断该极大值和极小值是否均属于载荷限定区间，若是则保留与该冲击力Fc对应的备选半周期，若两者之一不属于或两者均不属于载荷限定区间则剔除该备选半周期；
当第四步S4所得的是载荷下限值时，按第五步T4计算冲击力Fc的极小值，并判断该极小值是否大于或等于载荷下限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
当第四步S4所得的是载荷上限值时，按第五步T4计算冲击力Fc的极大值，并判断该极大值是否小于或等于载荷上限值，若是则保留与该冲击力Fc对应的备选半周期，若否则剔除该备选半周期；
T6、将保留的备选半周期作为蹦床运动员起跳动作的推荐跳跃半周期。
2.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第一步中，所述静载实验采用若干杠铃片加载于蹦床网面，各杠铃片中心对准蹦床网面中心；以蹦床网面中心点为原点，在蹦床网面建立直角坐标系；以杠铃片边缘与两坐标轴的四个交点、及蹦床网面中心点作为测量点；以杠铃片覆盖的蹦床网面区域作为蹦床网面中心区域。
3.根据权利要求2所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第一步中，静载实验的具体过程如下：
在蹦床网面中心区域加载预设初始质量的杠铃片，测量各测量点的纵向位移；继续多次加载预设质量的杠铃片，直至加载杠铃片的总质量达到预设质量上限，每次加载后测量各测量点的纵向位移；所述纵向位移在蹦床网面未加载杠铃片时为零。
4.根据权利要求3所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第一步中，静载实验计算蹦床网面中心点等效刚度k0的具体过程如下：
将加载杠铃片的总质量作为载荷Fy，并将其与蹦床网面中心点的纵向位移y的关系进行线性拟合，所得式为Fy＝Ay+B，A、B分别为常数；其中常数A即为蹦床网面中心点等效刚度k0。
5.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第二步中，动载实验的具体过程如下：
将预设质量的小球固定在蹦床网面中心点上，蹦床网面中心点即为小球与蹦床网面的接触点；对小球施加垂直向下的激励，使小球在垂直方向上振动；通过高速摄像机记录小球在震动过程中的最大振幅，并通过影像解析确定小球在每个衰减周期的振幅。
6.根据权利要求5所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第二步中，动载实验计算蹦床阻尼比ξ0的具体过程如下：
先利用下式计算对数减缩系数Λ：
其中，A1为最大振幅，Aj+1为振动开始后小球的第j+1个振幅，N为小球振幅的总数量；
再利用下式计算蹦床阻尼比ξ0：
7.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第三步中，建立有限元模型时，蹦床网面的结构采用空间铰接体系，结构边界为点支撑铰接或者线支撑铰接边界，且蹦床网面保持弹性状态。
8.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第四步S1中，利用下式计算温度载荷Δt：
其中，T为蹦床网面的初始预应力，E为蹦床网面的弹性模量，a为蹦床网面的膨胀系数，tn为蹦床网面的厚度。
9.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第五步T3中，先利用下式计算人体-蹦床系统的等效质量m1：其中，f为蹦床网面N阶固有频率；再以人体-蹦床系统的等效质量m1减去运动员质量，即得蹦床等效质量m0。
10.根据权利要求1所述的基于网面有限元建模的蹦床动态应力测试分析方法，其特征是，第五步T3中，利用下式计算蹦床阻尼比例系数c0：
CN201610148900.XA 2016-03-16 2016-03-16 基于网面有限元建模的蹦床动态应力测试分析方法 Expired - Fee Related CN105825010B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610148900.XA CN105825010B (zh) 2016-03-16 2016-03-16 基于网面有限元建模的蹦床动态应力测试分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610148900.XA CN105825010B (zh) 2016-03-16 2016-03-16 基于网面有限元建模的蹦床动态应力测试分析方法

Publications (2)

Publication Number Publication Date
CN105825010A CN105825010A (zh) 2016-08-03
CN105825010B true CN105825010B (zh) 2018-11-27

Family

ID=56987769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610148900.XA Expired - Fee Related CN105825010B (zh) 2016-03-16 2016-03-16 基于网面有限元建模的蹦床动态应力测试分析方法

Country Status (1)

Country Link
CN (1) CN105825010B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090298B (zh) * 2017-12-29 2021-06-18 青岛海聚仿真软件技术有限公司 一种基于有限元分析的蹦床跳布找形方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202410090U (zh) * 2011-11-28 2012-09-05 宋雅伟 蹦床运动动态监测系统
CA2849327A1 (en) * 2014-04-07 2015-10-07 Prohi Enterprise Co., Ltd. Improvement of trampoline net structure

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202410090U (zh) * 2011-11-28 2012-09-05 宋雅伟 蹦床运动动态监测系统
CA2849327A1 (en) * 2014-04-07 2015-10-07 Prohi Enterprise Co., Ltd. Improvement of trampoline net structure

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Competitive trampolining influence trabecular bone structure,;Lauren A. Burt, John D. Schipilow, Steven K. Boyd;《Journal of Sport and Health Science》;20150430(第4期);第469-475页 *
蹦床的有限元建模及最佳起跳周期的模拟计算;钱竞光 ， 吕超 ， 李广凯 ， 毛一玲 ， 李兆霞;《全国竞技体育科学论文报告会》;20140813;第445-446页 *
蹦床网面的载荷实验和有限元建模;钱竞光，吕超，毛一玲，李兆霞;《第16届全国运动生物力学学术交流大会》;20131022;第95-96页 *

Also Published As

Publication number Publication date
CN105825010A (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
CN107515972B (zh) 一种扬声器跌落过程的数值仿真分析方法
JP5787598B2 (ja) 複合ヘッドのモード減衰比を予測する方法
CN104866652A (zh) 一种基于abaqus的喷丸强化变形的有限元模拟方法
Huisman et al. Monte Carlo study of multiply crosslinked semiflexible polymer networks
CN110532727B (zh) 可用于常见非牛顿流体的数值模拟方法
CN106802969A (zh) 阻尼材料动态特性的验证系统及其验证方法
Yabuno et al. Nonlinear normal modes of a parametrically excited cantilever beam
Mao et al. Modelling and simulation of the straw-grain separation process based on a discrete element model with flexible hollow cylindrical bonds
CN110188417A (zh) 一种基于多级超单元的复杂结构有限元模型修正方法
CN110258435A (zh) 一种真实应力场下百米级高土石坝动力试验方法
CN105825010B (zh) 基于网面有限元建模的蹦床动态应力测试分析方法
JP5892257B2 (ja) シミュレーションプログラム、シミュレーション方法及びシミュレーション装置
Wang et al. Numerical study on the dynamic properties of wrinkled membranes
CN107657073A (zh) 一种三明治复合材料结构的参数识别方法
CN107862119A (zh) 一种振动可靠性测试方法
Esser et al. An extended finite element method applied to levitated droplet problems
Songsuwan et al. Nonlinear transient response of sandwich beams with functionally graded porous core under moving load
JP5241469B2 (ja) シミュレーション方法及びプログラム
Deng et al. Impact peak force measurement of potato
Dolfen et al. A multi-stage approach of simulating turbulence-induced vibrations in a wire-wrapped tube bundle for fretting wear prediction
US8229711B2 (en) Method of analysis for kinetic properties of golf club head and golf club therefor
Xu et al. Finite element analysis modeling research on the compression process of cotton fiber assembly
CN109186967A (zh) 一种基于bp神经网络的橡胶减振器性能预测及选型方法
CN105404702B (zh) 周视观察镜设备结构减振优化设计方法
Nasrollahi et al. A nondestructive evaluation approach to characterize tennis balls

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181127

CF01 Termination of patent right due to non-payment of annual fee