当前位置:首页 >> 机械/仪表 >>

ls-dynasph小水滴跌落


LS-dyna 小水滴跌落台阶 SPH sph 小水滴是在 lspp 中建立的,这点还是比较简单,其中 2.4 版本中第七张卡片有一个 sphgen 生成器,其中要填写的参数包括质心点的三方向坐标,密度,三方向的粒子个数等 等,基本大家都可以看懂,做完之后就导入到其他前处理中完成刚体建模,定义材料状态方 程,其中材料是使用 mat-null,配合状态方程 eos-gruneisen 完成水滴的材料,对于材料 参数,可以在论坛和 manual 上找到答案,在此不再详解,定义结构属性时,sph 粒子要定 义为 sph 粒子属性,接触采用通用的 a-n-t-s 的形式即可,然后就是控制参数,一般采用默 认值定义粘度沙漏时间步长能量等,control sph 内部有几个参数需要注意,此处只讲解流 程,不对具体定义进行阐述。至此模型建立完成,上图。

$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ LS-DYNA(970) DECK WAS WRITTEN BY: eta/VPG VERSION 3.2 $ $ ENGINEER: $ PROJECT: $ UNITS: M, KG, SEC, N $ DATE: Dec 19, 2011 at 14:40:19 $ $ NOTES: $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ VIEWING INFORMATION $ -.252462E+000.207475E+020.388527E+010.185637E+02 $ 0.996363E+00-.824454E-010.215457E-01 $ 0.592740E-020.319285E+000.947640E+00 $ -.850088E-01-.944065E+000.318613E+00 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *TITLE LS-DYNA keyword deck by LS-PRE $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8

*CONTROL_BULK_VISCOSITY $ q1 q2 type 1.500 0.06000 1 *CONTROL_CONTACT $ slsfac rwpnal islchk shlthk penopt thkchg orien enmass 0.10000 0.0 1 0 1 0 1 0 $ usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj 0 0 10 0 4.000 0 0 0 $ sfric dfric edc vfc th th_sf pen_sf 0.0 0.0 0.0 0.0 0.0 0.0 0.0 $ ignore frceng skiprwg outseg spotstp spotdel 0 0 0 0 0 0 *CONTROL_ENERGY $ hgen rwen slnten rylen 2 2 2 2 *CONTROL_HOURGLASS $ ihq qh 1 0.10000 *CONTROL_SPH $ ncbs boxid dt idim memory form start maxv 1 0 .1000E+21 3 500 0 0.0 .1000E+16 $ cont deriv 0 0 *CONTROL_TERMINATION $ endtim endcyc dtmin endeng endmas 0.05000 0 0.0 0.0 0.0 *CONTROL_TIMESTEP $ dtinit tssfac isdo tslimt dt2ms lctm erode ms1st 0.0 0.60000 0 0.0 0.0 0 0 0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *DATABASE_GLSTAT $ dt binary .5000E-3 1 *DATABASE_MATSUM $ dt binary .5000E-3 0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *DATABASE_BINARY_D3PLOT $ dt/cycl lcdt beam npltc .7000E-3 0 0 0 $ ioopt 0

$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *PART SphNode $ pid secid mid eosid hgid grav 2 1 1 1 0 0 *PART PART_1 1 2 2 0 0 0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *SECTION_SPH_TITLE $ secid cslh 1 1.200 *SECTION_SHELL_TITLE $ hmin 0.20000 hmax 2.000 sphini death 0.0 .1000E+21

adpopt 0

tmid 0

0

0

start 0.0

secid elform shrf nip propt qr/irid icomp 2 2 1.000 2.000 0.0 0.0 $ t1 t2 t3 t4 nloc marea 1.000 1.000 1.000 1.000 0 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *MAT_NULL_TITLE $ pr mid ro 0.0 pc 0.00100 mu 0.0 terod 0.0 cerod

setyp 0 1

ym 0.0

1 1000.000 *MAT_RIGID_TITLE $

0.0

mid ro e pr n couple 2 7830.000 .2070E+12 0.28000 0.0 0.0 $ cmo con1 con2 1.000 7 7 $ a1 a2 a3 v1 v2 v3 0.0 0.0 0.0 0.0 0.0 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *EOS_GRUNEISEN_TITLE $ eosid c s1 s2 s3 gama0 1 1484.000 1.979 0.0 0.0 0.11000 $ v0 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *SET_NODE_LIST_TITLE

m 0.0

alias 0.0

a 0.0 0.0

e0

