汽车碰撞仿真LS-DYNA控制卡片关键字.
汽车碰撞分析LS DYNA控制卡片设置共26页

谢谢你的阅读
❖ 知识就是财富 ❖ 丰富你的人生
71、既然我已经踏上这条道路,那么,任何东西都不应妨碍我沿着这条路走下去。——康德 72、家庭成为快乐的种子在外也不致成为障碍物但在旅行之际却是夜间的伴侣。——西塞罗 73、坚持意志伟大的事业需要始终不渝的精神。——伏尔泰 74、路漫漫其修道远,吾将上下而求索。——屈原 75、内外相应,言行相称。——韩非
汽车碰撞分析LS DYNA控制卡片设置
36、如果我们国家的法律中只有某种 神灵, 而不是 殚精竭 虑将神 灵揉进 宪法, 总体上 来说, 法律就 会更好 。—— 马克·吐 温 37、纲纪废弃之日,便是暴政兴起之 时。— —威·皮 物特
38、若是没有公众舆论的支持,法律 是丝毫 没有力 量的。 ——菲 力普斯 39、一个判例造出另一个判例,它们 迅速累 聚,进 而变成 法律。 ——朱 尼厄斯
LS-Dyna碰撞分析调试指南

LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有:a、质量增加百分比小于5%;b、总沙漏能小于5%;c、滑移界面能;d、检查各部件之间的连接、接触关系是否定义正确,检查模型的完整性;e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查:质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元,我们通常采取增加单元材料密度的方法来增大其时间步长,以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步,在Hyper view中调出glstat文件并检查mass scaling项(质量增加应该小于5%);调出matsum文件并检查各部件的质量增加情况,对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置(质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误,导致该单元的时间步长低于控制卡片中设置的最小时间步长,从而引起质量缩放)。
2、全过程检查。
调整模型使其符合初步检查的标准,计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中,最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查:沙漏能的出现是因为模型中采用了缩减积分引起的,我们常用的B-T单元采用的是面内单点积分,这种算法会引起沙漏效应(零能模式)。
具体介绍参见附录二。
检查:在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项,整个计算过程中沙漏能应小于总能量的5% 。
三、滑移界面能sliding interface energy的检查:滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能;未能检测到的穿透(undetected penetrations)常常会引起大的负值的滑移截面能。
Ls-dyna 高级教程 关键字讲解-_Lectures_1

Solver - Solution processing
DYNA3D developed at Lawrence Livermore National Laboratory by John Hallquist Low velocity impact of heavy, solid structures, military applications DYNA3D ported on Cray-1 Improved sliding interface Order of magnitude faster New material models - Explosive-structure, Soil-structure Impacts of penetration projectiles
The George Washington University
Pre-stress and Post-stress Interactive graphics Preprocessor - LS-PrePost Third party interfaces Postprocessor, LS-PrePost Other rigid body program coupling CAD d data t i interface t f
The George Washington University
Department of Civil and Environmental Engineering
Keyword Format Input File
*KEYWORD *TITLE SAMPLE INPUT FILE *CONTROL_TERMINATION 0.1000000 0 0.0000000 0 0.0000000 *DATABASE_BINARY_D3PLOT 1.00000-3 0 *DATABASE_BINARY_D3THDT 1.00000-3 *MAT_ELASTIC 1 7.89000-9 2.00000+5 0.3000000 *SECTION_SOLID 1 0 *SECTION_SHELL 1 2 1.0000000 1.0000000 1.0000000 1.0000000 0.0000000 *PART PART NAME 1 1 1 1 0 0 0 0 0 *NODE 1 0.000000000E+00 0.000000000E+00 0.000000000E+00 2 7.000000000E+00 0.000000000E+00 0.000000000E+00 3 0.000000000E+00 7.000000000E+00 0.000000000E+00 4 7.000000000E+00 7.000000000E+00 0.000000000E+00 5 0.000000000E+00 0.000000000E+00 7.000000000E+00 6 7.000000000E+00 0.000000000E+00 7.000000000E+00 7 0.000000000E+00 7.000000000E+00 7.000000000E+00 8 7.000000000E+00 7.000000000E+00 7.000000000E+00 *ELEMENT_SOLID 1 1 1 2 4 3 5 6 8 7 *PART PART NAME 2 2 2 2 0 0 0 0 0 *ELEMENT SHELL *ELEMENT_SHELL 1 2 1 2 4 3 *END
LS-Dyna碰撞分析调试指南

