LS-DYNA爆炸模拟资料

《LS-DYNA爆炸模拟资料》由会员分享，可在线阅读，更多相关《LS-DYNA爆炸模拟资料（11页珍藏版）》请在装配图网上搜索。

1、LS-DYNA爆炸模拟资料及技术问答汇总专帖i=s 本帖最后由 wyc412721 于 2010-10-30 22:38 编辑 /i:r t!o sqc9Z:zf,Dc8eYPw&F1| bsize=4/size/b.L$sObsize=4/size/b2G c0MF4U F; bsize=4color=red注：本帖属于资料及技术问题讨论帖，感谢辛苦之类的话就不要留言了。/color/size/b-?%Q6XtNbbsize=4/size/bG:J2l$aP6;C bsize=4/size/b6JFq.oBE+Vbsize=41、爆炸模拟的一般过程/size/bn bsize=4%8XX,+

2、K /size/b（1）首先要有一个总体思路，怎样font=宋体把爆炸程转化为程序模型（明白该怎样建立自己的模型，选取材料模型，采用哪儿种算法计算求解），这一步很重要，直接关系计算结果的合理性。/font9M)nca:a68x,y （2）前处理，包括：画图、建模、划分网格、定义边界及约束条件、设定分析控制选项、输出K文件等。%HBK2Mt/q$LO（3）有些前处理不支持LS-DYNA的K文件中关键字，需修改完善K文件（注：K文件的中关键字格式十分重要，需要注意）。;jg7Z Y%JC5RH bcolor=red/color/b（4）计算求解。如果有错误终止计算，根据错误提示，查找并修改出错的地

3、方，使计算正常结束。-Vl#P7UJ $I （5）后处理。结果分析，提取自己关心的数据，写报告。:i-yTIhynM d9j-_ DA3bsize=42、爆炸模拟的方法/size/b%w,PU9ofKCP w,u V1Eg（1）共用节点算法模拟爆炸yY:H!d6Q80O（2）接触耦合算法模拟爆炸.reQ0nYC F （3）流固耦合算法模拟爆炸 （详见白金泽LS-DYNA3D理论基础与实例分析）Ad h,*YMnF:vxyA9Y H size=4b3、推荐资料/b/size4U#Or)M86U3H /P （1）时党勇基于ANSYS/LS-DYNA8.1进行显式动力分析，论坛有电子版（不全），该书

4、重点给出了18个经典算例的操作求解过程，主要是通过GUI建模，非常适合于初学者。NG Ld y/M5d9C6mn8I a6G p（2）白金泽LS-DYNA3D理论基础与实例分析，第五章介绍了爆炸模拟常用的方法及建模过程，主要通过有命令流建模，并对比了三种方法的特点。Yp2| $? YtC#i q #n#NM4Qb!6 JV（3）赵海鸥LS-DYNA动力分析指南中介绍了流固耦合算法和爆炸模拟，侧重介绍关键字。9O#L2D0Zrr8Wet:d （4）尚晓江ANSYS_LS-DYNA动力分析方法与工程实例，有一章分别介绍关键字和工程爆破问题。 k0*?1xOO cqZ6C$t4tYov/p（5）LS

5、-DYNA关键字手册j-Y7g A )n;l*su)f8F color=redN$UHH:UYJBf p/_ 以上资料都可在公共邮箱中下载5LS(z,i1W /coloremail=ls_dyna2010color=redls_dyna2010/color/emailcolor=red，密码：12345607G/Y4c/l*?9W1c,a iG&j$!KXt%uaW/color（6）爆炸方面的理论知识，爆炸动力学及其应用、凝聚炸药起爆动力学、爆轰理论基础、水下爆炸、猛炸药爆轰学Q)i*R K_,D1o j4 T url=#v?Suv&J_ t/o.E?8U （7）应用DYNA3D进行爆炸分析9

