工程力学建模_杨冠声

合集下载

RRR-RRP平面六杆Ⅱ级机构的运动学仿真

RRR-RRP平面六杆Ⅱ级机构的运动学仿真
( . a cJ n x 3 0 6 , hn ; . a c a g U i r t , a c a g,i g i 3 0 1 C ia 1 N nh n a g o g U i r t N n h n ,i g i 3 0 3 C i 2 N n h n nv s y N n h n J n x 3 0 3 , hn ) v sy a a ei a
机构 的运动 分 析 , 要是 获 得 机 构 中某 些 构件 的动态 仿真 。采用 MA L B对 机构 进 行仿 真 , 乎 主 TA 几
的位 移 、 角速 度 和 加速 度 , 以及 某 些 点 的轨 迹 , 速度 所 有 的构件 运 动 参 数 都 在 仿 真 模 型 的数 据 线 上 传 和加 速度 。它是机 械设计 及评 价机 械运 动和动 力性 输 , 只要将该 数 据 线 上 的 信 息 引入 Sm u 模 块 , iot 就
柄 、R R RⅡ级 杆组 、 R R PⅡ级 杆 组 三 个 基 本 模 组 的 运 动 学 数 学 模 型 , 用 其 组 成 机 构 杆 组 并 搭 建 平 面 连 杆 机 构 的 运 动 学 仿 真 利
模型 , 充分利用 MA L B的 Sm l k TA i ui 仿真模型数据可视化的特点, n 观察和分析其运动参数的变化 。 [ 中图分类号 ]T 3 19 P 9 . [ 文献标志码 ]A [ 文章编号]10 4 2 (00)4— 0 6— 5 0 1— 9 6 2 1 0 04 0
能 的基础 , 是 分析 现 有 机械 优 化 综合 新 机 械 的基 可 以观察 到该 运 动 参 数 是 如 何 变 化 的及 相 应 的数 也
本手段 。

据, 并可 以 图形 的形 式直 观地 表现 出来 。 为 了利用 Ma a t b仿 真 软 件 包 的数 值 积 分 算 法 l

液压缸缸筒径向锻造成形工艺仿真分析

液压缸缸筒径向锻造成形工艺仿真分析
Copyright©博看网. All Rights Reserved.
#*#$ 8&0%!!!!!!!!!!!!!!! !!!重 型 机 械
+'++
*!前言
液压缸广泛应用于工程机械中" 是工程机械
中的关键部件之一" 其性能优劣直接影响工程机
械的使用寿命和生产安全." -#/ & 目前工程机械用
液压缸缸筒的生产方式主要有两种" 一种是切削
图 )!缸筒第一道次径向锻造
图 &!液压缸缸筒径向锻造有限元模型
该工艺采用室温下冷锻的方式" 分 $ 个道次 将外径 '$,* XX% 内径 '$#, XX的管坯逐步径 向锻造成大径段外径 '$)& XX内径 '$#, XX% 小径段外径 '$'" XX内径 '$#, XX的锻件& 图 ' 为液压缸筒坯三道次径向锻造工艺的成形过程 具体流 程 图& 先 将 外 径 '$,* XX% 内 径 '$#, XX% 长度为 '"* XX的管坯沿整个轴向径向锻 造成外径 '$)& XX% 内径 '$#, XX的管件" 然
"!液压缸缸筒径向锻造工艺原理
液压缸缸筒径向锻造工艺的三维模型图如 图 " 所示" 芯轴置于管坯中" 四个锤头沿着管 坯的圆周方向均布阵列" 机械手夹持管坯的一 端& 在径向锻造过程中" 四锤头沿管坯的径向 往复运动" 进行同步锻打" 使管坯产生塑性变 形# 在锻打间隙" 管坯在机械手的夹持下" 绕 自身轴线旋转# 机械手旋转的同时" 作轴向进 给运 动# 通 过 上 述 运 动 的 配 合" 实 现 管 坯 变 长% 壁厚变薄&

哈尔滨工业大学2004年第一次攻读博士学位研究生拟录取名单

哈尔滨工业大学2004年第一次攻读博士学位研究生拟录取名单
春季
仪器科学与技术
戴景民
张虎
秋季
仪器科学与技术
戴景民
辛春锁
秋季
仪器科学与技术
戴景民
宋扬
春季
仪器科学与技术
丁雪梅
陈洪芳
春季
仪器科学与技术
丁雪梅
陈浩
春季
仪器科学与技术
丁振良
佟庆彬
秋季
仪器科学与技术
丁振良
陈杰春
春季
仪器科学与技术
丁振良
李?
春季
仪器科学与技术
付平
刘雷
秋季
仪器科学与技术
彭喜元
付宁
秋季
仪器科学与技术
春季
应用化学
安茂忠
乐士儒
秋季
应用化学
傅宏刚
谢颖
秋季
应用化学
高云智
张一
秋季
化学工艺
胡信国
高昆
春季
化学工艺
胡信国
伊廷锋
秋季
化学工艺
胡信国
杨春巍
春季
应用化学
黄玉东
贾近
春季
应用化学
黄玉东
王哲
春季
应用化学
黄玉东
刘玉荣
春季
应用化学
黄玉东
张学忠
春季
应用化学
黄玉东
卢晓东
秋季
应用化学
黄玉东
徐志伟
秋季
化学工艺
姜兆华
崔瑞海
韩纪庆
林常志
秋季
计算机科学与技术
韩纪庆
张茹
秋季
计算机科学与技术
郝忠孝
马海涛
春季
计算机科学与技术
洪炳熔
杨晶东

第六章 混凝土简支梁桥的计算.

第六章 混凝土简支梁桥的计算.
计算时荷重 (桥面铺装、人行道、栏杆、灯柱) 均匀分配给各主梁。
M xg 2 lxgx2 xg 2 x(lx)
Qxg2lgxg 2(l2x)
5/30/2021
[例2-3-2]一座五梁式装配式钢筋混凝土简支梁桥的主梁和横 隔梁截面如下图所示,计算跨径l=19.5m,结构重要性系数 1.0。求边主梁的结构自重产生内力。(已知每侧的栏杆即人 行道构件重量的作用力为5kN/m)。
14
5/30/2021
(尺寸单位:cm)
第二节 主梁内力计算
一、结构自重效应计算 二、汽车、人群荷载内力计算
1.横向荷载分布的定义 2.横向荷载分布的计算 3.横向荷载分布系数m沿桥跨的变化 4.汽车、人群作用效应计算
三、主梁内力组合
5/30/2021
主梁内力计算中截面的选取:
对于跨径在10m内的简支梁,需计算: 跨中截面的最大弯矩 支点截面及跨中截面的剪力 其它各截面的剪力假设按线性规律变化,弯矩假设按二 次抛物线规律变化
d
d:为最外两荷载中心距离
5/30/2021
b.荷载位于支承边处
aa1ta22Ht3 l
c.荷载靠近支承边处
ax a'2xa
x—荷载离支承边缘的距离。 说明:荷载从支点处向跨中移 动时,相应的有效分布宽度可 近似地按45°线过渡。 按上述公式算得的所有分布宽 度,均不得大于板的全宽度。
5/30/2021
Ms MspMsg
5/30/2021
3)自由悬臂板
车轮荷载靠板的边缘布置
b1 b2 H
可变作用:
M sp(1)1 2pl0 2(1)4a P b1l0 2, (b1l0时 )
M sp(1)pb1(l0b 2 1)(1)2 P a(l0b 2 1), (b1l0时 )

板特征与杆特征在圆管自由振动计算中的统一

