冰雹冲击温室屋面玻璃荷载的数值模拟研究_龙彦文.pdf

第 3 4 卷 第 2 期 2 0 2 4 年 6 月 V o l 3 4 N o 2 J u n e 2 0 2 4 第 3 4 卷 第 2 期 2 0 2 4 年 6 月 V o l 3 4 N o 2 J u n e 2 0 2 4 湖 南 工 程 学 院 学 报 自 然 科 学 版 J o u r n a l o f H u n a n I n s t i t u t e o f E n g i n e e r i n g N a t u r a l S c i e n c e E d i t i o n 湖 南 工 程 学 院 学 报 自 然 科 学 版 J o u r n a l o f H u n a n I n s t i t u t e o f E n g i n e e r i n g N a t u r a l S c i e n c e E d i t i o n 冰 雹 冲 击 温 室 屋 面 玻 璃 荷 载 的 数 值 模 拟 研 究 龙 彦 文 戴 益 民 杨 伟 元 湖 南 科 技 大 学 土 木 工 程 学 院 湘 潭 4 1 1 2 0 1 收 稿 日 期 2 0 2 3 1 0 2 4 基 金 项 目 国 家 自 然 科 学 基 金 面 上 项 目 5 2 1 7 8 4 7 8 湖 南 省 自 然 科 学 基 金 面 上 项 目 2 0 2 2 J J 3 0 2 4 7 作 者 简 介 龙 彦 文 1 9 9 9 男 硕 士 研 究 生 研 究 方 向 结 构 风 致 及 风 雹 灾 害 等 通 信 作 者 戴 益 民 1 9 7 2 男 博 士 教 授 博 士 生 导 师 研 究 方 向 结 构 风 致 及 风 雹 灾 害 等 摘 要 为 揭 示 冰 雹 冲 击 玻 璃 温 室 屋 面 受 载 机 理 采 用 L S P r e P o s t 的 随 动 强 化 材 料 M A T 3 进 行 建 模 结 合 考 虑 冰 雹 下 坠 速 度 粒 径 大 小 及 屋 面 坡 角 因 素 开 展 自 然 条 件 下 冰 雹 低 速 冲 击 玻 璃 温 室 屋 面 荷 载 研 究 结 果 表 明 冰 雹 粒 径 对 屋 面 冲 击 能 量 影 响 明 显 冲 击 能 量 随 冰 雹 粒 径 增 加 而 非 线 性 增 加 随 屋 面 坡 角 增 加 冲 击 能 量 减 少 且 坡 角 在 2 0 2 5 之 间 有 效 应 力 达 到 最 大 值 研 究 结 果为量大面广的温室大棚屋面坡角设计及防灾提供了一定价值与参考 关 键 字 玻 璃 温 室 坡 角 冰 雹 冲 击 中 图 分 类 号 T U 5 0 2 6 文 献 标 志 码 A 文 章 编 号 1 6 7 1 1 1 9 X 2 0 2 4 0 2 0 0 9 0 0 5 0 引 言 农 业 高 质 量 发 展 是 国 家 农 业 供 给 侧 结 构 性 改 革 的 重 要 目 标 温 室 发 展 形 势 是 设 施 农 业 工 程 发 展 水 平 的 重 要 体 现 目 前 我 国 现 代 都 市 型 农 业 发 展 规 模 不 断 扩 大 以 荷 兰 V e n l o 型 玻 璃 温 室 为 代 表 的 现 代 化 农 业 温 室 已 经 成 为 使 用 面 积 最 大 数 量 最 多 应 用 最 广 的 温 室 结 构 1 近 年 来 强 对 流 天 气 频 发 全 球 范 围 内 的 风 雹 灾 害 对 人 类 造 成 巨 大 损 害 2 温 室 大 棚 结 构 也 因 大 风 或 风 雹 致 使 农 作 物 损 失 严 重 而 目 前 国 内 关 于 温 室 大 棚 安 全 性 的 研 究 主 要 集 中 在 风 雪 荷 载 对 冰 雹 冲 击 玻 璃 温 室 大 棚 的 相 关 研 究 基 本 为 空 白 玻 璃 材 料 在 冲 击 荷 载 下 呈 现 脆 性 破 坏 特 征 其 对 表 面 损 伤 非 常 敏 感 因 此 有 必 要 研 究 玻 璃 构 件 的 抗 冲 击 能 力 以 及 冲 击 作 用 对 构 件 损 伤 产 生 的 承 载 力 折 减 规 律 3 为 研 究 冰 雹 冲 击 的 规 律 特 征 张 晓 晴 等 4 使 用 A b a q u s 建 立 了 复 合 材 料 加 筋 结 构 有 限 元 模 型 从 多 方 面 预 测 了 冰 雹 冲 击 对 复 合 材 料 加 筋 板 的 冲 击 动 力 响 应 影 响 K i m 等 5 使 用 氮 气 炮 发 射 不 同 直 径 的 人 造 冰 雹 s i m u l a t e d h a i l i c e S H I 冲 击 碳 环 氧 复 合 材 料 板 得 到 了 复 合 材 料 板 的 损 伤 模 式 丁 敏 等 6 通 过 使 用 A N S Y S 的 S o l i d 1 6 4 单 元 建 立 冰 雹 冲 击 温 室 玻 璃 板 的 数 值 模 型 系 统 分 析 了 温 室 玻 璃 在 冰 雹 冲 击 下 的 破 坏 情 景 李 静 7 通 过 L S D Y N A 对 冰 雹 高 速 冲 击 铝 合 金 板 进 行 了 数 值 模 拟 取 得 了 理 想 的 效 果 王 星 尔 8 对 建 筑 夹 层 玻 璃 在 刚 体 冲 击 下 的 脆 性 破 坏 特 征 进 行 分 析 得 到 了 夹 层 玻 璃 的 破 坏 机 理 和 冲 击 吸 能 机 制 戴 益 民 等 9 进 行 了 风 雹 耦 合 机 理 试 验 得 到 了 不 同 粒 径 冰 雹 对 光 伏 结 构 冲 击 力 峰 值 规 律 并 优 化 建 立 了 G A B P 神 经 网 络 用 于 预 测 风 雹 下 单 颗 粒 冰 雹 冲 击 力 综 上 国 内 外 学 者 大 多 从 航 空 航 天 领 域 出 发 研 究 冰 雹 高 速 冲 击 构 件 的 破 坏 机 理 而 实 际 建 筑 结 构 遭 受 的 风 雹 灾 害 其 冰 雹 着 陆 速 度 一 般 不 超 过 5 0 m s 且 冰 雹 是 一 种 率 相 关 材 料 低 速 冲 击 下 的 材 料 性 质 与 高 速 冲 击 时 不 同 主 要 体 现 在 压 缩 屈 服 强 度 上 1 0 随 着 国 内 风 雹 灾 害 频 率 逐 年 递 增 趋 势 而 关 于 冰 雹 冲 击 温 室 大 棚 机 理 研 究 基 本 空 白 因 此 研 究 自 然 下 落 的 冰 雹 对 温 室 玻 璃 的 冲 击 作 用 有 重 要 意 义 本 文 综 合 冰 雹 粒 径 下 坠 速 度 及 屋 面 坡 角 DOI 10 15987 ki hgbjbz 2024 02 013第 2 期 影 响 因 素 建 立 了 不 同 直 径 的 冰 雹 以 对 应 的 着 陆 末 速 度 冲 击 不 同 坡 角 的 温 室 玻 璃 板 的 数 值 模 型 探 讨 坡 角 的 变 化 对 玻 璃 板 受 冲 击 能 力 的 影 响 以 及 定 量 分 析 不 同 直 径 冰 雹 的 冲 击 作 用 规 律 1 有 限 元 模 型 1 1 试 验 模 型 参 数 由 相 关 研 究 所 知 温 室 玻 璃 板 平 面 尺 寸 长 宽 比 在 1 5 左 右 是 比 较 理 想 的 数 值 范 围 6 因 此 本 文 玻 璃 板 平 板 尺 寸 选 择 1 5 m 1 m 厚 度 取 5 m m 自 然 界 冰 雹 形 状 各 异 一 般 以 球 形 和 椭 球 形 居 多 1 1 1 2 数 据 统 计 资 料 表 明 实 际 冰 雹 直 径 在 0 5 5 c m 之 间 而 直 径 过 小 的 冰 雹 带 来 的 危 害 有 限 6 因 此 本 次 试 验 的 冰 雹 直 径 选 择 2 5 5 c m 之 间 冰 雹 的 自 由 下 落 过 程 中 受 到 重 力 和 阻 力 的 共 同 作 用 两 者 平 衡 时 冰 雹 速 度 达 到 最 大 即 冰 雹 的 末 速 度 1 1 本 文 称 之 为 冰 雹 着 陆 速 度 冰 雹 下 落 时 的 阻 力 主 要 来 自 空 气 阻 力 阻 力 大 小 与 速 度 和 垂 直 下 落 方 向 的 最 大 截 面 积 有 关 计 算 公 式 如 式 1 所 示 冰 雹 的 重 力 与 其 体 积 大 小 有 关 而 实 际 冰 雹 一 般 以 球 形 和 椭 球 形 居 多 为 简 化 计 算 本 文 冰 雹 模 型 形 状 设 置 为 球 形 冰 雹 的 重 力 计 算 可 由 式 2 计 算 得 出 F C d S a 2 V 2 1 式 中 C d 为 空 气 阻 力 系 数 S 为 与 下 落 方 向 垂 直 的 冰 雹 最 大 截 面 积 a 为 空 气 的 密 度 V 为 冰 雹 相 对 于 空 气 的 速 度 G h d 3 g 6 2 式 中 h 为 冰 球 的 密 度 d 为 冰 球 的 直 径 g 为 重 力 加 速 度 由 相 关 研 究 1 2 1 3 得 到 球 形 冰 雹 的 末 速 度 v 0 计 算 公 式 如 式 3 所 示 v 0 1 8 1 0 2 d 0 4 0 7 5 h a g 0 5 d 5 0 1 5 2 1 0 2 1 0 0 1 8 d 4 2 h a g d 5 0 d 6 4 3 4 9 1 0 2 d 0 3 7 3 5 h a g 6 4 d 1 6 3 本 次 试 验 以 玻 璃 温 室 屋 面 的 坡 角 为 控 制 变 量 试 验 参 数 见 表 1 表 1 试 验 参 数 图 序 号 1 2 3 4 5 冰 雹 直 径 c m 2 5 3 0 3 5 4 0 5 0 末 速 度 m s 2 3 4 1 2 5 2 2 2 6 8 5 2 8 3 6 3 1 0 6 玻 璃 板 尺 寸 1 5 m 1 m 5 m m 1 5 m 1 m 5 m m 1 5 m 1 m 5 m m 1 5 m 1 m 5 m m 1 5 m 1 m 5 m m 坡 角 0 1 0 2 0 2 5 3 0 0 1 0 2 0 2 5 3 0 0 1 0 2 0 2 5 3 0 0 1 0 2 0 2 5 3 0 0 1 0 2 0 2 5 3 0 1 2 单 元 模 型 L S D Y N A 是 世 界 上 著 名 的 通 用 显 示 动 力 分 析 程 序 特 别 适 合 求 解 爆 炸 冲 击 问 题 后 被 A N S Y S 公 司 收 购 故 本 次 试 验 可 采 用 A N S Y S 内 置 的 L S D Y N A 求 解 器 进 行 计 算 建 模 在 L S P r e P o s t 中 进 行 玻 璃 平 板 因 其 中 两 个 维 度 尺 寸 远 远 大 于 厚 度 方 向 的 维 度 因 此 选 择 壳 单 元 类 型 壳 单 元 厚 度 选 择 玻 璃 板 实 际 厚 度 5 m m 对 于 壳 单 元 算 法 有 很 多 从 计 算 资 源 和 精 度 考 虑 本 次 算 法 选 择 默 认 的 B e l y t s c h k o T s a y 单 元 算 法 该 算 法 基 于 M i n d l i n R e i s s n e r 假 设 包 括 横 向 剪 切 力 易 收 敛 精 度 高 但 它 不 能 处 理 翘 曲 因 此 不 能 使 用 粗 网 格 1 3 本 次 试 验 网 格 尺 寸 为 1 m m 见 图 1 冰 雹 是 实 心 球 体 单 元 类 型 选 择 S o l i d 六 面 体 单 元 算 法 选 择 一 点 的 积 分 加 上 黏 性 沙 漏 控 制 本 次 试 验 单 元 数 量 为 7 0 0 0 见 图 2 图 1 玻 璃 平 板 图 2 冰 雹 龙 彦 文 等 冰 雹 冲 击 温 室 屋 面 玻 璃 荷 载 的 数 值 模 拟 研 究 9 12 0 2 4 年 湖 南 工 程 学 院 学 报 自 然 科 学 版 1 3 材 料 模 型 L S D Y N A 最 为 著 名 的 便 是 它 具 有 丰 富 的 材 料 库 能 够 足 够 真 实 地 模 拟 真 实 世 界 中 的 各 种 材 料 特 性 1 4 玻 璃 温 室 的 屋 面 覆 盖 玻 璃 多 采 用 浮 法 玻 璃 1 属 于 无 机 玻 璃 的 一 种 无 机 玻 璃 是 脆 性 材 料 在 静 态 荷 载 作 用 下 表 现 出 明 显 的 线 弹 性 力 学 行 为 但 在 动 态 荷 载 作 用 下 无 机 玻 璃 会 表 现 出 一 定 的 非 线 性 力 学 行 为 如 率 相 关 性 受 压 时 可 进 入 塑 性 阶 段 等 1 5 主 要 观 察 冰 雹 动 能 在 玻 璃 中 产 生 的 应 力 波 传 递 状 态 因 此 选 择 随 动 塑 性 材 料 即 L S D Y N A 材 料 库 的 M A T 3 冰 雹 是 一 种 应 变 率 相 关 材 料 在 高 速 冲 击 下 会 出 现 流 体 特 征 而 自 然 下 落 的 冰 雹 属 于 低 速 冲 击 因 此 冰 雹 材 料 也 选 择 随 动 塑 性 材 料 M A T 3 具 体 参 数 见 表 2 表 2 材 料 参 数 1 6 1 7 材 料 玻 璃 冰 雹 厚 度 m m 5 密 度 k g m 3 2 4 0 0 9 0 0 弹 性 模 量 G P a 7 0 9 3 8 泊 松 比 0 2 5 0 3 3 屈 服 应 力 M P a 5 8 8 5 2 1 4 控 制 参 数 温 室 的 覆 盖 玻 璃 采 用 铝 合 金 型 材 镶 嵌 式 安 装 且 由 于 玻 璃 自 身 提 供 一 定 刚 度 考 虑 玻 璃 蒙 皮 作 用 的 影 响 可 简 化 为 刚 接 形 式 采 用 L S D Y N A 关 键 字 B O U N D A R Y S P C S E T 约 束 玻 璃 板 四 周 节 点 的 全 部 自 由 度 1 冰 雹 冲 击 为 软 体 冲 击 事 件 1 0 属 于 柔 体 对 刚 体 这 一 类 接 触 问 题 试 验 接 触 方 式 选 择 面 面 接 触 A S T S 定 义 主 面 为 玻 璃 从 面 为 冰 雹 L S D Y N A 关 键 字 为 C O N T A C T A U T O M A T I C S U R F A C E T O S U R F A C E 本 次 L S D Y N A 试 验 分 析 的 单 位 制 为 k g m s 时 间 步 长 控 制 为 0 0 0 0 1 控 制 关 键 字 为 D A T A B A S E B I N A R Y D 3 P L O T 终 止 时 间 设 置 为 0 0 1 控 制 关 键 字 为 C O N T R O L T E R M I N A T I O N 2 结 果 分 析 2 1 冲 击 过 程 分 析 本 次 试 验 冰 雹 冲 击 点 位 置 为 玻 璃 几 何 中 心 位 置 以 3 5 c m 直 径 冰 雹 冲 击 0 屋 面 玻 璃 板 为 例 在 0 0 0 0 5 s 时 冰 雹 与 玻 璃 刚 接 触 时 接 触 点 位 置 开 始 产 生 有 效 应 力 2 5 9 5 M P a 随 即 在 0 0 0 0 1 s 后 达 到 本 次 冲 击 过 程 中 的 最 大 值 6 5 5 2 M P a 随 后 应 力 波 以 圆 形 的 形 式 从 局 部 扩 散 到 全 局 随 后 在 接 触 边 界 时 返 回 随 后 最 大 有 效 应 力 在 1 0 M P a 左 右 波 动 应 力 波 扩 散 云 图 见 图 3 a 0 0 0 0 5 s b 0 0 0 0 6 s c 0 0 0 0 7 s d 0 0 0 0 8 s e 0 0 0 0 9 s f 0 0 0 1 2 s g 0 0 0 1 5 s h 0 0 0 4 4 s 图 3 3 5 c m 冰 雹 冲 击 0 屋 面 玻 璃 板 应 力 波 扩 散 云 图 2 2 坡 角 影 响 分 析 冰 雹 竖 向 冲 击 不 同 坡 度 的 屋 面 覆 盖 玻 璃 会 产 生 能 量 传 递 图 4 a e 表 明 不 同 直 径 冰 雹 冲 击 下 横 坐 标 为 时 间 单 位 s 纵 坐 标 为 玻 璃 板 总 能 量 单 位 k g m 2 s 2 总 能 量 动 能 内 能 会 在 冲 击 作 用 后 短 暂 时 间 内 达 到 幅 值 后 缓 慢 衰 减 且 在 各 种 直 径 的 冰 雹 冲 击 下 屋 面 坡 角 在 0 时 能 量 的 幅 值 点 最 大 后 随 着 坡 角 的 增 加 幅 值 呈 现 非 线 性 下 降 趋 势 图 4 f 为 不 同 直 径 冰 雹 冲 击 不 同 坡 度 屋 面 玻 9 2第 2 期 龙 彦 文 等 冰 雹 冲 击 温 室 屋 面 玻 璃 荷 载 的 数 值 模 拟 研 究 璃 板 的 最 大 有 效 应 力 变 化 图 横 坐 标 为 坡 角 度 纵 坐 标 为 归 一 化 无 量 纲 最 大 有 效 应 力 值 如 图 4 所 示 当 冰 雹 直 径 为 2 5 c m 3 c m 3 5 c m 4 c m 时 曲 线 化 趋 势 一 致 随 着 坡 角 的 增 加 应 力 大 小 递 增 在 2 0 或 2 5 时 达 到 幅 值 而 后 下 降 在 5 c m 直 径 的 冰 雹 冲 击 下 坡 角 为 0 时 幅 值 最 大 而 后 的 曲 线 趋 势 与 前 者 相 似 也 在 2 0 时 达 到 幅 值 后 衰 减 0 2 5 0 2 0 0 1 5 0 1 0 0 0 5 0 0 0 玻 璃 板 总 能 量 k g m 2 s 2 时 间 s 0 0 0 0 0 0 0 2 0 0 0 4 0 0 0 6 0 0 0 8 0 0 1 0 0 1 0 2 0 2 5 3 0 时 间 s 0 0 0 0 0 0 0 2 0 0 0 4 0 0 0 6 0 0 0 8 0 0 1 0 0 8 0 6 0 5 0 4 0 3 0 2 0 1 0 0 玻 璃 板 总 能 量 k g m 2 s 2 0 1 0 2 0 2 5 3 0 a 2 5 c m 冰 雹 冲 击 玻 璃 平 板 b 3 0 c m 冰 雹 冲 击 玻 璃 平 板 0 1 0 2 0 2 5 3 0 时 间 s 0 0 0 0 0 0 0 2 0 0 0 4 0 0 0 6 0 0 0 8 0 0 1 0 1 7 5 1 5 0 1 2 5 1 0 0 0 7 5 0 5 0 0 2 5 0 0 0 玻 璃 板 总 能 量 k g m 2 s 2 0 0 0 0 0 0 0 2 0 0 0 4 0 0 0 6 0 0 0 8 0 0 1 0 时 间 s 3 5 3 0 2 5 2 0 1 5 1 0 0 5 0 0 玻 璃 板 总 能 量 k g m 2 s 2 0 1 0 2 0 2 5 3 0 c 3 5 c m 冰 雹 冲 击 玻 璃 平 板 d 4 0 c m 冰 雹 冲 击 玻 璃 平 板 时 间 s 0 0 0 0 0 0 0 2 0 0 0 4 0 0 0 6 0 0 0 8 0 0 1 0 1 0 8 6 4 2 0 玻 璃 板 总 能 量 k g m 2 s 2 0 1 0 2 0 2 5 3 0 有 效 应 力 归 一 化 指 标 屋 面 坡 角 度 0 2 4 0 2 2 0 2 0 0 1 8 0 1 6 0 5 1 0 1 5 2 0 2 5 3 0 2 5 c m 冰 雹 3 0 c m 冰 雹 3 5 c m 冰 雹 4 0 c m 冰 雹 5 0 c m 冰 雹 e 5 0 c m 冰 雹 冲 击 玻 璃 平 板 f 最 大 有 效 应 力 图 4 不 同 屋 面 坡 角 能 量 有 效 应 力 图 3 结 论 本 文 考 虑 了 自 然 界 冰 雹 的 着 陆 末 速 度 和 冰 雹 尺 寸 的 关 系 基 于 L S D Y N A 建 立 了 不 同 尺 寸 冰 雹 冲 击 不 同 坡 度 屋 面 覆 盖 玻 璃 的 有 限 元 模 型 结 合 冰 雹 和 玻 璃 材 料 的 应 变 率 敏 感 特 性 分 析 了 冰 雹 冲 击 过 程 中 玻 璃 内 部 应 力 的 发 展 规 律 得 到 了 以 下 结 论 1 冰 雹 冲 击 过 程 中 玻 璃 内 部 有 效 应 力 最 大 幅 值 出 现 在 接 触 之 后 的 瞬 间 呈 现 突 变 特 征 2 在 相 同 坡 度 下 玻 璃 板 受 到 冰 雹 冲 击 的 能 量 大 小 随 着 冰 雹 粒 径 增 加 而 增 加 3 在 相 同 直 径 的 冰 雹 冲 击 下 玻 璃 板 受 到 冰 雹 冲 击 的 能 量 大 小 随 着 坡 度 的 增 加 而 减 少 9 32 0 2 4 年 湖 南 工 程 学 院 学 报 自 然 科 学 版 4 玻 璃 板 受 到 冰 雹 冲 击 的 局 部 最 大 有 效 应 力 随 着 坡 度 的 增 加 呈 现 上 三 角 形 趋 势 在 2 0 或 2 5 时 达 到 峰 值 参 考 文 献 1 李 一 哲 骈 超 荷 兰 V e n l o 型 温 室 结 构 抗 震 性 能 分 析 J 农 业 工 程 学 报 2 0 2 1 3 7 1 9 2 4 3 2 4 9 2 戴 益 民 徐 瑛 李 怿 歆 等 工 程 结 构 抗 冰 雹 冲 击 研 究 进 展 及 展 望 J 湖 南 大 学 学 报 自 然 科 学 版 2 0 2 3 5 0 1 2 2 8 2 3 6 3 黄 小 坤 段 树 坤 刘 强 等 玻 璃 结 构 研 究 进 展 与 工 程 实 践 J 建 筑 结 构 学 报 2 0 2 0 4 1 6 1 2 0 4 张 晓 晴 丁 铁 龙 舒 畅 等 复 合 材 料 加 筋 壁 板 的 抗 冰 雹 冲 击 动 力 响 应 及 损 伤 预 测 J 华 南 理 工 大 学 学 报 自 然 科 学 版 2 0 1 7 4 5 5 1 2 0 1 2 8 5 K I M H W E L C H D A K E D W A R D K T E x p e r i m e n t a l i n v e s t i g a t i o n o f h i g h v e l o c i t y i c e i m p a c t s o n w o v e n c a r b o n e p o x y c o m p o s i t e p a n e l s J C o m p o s i t e s P a r t A A p p l i e d S c i e n c e a n d M a n u f a c t u r i n g 2 0 0 3 3 4 1 2 5 4 1 6 丁 敏 施 旭 栋 周 蕊 等 冰 雹 冲 击 下 温 室 玻 璃 破 损 情 景 数 值 模 拟 J 农 业 工 程 学 报 2 0 1 2 2 8 1 2 0 2 2 0 7 2 9 8 7 李 静 基 于 L S D Y N A 的 冰 雹 冲 击 铝 合 金 平 板 数 值 模 拟 D 南 京 东 南 大 学 2 0 1 5 8 王 星 尔 刚 体 冲 击 作 用 下 建 筑 夹 层 玻 璃 的 动 态 响 应 和 破 坏 行 为 研 究 D 上 海 上 海 交 通 大 学 2 0 1 8 9 戴 益 民 李 怿 歆 徐 瑛 等 基 于 G A B P 神 经 网 络 的 风 雹 耦 合 所 致 冰 雹 冲 击 力 预 测 O L 工 程 力 学 1 1 0 2 0 2 3 1 0 0 7 1 0 张 超 方 鑫 刘 建 春 复 合 材 料 层 板 冰 雹 高 速 冲 击 损 伤 预 测 及 失 效 分 析 J 北 京 航 空 航 天 大 学 学 报 2 0 2 2 4 8 4 6 9 8 7 0 7 1 1 徐 家 骝 冰 雹 阻 力 系 数 的 半 经 验 关 系 和 末 速 度 的 近 似 公 式 J 兰 州 大 学 学 报 1 9 7 8 1 1 9 0 1 0 3 1 2 王 凡 冰 雹 和 暴 雨 荷 载 作 用 下 正 交 异 性 鞍 形 膜 结 构 的 动 力 响 应 研 究 D 成 都 成 都 理 工 大 学 2 0 2 0 1 3 胡 仁 喜 康 士 廷 A N S Y S 工 程 应 用 系 列 A N S Y S 1 8 0 L S D Y N A 非 线 性 有 限 元 分 析 实 例 指 导 教 程 第 四 版 M 北 京 机 械 工 业 出 版 社 2 0 1 8 1 4 L S T C L s d y n a k e y w o r d u s e r s m a n u a l M L i v e r m o r e C a l i f o r n i a L i v e r m o r e S o f t w a r e T e c h n o l o g y A n A N S Y S C o m p a n y 2 0 2 1 1 5 胡 宗 文 袁 国 青 基 于 S P H 模 型 的 冰 雹 高 速 冲 击 复 合 材 料 结 构 响 应 规 律 研 究 J 力 学 季 刊 2 0 1 6 3 7 1 1 0 4 1 1 3 1 6 丁 敏 朱 丹 许 晶 等 风 雪 荷 载 作 用 下 V e n l o 型 温 室 结 构 整 体 性 能 研 究 J 中 国 农 业 大 学 学 报 2 0 1 7 2 2 1 1 2 0 1 2 8 1 7 T I P P M A N N J D K I M H R H Y M E R J D E x p e r i m e n t a l l y v a l i d a t e d s t r a i n r a t e d e p e n d e n t m a t e r i a l m o d e l f o r s p h e r i c a l i c e i m p a c t s i m u l a t i o n J I n t e r n a t i o n a l J o u r n a l o f I m p a c t E n g i n e e r i n g 2 0 1 3 5 7 1 4 3 5 4 N u m e r i c a l S i m u l a t i o n S t u d y o n t h e I m p a c t o f H a i l o n G r e e n h o u s e R o o f G l a s s L O N G Y a n w e n D A I Y i m i n Y A N G W e i y u a n S c h o o l o f C i v i l E n g i n e e r i n g H u n a n U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y X i a n g t a n 4 1 1 2 0 1 C h i n a A b s t r a c t I n o r d e r t o r e v e a l t h e l o a d b e a r i n g m e c h a n i s m o f h a i l i m p a c t o n t h e g l a s s g r e e n h o u s e r o o f t h i s p a p e r u s e s L S P r e P o s t s f o l l o w u p s t r e n g t h e n i n g m a t e r i a l M A T 3 f o r m o d e l i n g a n d c o n s i d e r s t h e f a c t o r s o f h a i l f a l l i n g s p e e d p a r t i c l e s i z e a n d r o o f s l o p e a n g l e a r e t a k e n i n t o c o n s i d e r t i o n A s t u d y o n t h e r o o f l o a d o f h a i l i m p a c t i n g g l a s s g r e e n h o u s e s a t l o w s p e e d u n d e r n a t u r a l c o n d i t i o n s i s c o n d u c t e d T h e r e s u l t s s h o w t h a t h a i l p a r t i c l e s i z e h a s a s i g n i f i c a n t i m p a c t o n r o o f i m p a c t e n e r g y T h e i m p a c t e n e r g y i n c r e a s e s n o n l i n e a r l y w i t h t h e i n c r e a s e o f h a i l p a r t i c l e s i z e A s t h e r o o f s l o p e a n g l e i n c r e a s e s t h e i m p a c t e n e r g y d e c r e a s e s W h e n t h e s l o p e a n g l e i s b e t w e e n 2 0 a n d 2 5 t h e e f f e c t i v e s t r e s s r e a c h e s m a x i m u m v a l u e T h e r e s e a r c h r e s u l t s p r o v i d e c e r t a i n v a l u e a n d r e f e r e n c e f o r t h e r o o f s l o p e a n g l e d e s i g n a n d d i s a s t e r p r e v e n t i o n o f n u m e r o u s g r e e n h o u s e s K e y w o r d s g r e e n h o u s e g l a s s s l o p e a n g l e h a i l i m p a c t 9 4