6、AMY!_%kurl=（8）一篇炮轰模拟培训的pptpBQ+j1Xk(N$ie- attach256994/attach1L(fd rpW6H6s?j*H8B2G （9）一篇模拟爆炸的英文文献SIMULATION OF LANDMINE EXPLOSIONUSING LS-DYNA3D SOFTWARE+KI1R IK9d W url= r Jx6U/wGG8Y(X（10）有关水下爆炸网格划分的一篇综合文章o d5dI;Bm url= sB/C |)Z （11）论文资料FV F4H)url= js4H j_7SiE!0Q(saG-G;G%c bsize=44、论坛中爆炸模拟的作品和算例（含K文

7、件）/size/bqe,ip5PSSm2&A 5SG)Zu J?（1）用ALE方法做的爆炸模型J?%s:bN5r ?)surl= eSJ:aY/K9z cL（2）ls-dyna 水下爆炸对结构的冲击经典算冽SE8UQ#P6v$U-y5VV url= ohpzVX:T7Q,w3c#mn*YA（3）柱壳内部爆炸的例子S RCFk2E- color=royalblueurl=l(IA6b9r,B（4）带水容器壁在爆炸荷载作用下的耦合W)cHencolor=royalblueurl=p(Ik:T(H#9l C4rQt z （5）爆炸成型弹丸侵彻钢靶的ALE 算法U/tFzKE#u3E url=gPwT

8、$R W&F0vx&N1G （6）爆炸模拟的算例 color=royalblueurl= Iea p7O*t6Y j color=#4169e1/color*m u(K9g5PEG color=royalblueurl=6L1l)yH%?wrR7la 8rA2m%Xcolor=royalblueurl=_W2W go ToZ#Au)l color=royalblueurl=*)a gM.color=royalblueurl=S4hE) color=royalblueurl=w5O&r Gp e9W.Bqurl=oqp4n4L391f%O +k (LQO%h0Yl eU3ri*Ebsize=45

9、、论坛中讨论帖（目前正在对讨论帖资料进行汇总，将个知识点的总结帖在楼下个层，欢迎广大会员积极参与，视情况将会给以积分奖励）/size/bZ#T;Q mU-X1x1m*o0H ZZ （1）也谈流固耦合.OC24Z q url=jb*D1y3Q BT%s/? J!m!c G-YT （2）原创共享DYNA爆炸流场计算算法概述&B-|$u KN(yurl= M4W#Ky2U;c3G$qhi/De （3）流固耦合分析爆炸过程中的各种能量问题讨论energy4g7iY3S8S?- url= w1bY$w3eoy,BY （4）【讨论】这两种爆炸方式有什么区别？;k7qSJ;R, +nK,rdurl=)*J0

10、a1JR 3l.Tv tN9qTo （5）流固耦合模拟爆炸问题:F*daB%G#Qz1Vurl= v*s3OH/U eA$HF J0? x |J*KbHnb（6）在关于岩石、混凝土爆炸（流固耦合）中，空气起什么作用？&Tj2ee url=!I;j3QP B;K 3Gt!Ru （7）LS_DYNA 能让炸药和空气，在爆炸一定时间后失效吗？$De+w2L9 YP H&a url= rCD6Z#?fs/K$kZ/Y0Up#mM（8）求助：空中爆炸冲击波模拟问题Q8E%fQ:Ymgvurl=W o0Zkbc&W HM MPu5Y2s （9）关于ANSYS/LS-DYNA8.1显示动力分析书上例4.2炸

11、药在土壤内部爆炸作用的一点不解#faoqy(yj,H OS url=f,DoQd( )iD:V$? （10）向高手请教爆炸计算问题OjXUZ6A url= XIw ud!j#L(V o*Bb C Y （11）炸药不爆炸EFjR%T C:Surl=#iG!u$g m zv1L&qVj %W6VCRu |&W MW7B.（12）空气爆炸边界条件怎么设定;x0C fC !c6vurl=wpet W:o （13）请教各位大虾：一个爆炸问题的K文件 +pW v l $VU:Nurl=CyHx;n YW7| (Z kg?W （14）爆炸中的经典问题求助！$YS+dt ix9z url=xE M:x /pS