板特征与杆特征在圆管自由振动计算中的统一
图 2 圆 管 的截 面 振 动 模 态
同长度 的梁单元 模型计 算结 果 中 , 出其 基频 , 取 可得
出圆管相 应 厚 度 下 梁 单 元 基频 随 长 径 比变 化 的 情
况 。 图 3显 示 的 是 板 厚 分 别 为 1 0 mm 和 4 mm 时 0
为 了 研 究 不 同 的 长 径 比 和 不 同 的 径 厚 比对 自振 模 态 的影 响 , 同 一 直 径 不 同 壁 厚 和 不 同 长 度 的 圆 对
轴线垂 直 , 而是 与 梁 的轴 线 呈一 定 的 角度 。该 理 论 中, 剪应变 沿梁厚 度方 向为 常数分 布 , 违反 了上下 表
面剪应 力 为 0的 自由边 界条 件 。为 消 除 这一 现 象 , 引入 与几何 形状 和 弹性 特 性 有关 的剪 切修 正 因子 。
此外 , i h n o梁 理 论 不 能 反 映 梁 横 截 面 的 翘 曲 Tmse k
管 进行分 析 。这里 认 为 A S S合 理 划分 网格 后 计 NY
长径 比与 基频 的关 系 , 厚 分 别 为 5 板 0 mm、5 m 2 m、
算 出的结 果在梁 和壳 分析 中是可信 的 。
1 2
1 m、0mm、 m 时得 出的长径 比与 频率之 间 的 6m 1 6m
和横 截面形 状 的变化 … 。
l 思路 与基 本 理 论
圆管 在 不加 任 何筋 板 和横 隔 板 的情 况下 , 振 其 动 呈现 出两 种基 本 的模 态 : 是 在管 的长度 方 向上 一 的振 动模态 , 即梁单 元振 动模 态 ; 二是 在 圆截 面上 的 振动模 态 , 即板壳振 动模 态 ( 图 1 图 2 。 见 、 )

三维机织圆管复合材料拉伸载荷下损伤过程的有限元模拟

三维机织圆管复合材料拉伸载荷下损伤过程的有限元模拟

mo es d l
三维 机织 复合 材料 作 为纺织 复 合材 料 的e td o h tan—te s c r e T u rc lmo e s v rfe r p ri s o s r a l r s n e n t e sr i sr s u v . he n me i a d li eiid. e
t nfr odig R M)e h i . odcr l i curdb tentecl lt nrsl n x r s ll ( T t n sG o or a o i aq i e e h a uai eut ade— a em n c c e tn s e w c o s
Ke r s:t e — i e so a v n o p st s u u a tucu e; m a e v l t n; e h nc l y wo d hr e d m n i n l wo e c m o ie ;t b lr sr t r da g e o ui o m c a ia
V0 . 5 1 3 NO .1 Fe 2 b. 01l
三 维 机 织 圆 管 复 合 材 料 拉 伸 载 荷 下 损 伤
过 程 的 有 限 元 模 拟
王 新 峰 刘 佳 , 光 明 。 一, 周
( 南京航空航 天大 学 1结 构强度研究所 ;. . 2 飞行器结构力学与控制教 育部重点实验室 , 江苏 南京 2 0 1 ) 10 6
( .ntueo t c rs& Srn t :. E K yL bo t c r Meh n sa dC nrl 1 Istt f r t e i Su u t gh 2 MO e a f r t e c a i n o t e Su u c o

工程力学建模

工程力学建模

与 应 用 脱 节 , 利 于 学 生 的技 能培 养 。在 教 学 实 际 中 , 以根 据 “ 靠 性 ” “ 济 性 ” 则 , 工程 力 学 建 模 分 解 为 “ 不 可 可 、经 原 将 主
体 的处 理 ” “ 的 处 理” “ 结 的 处 理 ” 个 基 本 过 程 。 、力 、联 三
工 程 力学 建 模
杨Байду номын сангаас冠 声
( 津现代 职 业技 术 学院 , 津 市 3 0 2 ) 天 天 0 2 2
摘 要 : 传 统 的 工 程 力 学教 学 , 接 就 力 学模 型 讲 授 理 论 的 应 用 , 乏 从 工程 1题 中抽 象 出力 学模 型 的 环 节 , 直 缺 " - 3
力 学 建 模 是 力 学 工 程 应 用 首 先 面 临 的 问 题 , 进 行 力 学 分 析 计 算 的基 础 。大 量 的 问 题 要 在 这 一 过 程 中完 成 , 是 包
括确定分析方案 、 分解 分 析项 目、 化 分 析 对 象 、 定 载荷 情 况 等 等 , 是 力 学 研 究 、 学 实 践 的重 要 环 节 。对 于 工 程 简 确 它 力
时也是对应用规律的描述 。建模过程可以使学生加 深对理 论本质 的理解 , 把握 理论 的核心要 素, 分清操作 中各 因素 的主次差别 , 领会在应用 中的关键 、 难点及一般解决方案中存在的问题 , 方便学生进一步研究与 分析深度 的加强 。
工 程 力 学 教 学 , 该 弱 化 理 论 推 演 , 强 应 用 环 节 , 模 过 程 是 达 成 这 一 目标 的重 要 手 段 。学 生 的认 识 过 程 , 喜 应 加 建 有
来 的一 个 必 须 的 环 节 。在 大 学课 程 教 学 中 , 该 加 强 这 一 能 力 的培 养 , 添 这 一 部 分 的 内容 。 应 增

211151054_横向蒙皮拉伸机运动学分析

211151054_横向蒙皮拉伸机运动学分析
)-$)
重 型 机 械!! !!!!!!!!!!!!!!!!%$%( 7+)%
横向蒙皮拉伸机运动学分析
杨!伟#%
#%
程军周

#%
张康武

#%
闫玉平