LS-DYNA 碰撞分析调试LS-DYNA碰撞计算模型的主要检查、调试项目有:a、质量增加百分比小于5%;b、总沙漏能小于5%;c、滑移界面能;d、检查各部件之间的连接、接触关系是否定义正确,检查模型的完整性;e、检查数值输出的稳定性。
一、质量缩放Mass scale的检查:质量缩放——对于时间步长小于控制卡片中设置的最小时间步长的单元,我们通常采取增加单元材料密度的方法来增大其时间步长,以减短模型的计算时间。
关于LS-DYNA中单元时间步长的计算方法请参见附录一。
1、初步检查。
让模型在dyna中运行2个时间步,在Hyper view中调出glstat 文件并检查mass scaling项(质量增加应该小于5%);调出matsum文件并检查各部件的质量增加情况,对于质量增加过大以及有快速增长趋势的部件应检查此部件的网格质量和材料参数设置(质量增加一般是由于单元的特征长度太小或者是材料参数E、ρ设置错误,导致该单元的时间步长低于控制卡片中设置的最小时间步长,从而引起质量缩放)。
2、全过程检查。
调整模型使其符合初步检查的标准,计算模型至其正常结束。
再按[初步检查]的要求检查调试整个模型直至达到要求。
一个计算收敛的模型在其整个计算过程中,最大质量缩放应小于总质量的5% 。
二、沙漏能Hourglass energy的检查:沙漏能的出现是因为模型中采用了缩减积分引起的,我们常用的B-T单元采用的是面内单点积分,这种算法会引起沙漏效应(零能模式)。
具体介绍参见附录二。
检查:在dyna中计算模型至其正常结束。
在Hyper view中调出glstat文件并检查energy的total energy 、Hourglass energy两项,整个计算过程中沙漏能应小三、滑移界面能sliding interface energy的检查:滑移界面能是由摩擦和阻尼所引起的。
剧烈的滑动摩擦会引起大的正值的滑移界面能;未能检测到的穿透(undetected penetrations)常常会引起大的负值的滑移截面能。
汽车碰撞分析LS_DYNA控制卡片设置

控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时,要检查所有的单元。出于稳定性原因,用0.9(缺省)来 减小时间步:Δt = 0.9 l/c ,特征长度l,和波的传播速度c,都与单元的类型有关。 DTINIT:初始时间步长,如为0.0,由DYNA自行决定初始步长; TSSFAC:时间步长缩放系数,用于确定新的时间步长。默认为0.9,当计算不稳定时,可以减小该值,但同时 增加计算时间; ISDO:计算4节点壳单元时间步长的(不同的值对应特征长度的不同算法,推荐使用2,因为此选项可以获得 最大的时间步长,但有三角形单元存在时会导致计算不稳定); TSLIMT:壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上;只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元,不建议使用该选项,因为使用DT2MS选项更好。 DT2MS:因质量缩放计算得到的时间步长。当设置为一个负值时,初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时,推荐用这种方法。用这种方法时 质量增加是有限的,过多的增加质量会导致计算终止。当设置为正值时,初始时间步长不会小于DT2MS。单 元质量会增件或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止,但更难以作出合理的解释。默认为0.0,不进行质量缩放; LCTM:限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系(可选择) ; ERODE:当计算时间步长小于TSMIN(最小时间步长)时体单元和t-shell被自动删除。
LS-DYNA 关键字用户手册说明书

Copyright © 2001, Dr Brian Gladman < brg@>, Worcester, UK. All rights reserved.
The free distribution and use of this software in both source and binary form is allowed (with or without changes) provided that:
LS-DYNA®, LS-OPT® and LS-PrePost® are registered trademarks of Livermore Software Technology, an Ansys company in the United States. All other trademarks, product names and brand names belong to their respective owners.
LST reserves the right to modify the material contained within Байду номын сангаасhis manual without prior notice.
The information and examples included herein are for illustrative purposes only and are not intended to be exhaustive or all-inclusive. LSTC assumes no liability or responsibility whatsoever for any direct or indirect damages or inaccuracies of any type or nature that could be deemed to have resulted from the use of this manual.
汽车碰撞精确分析LSDYNA控制卡片设置