12、 E6VQ U （15）请教大家无限水域中爆炸，怎么结果不合理呢？4kG/0sY.e2J,iAurl=#f#O6y kL(Rlr 3F4F9X2B G1u2V f（16）爆炸问题求助C7Q-a%qfm(|Wurl= T O bk A/wo*Q O A.u%U1bv/z（17）关于爆炸模拟的一些问题请教5WRG0y9l url=8i)J%a v-X+G （18）1公斤炸药在水下爆炸，传播1米后的压力还能剩多少？$O;rh/Gk r!+x_ Ourl=0jdBH3F(_ *4T OY(Z （19）请教：为什么我做的流固耦合爆炸被炸物没什么反应？M dUm url= kf(XX1W,u/ux4yCD

13、T0x%K （20）帮忙解释一下为什么爆炸波从岩石传播不到土壤层yt9f1D5%A3(iz url= gdyu_yu 发表于 2010-8-5 00:28我觉得可以先以wyc412721兄的这个爆炸模拟专帖为示范，摸索资料整理与归总的合理之路。同时我认为，按专业领域开若干个专帖进行各自的整理与归总是合理的。1dU X#Za *G WY+$P q$R1?$f&a:p3K i 下面提几点意见。未尽事宜，另请补充。 #N%u5obRP/u.j rLR TD （1）每一个专帖在主楼（1楼）总纲式给出以下方面：介绍该领域的基本情况或常识（扫盲之用）、该领域常用的前后处理软件（最好有简要的优缺点比较）、模

14、拟的一般流程及常用方法、推荐的知识来源（主要是有代表性的书籍、网站以及各自的链接地址）、本版中该领域相关的精华帖与讨论帖链接（挑选若干有选择性的帖子）、本版该领域模拟中的问题及解答总集（参考YUMINHUST2005斑竹的LS-dyna 常见问题汇总2.0的格式。这项工作持续时间最长，工作量最大也最重要。以WORD文档附件提供。阶段性的成果以1、2、3给出。）v,p*bJ2_eb5k （2）除主楼之外，占24层楼以备日后操作之用。7OQj5$W99ls y&S1y n2RO Z （3）开放权限，号召全体会员参与资料的整理与汇总。参与者先在所在专帖中报到，然后在所在的报到楼层中负责各自的整理工作

15、。f!$h9_$_ 9t%$i9Y+QnY （4）对于纯资料帖，只需要给出主题帖的名称及相应的链接地址（鼓励同时对该资料进行简要评价）；对于非资料帖，除了给出主题帖的名称和链接之外，要求对所提供的帖子内容进行问答整理与汇总。（在对应的楼层中上传整理的WORD文档附件。）5?Dy3f5#WwKr.s2I4H$dq+?5k （5）每个专帖指定专人负责。若楼主为本版管理员，则直接由该员负责该专帖资料的总编。若楼主为其他会员，则令指定1名本版管理员与其共同负责总编（需要协调好）。 pv3p&Cy0HI;u*u6R;S x （6）对所有参与资料整理与汇总的会员先予以50到100个贡献分鼓励。对于工作出色

16、的参与者，视其贡献大小再追加技术分，以示肯定。z,_-a C7GySW?cq#hxl+w（7）负责该专帖的本版管理员不定期对整理好的纯资料帖进行关闭主题操作。在日常的整理中，需要与参与的其他会员及时沟通，以便工作更高效、有序。 蓝牙 发表于 2010-8-5 17:41i=s 本帖最后由 蓝牙 于 2010-8-5 17:42 编辑 /i$XA1GSrv |F.j CB+A-K2U目前在爆炸方面;U7hy$S （1）共用节点算法模拟爆炸CO _4x0e5e-Z（2）接触耦合算法模拟爆炸v6tX4K5a/y4X X（3）流固耦合算法模拟爆炸，h8eIF5X流固耦合用的可能是最多的，它的实用范围可

17、能也是最广泛的，通过论坛答疑，反映最多的主要有如下几个问题k Ie-t 1.欧拉流体与拉格朗日SHELL和SOLID的耦合问题，我知道，在拉格朗日shell和流体的耦合方面可以通过*constraind_lagrange_In_solid_edge进行改进，是否有更好的方法，可以更好的控制渗漏等方面的问题呢？希望有人能跟进qFi8DnaL 2.爆炸的超压问题，目前验证计算超压是否正确，主要是和经验公式对比，在我的验证过程中发现网格的尺寸对超压的精确度有不小的影响，我目前只能不断细化网格然后通过验证两次结果的差别来校对，是否有达人总结过当量和网格的对应关系？sE1bek9E 3.LOAD_bla

18、st 可以用来替代欧拉流体，我总觉得其计算精度要高于FSI，同时网格规模大幅度减少，但是为什么用的人不是特别多，有些不解啊！ wyc412721 发表于 2010-9-3 19:19i=s 本帖最后由 wyc412721 于 2010-11-13 00:42 编辑 /iY+g.T-q$|2Q 关于起爆点设置的问题：1l e,h Cr D4CC LS-DYNA中是通过点来定义起爆的，即关键字*INITIAL_DETONTION，解释一下下面card的意思，第一个PID，定义炸药的ID，EQ:-1，考虑声学边界，可参考注释中的公式，详细可查看关键字手册；EQ:0，考虑所有炸药。下面三个数据定义点坐

19、标X,Y,Z的值，LT定义起爆时间。如果考虑声学边界，则要用到第二个card。BoJnE0q ga &r/H AAh下面纠正一个错误，一直以来很多人都认为设置两个起爆点即为线起爆，这种说法并不准确，k手册中并没有提及到线或面起爆，只是通过点来定义起爆的，因此，定义两个点就叫做线起爆是不对的，从另一个方面分析：线起爆只能通过多个点来近似模拟，并不是真正意义上的线起爆。T3X7Q Ju:c +Y 关于线起爆的说法，在时老师的书中提到过，在一个例子中，说是采用中心线起爆，设置了两个起爆点。这个问题我问过时老师，设置两个起爆点就说是线起爆是不准确的，线起爆可以通过多点近似模拟，并不能真正实现线起爆。:

20、YzMuzS#T所以ls_dyna中只能说是点起爆或者是多点起爆。3uV0lvg16gf补充一下，如果不设置起爆点，则系统会默认为考虑所有炸药同时起爆，下面是一个简单的例子，例子中未设置起爆点，K文件中包含3个part，part1、2为炸药，part3为空气，均采用共节点。为方便计算，采用“准二维”方式建模，即z轴只有一个单元的厚度。4C!w2S(I*K attach269941/attach 直面人生 发表于 2010-9-20 16:24对于超压的计算，在dyna中是不够理想，不过你可以通过autodye进行一下模拟如果要划分的网格较大的话最好采用二维模型，计算结果挺好的b url= i蓝

21、牙/i /b wyc412721 发表于 2010-10-20 19:22i=s 本帖最后由 wyc412721 于 2010-10-20 19:38 编辑 /iQA;O1I|n-hOs &y.BIqU3Iz9iJ1h论坛有很多讨论关于流固耦合的帖子,下面根据我的使用经验做以下总结,不足之处希望大家指教.Y zJP;C%W hT*s 流固耦合用以处理流体和结构相互作用的问题，说白了 ，就是ALE（或者Euler）和lagrange单元间的耦合。论坛中有中说法:流固耦合有三种方法,共节点法,接触法,和真正意义上的流固耦合,我觉得这种说法是不准确的,下面解释一下,流固耦合最明显的特征就是在建立模型

22、时,流体和结构(固体)之间必须有重合部分，一般是流体的网格包含结构的网格，举例说明一下，如时党勇一书中,炸药在土壤中的爆炸一例,炸药、土壤、空气、和混凝土板，前三个定义为ALE（流体）单元，混凝土板定义为LAG单元，流体网格包含结构网格，也就是流体和固体之间有重合的网格。G;O(U(t EY u;pD/j5u0S7g,v但是接触法和共节点法并不是流固耦合,接触法和共节点法也可以处理流体和固体之间的作用,但不是流固耦合.流固耦合可以包含接触法和共节点法两种方法,即在针对有多个流体或者多个固体,不同的流体之间（或者不同的固体之间）可以定义为接触,也可以根据具体情况，用共节点，但是流体和固体之间必须

23、用流固耦合，如果用接触或共节点的话，那就不叫流固耦合了。例如,在时党勇一书中，在钢板上爆炸的例子中，空气、炸药共节点，空气和钢板有重复网格，将钢板耦合在炸药和空气的流场中；在土壤中爆炸一例，炸药、土壤、空气是共节点，空气网格包含混凝土板网格，混凝土板耦合在空气和土壤的流场中。后面还用很多用到的例子，自己可以分析总结,活学活用,掌握主要特征就可以了。WX.Gl7?X流固耦合的资料，在论坛推荐的几本资料里都有介绍，比如赵海鸥、李裕春的书中都有介绍，还有就是关键字手册，这些都是最基础的，关于论文方面，中文论文期刊的可以参考一下，主要是应用；外文的看看读一下，最好看原始文献。k e!Kv$zC u#Q

24、!G5Eq7KHj 另外就是要多做些练习，积累些经验，分析总结。 wyc412721 发表于 2010-11-18 06:59接触定义中soft的使用小结!L1br9He Z n:YrXY nm v7pi)?3QxX 感谢new_little前辈（!E/W2ud!m 总结的已经很详细了。TX.w$w,H:s I n1Xv 2yE0MO:z9k%f8AZ 关于接触定义中soft选项的使用: j*yb!Lt NQ%iop jrm T8S6F dyna的接触定义的可选卡片A中有个soft选项，有0、1、2三个选项，很多人对该选项的使用不清楚，下面就把有关此选项的设置加以解释。 MKMQ2Bu q1x

25、&?2:kb PI实际上，当soft1时，与默认的罚函数接触算法并无本质区别，soft1除了在接触刚度上的确定方法有所不同之外，其他的和默认方法是一样的。soft1在计算接触刚度时考虑了时间步长以确保计算的稳定。换句话说，你可以将soft1时的接触看成一组简单的弹簧系统，每根弹簧都具有一个和实际计算使用的时步相匹配的Courant时步。在接触双方的弹性模量相差很大时，或是接触对之间的网格密度不一致时，用soft1选项可以使接触定义更为真实有效。 F%Ok2ART XA8L$Q*U 当soft1时，我们使用的是soft0与soft1中最大的刚度值，因此，当soft0时的刚度值更大时，减小SOFS

26、CL参数值是没有作用的。 Tn(aF-XG.W8lT&x1E A+U5H k = max(SLSFAC*SFS*k0, SOFSCL*k1) U4N0hcxl i4b IF 其中： $T.n0G(_.w& k罚刚度值； (G#r;Z,|-l_ dRSLSFAC*CONTROL_CONTACT中的用户输入值； X OZ ZiSFS*CONTACT第三个卡片中的输入值； lD |5/O:S&SOFSCL*CONTACT 可选卡片 A中的输入值； !r. N$Jk0由材料的体积模量和单元维数计算得出的刚度值； k1由节点质量和求解时步得出的刚度值。 G.Jb0N)LM:qV3 pd- a 需要注意：

27、对于双向接触类型，如*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE，需要用SFM代替SFS 参数。 Z c;LM5OP5?ODU sSFj当soft2时，接触算法为基于segment的接触算法，该方法基于实际的时间步长来计算接触刚度，此时在d3hsp文件中的报告的接触步长是无意义的。 2R+?)Wl#K r(Ywle*Ce!m1u2c 模型当中的初始穿透在计算过程中会被保留，并以此作为基线来探测额外的穿透并以额外穿透量为基准计算接触力。初始穿透不产生任何力。如果在计算过程中接触片分开到某一位置，使得当前穿透量小于初始穿透量，则程序以此较小得穿透量作为新的基线。这种算

28、法处理网格不匹配的接触问题非常有效，可以消除接触对间的互相穿透。 wyc412721 发表于 2010-11-18 11:35i=s 本帖最后由 wyc412721 于 2010-11-18 11:37 编辑 /ivV d%q0HMI%w ?-t kQHrA 关于LOAD_BLAST应用小结IoJ:_C(zN 1、感谢余版主提供的文献（ u2、load_blast不能应用于封闭场合，仅用在非常简单的场合，例如一个平板（没有任何遮挡）在爆炸载荷下的作用，其他的任何有限空间，坑道效应，拐角效应，曲面板等，应用该关键字会出现很大的偏差，不要迷信该关键字。它仅仅应用在最简单的场合。f$jy*xnr)L

29、)vl 5rNCuy b3、介绍一下*LOAD_BLAST关键字，他通常和*LOAD_SEGMENT，*LOAD_SEGMENT，或者*LOAD_SHELL。j9xt$eh7if7H$ WGT:TNT的等效质量；XBO,YBO,ZBO: X,Y,Z轴起爆点坐标；TBO:零点起爆；IUNIT:可以设置单位制；ISURF:起爆类型；CFM,CFL,CFT,CFP:分别为质量、长度、时间、压力单位的变换因数。/aan5 u?g 4、介绍一下LOAD_BLAST的用法，只使用于solid模型。FF*bAGlX5F （1）关键字手册中明确注明：A minimum of two load curves ,

30、 even if unreferenced , must be present in the model.可见两个曲线必须要定义，数据具体多少无关紧要，但必须有2条曲线,可以是任何形式。KJPE je_9Y j （2）定义*LOAD_BLAST关键字。8u&VL)m. R!k$b.k （3）定义*LOAD_SEGMENT_SETb6G+69I%d n$(Defines which shell to apply *LOAD_BLAST)0bGN _prBG4gQd_INef /其中SID是爆炸压力作用的面ID/CTb$e0_s_（4）定义*SET_SEGMENTs ytX1q5SaK!W*U _

31、3t9Vs/?2,P*LOAD_BLAST就是通过上面的步骤加载在SID的面上的，以上四个关键字缺一不可。)CPF U q1A5DG-关于*SET_SEGMENT如何产生，最简单的方法是将要加载的面定义为透射边界，然后在关键字中把透射关键字删除就剩下*SET_SEGMENT了，置于SID是如何定义，只要和*LOAD_SEGMENT_SET中的SID保持一致即可。 herovsdevil 发表于 2011-1-15 18:50万分感谢，正需求:) dt_lwd 发表于 2011-2-14 11:54真是万分感激，楼主辛苦了！ menlatin 发表于 2011-3-15 23:16非常感谢！辛苦

32、了！ jangking2010 发表于 2011-3-15 23:58我是新手，感谢楼主，感觉这些东西很有用！ wsdzyt 发表于 2011-4-1 09:57多谢楼主无私奉献 dong200712 发表于 2011-5-14 23:52这个资料好得没法说啊，直接收藏！ dong200712 发表于 2011-5-14 23:53这个资料好得没法说啊，直接收藏！我也是做爆炸模拟的！ NWLNLENWLNL0379 发表于 2011-5-16 15:27非常感觉，总结的挺好。关于流固耦合，进行三维流固耦合时，也是流体网格包括住固体网格吗？请教。重合这个词用的有点不精确，可以是包含，也可以是有公共部分啊，能否准确说明一下，谢谢 wgl0635 发表于 2011-5-18 22:18非常感谢楼主，我在模拟过程中遇到是问题是在爆破发生后很长时间后岩石材料的应力值仍然没有回到零点，是不是我的材料参数有问题啊，没有什么合适的参数可以参考呢？11