张!镭#%
#%
芦跃峰
!#H金属成形技术与重型装备全国重点实验室" 陕西 西安 "#$$#W& %H中国重型机械研究院股份公司" 陕西 西安 "#$$#W#
的剩余自由度" 即机构出现奇异位形%
通过拉伸机构位于不同铰链点布置时" 锁住
X 个驱动副后对拉伸机构约束螺旋系的线性相关
性进行分析" 得到了拉伸机构的 % 种奇异位形和
% 种非奇异位形" 并根据 =E;@@K;55 线几何原理
对这 % 种奇异位形和 % 种非奇异位形一一进行
说明%
由文献,#W- 可知" 拉伸机构中 bGb约束分
互平行的力线矢的最大线性无关数为 (" 即 #N$
%N$
(N$
N
X这
X
个约束力线矢线性相关"
所以锁
住 X 个驱动副后" 夹钳的自由度不为 $" 此时夹
钳剩余 # 个绕 E轴的转动自由度" 该机构处于奇
异位形%
如图 ( ! P# " 约束 力线 矢 #N$
%N$
(N$
N
X汇
交于空间一点 #Q" 根据 =E;@@K;55 线几何原理可
股份公司工程师" 研究方向( 蒙皮拉伸成形技术研 究及应用%
国 +G7公司 'c3型纵向蒙皮拉伸机$ 'c&型横

拉伸载荷下U形缺口薄板试件应力集中系数修正公式

拉伸载荷下U形缺口薄板试件应力集中系数修正公式

霹 = /s e
式 中 : 缺 口处 的真 实应 变 ; 为缺 口处 的名 义应 变。 为 e
B r t —N a 法 的表 达式 为… : ar a el a
() 2
K(823 )9 8 ) 0 =. +2 [ 0 6 00. 0 7 4 . ( 0
式 中 : 为缺 口深度 ; 为 缺 口半 径 ; t r D为试 件宽度 。
E —mal z a g h 1 6 @ 1 3 c r i:h n z p 9 2 6 .o n
第 1 期
张忠平 等 : 拉伸载荷下 U形 缺 口薄板试件 应力 集中系数修正公式
7 5
为试 件厚 度 , 为试 件 长度 。
表 1 各试件尺寸参数及应力集 中系数
T b 1 T e d me so s frv r u oc e p cme s a . h i n in a i sn th d s e i n o o
1 应力集 中系数 的定义及解析估算方法
应力集 中系数 的定义是 ¨ :
=' oS / () 1
式 中 : 为应 力集 中系数 ; 为缺 口处 的最 大真实应 力 ; K o r . s为缺 口处 的名义 应力 。
通常情 况下 , 应力 集 中系数 的估算方 法分 为数值计 算法 与解析 法 。数 值计算 法 主要采用 有限元 法 , 解析 法 主要利用 N ue 法 和 B r t—N a 法 , 中 N ue 法的表 达式 为 J ebr aa r a el 其 ebr :
第 l 卷第 1 1 期
21 00年 2月





ห้องสมุดไป่ตู้

基于ABAQUS自升式平台的桩腿弯曲数值模拟分析

基于ABAQUS自升式平台的桩腿弯曲数值模拟分析

斯托的垂直距离按规范确定 ;Cs 为构件的形状系数,其值可根据规范给出值确定 ;S
P= 1 ρ gH ch[k(z+ d)] cos(kx− ω t) + 3 H 2 ( ρ gh) × 1
2
ch(kd)
4 L sh(2 kd)
为平台处于平浮或倾斜状态时,受风构件在受风方向
单元,桩土接触面设置为摩尔库伦罚函数摩擦形式。 边界条件为 :土体侧面约束径向上的位移,底部约束 竖直方向的位移。分析步骤为 :第一步为地应力平衡 分析步,由于桩体和土体的重度不同,在设置时先假 设桩和土体的重度一致 ;然后通过第二步静力分析步 中给桩体加上桩土间重度的实际差值 ;最后在给桩设 置实际所受的载荷力 [8],见图 1。
Yang Yifei
(Department of offshore oil and gas engineering, School of petroleum engineering, Xi’an University of petroleum, Xi’an 710065, China)
Abstract: In order to analyze the reasons for the bending of the legs of a jack-up platform in the Bohai Sea when pulling out the piles, the ABAQUS finite element software was used to analyze the platform legs under the effects of wind and wave loads, focusing on the study of the legs under different wave parameters. The offset distance and stress of the leg center are small. According to the calculation results, the reasons for the offset distance and the force of the leg center under different conditions are analyzed, and the actual bending of the leg is compared. The results show that the possible causes of the leg bending when pulling out the pile are wind and wave loads and the inertial force of the platform.

圆板在冲击荷载作用下塑性动力响应的加权余量法解

圆板在冲击荷载作用下塑性动力响应的加权余量法解

圆板在冲击荷载作用下塑性动力响应的加权余量法解
刘福林
【期刊名称】《振动与冲击》
【年(卷),期】1991(010)001
【总页数】6页(P55-60)
【作者】刘福林
【作者单位】无
【正文语种】中文
【中图分类】O347.3
【相关文献】
1.简支圆板在冲击荷载作用下的塑性动力响应分析 [J], 王延斌;魏雪英;俞茂宏
2.冲击荷载作用下简支圆板的塑性动力响应统一解 [J], 王延斌;俞茂宏;肖耘;李林生
3.简支圆板在爆炸冲击荷载作用下的动力响应 [J], 刘永行;刘彦东;梁崇勇;延鑫年
4.圆板在冲击载荷下的大挠度塑性动力响应 [J], 于政文;张春雨
5.冲击载荷作用下阻尼介质中刚塑性圆板的动力响应 [J], 李庆明
因版权原因,仅展示原文概要,查看原文内容请购买。

研究生学位论文开题报告的查新技术

研究生学位论文开题报告的查新技术

模拟实验;(3)人机交互模拟实验。

目前一些高校已经开始用MATLAB软件为工程力学进行模拟计算和实验,可以借鉴这些先进成果进行转化应用。

3.增加力学建模考核内容。

目前我校仍然采用课后习题作业和闭卷考试作为考核学生的基本手段,这种方式考查学生记忆能力的较多,不利于发挥学生运用知识和动手实践的能力。

因此,有必要对现行考试方法和考试内容加以改革。

平时的考核可以采用课后习题作业、实验操作和读书报告的形式[6],课后习题是适量的基本概念的考核和理论计算;实验内容应贴近于工程实际,让学生综合运用理论知识,建立其力学模型,研究其理论解决方案;读书报告要求学生就一个工程实践中的力学问题,根据力学基本概念和定理对问题进行抽象简化并建立其力学模型,然后利用所学的理论知识确定计算方法,最后进行分析计算并给出解答。

期末考试内容不仅包括考核学生掌握力学基本概念的能力,还包括考核学生应用力学理论知识分析和解决问题的能力。

例如,给定一个简单的典型力学问题,让学生简化力学模型,给出受力分析,为各构件选取材料及截面形状和尺寸,并对该结构存在的问题谈自己的看法。

这种考核方式注重培养学生解决力学问题的能力和对所学工程力学知识进行归纳、总结的能力。

参考文献:[1]杨冠声.工程力学建模[J].天津职业院校联合学报,2007,(11):12-15.[2]黎杰松.浅析大学工程师教育存在的问题及对策[J].华北水利水电学院学报(社科版),2012,(4):155-157.[3]周丽珍,张涛.勘查技术与工程专业的“工程力学”课程教学方法探索[J].中国地质教育,2010,(增);[4]张涛,周丽珍.加强高校实验室建设培养学生创新思维[J].中国地质教育,2007,(增).[5]王彦生,侯中华.创建工程力学实验教学示范中心培养创新型人才[J].实验室研究与探索,2010,(8):75-77.[6]刘永寿,支希哲.工科理论力学考试改革的理论与实践[J].力学与实践,2004,(26):68-69.基金项目:中国地质大学(武汉)教学研究项目(2012B07)作者简介:周丽珍(1979-),女,在读博士研究生,讲师,主要从事工程力学方面的教学和科研工作。

工程力学杨庆生答案

工程力学杨庆生答案

工程力学杨庆生答案【篇一:北工大研究生复试科目】学院(研究院、所)名称:机械工程与应用电子技术学院注:初试科目选择813电工学的考试,复试科目只能选择机械设计基础或工程力学i。

学院(研究院、所)名称:电子信息与控制工程学院学院(研究院、所)名称:建筑工程学院学院(研究院、所)名称:建筑工程学院【篇二:ansys经典命令流资料裂缝】标题: re: wjc501请帮帮忙,有没有裂缝方面的ansys例题发信站: bbs 水木清华站 (mon jul 22 18:21:57 2002)供参考,不当之处请指正。

祝好运!(1) 如果你模拟结构体中裂缝扩展过程的模拟,在ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展,我曾用ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。

我模拟的过程相对来说比较简单,模拟过程中我们首先要知道裂缝的可能扩展方向,这样在裂缝可能扩展的带内进行网格加密处理,加密到什么程度依据计算的问题来确定。

(2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子,建模时只需在裂尖通过命令kscon生成奇异单元即可。

ansys模块中存在的断裂力学模块可以计算i、ii、iii型应力强度因子(线弹性断裂力学)和j积分(弹塑性断裂力学),在ansys中verification里面有一个计算i型应力强度因子的例子vm143,参见该例子就可以了。

(3) 如果通过断裂力学模拟裂缝的扩展过程,需要采用动态网格划分,这方面我没有做,通过ansys的宏命令流应该可以实现。

技术参考可参阅文献:杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报,1997,14(4).(4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律,我是通过位移插值得到不同时间点处的应力强度因子。

如果想这样做,可参阅理论参考中关于应力强度因子计算说明。

钢管混凝土系杆拱桥拱圈落架仿真分析

钢管混凝土系杆拱桥拱圈落架仿真分析

钢管混凝土系杆拱桥拱圈落架仿真分析马祥春【摘要】Combined with the double ballast track 96 m concrete-filled steel pipe arch bridge in Tianjin Nangang Railway engineering,a model is established with the Midas Civil.Through comparative analysis of fourkinds of arch ring landing structure schemes,scheme with minimum control cross section deformation and stress values in the process of arch ring landing is found out.It is concluded that scheme of successive symmetric landing from a quarter of the arch rib to the arch foot and vault is more favorable to arch ring deformation and stress distribution,and good results have been achieved in the construction process.And considering factors such as construction plant,environment,construction technology and condi-tion,the author puts forward three arch ring asymmetry landing structure construction schemes to provide theoretical refer-ence for monitoring in the process of construction and construction safety control.%结合天津南港铁路工程有砟轨道双线96 m 钢管混凝土系杆拱桥,利用 Midas Civil 建立模型,通过对4种拱圈对称落架方案的对比分析,得出在拱圈落架过程中控制截面的变形值和应力值最小的方案,提出了从1/4拱肋附近逐次对称向拱脚与拱顶均衡卸落的落架方案对拱圈变形和受力更为有利,并在施工过程中取得了较好的效果。

《工程力学(三)》军队优质课程建设探索与实践

《工程力学(三)》军队优质课程建设探索与实践

《工程力学(三)》军队优质课程建设探索与实践
方祥位;石少卿;徐守彬;张忠发;郭剑峰
【期刊名称】《训练与科技》
【年(卷),期】2006(027)003
【摘要】课程建设是院校教学的基础性工作,是提高教学质量的重要环节,院校的课程建设状态是教学工作水平的重要标志。

2003年,教育部做出决定,在全国高等学校(包括高职高专院校)中启动高等学校教学质量与教学改革工程精品课程建设工作。

近年来,精品课程建设受到各高等院校的广泛关注。

【总页数】6页(P44-49)
【作者】方祥位;石少卿;徐守彬;张忠发;郭剑峰
【作者单位】后勤工程学院军事建筑工程系
【正文语种】中文
【中图分类】TB12-4
【相关文献】
1.高分子物理与化学优质课程建设的探索与实践
2.《粉体工程》优质课程建设探索与实践
3.船舶流体力学优质课程建设的探索与实践
4.基于"三全育人"的《工程力学》课程教学探索与实践
5.基于网络平台研究生优质课程建设的探索与实践
因版权原因,仅展示原文概要,查看原文内容请购买。

金属材料动态损伤的微观数值模拟

金属材料动态损伤的微观数值模拟

第28卷第1期2007年1月江苏大学学报(自然科学版)Journal of J iangsu University (Natural Science Editi on )Vol .28No .1Jan .2007金属材料动态损伤的微观数值模拟叶建军,杨 健,郑津洋,贺世正(浙江大学化工机械研究所,浙江杭州310027)摘要:针对金属材料的动态损伤现象,采用分子动力学方法对存在微观表面缺陷的材料在瞬态强冲击作用下的损伤情况进行数值模拟.模拟对象为铜材料,表面缺陷为锲形、圆形和方形沟槽,采用目前较为合理的嵌入原子势(EAM ),并利用分子动力学程序(LAMMPS )进行动态数值模拟,有效地给出了存在不同微观表面缺陷的金属铜材料在瞬间强冲击作用下的损伤效应.结果表明,存在锲形缺陷的金属表面更容易发生微喷射现象,且微喷射速度随着锲形角增大而减小,随着冲击速率增大而增大.关键词:微喷射;数值模拟;分子动力学方法;E AM 势;动态损伤中图分类号:TG111;O347.3 文献标识码:A 文章编号:1671-7775(2007)01-0041-05M i cro 2scale nu mer i cal si m ul ati on on dynam i c damagefor met al materi a lsYE J ian 2jun,YAN G J ian,ZHEN G J in 2yang,HE S hi 2zheng(I nstitute of Che m ical Engineering Pr ocess and Machinery,Zhejiang University,Hangzhou,Zhejiang 310027,China )Abstract:A i m ed at vari ous surface defects of the material,the m icr o 2scale numerical si m ulati on f or dy 2na m ic da mage under instantaneous i m pact is conducted .The e mbedded at om s method (E AM )potential and the molecular dyna m ics (MD )code LAMMPS are adop ted .The shapes of surface defects are de 2scribed as rectangle,r ound and gr oove,res pectively .The da mage effect of the copper material with vari 2ous surface m icr o 2scale defects under instantaneous i m pact was si m ulated effectively .The results show that the surface with gr oove defect is easier t o induce the ejecti on of m icr o particles than the circular or square defects .The ejecti on vel ocity decreases with the increase of the gr oove angle,and increases with the increase of the i m pact s peed .Key words:m icr o 2particle ejecti on;numerical si m ulati on;molecular dyna m ics;EAM potential;dyna m ic da mage收稿日期:2006-09-18基金项目:国家自然科学基金资助项目(50475100);国家自然科学基金中澳科技合作特别基金资助项目(50610055)作者简介:叶建军(1981-),男,湖北鄂州人,博士研究生(salas_1981@t ),主要从事计算流体动力学及微尺度数值模拟研究.杨 健(1968-),男,浙江杭州人,副教授(通讯作者,zdhjkz@zju .edu .cn ),主要从事计算流体动力学及ME MS 研究. 当金属材料自由表面受到强冲击载荷作用时,有部分物质微粒以高于自由面的运动速度向外喷射,导致强烈的材料结构损伤,这一现象称为微喷射动态损伤.自从1953年W alsh 在试验中首次观察到此现象以来,金属材料在强冲击作用下的动态损伤效应及机理越来越为人们关注.不仅如此,相关物理研究表明[1],在瞬态强冲击载荷作用下,金属材料表面不同微观缺陷具有各自特殊的破坏效应,并且其动态过程十分复杂,迄今仍缺乏相应的材料微观动态损伤理论模型,而且想要进行微观特性试验研究,往往成本和难度都很高.因此,发展一种微观数值模拟方法[2,3],构筑一套微观模型与仿真算法,实现对金属材料在强冲击载荷下的动态损伤过程数值模拟试验,具有重要的学术价值和实际意义.自20世纪90年代以来,随着计算机软、硬件技术的发展与突破,基于微观的分子模拟技术更为人们从原子及分子尺度探讨材料的微观特性奠定了基础[4].最近,随着并行计算和网格技术的成熟,其模拟尺度和仿真能力正在时间和空间尺度上不断拓展[5],有可能突破经典分子动力学模拟的计算瓶颈,从而使得微观数值模拟真正步入应用时代.总结现有文献,A say等人在1969年利用统计力学原理对微喷射量和金属表面粗糙度的关系进行过研究,其又在1997年通过试验,研究了强冲击作用下的金属表面特性[6];在国内,陈军等人也对金属表面的微喷射现象进行了模拟研究[7],但其研究只限于锲形槽缺陷,没有进行过多种缺陷形式的动态损伤研究及对比分析.而且,目前关于不同缺陷形式下金属表面受到强冲击损伤的研究,在国际上也还未见报道.因此,文中采用分子动力学方法,在微观层面上对强冲击载荷作用下金属材料破坏形式进行动态特性的并行数值模拟,主要是通过对表面存在不同微观缺陷的金属材料,在瞬间强冲击载荷下的破坏现象进行模拟研究,揭示传统的研究手段所难以表达的金属材料的微观动态损伤机理,尝试建立一种有效的金属材料动态模拟微观数值方法.1 MD方法简介MD方法的基本原理是从微观尺度建立一个粒子系统来模拟所研究的对象,系统中各粒子之间的相互作用根据量子及分子动力学理论来确定.首先,通过求解粒子的动力学方程组可得到各粒子在相空间的运动规律和轨迹,然后,按统计物理原理得出该系统相应的宏观物理特性.由于分子动力学模拟具有沟通宏观特性与微观结构的作用,特别是可以捕捉到许多在理论分析和试验观察上都难以获得的现象,并能给出一定的微观解释,模拟一些极端条件下的微观现象,而且能应用广泛,因此被认为是理论和试验观察相联系的第三种科学手段[8].MD方法的模拟计算主要基于以下基本假设:(1)所有粒子之间的的运动都遵循经典牛顿力学定律;(2)粒子之间的相互作用能满足叠加原理.其具体的理论表述为,对于粒子数为N的物理体系,根据牛顿第二定律有F i(t)=m i a i(t)=m i52r i(t)5t2(1)其中,i为粒子序号,最大个数为N个,Fi(t)为i粒子在t时刻所受的力,mi为i粒子的质量,ai(t)为i粒子在t时刻的加速度,r i(t)为i粒子在t时刻的位移.根据势场理论中力与势能的关系,i粒子所受的力可由分子力学势能函数的负梯度求得F i(t)=-55r i U(r1,r2,…,r n)(2)其中,U为体系的总势能.将式(2)代入式(1),可以得到-55r i U(r1,r2,…,r n)=m i52r i(t)5t2(3)从以上推导过程可以看出,只要给定原子的初始位置和初始速度,并对很小的时间间隔作数值积分,便可得到该原子下一时刻的速度和原子位移即原子运动轨迹,再按统计物理原理可得出该系统相应的宏观物理特性.其中,要获得比较准确的模拟结果,最重要的是要给定粒子与粒子间作用势的数学量化表达式.2 嵌入原子势函数MD模拟的关键是作用势函数的选取,文中采用嵌入原子势(Embedded A t om s Method,简称为EAM)来表达其金属铜原子间的相互作用.E AM势已经被证明比较适合于金属表面研究,且能很好地描述金属原子之间的相互作用[9].文献查证[10,11], Da w和Baskes采用St ott和Zare mba的准原子思想,在准原子近似和有效介质理论的基础上根据密度泛函理论认为:某原子的原子核除了受到周围其他原子核的排斥作用外,还受到该原子的核外电子及其周围其他原子产生的背景电子的静电作用.于是,由其提出的E AM模型中,将N个原子组成的金属系统总能量表示为E tot=126Ni=16Nj=1≠iφij(rij)+6N i=1F i(ρi)(4)其中,φij是原子芯与原子芯之间的二体排斥势,Fi 是将原子i嵌入电子密度为ρi的位置处所必须的能量,ρi为系统中所有其他原子在原子i处产生的局域背景电子密度,因此ρi可以描述金属原子势能的多体性质.F i(ρi)仅依赖于嵌入原子的种类,与作用于ρi24 江苏大学学报(自然科学版) 第28卷的原子种类无关.因此,无论对于合金或纯金属都可采用相同的函数进行嵌入能的计算.显然,对于嵌入原子法,为了确定各物质的对势φij 、电子密度ρi以及嵌入能Fi,必须相应地确定一些参数.对此Daw 和B askes提供了在平衡状态下块体纯金属(单质金属)的弹性系数、相变潜热和空位形成能的试验值,确定了相应的对势参数和Fi[12].早期势函数的拟合过程比较复杂,后来Fo iles等人建立了一种比较简捷的拟合方法,即利用Ro se普适函数来表达:E(r)=-E c1+b rr e -1exp-brr e-1(5)其中,b=ΩeBe c1/2,Ωe为平衡原子体积,E c为结合能,B是体积弹性模量.令式(4)和式(5)相等,则当电子密度函数ρi 已知时,如果知道对势的具体形式,就可确定嵌入势;反之如果知道嵌入势,也可以唯一地确定对势.在对势部分,Da w和Foiles采用库仑力形式,Johns on 采用Born2Mayer势,Vet or等人采用Morse势.早期的研究中,电子密度函数ρi一般用Hartree2Fock理论计算得到的自由原子电子密度表示,最近的研究中则常常采用负指数形式.总之,E AM势函数能很好地描述金属原子之间的相互作用,是最新的金属体系模拟中比较贴切的作用势函数表达.3 MD模拟计算文中采用大规模并行分子动力学计算程序(LAMMPS),编程语言为C++.在最大CP U数为64的SGIOnyx3900超级计算机上进行并行计算,模拟粒子数量可达百万级[13-16].模拟研究对象以铜材料为代表,其面心晶格常数为0.344n m,势函数为金属铜原子的EAM势.在所构建的物理模型中,冲击主要沿x方向,y方向和z方向的表现性质基本相同,类似于二维问题.但实际计算中仍采用三维模拟,只是给出的结果分析及图示中主要表达的是x方向和y方向的变化.模拟的瞬间强冲击载荷由同种材料金属原子构成的飞片与靶相互碰撞造成.其中,被撞击的靶的长度为130 n m,高度为200n m,飞片的长度为60n m,高度为200nm,模拟的粒子数约为0.85×106个.垂直于冲击波传播方向的边界条件设置为周期性边界条件,计算时间步长设置为0.001p s,计算步数为10000步.模拟的金属表面缺陷为锲形、圆形和方形等沟槽,分别得到其不同微观表面缺陷在瞬间强冲击作用下的各自损伤状况.4 模拟结果与分析4.1 不同微观表面缺陷的动态损伤图1为方形、圆形和锲形等沟槽形式表面缺陷存在下的瞬间强冲击作用模拟结果,其中方形缺陷深50n m,宽20nm;圆形缺陷半径为30nm;锲形缺陷张开的锲角为60°,锲口宽度60nm.图1c所示的模拟计算结果表明锲形表面缺陷在瞬间强冲击作用下,锲形尖角处存在明显的粒子被激发后向外喷射的现象(以下简称微喷射).而这样的微喷射现象会直接影响金属材料的微观结构,从而局部发生急剧的崩塌损伤.图1a和图1b所示中,在同样的强冲击下,方形和圆形表面缺陷则没有这样的微喷射现象,甚至在图1a中还表现出会自我修复,在图1b中则表现出消融的现象.可见,不同微观表面缺陷在瞬间强冲击破坏下其损伤现象极为不同,而且材料表面存在锲形微观缺陷最易发生微喷射现象,会使得材料的抗冲击性能大大降低.因此,有必要重点研究锲形微观表面缺陷在强冲击加载过程中的具体破坏过程.(a)方形(b)圆形(c)锲形图1 微观表面缺陷受瞬间强冲击作用的模拟结果Fig.1 Si m ulati on result of surface defect underinstantaneous i m pact4.2 不同角度锲形缺陷的微喷射特性以上模拟结果表明,锲形微观表面缺陷受强冲击载荷作用易发生微喷射现象,显然其破坏效应最34第1期 叶建军等:金属材料动态损伤的微观数值模拟为严重.而且宏观研究表明,缺陷最易遭受破坏处为其尖端部分,如裂纹尖端,因为材料的缺陷尖端在冲击作用过程中受到的应力最大,破坏也往往从缺陷尖端开始.因此,需要重点研究锲形缺陷的尖端处粒子和锲形缺陷自由表面处粒子,各自在冲击作用过程中的速度和受力等的变化情况.如图2所示,不同角度下的锲形表面缺陷,其缺陷尖端的粒子和缺陷斜自由面粒子在冲击作用下有着不同的喷射速度.由图可见,随着锲形的张开角度逐渐增大,锲角顶端的粒子从锲形沟槽中的飞出速度增大,同时锲槽自由面的粒子沿冲击波方向的速度也增大,两者的速度差导致出现微喷射现象;而且由图2可知,两者的速度差值会随着锲形张开角的增大而减小,这可以解释材料表面微缺陷的锲角张开越大,微喷射的现象越不明显,说明更能承受较大的冲击;另一方面,材料表面锲形缺陷角度越大,说明材料本身的表面粗糙度越小,微观缺陷程度较小,这与实际宏观性质和现象分析结果相一致.图2 锲形缺陷的粒子速度随锲形张开角度的变化Fig .2 Relati onshi p bet w een gr oove angle and particles vel ocity4.3 不同冲击速率下的微喷射损伤效应为了给出不同冲击速率下的锲形缺陷存在的微喷射程度,取60°锲角的微观缺陷作为模拟对象,设飞片对铜材料的初始冲击速度变化为范围为5~20An /Ps,即冲击速度为500~2000m /s,分别计算在各个冲击过程中,锲形缺陷尖端及自由面的粒子速度变化.图3为锲形缺陷的尖端粒子喷射速度及自由面粒子速度随冲击速度的变化曲线.由图可见,当冲击速度增大时,材料内部应力明显增大,缺陷尖端粒子的微喷射速度增大,同时缺陷自由表面处粒子飞出速度也增大,但是显然尖端处粒子的喷射速度增加得更快,微喷射现象更明显.由此可以说明当冲击增大时,有锲形缺陷的金属材料更容易发生微喷射现象,也就是说有此类微观缺陷的材料更加容易发生结构破坏.图3 锲形缺陷的粒子速度随冲击速度的变化Fig .3 Relati onshi p bet w een i m pact vel ocity andparticles vel ocityy5 结 论针对金属材料的动态损伤现象,采用分子动力学方法,结合E AM 金属势,对存在微观表面缺陷的铜材料,在瞬间强冲击作用下的损伤情况进行数值模拟,研究表明:(1)在相同强度的冲击作用下,材料的锲形微观表面缺陷比方形和圆形更容易发生微喷射现象,说明若表面存在锲形缺陷,材料的抗冲击能力较差.因此,在涉及强冲击的场合,例如爆燃爆轰容器的材料应该尽量避免其表面出现锲形微观缺陷;(2)在相同强度的冲击作用下,锲角的变化对锲形缺陷本身的抗冲击性能有不同的效应,计算结果显示锲角越小越容易发生微喷射现象.因此,除了改进材料性能,提高材料表面加工精度和表观质量,也能直接提高其抗强冲击能力;(3)在锲角不变的情况下,随着冲击速度的增大,尖端粒子的速度比缺陷自由面速度增大更快,因此喷射现象也更为明显,这与微喷射试验结论相一致.以上数值模拟结果能定性地描述金属铜材料不同微观表面缺陷,在强冲击作用下的微观动态损伤破坏过程,这表明文中提出的针对金属材料动态损伤微观数值模拟的方法是有效的.在后续的研究中,还将重点通过模拟来研究金属材料在高温高压下的多种性能参数预测,如杨氏模量、泊松比和应力波等,以期获得进一步的定量模拟结果.参考文献(References)[1] Palaci os F J,I niguez M P .Molecular dyna m ics study ofcluster i m pact on the (001)and (110)surfaces of fee metals [J ].Co m putational M aterials Science ,2000,17(7):515-519.[2] Elber,Ron .Novel methods f or molecular dyna m ics44 江苏大学学报(自然科学版) 第28卷si m ulati ons[J].Current O pinion in S tructural B iology,1996(4):232-235.[3] Zhou Guohui,Gao Kehan.Molecular dyna m ics si m ula2ti on of m icr ocrack healing in alu m iniu m[J].Progress inN ature Science,2001,11(3):216-220.[4] 程晓农,戴起勋.材料计算设计及研究进展[J].江苏大学学报:自然科学版,2003,24(1):15-18.CHE NG Xiao2nong,DA IQ i2xun.Devel opment of mate2rials calculati on design[J].Journal of J iangsu U niversi2ty:N atural Science Edition,2003,24(1):15-18.(inChinese)[5] Cao Xiaolin,Mo Zeyao.Parallel computati on f or mole2cular dyna m ics si m ulati on based on cell2bl ock datastructures[J].Chinese Journal of Co m putationalPhysics,2004,21(5):377-385.[6] A sa y J R.The use of shock2structure methods f orevaluating high2p ressure material p r operties[J].Inter2national Journal of I m pact Engineering,1997,20(1):27-61.[7] 陈 军,经福谦,张景琳,等.冲击作用下金属表面微喷射的分子动力学模拟[J].物理学报,2002,51(10):2386-2392.CHE N Jun,J I N G Fu2qian,ZHANG J ing2lin,et al.Mo2lecular dyna m ics si m ulati on of m icr o particle ejecti onfr om a shock2i m pacted metal surface[J].A cta PhysicaS inica,2002,51(10):2386-2392.(in Chinese) [8] Frenkel,S m it.分子模拟———从算法到应用[M].汪文川译.北京:化学工业出版社,2002:51-92.[9] Hoagland R G,Da w M S.A t om ic model of crack ti p de2f or mati on in alu m inu m using an e mbedded at om potential[J].Journal of M aterials R esearch,1990,5(2):313-324.[10] 陈 强,曹红红,黄海波.分子动力学中势函数研究[J].天津理工学院学报,2004,20(2):101-105.CHE N Q iang,CAO Hong2hong,HUANG Hai2bo.A re2search on the interat om ic potential in moleculardyna m ics(MD)[J].Journal of T ianjin Institute ofTechnology,2004,20(2):101-105.(in Chinese) [11] 吴兴惠,项金钟.现代材料计算与设计教程[M].北京:电子工业出版社,2002:25-30.[12] Da w M S,BaskesM I.E mbedded at om method:deriva2ti on and app licati on t o i m purities,surfaces,and otherdefects in metals[J].Phys Rev B,1984,29(2):6443-6453.[13] 王 俊,丁 皓,朱 宇,等.并行计算在分子模拟中的应用[J].南京工业大学学报,2003,25(1):9-13.WANG Jun,D I N G Hao,ZHU Yu,et al.App licati on ofparallel computer in p r omoting devel opment of molecularsi m ulati on[J].Journal of N anjing U niversity of Tecnolo2gy,2003,25(1):9-13.(in Chinese)[14] Zhang Tao,Zhang J ing,W ang Hui,Hu Yuanzhong.A2t om ic stick2sli p fricti on bet w een commensurate self2as2se mbled monolayers[J].Progress in N ature Science,2002,12(7):501-505.[15] Refs on K,Moldy.A portable molecular dyna m ics si m u2lati on p r ogra m for serial and parallel computers[J].Co m puter Physics Co mm unications,2000,127(3):310-329.[16] Khol m ur odov K.MD si m ulati on of cluster2surface i m2pacts for metallic phases:Soft landing,dr op let s p rea2ding and i m p lantati on[J].Co m puter Physics Co mm uni2cations,2001,141(1):1-16.(责任编辑 徐红星)54第1期 叶建军等:金属材料动态损伤的微观数值模拟。

用三角形薄板广义协调元分析开孔矩形薄板的弹性稳定

用三角形薄板广义协调元分析开孔矩形薄板的弹性稳定

用三角形薄板广义协调元分析开孔矩形薄板的弹性稳定
李自林;刘兴业
【期刊名称】《工程力学》
【年(卷),期】1996()A01
【摘要】本文利用三角形薄板广义协调元分析了开孔矩形薄板的弹性稳定,其优点是自由度少,精度高,程序简便,文中算例计算了具有内方孔的四边简支方板的临界压力系数,其结果与文献(4)的结果符合良好。

【总页数】4页(P524-527)
【关键词】广义协调元;开孔薄板;临界压力;弹性稳定性
【作者】李自林;刘兴业
【作者单位】
【正文语种】中文
【中图分类】O343.9;TU311.2
【相关文献】
1.缓和开孔矩形薄板应力集中方法的弹性分析 [J], 赵卫平;王道锋;周世军
2.广义协调条件的薄板弯曲三角形元和矩形元 [J], 龙志飞
3.用广义协调元分析复合式多层无孔和开孔矩形薄板的振动 [J], 李自林;王荣霞;刘兴业
4.用三角形薄板广义协调元分析开孔矩形薄板的振动 [J], 李自林;宋德玉
5.两个高效稳定的厚薄板广义协调矩形元 [J], 岑松;龙志飞
因版权原因,仅展示原文概要,查看原文内容请购买。

基于两种梁理论对变幅锥形杆弯曲振动的特性分析及参数设计

基于两种梁理论对变幅锥形杆弯曲振动的特性分析及参数设计

基于两种梁理论对变幅锥形杆弯曲振动的特性分析及参数设计常婷婷;沈峰;鲍四元
【期刊名称】《振动与冲击》
【年(卷),期】2024(43)2
【摘要】为了研究圆锥形杆自由振动的特性,分别基于欧拉-贝努力梁理论和铁木辛柯梁理论,建立变截面杆自由振动的分析模型。

采用一种含三角函数的级数形式来表示欧拉-贝努利梁理论下杆的位移函数,以满足端部位移的条件;利用能量泛函极小化得到系数满足的线性方程组,进而获得不同边界条件下圆锥形杆在欧拉-贝努利梁理论下的若干阶固有频率;类似地,假设位移的级数形式并利用能量函数,建立锥形杆基于铁木辛柯梁理论的求解方法,可得各阶固有频率和模态;给出等截面杆在两种理论下固有频率的转化公式,并推广应用到圆锥形杆的固有频率近似转化。

算例分析锥形杆截面参数对结构固有频率的影响,并基于目标设计频率和若干限制条件对锥形杆的尺寸进行设计。

数值结果表明,在应用欧拉-贝努利梁理论和铁木辛柯梁理论时,所提方法都能够稳定收敛且计算效率高,具有较高的精确度。

该研究工作为超声工程中变幅杆的动力学特性提供了计算依据。

【总页数】9页(P114-122)
【作者】常婷婷;沈峰;鲍四元
【作者单位】苏州科技大学工程力学系
【正文语种】中文
【中图分类】O327
【相关文献】
1.一种圆锥形变幅杆弯曲振动固有频率的计算方法
2.基于Mindlin理论新型阶梯环型变幅器弯曲振动特性研究
3.A2弯曲振动变幅杆的设计
4.超声弯曲振动变幅杆的特性
5.大截面圆锥形超声变幅杆的设计及负载特性分析
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第9卷 第6期2007年11月天津职业院校联合学报Jou rnal of Tianjin Vocational Institutes NO.6Vol.9Nov.2007工程力学建模杨冠声(天津现代职业技术学院,天津市 300222)摘 要: 传统的工程力学教学,直接就力学模型讲授理论的应用,缺乏从工程问题中抽象出力学模型的环节,与应用脱节,不利于学生的技能培养。

在教学实际中,可以根据 可靠性 、 经济性 原则,将工程力学建模分解为 主体的处理 、 力的处理 、 联结的处理 三个基本过程。

关键词: 工程力学;建模;教学中图分类号:P731.22 文献标识码:A 文章编号:1673-582X(2007)06-0012-04收稿日期:2007-03-20作者简介:杨冠声(1969-),男,河南省人,天津现代职业技术学院副院长,副教授,博士,港口海岸及近海工程专业,从事力学教学与研究工作。

基金项目:天津市 十一五 规划课题 技能基础课培养学生应用能力和创新能力模式研究 资助项目一、工程力学教学中增加建模部分内容的必要性力学建模是力学工程应用首先面临的问题,是进行力学分析计算的基础。

大量的问题要在这一过程中完成,包括确定分析方案、分解分析项目、简化分析对象、确定载荷情况等等,它是力学研究、力学实践的重要环节。

对于工程力学而言,这一工作占有更加突出的位置,甚至会占用比力学计算更多的时间,直接决定结论的可信度。

不恰当的力学建模带来的分析误差,会将后续提高计算精度的努力冲减甚至否定。

工程力学是一般工程人员进行简单力学分析计算的工具,一般工程人员通常不具备十分精深的力学素养,达不到一看到实际结构就能在头脑中形成一个相应力学模型的境界,工程力学建模是他们将工程实际和力学理论结合起来的一个必须的环节。

在大学课程教学中,应该加强这一能力的培养,增添这一部分的内容。

工程力学建模是学生理解力学理论的必要过程。

理论是工程力学分析计算的工具,其应用对象是工程实际,同时也是对应用规律的描述。

建模过程可以使学生加深对理论本质的理解,把握理论的核心要素,分清操作中各因素的主次差别,领会在应用中的关键、难点及一般解决方案中存在的问题,方便学生进一步研究与分析深度的加强。

工程力学教学,应该弱化理论推演,加强应用环节,建模过程是达成这一目标的重要手段。

学生的认识过程,有喜欢追根求源、喜欢自己推演的特点,没有经过自己推导的东西即认为 没有理解 、 没有掌握 。

这是中学延续下来的学习习惯。

一旦进入工程应用,大量的理论、规范与材料,使得我们没有时间、没有精力,也不可能对每一信息逐一推演、考证,必须要相信一些东西,特别是列入规范为大家认同的理论。

工程力学具有很强的数学特征,正如一般所说 力学是数学的磨刀石 ,在学习中很容易陷入环环相扣的数学推演,脱离具体工程对象。

增加力学建模部分,有利于使学习过程始终围绕应用这个核心。

力学建模是培养力学方面工程习惯的重要环节。

与工程观念相对应的是工程习惯,应该形成这样一种直接反应,对于问题首先确定分析范畴,查找规范,进行分析,认定参数,计算并得出结论,形成报告。

对于工程力学分析范畴的问题,教材中的理论均符合规范的要求,分析过程首先是力学建模,然后根据力学理论,代入相关参数进行计算,最终得出结论,形成报告。

二、工程力学建模在现有教学体系中的状况现有工程力学教学体系中,从主流教材(包括本科用教材[1,2,3,4]和专科用教材[5,6,7])看,尚没有对力学建模的系统的讲解。

在教材中有部分相关内容,通常是在理论力学篇静力学部分的 约束与约束力 章节中,主要是结构联接处理的内容。

存在以下问题:1.没有明确提出力学建模的概念,没有将之确定为力学分析的一个必须过程,因而也没有系统的阐述。

在后续内容中,直接就力学模型进行分析[8],如图1所示;或就结构示意图进行分析[9],如图2所示。

在教学中由教师对力学简图的含义、代表的结构进行补充说明,或者对结构示意图中特定部分进行力学阐述,如说明图2b中的A点为固定端约束。

整体看来,是直接就力学模型讲授理论的应用,没有从工程问题中抽象出力学模型的环节。

图1 已知q、[ ]确定工字钢的型号图2 一端固定的管线受力分析图2.从约束表示的角度来讲述结构的表达,在处理实际问题时表现出知识的不完整。

例如,在讲述固定铰链的时候,给出实际结构图[8,9](如图3.a所示),然后说明这种结构称为固定铰链,可以用图3.b表示。

对于图4表示的齿轮结构,在齿轮轴的中心点同样可以处理为固定铰链,可以用相应的力学符号表示,学生在学习中通常会在这点发生困扰。

这也是列举法本身的弊病,不可能穷尽所有情况。

从力学建模的角度阐述,可以解决这个问题。

定义所有在特定点能够提供某类限制的结构即为某种约束,例如定义固定铰链为 限制研究主体在基础上特定点不能移动只能绕该点转动的结构 ,可以避免前述困扰。

图3 固定铰链图4 齿轮结构三、工程力学建模的教学探索在实践教学中可以增加相关内容,以解决现有工程力学教学与实践脱节的问题。

1.对工程力学建模的概念进行界定定义工程力学建模位从工程实际中抽象出力学问题的过程。

其基本内涵是从具体工程问题中抽象出力学相关因素,根据力学分析项目,用力学语言进行描述。

按照一般认同的对工程力学内容的界定,工程力学能够完成的分析项目主要有物体的静力分析、运动与动力分析,杆件的变形分析、破坏分析及稳定性分析等。

力学语言包括力学术语和力学符号,力学符号又包括文字符号和图形符号。

力学建模的基本要求,是要把实际问题中的非力学因素剔除,确定相关的基本力学因素。

2.确定工程力学建模的基本原则工程力学建模过程既要保证后期的分析计算结果不失真,也要保证后期的分析计算可以实施。

这就要定义所谓的 可靠性原则 和 经济性原则 。

可靠性原则是指对问题处理后能够反映实际情况,计算精度能够达到工程要求,通常也称为保守原则。

它应该是力学建模的基本原则,也应该是所有力学分析计算的基本原则,是力学在工程应用中得以认可的基础。

工程计算、工程鉴定,往往有连带的经济责任、职务责任和法律责任,工作可靠性是界定责任的重要因素。

可靠性界定的尺度,是工程规范。

要保障力学分析计算的可靠性首先要以工程规范为依据。

工程规范除了工作规范,还包括国家标准系列,这些是力学分析计算的基础。

通用的工程力学教材及权威性的工程手册,其内容都符合规范的要求,是可以信赖的。

如下图所示,图5带切槽的传动轴,在分析轴的破坏时,如果忽略切槽的存在将其简化为图6所示的等直截面轴,不符合可靠性原则,计算结果将不能保证应用的安全要求。

图5带切槽的传动轴图6 不当简化后的模型经济性原则是指适当简化力学模型以保证计算的经济性。

计算中考虑太多的技术细节必然会带来计算时间和成本的增加。

力学计算的成本主要来自人力成本,企业自有人员的项目费用相对较低,研究机构专业人员的使用费用相对较高。

工程初算、保守计算、定性分析、问题排除,对模型简化要求较低,考虑技术细节较少,是一般工程人员可以完成的;而对于关键零部件、精细计算与优化设计往往要由专家来进行。

一般方法是对安全有贡献的部分可以直接忽略,对安全有损害的部分简化后要补偿(以保证结果的可靠性)。

如图7所示工字钢,在应用中通常会在上面焊接一些其它结构(如角钢),它们对整体的强度有贡献。

在一般计算中通常将其忽略(如图8所示),既降低了工作量,又符合安全要求,兼顾了可靠性原则和经济性原则。

图7带筋板结构的工字钢图8 简化后的模型3.分解工程力学建模过程工程力学建模可以分解为 主体的处理 、 力的处理 、 联结的处理三个基本过程。

图9 基础的力学符号主体的处理指确定主体,并根据力学分析项目的不同,将之处理为质点、刚体或变形体,用简单的点、线、面和尺寸标注进行表示。

力学分析计算和工程报告是需要审阅、审批的,为了便于理解与交流,通常要表示出主体的基本特征,能够看得出来是哪种工程问题。

处理为质点和刚体的主体不涉及变形与破坏,通常可以不标示出其截面尺寸。

定义联接为主体与基础之间、各主体之间的结构。

联接对主体的运动有限制作用,通常也称为 约束 。

根据联接对主体的限制作用,用相应的力学符号进行指代。

沿用原有的符号,以保证教学内容与现有的工程习惯一致。

根据对主体限制的差异,可以分为基础与接触面、柔索、活动铰链、固定铰链、中间铰链等。

将基础作为一个单独的联接形式和约束方式进行阐述,定义 与主体相关、相对固定的结构都可以称为基础 。

其力学符号如图9所示,粗实线标示出基础与主体之间的基础面。

在基础的处理上,通常不画出其实际结构,只在主体与基础联接的地方,用基础的力学符号标示出基础的位置。

除直接搁置在基础上的主体外,主体与基础之间联接方式的力学符号中都包含有基础的力学符号。

在力学建模中,联接的力学符号代表提供特定限制的结构,只要具有特定的功能,不管联接的具体结构如何,都可以用相应的符号表示。

在研究中,如果用特定的符号表示不能提供某一方面限制的结构,只要对研究结果没有影响,也是可行的。

考虑一般工程应用的特点,将力的处理根据力作用范围的差异分为集中力和分布力的处理。

处理过程中主要确定力的作用位置、作用方向和在作用面上的大小分布,并用适当的方式(图示或文字表述)进行表达。

4.存在的问题对于实际问题,一般是三维问题,对于实际结构在三位情况下如何表现,还缺少公认的方法。

通常需要转换为多个二维问题,或者画出实际结构的简图。

只有固定端约束、空间铰链、向心推力轴承等少数几种结构有比较公认的表达符号。

这方面还需要进行进一步的研究、规范。

总之,工程力学建模方面的探索是针对现代高等教育重应用的要求,分析现有工程力学教学内容中的缺陷,形成的把实际问题和力学理论相结合的环节。

在实践中证明是比较受学生欢迎的,是实用的,加强了学生技能培养,收到了良好的效果。

相关内容制作的多媒体课件 工程力学建模 在教育部信息中心举办的第四届全国多媒体课件大赛中获得三等奖。

参考文献:[1]张定华.工程力学[M].北京:高等教育出版社,2000.[2]单辉祖.材料力学[M].北京:国防工业出版社,1997.[3]张秉荣.理论力学[M].北京:机械工业出版社,1991.[4]梁治明,等.材料力学[M].北京:高等教育出版社,1992.[5]刘鸿文.材料力学[M].北京:高等教育出版社,1986.[6]哈尔滨工业大学.理论力学[M].北京:人民教育出版社,1987.[7]罗远祥等.理论力学[M].北京:高等教育出版社,1984.[8]西北工业大学.理论力学[M].北京:人民教育出版社,1980.[9]吴福光等.振动理论[M].北京:高等教育出版社,1987.[10][苏]B. .里舍夫斯基.通俗力学[M].北京:科学普及出版社,1984.[11]中国教育部高等教育司.高职高专教育专业教学改革试点方案选编[M].北京:高等教育出版社,2003.[12]中国教育部高等教育司.高职高专教育改革与建设[M].北京:高等教育出版社,2002.[13]范印哲.教材设计导论[M].北京:高等教育出版社,2003.[责任编辑:冯金泉]Engineering M echanics M o delingYANG Guan-sheng(T ianj in M odern Vocational&T echnical I nstitute,T ianj in300222China)Ab s tr a c t:The traditio nal teaching of engineering mechanics instructs theories on mechanical mod-el directly.The theo ries are disarticulated from application,for there's no ne procedure abstracting me-chanical mo dels from engineering problems.The article presents two basic principles of reliability and eco nomy in teaching practices,and divides the modeling procedure into three steps of processing of main bodies , processing of forces and processing of co uplings .Ke y w o r ds:engineering mechanics;modeling;teaching。

相关文档
最新文档