ZYF

Month: 六月, 2007

通过Scheme进行Fluent计算结果自动处理

最近用Fluent计算非定常流场,保存了大量的非定常结果,各个时间步分别对应一个cas和一个dat文件。现在希望可以将各个时间步上的某个边 界上的压力分布输出。如果计算之前考虑到所需要的结果,可以在Fluent中设置按需输出结果。但是如果没有考虑周全,而需要事后处理的话,面对大量的 cas/dat文件,就有必要借助于Fluent中内建的scheme解释器来处理了。
下面的这个代码“output.scm“,应存放在当前运行目录下,在Fluent中运行方法: 提示符下输出 (load “output.scm”)
这个代码假设文件名的规律为 “bp10x-un-” + 5位小数(时间),从时刻“0.01147“开始。文件名中的时间以0.0001递增。共处理10组文件。输出ascii格式的边界“wall-bottom“和“wall-top“上的压力分布。
目前,在相应文件名处理上还有一点问题,scheme的 number->string 和 format 函数功能都有限,造成如果时间最后几位为0,这些0会被忽略掉。还在想办法处理。
我是scheme完全新手,写出这个也不容易啊。
—— output.scm ——
(define prefix “bp10-un-”)
(define tstart 0.01147)
(define timestep 0.0001)
(define steps 10)

;; make postfix part of filename
(define postfix
(lambda (n)
(let ((flonum-unparser-cutoff ‘(absolute 6 scientific)))
(number->string (+ tstart (* timestep n))))))

;; loop over all cas/dat [...]

夜未央

谭咏麟 – 夜未央
轻轻踏在月光里 好像走在你的心事里
那年黯然离别后 再也没有人与我同饮
飞花轻似雾奈何风吹起 终究如烟纷飞东西
细雨细如愁忘了看个清楚你眼中默默深情
雨中路遥遥梦里风萧萧 仿佛中你在微笑
漫漫长夜里梦醒的太早
想起我轻狂的年少
无声又无息花落了满地 只留下纷芳依稀
蓦然再回首梦还是一样
为你等在夜未央 不知风寒

实用断裂力学

David Broek. The Practical Use of Fracture Mechanics. Kluwer. 1988.
这 本书从工程使用的角度出发介绍断裂力学的使用,书中没有涉及过多的理论知识,大部分的篇幅是讲解如何在实际分析中使用断裂力学,并解释一些容易混淆的概 念。因为不涉及数学公式等复杂的理论基础,对于工程师学习和使用断裂力学有很大帮助。另外,这本书也可以作为参考书使用。
刚刚翻看了第一章的基本介绍。裂纹(crack)、断裂(fracture)
任 何结构都不可避免的具有一定缺陷,这些缺陷将成为机械运行时问题的根源。细小的裂纹在运行中受到几种载荷形式作用,逐渐增大,当裂纹增长到一定阶段,其 residual strength 降低到与 permissable strength 同等水平的时候,就会发生断裂。
当然 crack growth 具有不同的模式,fracture 也具有不同的模式。
crack tip stresss 来源于三种 loading:

opening mode, tension mode
shear mode, in-plane
tearing mode, out-of-plane

最主要起作用的是第一种 opening mode,另外两种也会产生,大部分是和 opening mode 同时作用的。而基本上 crack 的发展方向也是 opening mode 的方向。这三种模式具有相同的作用方程,可以举一反三。
通过 Damage tolerance analysis 获得 residual strength diagram (fracture due to cracks) [...]