按天归档:16/11/2011

[转帖]win32下lsdyna mpich并行计算试验

win32下lsdyna mpich并行计算试验

今天尝试了lsdyna的并行计算,过程如下: 
1. 安装virtualbox虚拟机到自己的笔记本,虚拟机中安装xp,克隆成两个,与笔记本构成一个小型PC集群,主机名和IP分别为host1:192.168.1.1(笔记本) ;host2:192.168.1.2(虚拟的xp1);host3:192.168.1.3(虚拟的xp2)。在这三个主机中创建相同的管理员帐户如 admin , 密码 123 
2. 在三个主机中安装mpich1.2.5,安装好后打开“任务管理器”中的“进程”选项卡,查看是否有一个mpd.exe的进程。如果有的话说明安装成功。 
3. 安装好mpich之后还必须对每台主机进行注册才能使用:三个主机分别执行程序目录(默认为C:\Program Files\MPICH\mpd\bin\)下的mpiregister.exe,首先会提示输入用户账号(第一步中创建的admin),然后会提示输入两遍密码,之后会问你是否保持上面的设定。如果选择是,则上面的信息将写入硬盘,否则保存在内存中,再重新启动之后就不存在了,这样mpich才能在网络环境中访问每台主机。 
4. 下载mpp971_s_R5.1.1_Intelsse_win32_mpich125.zip,解压出来的mpp971.exe是求解器文件,I2a.exe用于转换结果文件,将求解器放置到三个主机中无空格无中文的目录如:C:\LSDYNA\mpp\中,设置系统环境变量LSTC_FILE=license的路径 
5. 为方便操作,将mpich程序目录(C:\Program Files\MPICH\mpd\bin\)加入系统环境变量path中;在3个主机中的C:\LSDYNA\目录下放一个算例testcase.k,执行以下命令进行单机测试[code]mpirun -np 2 C:\LSDYNA\mpp\mpp971.exe i=C:\LSDYNA\testcase.k[/code],normal termination的话再进行多机并行计算,2表示cpu数

6. 执行并行计算有以下3种方法,前两种必须保证三个主机中的求解器路径相同,k文件路径也相同,如求解器都放在C:\LSDYNA\mpp\mpp971.exe ,k文件都放在C:\LSDYNA\testcase.k: 
    方法一:用MPIConfig.exe建立主机列表,执行[code]mpirun -np 6 C:\LSDYNA\mpp\mpp971.exe i=C:\LSDYNA\testcase.k[/code]这种方法最方便,MPIConfig的设置方法参见附图。 
    方法二:直接执行[code]mpirun -hosts 3 host1 2 host2 2 host3 2 C:\LSDYNA\mpp\mpp971.exe i=C:\LSDYNA\testcase.k[/code](第一个3表示host数;后面三个2表示cpu数;host1(2,3)可用ip替换,效果一样的) 
    方法三:使用配置文件,这种方法最灵活,可设置的参数也较多(请参考mpich的user guide),我的设置如下:[code] exe C:\LSDYNA\mpp\mpp971.exe i=C:\LSDYNA\testcase.k 