THKCHG:在单面接触时考虑壳厚度的改变(默认时不考虑)。
ORIEN:在初始化时可选择性的对接触面部分自动再定位。
控制卡片参数说明
ENMASS:接触单元被腐蚀的质量处理。0-节点被移除,1-体单元节点被保留,2-体单元壳单元节点被保留。 USRSTR:每个接触面分配的存储空间,针对用户提供的接触控制子程序。 USRFRC:每个接触面分配的存储空间,针对用户提供的接触摩擦子程序。 NSBCS:接触搜寻的循环数(使用三维Bucket分类搜索),推荐使用默认项。 INTERM:间歇搜寻主面和从面接触次数。 XPENE:接触面穿透检查最大乘数,默认4.0。 SSTHK:在单面接触中是否使用真实壳单元厚度,默认0,不使用真实厚度。 ECDT:时间步长内忽略腐蚀接触。
16. DATABASE_BINARY_RUNRSF 设置如下:
控制卡片参数设置
17. DATABASE_BINARY_RUNRSF 设置如下:
控制卡片参数说明
*CONTROL_TERMINATION
$ ENDTIM ENDCYC DTMIN ENDENG ENDMAS
150
0
0.0
0.0
0.0
SLSFAC:滑动接触惩罚系数 ,默认为0.1。当发现穿透量过大时,可以调整该参数;
RWPNAL: 刚体作用于固定刚性墙时,刚性墙罚函数因子系数,为0.0时,不考虑刚体与刚性墙的作用,>0时, 刚体作用于固定的刚性墙,建议选择1.0;
ISLCHK:接触面初始穿透检查,为0或1(默认)时,不检查。为2时,检查。
后面将逐一介绍碰撞分析中经常用到的控制卡片,并对每个卡片的作 用进行说明。
控制卡片使用规则
卡片相应的使用规则如下:
碰撞模拟中的控制卡片设置及意义.

1.指定输出文件【ABSTAT】——气囊统计表。
输出体积、压强、内能、气体质量流入率、气体质量流出率、质量、温度、密度。
【AVSFLT】——A VS数据【BNDOUT】——边界环境的力和能量。
输出三个方向的力。
【DEFGEO】——变形的几何体的文件【DEFORC】——离散单元。
输出三个方向的力。
【ELOUT】——单元数据。
(见DATABASE_HISTORY_OPTION)梁单元平面应力块平面应变轴向合力xx,yy,zz 应力xx,yy,zz 应力xx,yy,zz 应变S方向剪切合力xy,yz,zx 应力xy,yz,zx 应力xy,yz,zx 应变T方向剪切合力塑性应变有效应力下表面应变S方向合力矩屈服函数上表面应变T方向合力矩扭力合力【GCEOUT】——几何接触实体。
包含三个方向力和力矩。
【GLSTAT】——总体数据。
【JNTFORC】——运动副力文件【NATSUM】——材料能量。
GLSTAT JNTFORC MATSUM动能x,y,z三方向的力动能内能x,y,z三方向的力矩内能总能量沙漏能比率x,y,z三方向的动量刚性墙能量x,y,z三方向的刚体速度弹簧和阻尼能量总动能沙漏能总内能阻尼能总沙漏能滑移面能量外功x,y,z三方向速度时间步单元ID号控制的时间步【MOVIE】——【MPGS】——【NCFORC】——接触面节点力【NODFOR】——节点力组【NODOUT】——节点数据NCFORC NODOUT NODFORx方向力位移x,y,z三方向力y方向力速度z方向力加速度转动量角速度角加速度【RBDOUT】——刚体数据【RCFORC】——接触面合成力【RWFORC】——刚性墙所受的力RBDOUT RCFORC RWFORC三方向合位移三方向合力法向力三方向合速度三方向合力三方向合加速度【SBTOUT】——安全带输出文件【SECFORC】——横截面通过的力(见DA TABASE_CROSS_SECTION_OPTION)【SLEOUT】——滑移面的能量【SPCFORC】——单点约束的反作用力【SPHOUT】——SPH数据(见DATABASE_HISTORY_OPTION)【SSSTAT】——子系统数据【SWFORC】——节点约束反力(焊点和铆钉)SECFORC SLEOUT SPCFORC SWFORC x,y,z三方向力Slave能量x,y,z三方向力轴向力x,y,z三方向力矩Master能量x,y,z三方向力矩剪切力x,y,z三方向中心面积合力【TPRINT】——结构对的热量输出【TRHIST】——追踪质点时间历程信息2.时间步控制【DTNIT】——初始时间步长。
第一章 LS-DYNA关键字简介讲解

