ABAQUS地震波资料.doc

提供常用的P50%10（50年超越概率10%），一般的工程设计地震常用这个,时间增量0.02秒。*Amplitude, name=HAMPX 0.02, 0.014, 0.04, 0.014, 0.06, 0.014, 0.08, 0.065 0.1, 0.014, 0.12, 0.016, 0.14, 0.219, 0.16, 0.13 0.18, 0.082, 0.2, 0.3, 0.22, 0.583, 0.24, 0.129 0.26, -0.263, 0.28, -0.948, 0.3, -0.105, 0.32, -0.524 0.34, -0.952, 0.36, 0.088, 0.38, 0.843, 0.4, 1.152 0.42, 1.716, 0.44, 2.523, 0.46, 0.07, 0.48, -1.69 0.5, -0.708, 0.52, -1.42, 0.54, -1.807, 0.56, -1.091 0.58, -1.674, 0.6, -2.547, 0.62, -1.639, 0.64, -2.514 0.66, -5.463, 0.68, -5.08, 0.7, -5.128, 0.72, -6.955 0.74, -7.118, 0.76, -5.805, 0.78, -3.695, 0.8, -1.871 0.82, 3.558, 0.84, 6.373, 0.86, 4.406, 0.88, 5.769 0.9, 10.47, 0.92, 11.534, 0.94, 10.337, 0.96, 12.44 0.98, 6.454, 1., 8.596, 1.02, 5.458, 1.04, 6.403 1.06, 0.05, 1.08, 1.007, 1.1, -5.859, 1.12, -9.448 1.14, -6.851, 1.16, -8.897, 1.18, 12.645, 1.2, 16.45 1.22, 13.529, 1.24, 12.146, 1.26, 16.093, 1.28, 10.12 1.3, -9.287, 1.32, 24.022, 1.34, 22.118, 1.36, 21.657 1.38, 17.831, 1.4, -1.39, 1.42, 10.005, 1.44, 8.174 1.46, 4.502, 1.48, -2.972, 1.5, -7.108, 1.52, -8.6354 1 U! U% y/ S- k, A& S Y# G0 h( i: ?3 为方便大家使用，已经将其转化为标准的ABAQUS 输入格式，数据文件是加速度，加速度单位是cm，请在加界中按0.01缩放！在INP中加入以下字段：; P, zC: b) V! Ce2 y3 D( 1 g( m- k; D/ Y8 t *AMPLITUDE, NAME=HAMPx, INPUT=X.inp8 j& N6 F4 I$ t& q. + y/ k+ p% I *AMPLITUDE, NAME=VAMPy, INPUT=Y.inp$ r. X3 d8 j L/ x6 c5 c *AMPLITUDE, NAME=VAMP, INPUT=Z.inp7 KH+ b, h P X7 k5 c0 h-8 C5 K( s) f- # s6 x$ L1 V A + k) 8 Q, e/ X *Boundary, op=NEW, amplitude=HAMPx, type=ACCELERATION) l0 |9 F$ K$ f2 C “定义的约束集合名”, 1, 1,0.01（红字） D+ y! e$ e. j- N0 s; T场地土层反应计算采用的输入地震波是以地震危险性分析结果得到的基岩加速度峰值和基岩加速度反应谱¾基岩地震相关反应谱作为目标谱，用人工数值模拟方法合成得到的，并以此作为场地地震反应计算的输入地震波。; y! |4 E( R6 3 p+ 8 C按照本章6.2.1的技术思路，利用计算机自动产生的不同随机相位谱，按50年超越概率10%、5%和100年超越概率2%三个概率水准合成基岩加速度时程，其中每一个概率水准合成三条时程，分别对应于三个不同的随机相位，时程采样步长为0.02秒，目标谱与合成时程的反应谱(计算谱)之间的相对误差均小于5%。目标谱与计算谱的比较、拟合精度、基岩地震波图见图。 N l L4 i O8 F) h) yg8 J- - t1 2 i8 j+ D8 u1 s 2010-12-9 16:35:10 上传下载附件 (70.51 KB) * b: t: P6 H) v8 q6 K 2010-12-9 16:35:10 上传下载附件 (73.42 KB) , N( O$ P# g2 M; v$ k9 7 L! s9 K离散塔身，对进口段使用三维实体单元整体建模，模型示意图见插图X-7，主要参数：顺水流方向（X轴）长30.0m，上游边界23.4m，下游边界23.4m，塔体总高86.0m，模型节点总数11426个，C3D8R（为获得时间与精度的平衡而选用）单元总数8912个。) b; T/ J) u/ r% A$ X) g A2011-1-18 10:15:25 上传下载附件 (23.59 KB) % S) C) I9 z, 1) 计算荷载和计算工况6 B y i5 X2 i1 i N( / mz, R% H 主要考虑的荷载有：自重荷载、正常蓄水位水压力、地震加速度。5 G, z( 7 a6 h, u* e 计算工况与分析时序相同，共三个工况：工程完建工况(自重荷载施加)正常蓄水工况(静水压力施加)正常蓄水+地震工况(动水压力施加和地震加速度施加)。! T7 e* x% k. |% F) l; 5 j2) 地震波输入7 C1 j( R! u+ D% I* ) z计算地震波输入采用地震安评报告所提供的闸址(基岩)场地谱人工合成的地震波，概率水平为50年超越概率10%，地震波时程见插图X-X。; D, k- b1 T# W# w6 a! l计算时，在模拟岸塔式进口施工完成及蓄水过程的基础上，假定运行期某一时刻发生地震，同时输入顺水流向、垂直水流向和竖直向地震曲线，本次计算将No.1、No.2分别作为顺水流向和垂直水流向输入，No.3作为竖向输入。并将竖直向峰值加速度调整为水平向的2/3倍，即顺水流向、垂直水流向、竖向的峰值加速度分别为100cm/s2、100cm/s2、66.7cm/s2。1 8 D: P; l% w+ U4) 计算成果及分析0 r0 j. s1 9 N5 O. U 位移成果：3 X1 a. h+ I3 i* T9 v蓄水期闸顶(相对于闸基岩)顺水流方向的水平位移为-1.1mm，蓄水+地震工况下最大水平位移达15.0mm，见插图X-8。蓄水期闸顶沉降2.9mm，蓄水+地震工况下沉4.8mm，增加沉降约1.9mm，见插图X-9。: c3 w# A3 u4 j2011-1-18 10:15:25 上传下载附件 (96.36 KB) * D( q, ( n5 Q# F3 p/ e4 K) K* o0 T（一个完整的位移时程如图所示，从小到大，再趋近于0（永久大变形除外）0 H f, F. ( f 2011-1-18 10:15:25 上传下载附件 (80.56 KB) 应力成果：: b* r! V4 z9 N6 8 泄洪放空洞进口竣工期：拉应力主要出现在胸墙底部，最大值约0.80MPa(见插图X-10)，压应力主要分布于进口两侧边墙与底板底部，最大约2.68MPa；蓄水期：拉应力极值出现在工作闸门槽两侧边墙处，约1.07MPa(见插图X-11)，压应力主要分布于底板底部，极值出现底板末端靠山侧，约为3.71MPa；蓄水期+地震：拉应力极值出现在底板前端靠进口处，约1.28MPa(见插图X-12)，压应力极值仍然出现在底板末端靠山侧，约5.89MPa。# 9 S. i P$ L* s2011-1-18 10:15:25 上传下载附件 (70.4 KB) 3 d7 r/ n! d4 d# e9 d! y: t4 f; T0 q# V0 _ R2 U+ K t; A+ v$ M7 R% q& v- X需要说明的是，在地震时程步中水平加速度未计入Y方向（实际意义不大）。% q! A$ J0 P; e) eFGC-30S4.rar (221.68 KB, 下载次数: 118, 售价: 20 仿真币) 2011-1-18 10:15:25 上传下载次数: 118 售价: 20 仿真币 记录 购买 该INP中已包含马尔康加速度，时间步等可自行调整。完整计算时间约为6个小时（电脑E5200，4GB）。2 u2 X6 T- V/ - x& p粗看了一下inp文件，觉得有以下几个方面可以予以考虑： B9 |( J7 M. s$ fJ9 h（1）模型中未见阻尼设置选项，是否忽略了阻尼；4 Q/ h2 p ) T1 y* t+ ( O（2）最后的dynamic分析是否应该采用自动步长比较好，即增加haftol选项: j S) 6 Fa6 I*Step, name=Dynamic, nlgeom=YES, inc=30008 JIA* u1 / s*Dynamic- A0 l) i( 7 t3 EL, W+ k0.1,30.,0.00031 K( L& Y4 g+ HSuggested initial time increment. For implicit integration, this same time increment will be used throughout the step unless contact impacts or releases occur or the automatic time incrementation scheme is used. D7 , HC8 A/ q时程步长为0.02s5 u& ?! i2 s; d1 y- V% z+ c8 _dynamic的参数如下是否好些：* ?, S4 Z: N0.02,33,1e-3,0.02- Y P$ E$ $ , t x6 r/ |（3）dynamic分析中是否应该释放掉u1,u2,u3, L- b! f! j# s) b, F2 o8 M. a这一点很奇怪，呵呵，不知道不释放的话算出来是什么结果（相对的？绝对的？）7 R9 Y4 h- V/ h7 I! L0 _a问题2关于dynamic步中约束的释放，ABAQUS自动给完成了，这里我做过一些试验可以有这个结论：$ H( k( l& I# X_& I为解决“震飞”现象，将模型以INP形式导入以便直观查看，导入后可以发现ABAQUS6.10自动修改在动力分析步上，自动释放了前一步的约束，从计算结果文件可以看出，动态分析步即使加约束与不加约束，所得到的结果完全一致。; l( 9 C$ I R因此可以总结为：对于三维空间模型，选择acceleration/angular acceleration类型输入地震动时，采用CAE模型接交作业时，有同方向（同作用点）的位移约束将报错；采用INP提交时不管有无约束均能计算，其成果完全一致，可以认为程序自动忽略加速度方向的位移约束。故建议在动力分析步解除地震动方向的约束，以避免不必要的麻烦。! 1 L2 H! h& X f/ |2 X关于地震时程分析模型位移过大问题（模型飞走），不知有没有解决了的。老问题再谈。我是一名设计院工作者，面对的是实际工程的地震时程分析之前也做过不少闸室之类的分析。! O6 m% v F6 3 f) k! M% u+ % u0 w. ?7 k2 G! $ R当地震输入数据是加速度时，地震动的输入的前提是要释放该方向的约束，如不释放上一个step该边界的约束，该动力分析步会提示过约束的错误，导致计算中断：7 b4 D) k- M8 C- Z; x0 q; H, Y9 ?( v6 v: y6 Z* t从导入INP后ABAQUS6.10自动修改后的文件可以看出在动力分析步上，ABAQUS自动释放了以前的约束：7 f- b : 1 j% K- a+ j*AMPLITUDE, NAME=HAMPx, INPUT=H.inp0 G) J W7 v?*AMPLITUDE, NAME=VAMP, INPUT=V.inp+ w8 C * u* U* B7 ?% D8 X(数据格式一排四组：时间,加速度时间,加速度时间,加速度时间,加速度)( l% h1 j+ n7 M) B% u% M* yw( R w 3 e1 3 K/ F* STEP: Dynamic8 E$ s L+ 6 XP5 v*Step, name=Dynamic, nlgeom=YES, inc=3000b 6 E5 U2 v! J, c: D*Dynamic3 4 7 O1 Q7 b) Y( Q0.02,30.,0.0003 I/ jA# L7 p3 0 8 I( ? V*Boundary, op=NEW, amplitude=HAMP, type=ACCELERATION3 x; TF. OA! C, U* p2 K$ W! L* W_PickedSet94, 1, 1,0.01s0 E( x4 Y& i（这个0.01是因为我的地震数据单位是cm，将基调整为m）, g! pf9 B1 # _1 y- k. o2 *_PickedSet94, 2, 2,0.01/ ! S# q5 z7 f) 6 H! o（为实验，专门去掉了2方向的acceleration，因此下图中为一条直线）$ B0 V: E0 Y9 : . a B*Boundary, op=NEW, amplitude=VAMP, type=ACCELERATION# T9 G6 z9 _# D- t2 f( y0 _PickedSet94, 3, 3,0.01; p9 x5 c9 F% j$ ?; F*Boundary, op=NEW/ ( k4 I- A1 y6 X: h1 c8 o6 （从这里可以看出，ABAQUS自动释放了以前的约束）1 d D8 0 T5 H$ p( $ N# 4 A: n0 g- s. P. Y! D( C) i+ g, W这里问题出现了，当地震时程很长时，模型的位移会随增量步的增加变得越来越大，从下图可以看出一个电站建在基础上的取水口直接跑了12m，这与事实是完全不相符的（即使震毁也应该是原地破坏）。一般做法是用特征点（自定）减去基础点的位移值来获得结构在地震动作用下的变形，如图2。2 u; H5 E. # d0 Q: U g8 |- h5 7 y4 j* j: J5 l Q( F) n7 f; s引用shanhuimin923 版主的原话：& i/ c. V: R0 ?7 U$ “对于三维空间模型，选择acceleration/angular acceleration类型输入地震动时，建议尽量解除地震动方向的约束，尽管有时不解除约束计算也有可能进行，但更多的时候程序会报错。无论是否解除该方向的约束，只要能进行计算，结构都是正确的，对结果没有影响。”( b f3 hH+ C& Dp& : % w3 b: B) i# b1 U我现在很想弄出在原地震动的结果，即一个点无需操做XY数据（减操做），直接输出波形文件，这样感觉更合情合理。3 M3 F5 q; _4 A8 F3 D* q O, p/ B : ! o6 I5 k- k7 N7 N/ U9 S# U- U) ! o3 p8 4 g下图1为特征点历时（合计值）的变位，可以看出飞了12米。$ E4 p: qC( T ; M, t$ x# I. X8 t- N$ e, J4 _* y7 k下图2为特征点减底部的历时变位（看出很合理）9 _4 L& W G q E- u; g) P. p. i; d6 ) L下图3为震飞的样子。下图3为震飞的样子。未命名.jpg (32.22 KB, 下载次数: 2) 2010-12-9 12:51:03 上传下载次数: 2未命名2.jpg (59.88 KB, 下载次数: 4) 2010-12-9 12:51:03 上传下载次数: 4未命名3.jpg (35.17 KB, 下载次数: 1) 2010-12-9 12:53:21 上传下载次数: 1CEL-水滴入水的例子Abaqus Release Notes.jpg (31.25 KB, 下载次数: 53) 2008-12-15 07:41:29 上传下载次数: 53waterdrop.gif (121.12 KB, 下载次数: 45) 2008-12-15 07:41:29 上传下载次数: 45part與assembly module:6 j+ f$ O G5 | A) M& h由於是水滴入水,所以全是eulerian element,並edit為eulerian part;如果有Lagrangian element,需要另外建立網格,asembly時可以接受干涉.+ I+ P4 D7 s- Kmaterial module:V) w: 1 X( 1 X5 X. ?水的材料參數中,其中C0為聲波在水中之速度,s與Gamma0小弟並不了解,只是參考作設定,如有高手了解,請幫忙作解釋.8 U$ u. a: Dsuboptions中的Eos shear中的Viscosity設定似乎並無定數,估計需依題目調整. X ?- R7 A7 G# M+ vstep module:) C( U9 Q% o A+ Q6 A在field output中選擇EVF才能觀察eulerian element含material的部份.M( kI+ l! n/ e( z: y注意eulerian element無法使用mass scaling.7 Q* R8 c: n) d f; rInteraction module:( E) M) z* I, x如全為eulerian element,則不需做contact interaction;如同時有Lagrangian element,設需設定general contact, contact property的friction,原廠洗衣機的例子是設定為rough,請參考.) K- y& Qa. 0 oa1 T4 r4 PLoad module:# & k; E* M g k# M在eulerian element的邊界需設定velocity為0以防止material流出., H( P% Y* L, i8 P: J) W. E* Zpredifined fieid中需設定material assignment指定eulerian element的volid及material區域0 G x7 q; w& H: A# R; T1 F8 H在例子中只讓水滴落下的速度為100 mm/s,各位可以自行調試看看是否能否激起水花. 2 I) w/ F7 y! oMesh module:4 t0 wn& B. U+ R此處調低了relax stiffness使element較為柔軟,並調低了bulk visocity使應力波較為尖銳,各位可以自行調試看看是否能否激起水花.2 7 g1 C8 n4 o+ & $ - MVisualization module:% + 9 N& c* a0 X6 D, p在results/options/改勾選compute scalars after averaging.; Z: K/ A0 v; i y; 0 L在tools/viewcult/manager中勾選EVF_VOID即可顯示eulerian element含material的部份.predifined fieid中需設定material assignment指定eulerian element的volid及material區域4 9 / I7 X# f2 x/ L. . v* B, i/ C应用abaqus的SPH技术分析【弹体穿靶模型】以及教程1：单元的转换问题- ; ? a. y+ ; Jc2：接触设置问题: 4 ?1 u, Z& # U5 u/ k3：同一物体不同区域（SPH区与非SPH区）的衔接问题。9 B) D: W7 b+ F- S j6 如下是模型计算结果：/ y3 L8 ) w. F2 . F8 I/ i2011-9-10 04:40:47 上传下载附件 (38.54 KB) 2011-9-10 04:41:02 上传下载附件 (49.94 KB) 2011-9-10 04:41:15 上传下载附件 (51.59 KB) 2011-9-10 04:41:26 上传下载附件 (50.77 KB) 2011-9-10 04:41:50 上传下载附件 (56.72 KB) ! b, y z4 v* b B K* b: p; Y3 ) c6 c2 f& J t+ P6 Z. N) y$ 4 % I: l9 U 2011-9-10 04:40:08 上传下载附件 (545.14 KB) 2011-9-10 04:36:22 上传下载附件 (347.88 KB) 9 i, g+ H2 w3 # $ xP( E7 X/ v( h: r / S- xU* t# Y0 O V( 如下是弹体的速度与加速度变化曲线：& t2 o6 S& r: ; D7 m2011-9-10 04:44:50 上传下载附件 (35.1 KB) 有限元网格准备，本模型的网格是在hm中导入，见附件的HM文件。 sph_hm.rar (676.16 KB, 下载次数: 94) 2011-9-10 04:49:48 上传下载次数: 94 ! Y- L Q* y# T4 S7 z5 X1 d2：从hm导出的文件已在上面的文件包里面给出，这个inp文件将是我们后面编辑inp文件的主体文件，可以先保存为sph_input.inp,3 w- c5 m. f- t8 d# M+ ?& u9 R$ O 然后，通过PYTHON脚本进行处理，详细见贴：/forum.php?mod=viewthread&tid=1002689&highlight=sph+ i/ W6 4 W- E2 X( e$ s 需要注意的是，因为含有多个部分，所以输入命令改为：abaqus python solidtosph.py -inp - s/ v- N* E. o/ a0 v) t) U4 Y3：生成新的inp文件，打开这个inp文件会发现，里面的C3D8/C3D8R/C3D8I单元全被转换为了PC3D单元，因为我在hm中对靶体的sph部分与8 c0 D* w$ M, d0 N; i. L3 g 非sph部分以及弹体分别建立了单元集，所以这个新生成的inp文件里面会有3个PC3D单元集，将其中代表SPH区域的单元集复制过来，覆盖1 t& E% ?( S; + B7 q sph_input.inp中原为C3D8R单元集的SPH部分。4 l, G w x V* g) i2 |( O- M: b4：建立面集：首先将SPH区域基于单元集建立节点集，再基于节点集建立面，这个面用于后面与弹体的接触以及与非sph区域的衔接，相应! F1 l% B9 F/ + n c 命令行：*nset,nset=nplate,elset=sph（SPH区域单元集）8 k1 h I8 T# _+ E1 K *- m* R9 |4 w( I7 m# h7 ; a *surface,name=nodes,type=node2 PO1 O( l; 3 0 I2 D nplate3 h* v/ f$ Z# l1 L- T/ D* F 其次将非SPH区域基于单元集建立face，用来与SPH区域建立衔接。$ 6 r. b4 s6 F( o 第三，将弹体基于单元集建立face,用来定义与SPH区域的接触。. Z) 3 q, p1 Qk+ Z1 U5：建立接触，将前面建立的SPH区域与非SPH区域的面，建立面面接触，相应命令4 b3 I2 l7 1 R5 Y6 n1 h8 k *contact9 I9 v/ n( W9 A) n: E *contact inclusions$ t: Z9 R- v1 A6 E1 n nodes, project/ J3 _# J. W2 I *Contact Property Assignment5 q E2 y C8 a6 U ,frict/ M; G$ n, q9 |+ s5 N ( Q3 a6：建立SPH区域与非SPH区域的衔接，这个是靠“tie”来实现的，* J2 x% y& 3 b O* x o* c *tie,name=plates,position tolerance=0.5,adjust=no, Q7 b4 L2 I: K1 L5 Z nodes,plate2& j: 6 H/ m g- p% z2 x/ M剩下的就是建立截面属性（SPH的按默认的设置即可，也可以按照手册上的介绍修改参数），材料，建立分析步，输出选项等等，见： sph_input.rar (240.15 KB, 下载次数: 77) 2011-9-10 05:12:42 上传下载次数: 77 , i L7 Y3 1 & : e. N?$ t2 W/ |/ Z. d, k+ R&3 V d3 i1 A3 r, |9 $ / g) y3 M4 I1 _2 o6 a说明：对直接建立PC3D单元，我觉得也较简单，我查看了一下转换得到的PC3D单元，无非是原模型中SPH单元的节点号，一个PC3D单元对应一个节点号，将这个: Q# % X# k) L: B9 z c8 k4 W6 节点号写两次，就是对应的PC3D单元。; f5 R7 l- 1 x1 z: I5 r( M! j3 1 D例如：*ELEMENT,TYPE=PC3D,ELSET=sph9 X3 H$ C4 k, f1 p 186, 1869 e& B0 d$ 6 |; H0 E5 f5 b 187, 187 x z: n$ U+ d* 3 W2 d 188, 188& q% o; 0 s3 E! A+ Z& 。+ H: d# & 5 t N/ v 11084, 110847 ( U9 E7 D- w$ A 11085, 11085- ( i* t5 _y. u1 Y0 F* 1 11086, 110869 D7 e0 r9 n: E5 c9 K$ ? 11087, 11087# 6 o# O/ H- s 11088, 110888 5 h$ d( p$ 2 Y; 8 N; b/ O这是转换出来的SPH区域的单元，如下是原单元的节点：4 u7 |5 L! N3 s w3 b2 NA+ M/ U/ S186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201$ i, J/ m6 _( u8 K) 8 M5 ) I8 N 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 2173 V9 r L/ w# v$ t8 G; $ T* s 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233 J9 W% w, t! 1 a* y) ; O1 t 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249$ T; z3 KS$ v) e2 p。( ; U& D0 F9 N( j( G11052, 11053, 11054, 11055, 11056, 11057, 11058, 11059, 11060, 11061, 11062, 11063, 11064, 11065, 11066, 11067+ a5 J; R5 K+ 7 o0 E6 QT& I1 q11068, 11069, 11070, 11071, 11072, 11073, 11074, 11075, 11076, 11077, 11078, 11079, 11080, 11081, 11082, 11083- F6 V8 L7 ?/ f% M2 K11084, 11085, 11086, 11087, 110888 D7 L. V_+ d是不是都对上了？7 a6 t, K& X5 |; M知道了这个规律，手写PC3D单元应该也不难。: U6 y- u; S1 ?b; l4 d5 1 e# W2 D$ N% L$ Y) h* L2011-9-10 04:45:03 上传下载附件 (32.71 KB) + O0 Q3 v* v k2 HLamb波的问题一块450mm*450mm*1.5mm的铝板，在其上激励lamb波，激励源是半径为1mm的圆，在距离30cm处放置两个传感器（一个与激励在同一水平面，另一个偏移30度），接收波形，lamb波中心频率为200k。, g9 2 A _6 j0 U- m( |0 Q+ k! e2 K9 y) G3 G% G用ABAQUS仿真过程：% ) : W( H! 5 E0 a! O 1、建立3维板模型，尺寸如上，确定激励和传感的位置，在激励处利用拆分草图及拆分几何元素选项作圆，然后将板拆分，如图：; r( Q6 q0 O- O9 P: B8 c 2011-12-15 07:54:09 上传下载附件 (32.1 KB) / r0 L. |; Y3 + c2 k$ W9 P$ 5 j+ j M3 ? 2、指定材料属性，密度2780kg/m3,杨氏模量73.1GPA，泊松比0.33，创建截面，指派截面，装配。. V9 Z. I6 / H/ X- C+ t6 S 3、划分网格，网格大小为1mm。3 D/ 9 m2 C- q* h- L 4、创建动力、显式的分析步，时间为0.001s，最大时间增量5e-007，分别在两个传感处建立集，每个集包含一个单元的三个节点，如图：- ZW p 4 % e 2011-12-15 08:02:51 上传下载附件 (3.45 KB) ) H: n9 X3 / H 0 V j$ V, n% _4 g R; t4 iB- y8 l* B0 S& 5、建立场输出和历程输出，场输出为整个模型，输出变量为位移U，采样时间为5e-007；历程输出为两个集，输出变量为位移U1、U2、U3，采样时间为5e-007。N7 J* b9 i1 Q 6、创建作业，选择写入输入文件选项，已得到inp文件，方便为载荷选点。6 p% R L1 BR7 m 7、导入inp文件，创建载荷，幅值为预先由MATLAB得到的数据，载荷为集中力，加载方法如图：& n$ Q. ?, 9 9 T 2011-12-15 08:08:53 上传下载附件 (39.15 KB) * g; h+ K6 - r: U8 i# P3 P0 m/ K$ y0 M, N U, T0 k; p 8、创建作业，提交，得结果5 z9 i8 d# D# F5 y3 H+ |5 5 o) g- D2 G) k n8 P& y5 Q9 l遇到问题：5 Z& a! M3 T3 L8 u% r& # l8 s* J. f& i 1、按照LAMB波在铝板中的传播理论，在100us附近，X轴上的集应该可以观察到S0模式的反射，而仿真结果显式，两个相邻节点处的U1只有一个有反射，另一个却没有，两个节点相距仅1mm而已，如下图：; z0 x0 |3 Y4 o% c/ e5 pu 2011-12-15 08:15:18 上传下载附件 (85.29 KB) % R* o: L) w. m$ % c0 ?- X m# G& : K5 U; B+ d4 f$ j6 - E& Z& s& v( e 2、因为两个传感位于同一波阵面，其U3波形应该一致，而实际结果相差较大，如下图：, N! 0 d# on! e 2011-12-15 08:18:25 上传下载附件 (89.7 KB) 波形2 8 & 7 u( Q) l 2 I; R6 d8 u* n9 Z% G8 p f3 Q1 b5 w5 c$ v, Y7 ! W+ I; G) d, k我实际想观察的量为应变，但参考相关文献，输出的量一般都为位移，如何得到应变量，由位移差直接相加可以吗？对于以上过程，麻烦各位高手帮忙看下，问题出在哪里？或者是我理解不到位，关于这个问题该怎么由ABAQUS实现，请指点！万分感谢！$ j$ g2 TB% |8 Q$ e- dlamb不是波的名字, 是一类问题的总称, lambs problem, 其分析的关键有下面两点/ x. B) Y) o# x& 4 A1.找到正确可用的impuls点荷载所对应的加载函数0 s1 x7 c2 I2.确认波的传播速度, 按照其传播速度确定边界是否采用无限元或是visco elastic 边界振动 复合材料板中超声波的模拟 首先要先了解超声波，超声波和振动是属于同一个类型，只不过一个高频一个低频。超声波在传播的时候会出现频散的现象。这样就造成了一个波包不同的频率具有不同的速度，波包在传递一段时间以后就会被拉长。按照质点振动的情况我们把超声波的模态分成对称模态S和反对称模态A。例子中给出的是模拟对称模态也就是S的情况。对于固定的材料无论是各项同性的还是像复合材料这种正交各项异性的材料，只要材料的参数给定就可以计算出每个频率所对应的速度。S/ i* g9 - y 我们现在来讨论abaqus中的模拟，在这里我讨论三个问题，首先对于复合材料的设定，这里我就不多说了，论坛中很多这类的资料说的都非常的详细。我在例子中为了简化起见用的是单层的正交材料。第二个问题，单元的尺寸，对于超声波的模拟单元的尺寸非常重要，一般单元的尺寸要比波长小一个数量级，比方说我例子中的波长大概是20mm左右，那么我选用的单元尺寸为2mm。以前曾经实验过，大概小于波长的1/5都还算是准确的，当然最好是在1/10左右。第三个问题就是时间，一般选用的时间跟频率有关系，可以选用1/(20*fmax)，关于这个问题还请大家多多去探究。最后的两个问题是控制准确性的非常重要的指标。大家可以多去实验。- P* ( O! G1 Q8 I! w红色和黑体并没有特别的意义，只不过同一种字体读起来比较单调。o% q4 C# k+ r( ! P 2011-8-17 14:51:38 上传下载附件 (51.89 KB) 2 s! p% D% z% $ n5 t 超声波在变截面棒中传播模拟的一个问题4 R0 5 7 u# j m) ! C O4 Z最近在做一个超声波在棒中的传播（inp文件在附件中）。根据理论计算，波在变截面杆中传播时（如图），假设大端直径为D1，面积为S1，小端直径为D2,面积为S2，那么波在小端和大端中引起的质点振动位移，振动速度之比V1/V2=位移s1/位移s2=面积S2/S1. 但是我模拟出来的结果显示，在小端振动有所增强，但是不满足上诉关系。请各位大侠指点，是我的模拟设置不正确还是其他问题呢？F