hosts 
192.168.1.1 2 
192.168.1.2 2 
192.168.1.3 2 
[/code]保存成一个文本文件,如cfg.txt 
执行[code] mpirun cfg.txt [/code]即可 
当然,以上三个ip也可以用host1,host2,host3替换 
计算完成后在host1中会生成主要的结果文件,host2,host3中也会生成一些文件。 
7.若求解器或k文件在三个主机中的安装位置不同,可采用上一步中第三种方法,蛋疼的假设: 
host1中求解器路径C:\LSDYNA\mpp\mpp971.exe ,k文件路径C:\LSDYNA\testcase.k; 
host2中求解器路径C:\1\mpp\mpp971.exe ,k文件路径C:\1\testcase.k; 
host3中求解器路径C:\a\mpp\mpp971.exe ,k文件路径C:\a\testcase.k; 
配置文件则设置为:[code] exe C:\LSDYNA\mpp\mpp971.exe i=C:\LSDYNA\testcase.k 
hosts 
host1 2 
host2 2 C:\1\mpp\mpp971.exe i=C:\1\testcase.k 
host3 2 C:\a\mpp\mpp971.exe i=C:\a\testcase.k 
[/code]host1无求解器路径和k文件路径,则使用第一行”exe “后面的路径。 
host2,3则使用指定路径求解器和文件 
结果文件的生成位置:在host1中执行mpirun cfg.txt命令时,其“当前目录”为结果文件的生成位置,若其他两个主机中不存在该目录,host2和host3中的结果文件会跑到C:\windows\system32里去了,k文件手册中说可设置p=pfile参数改变输出结果的位置,方法3中设置pfile似乎行不通。 
当前目录即命令提示符(cmd)中“>”前面的那串字符表示的目录,如 
C:\Documents and Settings\Administrator> 表示当前目录为“C:\Documents and Settings\Administrator\” 
mpirun的运行参数,mpiregister及mpiconfig的设置可以看mpich1.2.5程序自带的user guide的第五节Tools; 
以上过程均在win32上进行,win32的mpp_dyna不知道支持不支持mpich2,所以安装的是mpich1.2.5,64位的mpp_dyna需要mpich2,所以操作会有所不同,但总体过程相似的。 
===========

07.09 补充1:win32版mpp_dyna可支持mpich2,安装mpich1的dll文件,再安装mpich2即可。mpich2的执行程序换成了mpiexec,配置文件设置方式有所不同(见mpich2 user guide p.5-6) 
07.09 补充2:对比了下mpp版本与smp版本在单机4核8线程(伪8核)64位机器上的执行效率,计算一个DYNA自带算例airbag_deploy.k,mpp版选8核用了15s,选6核用了17s;smp版选8核用了28s,选6核用了27s。结论是1.mpp版本多核计算的效率要明显高于smp版;2.smp版不一定选最高核数时计算效率最高。

07.09 补充3:mpp计算时不能Ctrl+C,sw*,因为Ctrl+C会导致程序终止。解决方法是:新开一个cmd窗口,cd进入工作目录,键入命令[code] echo sw2 > d3kil[/code]其他命令类似,原理可见《常见问题2.0》第36问。

菊子曰 这就是菊子曰啦!

【转帖】在linux下安装和使用ls-dyna的步骤

在linux下安装和使用ls-dyna的步骤

今天安装了64位linux系统(ubuntu 11.04 server),无Xwindow界面,决定尝试下linux下的dyna计算,下面把过程叙述一遍: 
1.用另一台电脑从lstc官方网站下载64位linux求解器:ls971_s_R5_1_1_intel64_redhat54.gz,其他版本也可以,反正里面没ubuntu版的; 
2.为方便操作,将解压出来的求解器改名为ls971,通过U盘拷贝到linux机器下的/home/username(你的帐户名)/lsdyna/ 文件夹下; 
3.在/home/username/下创建keyfiles文件夹,通过U盘拷贝测试用的testcase.k文件到keyfiles下,这一步可与上一步一起完成; 
4.进入lsdyna文件夹,输入./ls971 info,./ 表示执行程序,若程序不在当前文件夹下,则程序名需要带上路径(见第7步),info参数用于生成注册文件lstc.log,但如果你安装的linux也是ubuntu,将会报错,错误为:缺少libnuma.so.1,网上可下载,放入/usr/lib64/文件夹下即可; 
5.重新执行./ls971 info,同文件夹下生成lstc.log,联系LSTC公司或者其代理商获取LICNESE文件 
6.讲得到的LICENSE文件LSTC_FILE放到/usr/local/lstc/文件夹下,并设置全局环境变量 LSTC_FILE=”/usr/local/lstc/LSTC_FILE” 
7.进入/home/username/keyfiles/文件夹,执行算例:  ./home/username/lsdyna/ls971 i=testcase.k 当然还可以设置memory、cpu等参数,与win下相同;计算结果文件会输出到当前文件夹下,应该可以通过指定o=参数更改 
8.程序正常执行,开始计算。 
若要后台计算,可Ctrl+Z挂起程序,再输入bg进行后台运行,或者在执行dyna计算时命令后面加入 & 符号可以直接后台运行。 
若想与执行ls等命令一样只输入一个单词,而不是通过 ./home/username/lsdyna/ls971 这种形式执行dyna,可将dyna程序放到/bin/目录下。