•1986 梁、壳、刚体
单一面接触
多平台支持
•1988
支持汽车行业应用
LS-DYNA3D
•1989
LSTC Full-Time 更多的选项和特性
ASME WAM’89 - San Francisco
-交通运输系统中的防撞性和乘
客安全
登录协同仿真时代
•1993 关键字格式
新的自动单一面接触
第一届LS-DYNA3D用户 年会
*COMPONENT假人组件 刚体假人. Integrated implicitly - separately from FE model. 由接触面控制相互作用. GEBOD 男性,女性,小孩 Hybrid III 5th, 50th and 95th 百分位成年人
登录协同仿真时代
*CONSTRAINED 结构体内和结构体之间
后处理: • »LS-POST, LS-TAURUS • »GLview, ANSYS,
HyperMesh, Fem-B
对于同一个物理系统总是有 多种方法构造有限元分析模 型 所有的模型不能完全重现物 理系统
登录协同仿真时代
输出文件
二进制文件 » d3dump01 » d3plot » d3plot01 » d3plot02 » …. » d3thdt » D3thdt01
登录协同仿真时代
*INCLUDE 把文件分成子文件 把子文件分成子-子文件 专门应用:冲压 part导入和转换(缩放、转动和平移)
登录协同仿真时代
*INITIAL 点火和动量 初始应力/应变 初始温度 初始汽车运动关系 初始速度
*INTERFACE • 研究大型结构中小部件的详细响应 •隔离关键部件 记录界面/边界的运动: »面 » 节点线 » 节点点
汽车碰撞分析LS_DYNA控制卡片设置

控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时,要检查所有的单元。出于稳定性原因,用0.9(缺省)来 减小时间步:Δt = 0.9 l/c ,特征长度l,和波的传播速度c,都与单元的类型有关。 DTINIT:初始时间步长,如为0.0,由DYNA自行决定初始步长; TSSFAC:时间步长缩放系数,用于确定新的时间步长。默认为0.9,当计算不稳定时,可以减小该值,但同时 增加计算时间; ISDO:计算4节点壳单元时间步长的(不同的值对应特征长度的不同算法,推荐使用2,因为此选项可以获得 最大的时间步长,但有三角形单元存在时会导致计算不稳定); TSLIMT:壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上;只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元,不建议使用该选项,因为使用DT2MS选项更好。 DT2MS:因质量缩放计算得到的时间步长。当设置为一个负值时,初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时,推荐用这种方法。用这种方法时 质量增加是有限的,过多的增加质量会导致计算终止。当设置为正值时,初始时间步长不会小于DT2MS。单 元质量会增件或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止,但更难以作出合理的解释。默认为0.0,不进行质量缩放; LCTM:限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系(可选择) ; ERODE:当计算时间步长小于TSMIN(最小时间步长)时体单元和t-shell被自动删除。
ls-dyna关键字命令介绍

LIVERMORE SOFTWARE TECHNOLOGY CORPORATION
*DATABASE_option 所谓的 ‘ASCII’ 输出文件
• 指定的时间历程输出 • 后续的幻灯片给出输出内容 • 输出格式的选择 (任何一种都可以用 LS-Prepost 处
理) • 真实的 ASCII ( 文本 ) 输出 • 压缩二进制输出 (binout)
在输出数据中选择另外的名字如在命令行中使用gfilename输出请求可能随分析的不同而不同通常50到100个状态的输出是足够的了对大模型来说使用d3thdt数据输出得到有选择的节点或单元的时间历程数据livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseextentbinarydatabaseextentbinary对写入到二进制输出文件中的数据进行控制默认输出3个积分点midbottomtop的数据默认情况下lsdyna计算的弹性应变是不输出的各向异性正交异性的材料其应力可以在全局坐标系中输出也可以在材料坐标系中输出matorthotropicelastic除外livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseextentbinarycontd控制每一个输出状态创建一个单独的d3plot文件选项ieverp梁的积分点输出数目beamip如果打开该选项可以对壳材料的沙漏能密度进行云图显livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationdatabaseoption所谓的ascii输出文件文本输出压缩二进制输出binoutlivermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationascii输出文件airbagstatisticsdatabaseabstatinternalenergymassflowoutputmassflowrateboundarynodalforcesdatabasebndout对应于节点和刚体边界条件的节点力livermoresoftwaretechnologycorporationlivermoresoftwaretechnologycorporationascii输出文件discreteelementdatadatabasedeforc离散单元的力和力矩
第一章 LS-DYNA关键字简介讲解

