CN104091055B - 一种在欧拉坐标系下计算二维理想弹塑性固体的技术 - Google Patents
一种在欧拉坐标系下计算二维理想弹塑性固体的技术 Download PDFInfo
- Publication number
- CN104091055B CN104091055B CN201410300869.8A CN201410300869A CN104091055B CN 104091055 B CN104091055 B CN 104091055B CN 201410300869 A CN201410300869 A CN 201410300869A CN 104091055 B CN104091055 B CN 104091055B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- plastic solid
- solid
- ideal elastic
- elastic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Toys (AREA)
Abstract
本发明提出了一种在欧拉坐标系下计算二维理想弹塑性固体的技术,其发明的主要内容为一套完整的在欧拉坐标系下利用二维理想弹塑性固体的模型计算其相关物理量的技术。本发明的创新点主要体现在二维Hooke定律中的物质导数在欧拉坐标系下的计算方式。与一维理想弹塑性固体的计算不同的是,二维情况下模型方程均发生变化,导致其计算更加困难、计算过程更加复杂。本发明的提出,可直接用来计算二维理想弹塑性固体,并在二维理想弹塑性固体受外力作用、与其他介质耦合作用等实际工程应用中具有重要意义。
Description
技术领域
本发明涉及一种计算二维理想弹塑性固体的技术,具体涉及一种在欧拉坐标系下计算二维理想弹塑性固体的技术。
背景技术
二维理想弹塑性固体模型可以较准确描述固体(如铝、钢等金属)在受一般强度的外力作用下的各物理量变化情况。因此,研究理想弹塑性固体的计算技术,在实际工程中具有重要的应用价值和广泛的应用前景。
目前,虽然已经存在一些计算理想弹塑性固体的技术,但均与本发明提出的技术不同。比如,M.L.Wilkins在1964年在提出理想弹塑性固体的模型之后,采用有限差分法对该模型进行求解,其中使用了复杂的全离散形式。再如,B.P.Howell在2000年采用FreeLagrange方法对理想弹塑性固体进行计算。该方法是在拉格朗日坐标下进行计算的,虽然在计算一些变量(如偏应力)上可以得到简化,但推广到高维时变得十分复杂。为了使理想弹塑性固体的计算既简单有准确,本发明直接在欧拉坐标系下进行计算,只需将Hooke定律中的导数作为物质导数进行处理。值得一提的是,本发明是在1993年M.B.Tyndall 的工作上受到启发而提出的。然而,M.B.Tyndall的计算方法中却存在一些错误。首先,在计算欧拉坐标下每个网格点在上一时间步的位置时,他在固定的网格点上对速度采取了时间平均。这种计算方法在拉格朗日坐标下是正确的,在欧拉坐标下却是错误的,与他建立的欧拉坐标下的控制方程相矛盾。其次,他在计算相关物理量在每个网格点上一时间步处的值的时候,采用了该点两侧网格点的抛物插值(二次函数插值)。可是,理想弹塑性固体的控制方程是双曲方程,传播具有方向性,采用抛物插值会造成计算的不准确,甚至造成不稳定并产生错误。实际的数值计算也验证了他的方法确实存在一些错误。对于这一问题,本发明则直接采取迎风线性插值。总之,本发明提出的一维理想弹塑性固体的计算技术兼顾了方法的简单性与正确性。
同时,与一维理想弹塑性固体的计算不同的是,二维情况下模型方程均发生变化,导致其计算更加困难、计算过程更加复杂。这些都需要重新进行考虑。
发明内容
本发明提出的的计算二维理想弹塑性固体的技术,其发明内容主要体现在在欧拉坐标系下的一套完整的计算二维理想弹塑性固体的技术,其创新点主要体现在二维Hooke定律中的物质导数在欧拉坐标系下的计算方式。
对于二维情况,理想弹塑性固体在欧拉坐标系下的控制方程为
此处,ρ是密度,u是x方向速度,v是y方向速度,p是压力,E是总能,σx是x方向的总应力,σy是y方向的总应力,sxy是剪切偏应力。此外,对于理想弹塑性固体,其总应力和压力还满足下面关系:
σx=-p+sx,σy=-p+sy
其中,sx是x方向的偏应力,sy是y方向的偏应力。当理想弹塑性固体处于弹性状态,有
其中K是体积模量,μ是剪切模量,εx、εy、εxy分别是相对应的应变,并且有
当理想弹塑性固体处于塑性状态,有
其中c0,ρ0,γs均为与具体固体有关的常数,Y0是屈服强度。当理想弹塑性固体满足下面方程时为弹性状态
sx 2+sy 2+2sxy 2≤2Y0 2/3
当上述不等式不成立时,固体处于塑性状态。
本发明的具体发明内容可以归结为如下计算技术。假设已知一维理想弹塑性固体在第n个时间步的各变量值需要将这些变量值推进到第n+1个时间步,得到 其计算技术通过以下六个步骤来实现:
1.求解控制方程(1),将控制方程中第n个时间步的各变量值更新至第n+1个时间步,得到
2.计算欧拉坐标下每个网格点在第n个时间步的位置,记为(xold,yold),有
3.采用迎风线性插值,计算ρ,p,sx,sy,sxy在(xold,yold)处的值,记作ρold,pold,sxold,syold, sxyold如
其中
pold,sxold,syold,sxyold也可通过类似方式进行计算。
4.利用Hooke定律及迎风线性插值,得到初步的
类似地,也可得到
5.通过von Mises屈服条件判断在每一欧拉网格点处的弹塑性状态,并更新压力值至如果某一网格点处满足vonMises屈服条件,即sx 2+sy 2+2sxy 2≤2Y0 2/3,则固体处于弹性状态,并且压力通过Hooke定律计算
若在该网格点不满足vonMises屈服条件,即sx 2+sy 2+2sxy 2>2Y0 2/3,则固体处于塑形状态,并且压力通过状态方程计算
同时,使偏应力满足理想塑性条件
6.返回步骤1直到达到设定的时间迭代要求。
附图说明
图1是本发明在欧拉坐标系下计算二维理想弹塑性固体的流程图;
图2至图4是本发明在欧拉坐标系下计算二维理想弹塑性固体的一个算例结果。
具体实施方式
为了说明本发明的具体实施方式,下面将对一个算例进行演示。考虑在二维平面中的问题,一个半无限长的铝棒冲击一个半无限大的铝块,其中左侧铝棒的无量纲初值为uL=20.0,pL=0.0,ρL=2.7, sxL=0.0,syL=0.0,sxyL=0.0。右侧铝块的无量纲初值为 uR=0.0,pR=1.0,ρR=2.7,sxR=0.0,syR=0.0,sxyR=0.0。该问题的无量纲的求解区域为 x×y∈[-0.02,0.02]×[-0.02,0.02],其中铝棒的初始区域为x×y∈[-0.02,0.00]×[-0.006,0.006],铝块的初始区域为x×y∈[0.00,0.02]×[-0.02,0.02]。同时,铝的理想弹塑性模型的相关无量纲参数分别为ρ0=2.71,c0=538.0,γs=2.71,K=740000.0,μ=265000.0,Y0=3000.0。
该问题将在界面左右两侧同时产生弹性波和塑性波。取时间步长Δt=0.0000001,x方向和y方向的空间步长均为0.00001。采用二维Lax-Friedrich格式计算,得到在时间分别为0.00005,0.0001, 0.00015时铝中的负的x方向总应力如图2至图4所示。
Claims (2)
1.一种在欧拉坐标系下计算二维理想弹塑性固体的方法,其特征在于,该方法是一套完整的在欧拉坐标系下计算二维理想弹塑性固体相关物理量的方法,所述相关物理量包括密度速度压力偏应力所述二维理想弹塑性固体为铝棒和铝块;
Hooke定律中的导数在欧拉坐标系下需要转换成物质导数计算,具体离散形式表示为
及
其中,ρold、pold、sxold分别表示欧拉坐标系下的网格点在上一时间步位置的密度值、压力值、偏应力值,它们均通过类似如下形式的迎风线性插值得到
其中
pold,sxold,syold,sxyold也通过类似方式进行计算,同时,可按上述的计算过程计算出
2.如权利要求1所述的在欧拉坐标系下计算二维理想弹塑性固体的方法,其特征在于,当其处于弹性状态,满足Hooke定律
其中K是体积模量,μ是剪切模量,εx、εy、εxy分别是相对应的应变;当其处于塑性状态,满足关系
其中c0,ρ0,γs均为与具体固体有关的常数,Y0是屈服强度,并且该固体的屈服条件为von Mises屈服条件,即当其偏应力满足
sx 2+sy 2+2sxy 2≤2Y0 2/3
时,固体处于弹性状态,当上述不等式不成立时,固体处于理想塑性状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300869.8A CN104091055B (zh) | 2014-06-27 | 2014-06-27 | 一种在欧拉坐标系下计算二维理想弹塑性固体的技术 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300869.8A CN104091055B (zh) | 2014-06-27 | 2014-06-27 | 一种在欧拉坐标系下计算二维理想弹塑性固体的技术 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104091055A CN104091055A (zh) | 2014-10-08 |
CN104091055B true CN104091055B (zh) | 2018-01-02 |
Family
ID=51638771
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410300869.8A Active CN104091055B (zh) | 2014-06-27 | 2014-06-27 | 一种在欧拉坐标系下计算二维理想弹塑性固体的技术 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104091055B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108073731B (zh) * | 2016-11-10 | 2021-02-19 | 中国石油化工股份有限公司 | 一种地震波数值模拟的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1291225C (zh) * | 2004-08-13 | 2006-12-20 | 大庆油田有限责任公司 | 应力性套损预测方法 |
CN102819633A (zh) * | 2012-07-27 | 2012-12-12 | 哈尔滨工业大学 | 焊接热循环温度与热变形历史材料本构关系建立方法及msc.marc二次开发 |
-
2014
- 2014-06-27 CN CN201410300869.8A patent/CN104091055B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104091055A (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018142394A3 (en) | Computer aided driving | |
Farhangi et al. | Classic grounded theory; description of generating gravity center of Iranian health tourism brand identity theory | |
CN105608708A (zh) | 基于分块平面拟合的图像二值化方法及装置 | |
CN104182585A (zh) | 硅胶弹性体手感分析及疲劳寿命预测的有限元模拟方法 | |
MX2019003688A (es) | Dispositivo de determinacion de fracturas, programa de determinacion de fracturas y metodo de la misma. | |
CN104091055B (zh) | 一种在欧拉坐标系下计算二维理想弹塑性固体的技术 | |
Vorkov et al. | The multi-breakage phenomenon in air bending process | |
WO2016027056A3 (en) | Deriving contact stress or contact load using ultrasound data | |
CN104036150B (zh) | 一种在欧拉坐标系下计算一维理想弹塑性固体的技术 | |
US20190065646A1 (en) | Methods And Systems For Conducting A Time-Marching Numerical Simulation Of A Deep Drawing Metal Forming Process For Manufacturing A Product or Part | |
Soheili et al. | A moving mesh method with variable mesh relaxation time | |
CN102819237A (zh) | 一种实体法铣削仿真过程中仿真目标毛坯的产生方法 | |
Li et al. | Numerical simulation method for ship-bridge collision considering fluid effect | |
CN106697328A (zh) | 一种飞行器薄壁结构传载特性同材料模型试验设计方法 | |
Ali et al. | Prediction of high cutting speed parameters for Ti-6Al-4V by using Finite Element Modeling | |
Farias et al. | Grid and time discretization issues affecting the application of the generalized material point method (GIMP) to simulate wedge penetration in soft soil | |
Bermudo Gamboa et al. | Selection of the optimal distribution for the upper bound theorem in indentation processes | |
蒋磊 et al. | Research on full process forming simulation of side outer panel based on die surface engineering | |
Walters et al. | Consequences of using the plane stress assumption for damage calculations in crash analyses | |
Dal Santo et al. | Cold forming by stretching of aeronautic sheet metal parts | |
Guo et al. | Explicit finite element simulation of oblique cutting process | |
BR112022006137A2 (pt) | Método implementado por computador, dispositivo eletrônico e programa de computador | |
Chiorescu et al. | The Influence of Dimensional-Type Mathematical Model for Cylindrical Deep Drawing on the Distribution of Extreme Unitary Stresses | |
Jajimoggala | Selection of Optimal Hot Extrusion Process Parameters for AA6061-Fly Ash Composites: Using Simulation | |
Zhang et al. | The Research of Center Shaft Bracket Stamping Simulation Technology Based on Dynaform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210618 Address after: Room 501, 5 / F, building 4, 1418-66 Moganshan Road, Hangzhou, Zhejiang 310000 Patentee after: Sanduo (Hangzhou) Technology Co.,Ltd. Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road Patentee before: BEIHANG University |
|
TR01 | Transfer of patent right |