菊子曰 这就是菊子曰啦!

【转帖】LS-DYNA各版本下载及破解

LS-DYNA各版本下载及破解

首先是下载,可以进入官方服务器下载各版本: 
http://ftp.lstc.com/user 用户名:user 密码:computer 
里面可以下载到ls-dyna各种版本及相关的一些资料,很多都是免费的· 
若支持正版,请找中国代理商购买许可证,若支持盗版,请接着往下看:

用hexworks等十六进制编辑软件操作 
ls971_s_r31_win32_p_crack 
查找hex 
E878520000  83C4108945A08B450883380274088B4D088339037507    
更换成 
9090909090  83C4108945A08B450883380274088B4D088339037507  
注意,这里“E878520000 ”及“9090909090”后的空格是我加的,查找的时候不要加。下同

 

查找hex 
E864100000  83C4108945F8837DF8007405E99F000000 
更换成 
9090909090  83C4108945F8837DF8007405E99F000000

ls971_d_r31_win32_p_crack 
查找hex 
E811570000 83C410998945908955948B4508898510FDFFFF8B8D10FDFFFF    替换成       
9090909090 83C410998945908955948B4508898510FDFFFF8B8D10FDFFFF

查找hex 
E864100000 83C4108945F8837DF8007405E99F0000006 
A00E859DC350083C4048945C0 
替换成 
9090909090 83C4108945F8837DF8007405E99F0000006 
A00E859DC350083C4048945C0        
                        
mpp971_s_7600.1224.exe 
查找HEX 
E8EA8500008 98424AC000000488B842450030000833 
802740D488B842450030000833803750B 
替换为 
9090909090 898424AC000000488B842450030000833 
802740D488B842450030000833803750B

查找HEX 
E802150000 89442434837C2434007405E9C90000003 
3C9E8CBFF44004889842488000000488D8C2488000000E876004500    
替换为 
9090909090 89442434837C2434007405E9C90000003 
3C9E8CBFF44004889842488000000488D8C2488000000E876004500

ls971_d_R4.2_win32_p 
E832A20000 83C4 10998945908955948B4508898548FDFF 
改成 
9090909090 83C4 10998945908955948B4508898548FDFF

E8D40B0000 83C4108945F8837DF8007405E982000000C745F8 
改成 
90909090908 3C4108945F8837DF8007405E982000000C745F8

ls971_s_R4.2_win32_p 
E8 F0 9C 00 00 83 C4 10 89 45 A0 8B 45 08 83 38 02 74 08 8B 4D 08 83 39 03 75 07 C7 45 A0 00 00 00 00 68 98 D7 
改成 
90 90 90 90 90 83 C4 10 89 45 A0 8B 45 08 83 38 02 74 08 8B 4D 08 83 39 03 75 07 C7 45 A0 00 00 00 00 68 98 D7

E8 D4 0B 00 00 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 82 00 00 00 C7 45 F8 
改成 
90 90 90 90 90 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 82 00 00 00 C7 45 F8 

ls971_d_R3.2.1_win32_p 
E8 D3 5C 00 00 83 C4 10 99 89 45 90 89 55 94 8B 45 08 89 85 10 FD FF FF 
改成 
90 90 90 90 90 83 C4 10 99 89 45 90 89 55 94 8B 45 08 89 85 10 FD FF FF

E8 6A 11 00 00 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 A2 00 00 00 6A 00 E8 D2 00 00 00 83 C4 
改成 
90 90 90 90 90 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 A2 00 00 00 6A 00 E8 D2 00 00 00 83 C4

ls971_s_R3.2.1_win32_p 
E8 7B 58 00 00 83 C4 10 89 45 A0 8B 45 08 83 38 02 74 08 8B 4D 08 83 39 03 75 07 C7 45 A0 00 00 
改成 
90 90 90 90 90 83 C4 10 89 45 A0 8B 45 08 83 38 02 74 08 8B 4D 08 83 39 03 75 07 C7 45 A0 00 00