登录协同仿真时代
谢 谢!
Welcome to
ANSYS
China@
登录协同仿真时代
•1995 热分析
建立培训 LAB
交通运输系统中的防撞性 和乘 客安全
第三届LS-DYNA3D用户 年会
•1996 PC 版本(KBS2) 整合ANSYS •1997 •LS-DYNA(940版) •Eulerian/Lagrangian耦合 •自适应网格 •1999 •LS-POST •» LS-OPT
DYNA高级应用培训
1、关键字简介
发展历史
•1976 DYNA3D:Lawrence Livermore
国家实验室 John O. Hallquist 固体结构低速冲击 •1979 Cray-1 上运行 改进的滑动界面 Order of Magnitude faster •1981 新的材料模式 -炸药结构、土壤结构 射弹冲击
*CONTROL •终止:时间,能量或质量的变化 •时间步 •全局参数 » 沙漏类型 » 接触罚函数缩放因子 » 壳单元公式 » 阻尼 • 自适应网格, ALE, CFD,隐式
登录协同仿真时代
*DATABASE 控制输出 »ASCII 文件 » Binary文件
*DEFINE 与其他方式联动的输入: 曲线 限制几何范围的Boxes 局部坐标系 矢量 转换 (缩放,旋转,平移) » 必须在 *INCLUDE_TRANSFORM前 定义
*COMPONENT假人组件 刚体假人. Integrated implicitly - separately from FE model. 由接触面控制相互作用. GEBOD 男性,女性,小孩 Hybrid III 5th, 50th and 95th 百分位成年人
汽车碰撞仿真LS-DYNA控制卡片关键字.答案

控制卡片使用规则
卡片相应的使用规则如下:
•������ 大部分的命令是由下划线分开的字符串, 如*control_hourglass
•������ 字符可以是大写或小写 •������ 在输入文件中,命令的顺序是不重要的(除了*keyword 和*define_table)
•������ 关键字命令必须左对齐,以*号开始
8. CONTROL_OUTPUT 用于设置输出参数。设置如下,按红色线框内设置,其余默认值。
控制卡片参数设置
8. CONTROL_PARALLEL 并行计算控制,最大可以使用4个CPU。
9. CONTROL_SHELL 壳单元控制
10. CONTROL_TERMINATION 计算终止控制卡片,控制计算终止时间,设置如下:
ISTUPD:单元厚度改变选项。该选项对所有壳单元变形有影响。
THEORY:壳单元使用的理论。(默认的是Belytschko-Tsay,面内单点积分,计算速度很快,采用Co-rotaional 应力更新,单元坐标系统置于单元中心,基于平面单元假定,建议在大多数分析中使用) BWC:针对Belytschko-Tsay单元的翘曲刚度。 MITER:平面应力塑性选项,默认为1。(运用于材料3,18,19和24)。 PROJ:在Belytschko-Tsay和Belytschko-Wong-Chiang单元中翘曲刚度投影方法。这个方法主要运用于显示分 析,如果是隐式分析,那此项无效 。 默认为0.
ORIEN:在初始化时可选择性的对接触面部分自动再定位。
控制卡片参数说明
ENMASS:接触单元被腐蚀的质量处理。0-节点被移除,1-体单元节点被保留,2-体单元壳单元节点被保留。 USRSTR:每个接触面分配的存储空间,针对用户提供的接触控制子程序。 USRFRC:每个接触面分配的存储空间,针对用户提供的接触摩擦子程序。
利用LS-DYNA进行成形仿真的输入控制参数

