- 积分
- 58
- 注册时间
- 2010-4-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2010-11-18 18:36:21
|
显示全部楼层
来自 天津
本帖最后由 wyc412721 于 2010-11-20 13:45 编辑
采纳lojade兄的意见,有愿意者可以跟帖报名,也可以申请自己愿意负责的部分。并将自己整理后的资料编辑在所在楼层。
我在爆炸模拟帖里整理了关于起爆点设置、流固耦合、接触设置中SOFT的使用和LOAD_BLAST的用法。
一、关于起爆点设置的问题:
LS-DYNA中是通过点来定义起爆的,即关键字*INITIAL_DETONTION,解释一下下面card的意思,第一个PID,定义炸药的ID,EQ:-1,考虑声学边界,可参考注释中的公式,详细可查看关键字手册;EQ:0,考虑所有炸药。下面三个数据定义点坐标X,Y,Z的值,LT定义起爆时间。如果考虑声学边界,则要用到第二个card。
下面纠正一个错误,一直以来很多人都认为设置两个起爆点即为线起爆,这种说法并不准确,k手册中并没有提及到线或面起爆,只是通过点来定义起爆的,因此,定义两个点就叫做线起爆是不对的,从另一个方面分析:线起爆只能通过多个点来近似模拟,并不是真正意义上的线起爆。
关于线起爆的说法,在时老师的书中提到过,在一个例子中,说是采用中心线起爆,设置了两个起爆点。这个问题我问过时老师,设置两个起爆点就说是线起爆是不准确的,线起爆可以通过多点近似模拟,并不能真正实现线起爆。
所以ls_dyna中只能说是点起爆或者是多点起爆。
补充一下,如果不设置起爆点,则系统会默认为考虑所有炸药同时起爆,下面是一个简单的例子,例子中未设置起爆点,K文件中包含3个part,part1、2为炸药,part3为空气,均采用共节点。为方便计算,采用“准二维”方式建模,即z轴只有一个单元的厚度。
二、流固耦合
流固耦合用以处理流体和结构相互作用的问题,说白了 ,就是ALE(或者Euler)和lagrange单元间的耦合。论坛中有中说法:流固耦合有三种方法,共节点法,接触法,和真正意义上的流固耦合,我觉得这种说法是不准确的,下面解释一下,流固耦合最明显的特征就是在建立模型时,流体和结构(固体)之间必须有重合部分,一般是流体的网格包含结构的网格,举例说明一下,如时党勇一书中,炸药在土壤中的爆炸一例,炸药、土壤、空气、和混凝土板,前三个定义为ALE(流体)单元,混凝土板定义为LAG单元,流体网格包含结构网格,也就是流体和固体之间有重合的网格。
但是接触法和共节点法并不是流固耦合,接触法和共节点法也可以处理流体和固体之间的作用,但不是流固耦合.流固耦合可以包含接触法和共节点法两种方法,即在针对有多个流体或者多个固体,不同的流体之间(或者不同的固体之间)可以定义为接触,也可以根据具体情况,用共节点,但是流体和固体之间必须用流固耦合,如果用接触或共节点的话,那就不叫流固耦合了。例如,在时党勇一书中,在钢板上爆炸的例子中,空气、炸药共节点,空气和钢板有重复网格,将钢板耦合在炸药和空气的流场中;在土壤中爆炸一例,炸药、土壤、空气是共节点,空气网格包含混凝土板网格,混凝土板耦合在空气和土壤的流场中。后面还用很多用到的例子,自己可以分析总结,活学活用,掌握主要特征就可以了。
流固耦合的资料,在论坛推荐的几本资料里都有介绍,比如赵海鸥、李裕春的书中都有介绍,还有就是关键字手册,这些都是最基础的,关于论文方面,中文论文期刊的可以参考一下,主要是应用;外文的看看读一下,最好看原始文献。
另外就是要多做些练习,积累些经验,分析总结。
三、接触定义中soft的使用小结
感谢new_little前辈(http://forum.simwe.com/viewthread.php?tid=127229&highlight=%BD%D3%B4%A5%B6%A8%D2%E5)
总结的已经很详细了。
关于接触定义中soft选项的使用:
dyna的接触定义的可选卡片A中有个soft选项,有0、1、2三个选项,很多人对该选项的使用不清楚,下面就把有关此选项的设置加以解释。
实际上,当soft=1时,与默认的罚函数接触算法并无本质区别,soft=1除了在接触刚度上的确定方法有所不同之外,其他的和默认方法是一样的。soft=1在计算接触刚度时考虑了时间步长以确保计算的稳定。换句话说,你可以将soft=1时的接触看成一组简单的弹簧系统,每根弹簧都具有一个和实际计算使用的时步相匹配的Courant时步。在接触双方的弹性模量相差很大时,或是接触对之间的网格密度不一致时,用soft=1选项可以使接触定义更为真实有效。
当soft=1时,我们使用的是soft=0与soft=1中最大的刚度值,因此,当soft=0时的刚度值更大时,减小SOFSCL参数值是没有作用的。
k = max(SLSFAC*SFS*k0, SOFSCL*k1)
其中:
k-罚刚度值;
SLSFAC-*CONTROL_CONTACT中的用户输入值;
SFS-*CONTACT第三个卡片中的输入值;
SOFSCL-*CONTACT 可选卡片 A中的输入值;
k0-由材料的体积模量和单元维数计算得出的刚度值;
k1-由节点质量和求解时步得出的刚度值。
需要注意:对于双向接触类型,如*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE,需要用SFM代替SFS 参数。
当soft=2时,接触算法为基于segment的接触算法,该方法基于实际的时间步长来计算接触刚度,此时在d3hsp文件中的报告的接触步长是无意义的。
模型当中的初始穿透在计算过程中会被保留,并以此作为基线来探测额外的穿透并以额外穿透量为基准计算接触力。初始穿透不产生任何力。如果在计算过程中接触片分开到某一位置,使得当前穿透量小于初始穿透量,则程序以此较小得穿透量作为新的基线。这种算法处理网格不匹配的接触问题非常有效,可以消除接触对间的互相穿透。
四、关于LOAD_BLAST应用小结
1、感谢余版主提供的文献(http://forum.simwe.com/viewthread.php?tid=833049&highlight=)和tianhu、lovehf的讨论(http://forum.simwe.com/viewthread.php?tid=833208&highlight=LOAD%2BBLAST)
2、load_blast不能应用于封闭场合,仅用在非常简单的场合,例如一个平板(没有任何遮挡)在爆炸载荷下的作用,其他的任何有限空间,坑道效应,拐角效应,曲面板等,应用该关键字会出现很大的偏差,不要迷信该关键字。它仅仅应用在最简单的场合。
3、介绍一下*LOAD_BLAST关键字,他通常和*LOAD_SEGMENT,*LOAD_SEGMENT,或者*LOAD_SHELL。
WGT:TNT的等效质量;XBO,YBO,ZBO: X,Y,Z轴起爆点坐标;TBO:零点起爆;IUNIT:可以设置单位制;ISURF:起爆类型;CFM,CFL,CFT,CFP:分别为质量、长度、时间、压力单位的变换因数。
4、介绍一下LOAD_BLAST的用法,只使用于solid模型。
(1)关键字手册中明确注明:A minimum of two load curves , even if unreferenced , must be present in the model.可见两个曲线必须要定义,数据具体多少无关紧要,但必须有2条曲线,可以是任何形式。
(2)定义*LOAD_BLAST关键字。
(3)定义*LOAD_SEGMENT_SET
$(Defines which shell to apply *LOAD_BLAST)
//其中SID是爆炸压力作用的面ID//
(4)定义*SET_SEGMENT
*LOAD_BLAST就是通过上面的步骤加载在SID的面上的,以上四个关键字缺一不可。
关于*SET_SEGMENT如何产生,最简单的方法是将要加载的面定义为透射边界,然后在关键字中把透射关键字删除就剩下*SET_SEGMENT了,置于SID是如何定义,只要和*LOAD_SEGMENT_SET中的SID保持一致即可。 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|