E8 6A 11 00 00 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 A2 00 00 00 6A 00 E8 D2 00 00 00 83 C4 
改成 
90 90 90 90 90 83 C4 10 89 45 F8 83 7D F8 00 74 05 E9 A2 00 00 00 6A 00 E8 D2 00 00 00 83 C4

ls971_s_R3_1_winx64_p-crack 
E8413A0000 8BD88B450083F802740583F8037502 
9090909090 8BD88B450083F802740583F8037502

E8B1E9FFFF  85C00F858400000033C9E8826E4C00 
9090909090 85C00F858400000033C9E8826E4C00

ls971_s_R4.2.1_win32_p 
ls971_d_R4.2.1_win32_p 
是我自己修改的,提供修改好文件的下载地址: 
http://www.brsbox.com/filebox/down/fc/c938204e2682d943c68743abef6f58b5

license生成,做法如下: 
1.在LS-DYNA Progrm Manager 中点License/Registration,按照要求填写公司名字,email,传真,电话,联系人(以上任意填),然后点save lstc.log file 。之后它会问你是否发email,不发。 
2.利用生成的lstc.log文件做一个名为lstc_file的文件。将lstc.log中有用的信息(如下例)提取出来: 
nihao 
    EMAIL: woaini@163.com 
      FAX: 58413142  
TELEPHONE: 58413142 (nihao) 
# 
PC-DYNA_970    05112005 
# 
MACHINE: E6 8D 27 41 
MACHINE: B9 62 85 2C 
上面就是一些有用的东西,程序会一个个检查,缺一不可。还有一个重要的东西要加上,就是:KEY 
以上面的例子做一个lstc_file文件如下: 
nihao 
    EMAIL: woaini@163.com 
      FAX: 58413142  
TELEPHONE: 58413142 (nihao) 
# 
PC-DYNA_970    05112005 
# 
MACHINE: E6 8D 27 41 
MACHINE: B9 62 85 2C 
# 
KEY:00 11 22 33 44 55 66 77 88 99 
上面就是lstc_file的制作过程。 
3.再次打开LS-DYNA Program Manager,点license/Import license file  然后指向你lstc_file,点打开,就OK了!

菊子曰 这就是菊子曰啦!

大胆举手发言,苦练杀敌本领!!

【考试心理学:怀疑自己答错了,改不改答案?】做完答卷,检查的时候发现两个答案不知道选哪个,这时候该不该改答案?有调查显示:55%的人认为改答案会损害成绩,只有15.5%认为会改善;但实际的情况是:58%把答案改对了,只有20%把对改错了,还有22%把错的改错了。因此,下碰到这种情况:改!

菊子曰 我用菊子曰写博客,你呢?

SecurityKISS––来自欧洲的免费多国VPN代理【转帖】

SecurityKISS是一家英国的VPN提供商,使用的是OpenVPN的网络连接,它不但能够隐藏你的真实IP地址,而且还加密你的数据连接。不象通常的PPTP VPN服务商,SecurityKISS 已经把他们的服务制作成一个自己的软件,结合了OpenVPN客户端,所以你只需要下载和安装这个软件。之后,打开软件进行连接。并不需要进行注册,也没有用户名和密码,直接就可以连接他们的服务器。

我对这个VPN做了一下测试,可以正常连接,速度也相当快。SecurityKISS 现在是免费提供,任何人都可以下载和使用。但是它的网站并没有说他们的免费VPN服务是永久免费还是暂时免费。当然,他们也是提供有收费的VPN的。

它与另一个VPN服务非常类似 – Hotspot Shield,但是SecurityKISS 目前并没有任何广告,这是最大的优点。

SecurityKISS每天的流量限制是50M,如果是Happy Days每天有300M流量。Happy Days应该是周末和节假日吧,具体也不清楚,不过软件页面都有显示可以使用的流量 。

菊子曰 我用菊子曰写博客,你呢?

周日没事,拍个视频玩玩,大家看我拍的肿么样啊


慢速摄影,用偶的卡片机拍滴。 
ps:可以调一下画质啊,偶上传的是720P的。

菊子曰 我用菊子曰写博客,你呢?

ANSYS V12安装问题1:Invalid license file:the FLexlm ID in the file does not mat …