利用LS-DYNA进行成形仿真的输入控制参数Translated by SunnyWinterLS-DYNA已广泛用于汽车碰撞分析。
默认的输入参数一般能给出有效,精确的碰撞模拟结果。
但是,这些默认值对于成形仿真分析并不一定理想。
下面是一个标准的金属成形过程。
为及时参考,推荐输入参数用黑体字标识,并包含在盒状关键字输入框中。
模型明确要求的数据,如终止时间等参数,输入黑体的0值。
一般问题设定在显式成形仿真中,利用质量比例缩放和(或者)人为的高工具速度,运行时间可以大大缩减。
这两种方法都会引入人为的动力学影响,因此必须将其减小到在工程意义上合理的水平。
一个单独的描述人为动力影响的参数是:工具每运动1毫米所采用的显式时间步进值(或周期)数目。
当成形过程允许大的无限的板料运动,比如冲击成形,需要更多的毫米周期数。
当板料被压边圈和冲模支撑较强的约束住时,较少的毫米周期数是必要的。
对大多数的仿真来说,100到1000之间的毫米周期数能产生合理的结果。
如果可能,或者有必要重复一个仿真,可利用两个不同的毫米周期值并比较分析结果去估计其对人为动力学影响的敏感性。
推荐选择的一个最大工具速度是2.0mm/ms,起始和结束速度为0。
可以使用简单的梯形速度轮廓(如图1)。
利用大的时间缩放步参数dt2ms获得要求的毫米周期数,可参考下面的公式:时间步大小=1.0/(最大工具速度*毫米周期数)工具速度,时间步大小和结束时间必须在协调的参照系中选择。
如果所有的工具运动给定,可用下面的步骤设置模拟参数:已知:工具全部行程(mm):D最大工具速度(mm/ms):2.0速度轮廓: 2.0毫秒上升和2.0毫秒下降的梯形(如图1)选择:毫米周期数:ncpm计算:结束时间(ms): T=2.0+D/2速度数据点:(0.0,0.0)(2.0,2.0)(T-2.0,2.0) (T,0.0)时间步大小(ms): dt2ms=1/(2*ncpm)上面的运算提供速度轮廓数据点用于下面的工具运动部分。
LS-DYNA关键字手册

