有限元分析中的单位问题..doc

有限元分析中的单位问题 葛颂（浙江大学化机所） 2004年11月版 摘要： 本文在前人基础上对使用有限元软件分析工程问题时的材料性能单位问题作了一些探讨，通过实例说明了如何统一各物理量的单位，以保证分析结果的正确。 ? 关键词：有限元、单位 ? 大多数有限元计算程序都不规定所使用的物理量的单位，不同问题可以使用不同的单位，只要在一个问题中各物理量的单位统一就可以。但是，由于在实际工程问题中可能用到多种不同单位的物理量，如果只是按照习惯采用常用的单位，表面上看单位是统一的，实际上单位却不统一，从而导致错误的计算结果。 比如，在结构分析中分别用如下单位：长度 – m；时间 – s；质量 – kg；力 - N；压力、应力、弹性模量等 – Pa，此时单位是统一的。但是如果将压力单位改为 MPa，保持其余单位不变，单位就是不统一的；或者同时将长度单位改为 mm，压力单位改为 MPa，保持其余单位不变，单位也是不统一的。由此可见，对于实际工程问题，我们不能按照手工计算时的习惯来选择各物理量的单位，而是必须遵循一定的原则。 物理量的单位与所采用的单位制有关。所有物理量可分为基本物理量和导出物理量，在结构和热计算中的基本物理量有：质量、长度、时间和温度。导出物理量的种类很多，如面积、体积、速度、加速度、弹性模量、压力、应力、导热率、比热、热交换系数、能量、热量、功等等，都与基本物理量之间有确定的关系。基本物理量的单位确定了所用的单位制，然后可根据相应的公式得到各导出物理量的单位。具体做法是：首先确定各物理量的量纲，再根据基本物理量单位制的不同得到各物理量的具体单位。 基本物理量及其量纲： ????????? 质量 m; ????????? 长度 L; ????????? 时间 t; ????????? 温度 T。 导出物理量及其量纲： ??????? 速度：v = L / t; ??????? 加速度： a = L / t 2; ??????? 面积： A = L 2; ??????? 体积： V = L 3; ??????? 密度：ρ= m / L 3; ??????? 力： f = m · a = m · L / t 2; ??????? 力矩、能量、热量、焓等： e = f · L = m · L 2 / t 2; ??????? 压力、应力、弹性模量等： p = f / A = m / (t 2 · L) ; ??????? 热流量、功率：ψ= e / t = m · L 2 / t 3; ??????? 导热率： k =ψ/ (L · T) = m · L / (t 3 · T); ??????? 比热： c = e / (m · T) = L 2 / (t 2 · T); ??????? 热交换系数： Cv = e / (L 2 · T · t) = m / (t 3 · T) ??????? 粘性系数： Kv = p · t = m / (t · L) ; ??????? 熵： S = e / T = m · l 2 / (t 2 · T); ??????? 质量熵、比熵： s = S / m = l 2 / (t 2 · T); 在选定基本物理量的单位后，可导出其余物理量的单位，可以选用的单位制很多，下面举两个常用的例子。 1 基本物理量采用如下单位制： ????????? 质量 m – kg； （应该采用Mg单位才统一，具体可以参考MSC.MARC中的材料库统一单位，推导） ????????? 长度 L – mm； ????????? 时间 – S； ????????? 温度 – K (温度 K与 C 等价)。 各导出物理量的单位可推导如下，同时还列出了与 kg-m-S 单位制或一些常用单位的关系： ??????? 速度： v = L / t = mm / S = 10-3 m / S; ??????? 加速度： a = L / t 2 = mm / S 2 = 10-3 m / S 2; ??????? 面积：A = L 2 = mm 2 = 10-6 m 2; ??????? 体积：V = L 3 = mm 3 = 10-9 m 3; ??????? 密度：ρ= m / L 3 = kg / mm 3 = 10-9 kg / m 3 = 10-6 g / cm 3; ??????? 力： f = m · L / t 2 = kg · mm / S 2 = 10 –3 kg · m / S 2 = mN (牛); ??????? 力矩、能量、热