在安装ansys V12的过程中,碰到了这个问题。在最后加载license时,弹出Invalid license file:the FLexlm ID in the file does not match the system Flexlm ID. 出现这个问题是因为ansys自动找到的网卡地址,与实际的网卡地址不符,一般是因为电脑禁止了网卡所致。

你可以‘运行’——cmd——ipconfig/all 看能不能找到网卡ID,也就是常说的MAC地址,如果找不到,就肯定是这个原因了。

解决方法:

右击“我的电脑”——选择“属性”——选择“硬件”——“设备管理器”——查看网卡是否有问题(打叉或者问号)——右击启用网卡,即可。

再重新生成license文件,放在指定文件夹下,重新安装加载license,问题就解决了。


PS:如果上述方法不能解决的话,也可能是在安装时没有关闭笔记本的无线网卡导致。

菊子曰 我用菊子曰写博客,你呢?

云 Alpha 1

优酷上的视频还是不清楚,清晰版的传到115网盘上了,需要的可以下。
清晰版:下载
菊子曰 我用菊子曰写博客,你呢?

ANSYS/LS-DYNA模拟冲压、锻压和铸造ANSYS/LS-DYNA模拟冲压、锻压和铸造 …

下面仅就LS-DYNA在模拟冲压、锻压和铸造等工艺过程的功能和特色进行说明: 
1. 冲压 薄板冲压过程的物理描述是:在模具各部件(通常是凸模、凹模和压料板)的共同作用下,板料发生大变形,板料成形的变形能来自强迫模具部件运动外功,而能量的传递完全靠模具与板料的接触和摩擦。由此可见,对于成形过程的模拟,软件的接触(contact)算法的理论和精度决定程序的可靠性,除此之外,由于板料的位移和变形很大,用来模拟板料的单元类型应满足这一要求。进行一定的假设:模具为刚体,模具的运动可直接作为冲压系统的位移边界条件。将冲压过程的物理模型转化为力学模型,即动量方程、边界条件、初始条件。可描述为:在给定的模具位移条件下,求得板料的位移函数,并在任意时刻同时满足动量方程、边界条件和初始条件。这已经是一般性的力学问题,可采用有限元的方法进行求解。 LS-DYNA在分析冲压时模具定义为刚体,因此板料和模具都应用壳单元进行离散。LS-DYNA的单元都采用Lagrange增量方法进行描述。其壳单元算法共有16种,可用于板成形分析使用的单元有11中,可分类为四节点和三节点单元;单点积分、多点积分单元和缩减积分(select-reduced)单元。单元采用co-rotational坐标系统分离单元运动中的变形和刚体运动。使用单点积分的求解速度很快,一般都可得到可靠的结果。当单元的翘曲和弯曲变形较大时,可通过增加沿壳厚度方向的积分点数目保证精度。用于板料成形的材料模式是各种弹塑性材料,可考虑各向异性、强化特征。强化类型包括指数强化、随动强化、等向强化、混合强化以及应变率对材料强化的影响。应变率的影响归结为两种方式,1.采用Cowper-Symonds模型;2.以表格方式给定任意应变率下的应力-应变曲线。部分材料模式引用Hill或Barlat的各向异性屈服假设,并假定壳单元的平面应力状态,因此几乎专用于板成形模拟。并且还能够通过给定材料的FLD(flow limit dia.)判断板料在拉延过程中局部开裂现象。 LS-DYNA目前的接触类型有30余种,适于板成形分析的有12种,都采用罚函数方法(penalty),在接触计算过程中考虑壳单元厚度及其变化。值得说明的是:1.拉延筋与板料接触(contact-drawbead),可认为是非线性弹簧算法,需给定单位长度拉延筋的对板料的阻力变化曲线。2.LS-DYNA 新增加三种接触类型(forming类型接触)专用于板成形模拟,这些接触类型降低了对模具网格的连续性要求,并且计算速度更快。 LS-DYNA进行板成形分析时可选择使用3D adaptive mesh功能,可在计算过程中对板料网格进行局部加密,网格加密的准则可选择为:1.板厚变化;2.曲率变化;3.单步长接触穿透深度值。 
2. 锻压 锻压过程是金属体积成形过程,与板成形相比,其物理描述和力学模型中相同,但单元、材料、模具定义不同。在锻压过程中往往考虑模具的变形,单元采用实体单元,材料在多数情况下经历较大的温度变化,为热塑性材料。LS-DYNA的实体单元可分为三大类:1.结构单元;2.ALE单元(包括Euler流体单元);3.声单元。进行锻压分析时要采用结构实体单元,这些单元可分为单点积分、多点积分和缩减积分(select-reduced)单元;节点带旋转自由度(nodal rotations)和不带旋转自由度单元。单元采用co-rotational坐标系统分离单元运动中的变形和刚体运动,并在应力更新中采用Jaumann应力率,避免因刚体运动产生应力。在剪切变形较大时,可选择使用Green-Naphdi应力率。变形结构单元为8节点6面体,可退化为6节点5面体或4节点4面体。 LS-DYNA的热塑性材料通过列表给定不同温度下的材料性质,例如常用的一种各向同性热塑性材料可将整个温度范围分成7段,每个温度段内可定义不同弹性模量、泊松比、屈服应力、硬化模量、热膨胀系数等参数,这种材料采用线性硬化模式。材料的热性能(比热、导热系数等)可为各向同性或各向异性。 在LS-DYNA中结构材料和热材料的定义是分开的,并且在接触传热分析中定义相关热接触界面,因此可进行结构和热场的耦合分析。 在多数锻压分析中,随着金属件成形过程的继续,初始网格的变形逐渐加大,将导致单元精度降低甚至发生畸变,因此必须使用网格重新划分功能(remeshing)。网格重划分包括以下几个步骤:1.检查网格的变形程度,若超过规定的变形度停止计算,保存结果;2.检查需要改变位置的节点,调整节点位置,保证材料边界不变,材料内部节点可自由移动。3.将保存的结果映射到新的网格上。4.重新对网格初始化并进行计算。LS-DYNA对于二维与三维网格,皆提供重划分网格的功能。另外,LS-DYNA早已采用一种更为先进的网格ALE,即任意拉格朗日-欧拉网格。ALE网格进行Rezoning的目的和过程与Remeshing基本相同,但两者的网格描述存在本质差异(后者是拉格朗日网格)。ALE结合拉格朗日和欧拉网格各自的优点,已广泛用于结构材料的极度变形。有关ALE的技术在下面详细说明。 
3.浇注 前面已经提到,结构单元运动描述采用Lagrange方法,这是因为Lagrange描述中始终以初始构形为求解的参考构形,由材料点(material point,在Total Lagrange中是初始构形的X0 ,在Updated Lagrange描述中参考构形是上一个积分步的构形,即X n-1)来确定动量方程、运动-应变关系、应变-应力关系。由此可见,任一单元的积分点在整个过程当中可以保持不变,即为同一材料点,这对于求解历史相关的变形问题是极为重要的,因为对于固体结构材料而言,正是如此。对流体介质,LS-DYNA采用Euler描述,即以当前构形(通常记为spatial point x),来确定动量方程、变形-应变关系、应变-应力关系,因此不同时步单元积分点不是相同的材料点,即物质可以在Euler网格间输运,由物质的运动导致压力和能量在Euler区传播。 Lagrange和Euler 是对物质运动的两种表述,这两种方法本质上是一样的,但由于采用的自变量(其自变量分别为X,t和x,t)不同,各自具备特点。在形式上,前者的网格节点与材料点在物质运动过程中一一对应;而后者节点不动,材料点在Euler网格中移动。在前者,任意材料点的历程都可以得到,在后者,只能得到t时刻Euler节点处材料点的特性,在下一时刻是另外一个材料点的特性,而处于此节点处的材料从何处来到那里去难以确定。此外,在物理描述上,Lagrange和Euler在确定动量方程、质量方程、和能量方程时存在较大的不同,通常Euler方程采用保守形式,而Lagrange方程往往应用工程假设采用简化格式,这一点在质量方程的表达上尤为明显。 
总之,Lagrange和 Euler是对连续介质运动的两种描述,由于参考构形的选择不同(或者说是观测者位置不同),导致对物质运动的观测和描述的侧重点存在差异。 ALE方法最早是由Noh(1964)以耦合欧拉-拉格朗日的术语提出的,至80年代末90年代初才形成成熟理论并在少数分析程序中出现。在ALE描述中,网格点可以随物质点一起运动,但也可以在空间中固定不动,甚至网格点可以在一个方向上固定,而在另一个方向上随物体一起运动。ALE中,有限单元的剖分是对参考构形进行的,网格点就是参考点,网格是独立于物体和空间运动的,亦即参考构形是已知的,初始构形和现时构形是待求解的。 由于任意拉格朗日-欧拉(ALE)方法综合了纯拉格朗日和纯欧拉描述的优点,克服了各自的缺点,成为目前非线性连续介质力学中大变形分析的十分先进有效的方法。早在91年,DYNA程序中就成功地引入ALE算法,在流体动力学、流体-结构相互作用、加工成型、碰撞、爆炸冲击、接触等大变形问题中得到了广泛的应用,如海啸、坝的决口、容器中流体的大幅度晃动和液体泄露、液体中高压气泡的扩展、水下爆炸、超高速碰撞、成型装药、鸟撞飞机、锻压等等。 ANSYS/LS-DYNA的算法除拉格朗日和ALE外,还包括欧拉和多物质流体求解。欧拉构形主要有三种:一阶精度的Donor Cell;二阶精度的Van Leer;二阶精度的Van Leer +Half Index Shift。多物质流体的单元构形主要有二种:流体+空材料和全空材料;多种材料的混合单元(压力平衡)。 这些模型都可以和通用的固体结构单元如solid、shell、brick和beam等单元自动耦合,不需要滑移界面。同时,此类求解器的加入,使ANSYS/LS-DYNA具有了可压缩流体流动分析的能力,可求解如自由界面流动、波浪破碎、任意管道流动、流体混合、复合材料等的注塑成型、金属构件浇注成型、高速高压气体注入等复杂的流体和流体-结构耦合问题。 LS-DYNA在进行浇注模拟时,模具的空腔定义为Euler区,并将其材料定义成空(void)或任何物质(如空气),浇口处单元定义为Euler源(Euler ambient),即物质由此进入Euler区,物质运动的动力是压力和(或)重力。 LS-DYNA的流体介质定义为流体动力材料,其性质主要包括密度和粘性,单元的压力以及可压缩性由附带的状态方程决定(状态方程即压力方程,其自变量包括密度、温度、内能)。 随着物质由浇口流入Euler区,空腔和浇口的压力差逐渐降低,最终达到平衡,模拟即可终止。 在浇注分析中可考虑热扩散,LS-DYNA中可方便施加温度边界条件和热生成。 总之,LS-DYNA时间积分器采用中心差分格式,对未知量显式求解。由于质量矩阵进行对角化处理,可进一步加快求解速度。例如一般的冲压、锻压、铸造等问题合理控制有限元规模,在PC机上运行5-20小时能得到理想结果,这样的效率是其它程序难以相比的。
菊子曰 我用菊子曰写博客,你呢?

奇异性,要死的刚度奇异性!!

下述条件会引起求解过程出现奇异性:(The following conditions may cause singularities in the solution process: )
·约束条件不足
·模型中有非线性单元:如间隙元、滑动元、铰链元、缆束员等。结构的一部分可能已经塌陷或分散了
·材料特性为负:如在瞬态热分析中规定的密度或温度
·连接点无约束,单元排列可能会引起奇异性。例如:两个水平梁单元在连接点的垂直方向存在无约束自由度,在线性分析中,将会忽略加在该连接点的垂直载荷。另外,考虑一个与梁单元或管单元垂直相连的无板面内的旋转刚度的壳单元,在连接点处不存在板面内的旋转刚度。在线性分析中,将会忽略加在该点处的板面内力矩。
·屈曲。当应力刚化效果为负(压缩)时,结构受载后变弱。若结构变弱到刚度减小到零或为负值,就会出现奇异解,且结构已经屈曲。会打印出“主元值为负”的消息。
~~~~~~~~~~~~~~~~~~~~~~~~
非分散缝隙又是什么!!!

菊子曰 我用菊子曰写博客,你呢?

无觅相关文章插件,快速提升流量