LS-DYNA关键字手册各种前处理软件得到的k 文件往往不能满足使用要求,或者存在一些错误,这时就要自己修改、添加关键字面是自己总结的几条:0. k 文件格式分为标准格式和自由格式(数据之间用逗号隔开)两种,在一个k 文件中,两种方式可以并存,但是在一个数据卡中,只能选择一种方式1. 如果选择标准格式:k文件中除了节点(node)和单元(element )关键字外,通常每一行总共占80个字符长度,每个数据占10个字符长度,修改时千万不要超越这10 个字符长度的位置,也不要跑到别的数据的10 个字符位置2. 如果关键字手册里的card 介绍中没有提到optional ,那么每一行card 都不能省略,哪怕它们都是03. 为了方便查看10个字符长度,可以用ultraedit 软件4. 每一个关键字必须以* 开头,并且必须顶格写5. 在k 文件中$后面的是注释,求解时不考虑6. 为了查找和发现具体是那一关键字出错,可以使用lspost 打开k 文件,然后选择view 选项查看其实原理和定义材料曲线类似可以分三步:1,先将你的加速度曲线离散成一系列的点,找出这些点的横(时间) 纵(加速度)坐标值定义两个数组,一个用来存放横坐标值;另一个存放纵坐标值如:*dim,time,,5time(1)=0,.025,.05,.075,.1*dim,ACL,,5ACL(1)=0,1,5,8,9 根据你的曲线复杂程度,可适量增减节点2,定义和这条曲线相关的ID 号,利用上面确定好的数组进行定义如:edcurve,add,1,time,ACL这里面的 1 是id 号,可根据自己需要选择合适的号码3,将你定义好的加速度曲线施加到需要的物体上如:edload,add,fy,,comp,,,,1,1.0这里的comp 是你要施加作用的物体的组件名注意:应用此命令时,不要再加入别的数组了说完全不同也是不太正确的,ansys自己已经建立了对应的Isdyna如下:ANS YS Comma nd Correspo nding LS-D YNA Comma nd EDADAPT *PARTEDALE *CONTROL_ALEEDASMP *SET_PART_LISTEDBOUND *BOUNDARY_SLIDING_PLANE*BOUNDARY_C YCLICEDBVIS *CONTROL_BULK_VISCOSITYEDBX *DEFINE_BOXEDCADAPT *CONTROL_ADAPTIVEEDCGEN *CONTACTEDCMORE *CONTACTEDCNSTR *CONSTRAINED_EXTRA_NODES_SET*CONSTRAINED_NODAL_RIGID_BOD Y*CONSTRAINED_SHELL_TO_SOLID*CONSTRAINED_RIVETEDCONTACT *CONTROL_CONTACTEDCPU *CONTROL_CPUEDCRB *CONSTRAINED_RIGID_BODIESEDCSC *CONTROL_SUBC YCLEEDCTS *CONTROL_TIMESTEPEDCURVE *DEFINE_CURVEEDDAMP *DAMPING_PART_MASS*DAMPING_PART_STIFFNESSEDDC *DELETE_CONTACTEDDRELAX *CONTROL_D YN AMIC_RELAXATIONEDDUMP *DATABASE_BINARY_D3DUMPEDENERG Y *CONTROL_ENERG YEDGCALE *CONTROL_ALEEDHGLS *CONTROL_HOURGLASSEDHIST *DATABASE_HISTORY_NODE *DATABASE_HISTORY EDHTIME *DATABASE*DATABASE_BINARY_D3THDTEDINT *INTEGRATION_BEAM*INTEGRATION_SHELLEDIPART *PART_INERTIAEDIS *STRESS_INITIALIZATION*STRESS_INITIALIZATION_DISCRETEEDLCS *DEFINE_COORDINATE_VECTOR*DEFINE_COORDINATE_S YSTEM*DEFINE_VECTOREDLOAD, , (FX, FY, FZ, MX, MY, MZ) *LOAD_NODE_SET EDLOAD, , (UX, UY, UZ, ROTX, ROTY, ROTZ, VX, VY, VZ, AX, AY, AZ) *BOUNDARY_PRESCRIBED_MOTION_SET EDLOAD, , (ACLX, ACL Y, ACLZ, OMGX, OMG Y, OMGZ)*LOAD_BODY_GENERALIZEDEDLOAD, ,TEMP *LOAD_THERMAL_VARIABLEEDLOAD, , PRESS *LOAD_SEGMENT*LOAD_SEGMENT_SET*LOAD_SHELL_SETEDLOAD, , (RBUX, RBU Y, RBUZ, +RBRX, RBR Y, RBRZ, RBVX, RBV Y, RBVZ,RBOX, RBO Y, RBOZ)*BOUNDARY_PRESCRIBED_MOTION_RIGIDEDLOAD, , (RBFX, RBF Y, RBFZ, RBMX, RBM Y, RBMZ)*LOAD_RIGID_BODYEDMP,HGLS *HOURGLASSEDNB *BOUNDARY_NON_REFLECTINGEDNROT *BOUNDARY_SPC_SETEDOPT *DATABASE_FORMATEDOUT *DATABASE_OPTIONEDPART*PARTEDPVEL *SET_NODE*INITIAL_VELOCITY*INITIAL_VELOCITY_GENERATION*CHANGE_VELOCITY*CHANGE_VELOCITY_ZEROEDRC *RIGID_DEFORMABLE_CONTROLEDRD *DEFORMABLE_TO_RIGID*RIGID_DEF0RMABLE_D2R*RIGID_DEF0RMABLE_R2DEDRI *DEFORMABLE_TO_RIGID_INERTIAEDRST *DATABASE_BINARY_D3PL0TEDSHELL *CONTROL_SHELLEDSP *CHANGE_SMALL_PENETRATIONEDSTART r = d3dump nn Isd yna comma nd line optio n EDTERM *TERMINATION_NODE*TERMINATION_BOD YEDVEL *SET_NODE*INITIAL VELOCITYINITIAL VELOCITY GENERATION*CHANGE_VELOCITY*CHANGE_VELOCITY_ZEROEDWELD *CONSTRAINED_SPOTWELD。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
7. CONTROL_HOURGLASS 用于沙漏的控制,通过施加附加力来阻止沙漏变形。IHQ设置为5,QH默认为0.1。
8. CONTROL_OUTPUT 用于设置输出参数。设置如下,按红色线框内设置,其余默认值。
•������ 第一列的“$”表示该行是注释行
•������ 输入的参数可以是固定格式或者用逗号分开 •������ 空格或者0 参数������ 使用该参数的默认值
控制卡片的建立
控制卡片可通过以下方式建立:
•用hypermesh在LS-DYNA模板下,选择Analysis面板点击 control cards,选择相应卡片。 •直接在key文件中输入
汽车碰撞分析LS-DYNA 控制卡片的设置
作者:张远岭
2011-4-14
控制卡片
碰撞分析控制卡片包括求解控制和结果输出控制,其中KEYWORD、 CONTROL_TERMINATION、 DATABASE_BINARY_D3PLOT是必不可少的。其他一 些控制卡片如沙漏能控制、时间步控制、接触控制等则对计算过程进行控 制,以便在发现模型中存在错误时及时的终止程序。 后面将逐一介绍碰撞分析中经常用到的控制卡片,并对每个卡片的作 用进行说明。
下面介绍在hypermesh中给出碰撞分析中经常使用的卡片的参
数设置
控制卡片参数设置
1. KEYWORD KEY文件起始关键字。该卡片可不作任何设置。
2. TITLE 输入标题名称xxx。
3. CONTROL_ACCURACY 提高计算精度的控制卡片。设置INN值为2,其余默认,不起作用。
控制卡片参数设置
4. CONTROL_BULK_VISCOSITY
体积黏性项q,人工施加压力,用于处理应力波传播突变引起的不连续。Q1,Q2设置为默认,IBQ设为-1。
5. CONTROL_CONTACT 用于接触面的计算修改。控制摩擦惩罚系数、接触初始穿透检查、接触厚度等,设置见红框,其余默认。
控制卡片参数设置
12. DATABASE_BINARY_D3DUMP 设置如下:
13. DATABASE_BINARY_D. DATABASE_BINARY_D3THDT 设置如下:
15. DATABASE_BINARY_RUNRSF 设置如下:
16. DATABASE_EXTENT_BINARY 设置如下:
控制卡片参数说明
*CONTROL_TIMESTEP(时间步长控制卡片) $ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MSIST 0.0 0.9 2 0.0 -0.001 0 1 1 $ DT2MSF DT2MSLC 计算所需时间步长时,要检查所有的单元。出于稳定性原因,用0.9(缺省)来 减小时间步:Δt = 0.9 l/c ,特征长度l,和波的传播速度c,都与单元的类型有关。 DTINIT:初始时间步长,如为0.0,由DYNA自行决定初始步长; TSSFAC:时间步长缩放系数,用于确定新的时间步长。默认为0.9,当计算不稳定时,可以减小该值,但同时 增加计算时间; ISDO:计算4节点壳单元时间步长的(不同的值对应特征长度的不同算法,推荐使用2,因为此选项可以获得 最大的时间步长,但有三角形单元存在时会导致计算不稳定); TSLIMT:壳单元最小时间步分配 ,使单元的时间步长控制在最小时间步长之上;只适用于使用 *mat_plastic_kinematic,*mat_power_law_plasticity*mat_strain_rate_dependent_plasticity,*mat_piecewise_linear_pla sticity等材料模型的壳单元,不建议使用该选项,因为使用DT2MS选项更好。 DT2MS:因质量缩放计算得到的时间步长。当设置为一个负值时,初始时间将不会小于TSSFAC*|DT2MS|。质 量只是增加到时间步小于TSSAFC*|DT2MS|的单元上。当质量缩放可接受时,推荐用这种方法。用这种方法时 质量增加是有限的,过多的增加质量会导致计算终止。当设置为正值时,初始时间步长不会小于DT2MS。单 元质量会增加或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终 止,但更难以作出合理的解释。默认为0.0,不进行质量缩放; LCTM:限制最大时间步长的Load-curve,该曲线定义最大允许时间步长和时间的关系(可选择) ; ERODE:当计算时间步长小于TSMIN(最小时间步长)时体单元和t-shell被自动删除。
控制卡片参数设置
17. DATABASE_BINARY_RUNRSF
设置如下:
控制卡片参数说明
*CONTROL_TERMINATION $ ENDTIM ENDCYC DTMIN 150 0 0.0
ENDENG 0.0
ENDMAS 0.0
ENDTIM:强制终止计算时间,必选,默认0.0; ENDCYC:终止循环。在计算终止时间之前,程序达到指定循环次数即终止计算。 DTMIN:确定最小时间步长TSMIN的因子。TSMIN=DTMIN * DTSTART,其中DTSTART为程序自动确定的 初始步长。当迭代步长小于TSMIN时,程序终止。 ENDENG:能量改变百分比,超过设定值则终止计算。默认0.0,不起作用; ENDMAS:质量变化百分比,超过设定值则终止计算。仅用于质量缩放DT2MS被使用时。默认0.0 ,不起作 用。
控制卡片参数设置
8. CONTROL_PARALLEL 并行计算控制,最大可以使用4个CPU。
9. CONTROL_SHELL
壳单元控制
10. CONTROL_TERMINATION 计算终止控制卡片,控制计算终止时间,设置如下:
控制卡片参数设置
11. CONTROL_TIMESTEP 计算时间步长相关参数
控制卡片使用规则
卡片相应的使用规则如下:
•������ 大部分的命令是由下划线分开的字符串, 如*control_hourglass •������ 字符可以是大写或小写 •������ 在输入文件中,命令的顺序是不重要的(除了*keyword 和*define_table) •������ 关键字命令必须左对齐,以*号开始