理论力学、材料力学和结构力学等力学课程是土木工程及相关专业的重要基础课，而随着新课程的不断开发，力学课程的课时一再降低.力学课程涉及到的大量概念和计算并不容易掌握，在有限的课时内要达到理想的教学效果是任课教师面临的挑战.随着计算机硬件的高速发展以及计算算法的不断进步，计算软件也越来越丰富，很多用于工程计算和科学研究的有限元分析软件，可以方便地用于力学教学中.正如钱学森所说 “今日力学是一门用计算机计算去回答一切宏观的实际科学技术问题，计算方法非常重要；另一个辅助手段是巧妙设计的实验”[1].CAE软件在力学教学中的应用就印证了这句话.CAE软件作为工具可以巧妙地设计力学教学中的实验 [2]，节约了实验资源和实验时间.应用广泛的 CAE软件有ANSYS，ABAQUS，LS-DYNA，ADINA，COMSOL Multiphysics，FLUENT，ADMAS，NASTRAN 等.李丽君等 [3]将 ANSYS的屈曲分析应用到材料力学压杆稳定教学中，使抽象的欧拉公式更加容易理解.程学晶等[4]利用Ansys Workbench软件得出位移和应力的分布云图，展示了材料力学的相关概念和定理.王俊奇等[5]把ABAQUS引入土力学教学中，结合几个案例，提高了学生分析问题的能力，并培养了创造性思维.徐耀玲等[6]在LS-DYNA平台上对小球自由落体冲击简支梁进行了数值模拟，得到了梁的位移和应力的时程曲线，取得了良好的教学效果.

LS-DYNA软件以 Lagrange算法为主，兼有ALE和 Euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体--结构耦合功能；以非线性动力分析为主，兼有静力分析功能；它的显式算法在动力学计算方面具有优势[7].LS-PREPOST是专门为LS-DYNA定制的前后处理软件，完全支持LS-DYNA的全部关键字.对力学中的简单模型可以方便地建模，并能够形象地演示分析结果以及提取数据进行结果处理.本文以两个基于LS-PREPOST和LS-DYNA建模和计算的案例说明这两个软件在力学实验设计和科研成果转化经典案例两个方面的探索应用.

1 力学实验--牛顿摆球实验

1.1 计算模型的建立

牛顿摆是5个质量相同的球体由吊绳固定，彼此紧密排列，可以验证动量和能量的守恒.牛顿摆是由法国物理学家伊丹·马略特最早于1676年提出的.摆动最右侧的球在回摆时碰撞紧密排列的另外4个球，最左边的球将被弹出，并仅有最左边的球被弹出.这个经典实验经常出现在理论力学的教学中.

就矿床的成矿空间与成因来看，这类矿床一般都与深层的白岗质花岗斑岩岩体有关。大部分的矿化富集部位都在花岗斑岩岩株的隆起部位以及岩体的外接触带，最远可达1 000～1 500 m，但矿体的厚度很少超过500～700 m。

  

图1 牛顿摆球实验原型

图1显示了牛顿摆球的实验原型.图 2为模拟实验，在 LS-PREPOST软件中建立 5个直径100mm的钢球，每个钢球由夹角为60°长500mm的钢丝线悬挂，钢丝线直径 5mm，最右边的钢球与两根悬线组成的平面与垂直于纸面的竖直面夹角为 30°，在重力 (关键字为 LOADBODY Z)的作用下，最右边的钢球静止释放后开始运动.为了提高计算效率，钢球和钢丝悬线均采用刚体材料(MAT 020/MATRIGID)，悬线与钢球的刚体连接采用关键字 CONSTRAINEDRIGIDBODIES，模拟悬线与钢球的焊接连接方式.悬线上端由两根水平钢梁固定，水平梁的材料为线弹性材料(MAT 001/MATELASTIC)，直径为 5mm.钢球采用实体单元 (ELEMENTSOLID),如果采用壳单元模拟钢球，需要考虑壳单元的厚度，以防产生单元的初始侵彻，悬线和水平钢梁均采用梁单元(ELEMENT BEAM).碰撞类型采用关键字 CONTACT AUTOMATIC SINGLESURFACE.边界条件为：约束悬线顶端与两根水平梁交点处节点的3个坐标轴方向的位移，对应的关键字为BOUNDARY SPCSET.为了计算的稳定性，施加了数值为0.001的系统阻尼(DAMPINGGLOBAL).计算时间为1.5s.计算模型如图2所示，为方便叙述，小球从右往左依次编号为1,2,3,4,5号.

  

图2 牛顿摆球计算模型

1.2 牛顿摆球碰撞结果分析

牛顿摆球的碰撞过程如图3所示，从1号球释放开始，0.36s碰撞静止的2号小球，1～4号小球均保持静止，5号小球获得初速度开始运动，5号小球0.7s达到最高点，速度为零，开始回落，到1.06s碰撞4号小球，随即5号小球停止运动，与2,3,4号小球保持静止，1号小球开始运动，1.42s达到初始位移，完成一个运动周期.由于数值计算的误差和微小阻尼的作用，小球运动状态并非全完理想的状态，能够观察到静止的小球有轻微的水平摆动.图4显示了5个小球的速度时程曲线，图5显示了1号和5号球的竖向位移曲线(其他小球的竖向位移均为零).从图4和图5可以看出，初始状态1号摆球从最高点竖向位移63.8mm开始释放(以静止小球的质心位置为竖向坐标原点)，1号摆球0.35s到竖向坐标原点，速度达到最大值1.11m/s.此时5号摆球获得初速度1.08m/s，由于2,3,4号摆球的轻微震荡，5号摆球速度略有所损失，5号摆球0.7s达到最高点59.7mm，小于1号摆球的初始位移.5号摆球1.06s碰撞4号小球，1号小球同时获得速度1.13m/s，小于5号摆球的初速度.图 6显示了系统动能、势能以及动能与势能之和的时程曲线，动能的峰值处恰好势能为零，反之亦然，动能和势能之和基本保持恒定.

以桥梁上部结构碰撞系数研究[8]为例,说明科研成果在力学教学中的应用，研究过程中进行的桥梁模态分析涉及到《结构动力学》包含的与模态有关的相关概念.桥梁的护栏、桥面板和支撑桥面板的大梁均采用分层壳单元进行模拟，分别采用不同厚度的壳单元截面类型，其中大梁和板的单元类型交替布置.材料类型为线弹性材料.支座模型采用LS-DYNA软件中的离散梁单元模拟.这种单元类型可以方便地定义3个坐标轴方向的拉伸刚度和旋转刚度.桥梁模型的边界约束条件施加在支座底部节点上，并且忽略支座与大梁和墩台的摩擦.材料常数的选择在文献[8]有详细的描述.图7只显示了与实验结果相关的3个模态形式和对应的频率.在LSPREPOST平台上可以演示模态分析的全部结果，展示不同平面的位移模态动画，给学生更清晰的概念.同时传达学生本课题研究的背景和目的以及思路和方法，在课堂上展开讨论，引发学生更进一步的思考.

  

图3 牛顿摆球的运动过程

  

图4 5个摆球的速度时程曲线

  

图5 1号摆球和5号摆球的竖向位移时程曲线

  

图6 系统动能、势能和能量和的时程曲线

2 科研成果转化成力学教学中的经典案例——模态分析

LS-DYNA在力学教学中可以方便地建模、分析和演示.不仅可以模拟碰撞问题，还可以进行压杆稳定的屈曲分析、结构内力分析，不仅可以进行动力分析，也可以进行静力分析.再现力学教材中的实验和案例能够给学生更直观的印象，加深记忆和理解.另一方面，教师科研成果可以转化成教学中的经典案例，促使学生更深入地认识问题、思考问题，激发学生学习力学的热情和从事科研的好奇心.把 LSDYNA等有限元软件和前后处理软件作为教学辅助工具引入力学教学中，能够节约时间和实验资源，并能取得良好的教学效果，值得借鉴和推广.

这个数值算例可以方便地扩展到更多工况，比如两球、三球、四球同时在一定高度释放，碰撞其他小球，或者模拟总数少于五球的碰撞实验.通过牛顿摆球的力学实验模拟，能够使学生获得更加清晰的动量和能量守恒的力学概念，并且增加了教学的趣味性.

国家局提出数据中心建设，归根到底是要为行业发展和企业决策提供及时、准确、系统的信息支撑，因此，对数据进行统一管理显得尤为重要。一是要保证信息数据源头的质量，即各应用系统的数据是严格按照相应的规范采集的，是依法、依规、依标准的可加工和交换的数据；二是要建立数据质量的控制体系，数据的加工、存储、交换的过程严谨有序，不出差错，以保证数据中心的数据质量；三是要按照“授权访问，各取所需”的规则发布数据，避免信息发布的失控，做到“数入一库，数出一门”；四是要采用专业的统计分析技术，对数据中心海量的数据信息进行深度的挖掘和分析，将数据资源提炼为决策资源，以提高信息资源的增值服务。

  

图7 桥梁模态和频率的有限元结果和实验结果的对比

3 小结

LS-DYNA不仅可以重现力学书中的实验，也可以把与力学概念相关的科研成果展示给学生.把科研成果转化成为教学中的经典案例，也是科研成果价值的体现.高等院校任教的教师，工作内容不仅包含教学，同时也进行课题研究，一些与教学相关的研究思路、研究过程和研究结果可以作为课堂授课的补充内容传达给学生.不仅开阔学生的视野、提高对本专业的学习热情，也激发对科学的热爱、向往和好奇心，为将来进一步深造奠定基础.

调查显示，90%左右的大学生对古建筑具有浓厚的兴趣，其了解古建筑的方式依次为学校课堂学习（76.25%），互联网、社交媒体等新媒体（63.71%），旅游活动（61.37%），学校文化宣传交流活动（58%、86%），课外书刊等传统媒体（54.18%）,社会实践交流活动（50%），家庭长辈亲朋的影响（40.47%）以及历史文化氛围、电视节目等。反映了对大学生进行传统文化知识教育和文化自信培养，需要采取校内外、线上线下、整合各种力量共同进行。

参考文献

1 钱学森.我对今日力学的认识.力学与实践,1995,17(4):1

2 宋少云.将 CAE引入力学教学.力学与实践，2006,28(1):74-75

3 李丽君,沈玉凤.ANSYS对于压杆稳定的辅助教学.实验科学与技术,2008,7(1):121-122,125

4 程学晶,杨家军,谭丽辉.Ansys Workbench软件在材料力学教学中的探索与应用.科技教育,2015(24):171-172

5 王俊奇,颜月霞.土力学课程教学内涵及依托ABAQUS软件的研究性教学改革.中国冶金教育,2016,177(6):19-25

6 徐耀玲,谭文峰.用 LS DYNA辅助《材料力学》冲击问题的教学.力学与实践,2008,30(6):87-88

7 尚晓江,苏建宇,王化峰.ANSYS/LS-DYNA动力分析方法与工程实例,第2版.北京:水利水电出版社,2008

8 王娟,钱江.重型车辆与桥梁耦合振动及桥梁冲击系数的数值模拟研究.结构工程师,2014,30(5):35-41