NODE SET_1 $ sid da1 da2 da3 da4 1 $ nid1 nid2 nid3 nid4 nid5 nid6 nid7 nid8 296 297 298 299 300 301 302 *INITIAL_VELOCITY $ nsid nsidex boxid irigid 1 0 0 0 $ vx vy vz vxr vyr vzr 0.0 0.0 -200.000 0.0 0.0 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *CONTACT_AUTOMATIC_NODES_TO_SURFACE_ID $ cid 1 $ ssid msid sstyp mstyp sboxid mboxid spr mpr 1 1 4 3 0 0 0 0 $ fs fd dc vc vdc penchk bt dt 0.20000 0.20000 0.0 0.0 0.0 0 0.0 .1000E+21 $ sfs sfm sst mst sfst sfmt fsf vsf 1.000 1.000 0.0 0.0 1.000 1.000 1.000 1.000 $ soft sofscl lcidab maxpar sbopt depth bsort frcfrq 0 0.10000 0 1.025 0.0 2 10 1 $ penmax thkopt shlthk snlog isym i2d3d sldthk sldstf 0.0 0 0 0 0 0 0.0 0.0 $ igap ignore 1 0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *LOAD_BODY_Z $ lcid sf lciddr xc yc zc 1 1.000 0 0.0 0.0 0.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *DEFINE_CURVE_TITLE $ sfa sfo offa offo dattyp 1 0 1.000 1.000 0.0 0.0 $ abscissa (time) ordinate (value) 0.000000E+00 0.000000E+00 0.100000E-02 0.980000E+01 0.100000E+02 0.980000E+01 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *NODE $ nid x y z tc 1 1.000000000E+01 0.000000000E+00 1.000000000E+01 lcid sidr

0

rc

2 1.000000000E+01 1.000000000E+00 1.000000000E+01 33845 2 1.95312E+00 33846 2 1.95312E+00 33847 2 1.95312E+00 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *END


赞助商链接
相关文章:
DYNA基础,,实例
LS-DYNA 的最新版本是 2004 年 8 月推出的 970 ...4 9.SPH 算法 SPH(Smoothed Particle Hydrodynamics)...等 电子领域 跌落分析 包装设计 热分析 电子封装 ...
LS-dyna-常见问题汇总2.0
LS-dyna-常见问题汇总2.0_调查/报告_表格/模板_...36 请问 LS-Pre/Post 有没有办法改变 SPH 颗粒的...跌落物体一个初速度,把 它放在离地面一个非常小的...
LS-DYNA简介
LS-DYNA简介_IT/计算机_专业资料。LS-DYNA 是世界...SPH 单元 4.接触方式(50 多种) 柔体对柔体接触...等 电子领域 跌落分析 2 3 4 5 6 7 8 包装...
LS-DYNA问答总结
LS-DYNA(ANSYS)中文培训教... 40页 免费如要投诉违规内容,请到百度文库投诉中心;如要提出功能问题或意见建议,请点击此处进行反馈。 ...
SPH粒子教程
SPH粒子教程_交通运输_工程科技_专业资料。Ls-dyna 中无网格法的一点心得 前处理: 一. 采用 ls-prepost 在 ls-prepost 中第 7 页有 sphgen 的功能,利用...
LS-DYNA 简介
LS-DYNA 简介 LS-DYNA 是世界上最著名的通用显式...SPH 单元 4.接触方式(50 多种) ¨¨¨ ¨ 柔体...等 电子领域 ¨ 跌落分析 ¨ 包装设计 ¨ 热分析...
LS-DYNA精编教材
LS-DYNA精编教材_工学_高等教育_教育专区。LS-DYNA...SPH 算法适用于超高速碰撞、靶板贯穿等过程的计算...等 电子领域 跌落分析 包装设计 热分析 电子封装 ...
Ls-Dyna keywords
Ls-Dyna keywords_机械/仪表_工程科技_专业资料。Ls-Dyna部分重要keywords ...For SPH elements, particles will bedeactivated when TSMIN is reached For ...
LS-DYNA常见问题简化中文版
LS-DYNA常见问题简化中文版_机械/仪表_工程科技_专业...请问 LS-Pre/Post 有没有办法改变 SPH 颗粒的大小...话说回来, 「破裂」现像还有许多判别的方式,利用小...
LS-DYNA
解答: LS-DYNA 不只能分析实体元素(Tetra、Hexa)、薄/厚壳元素(Quad.、Tri),同时还有 梁元素、质点、 Spring/Damper、spotweld 等元素型态,另外还有 SPH ...
更多相关标签: