CN113343521A - 基于comsol预测选区激光熔化过程中层间热应力分布的方法 - Google Patents
基于comsol预测选区激光熔化过程中层间热应力分布的方法 Download PDFInfo
- Publication number
- CN113343521A CN113343521A CN202110584593.0A CN202110584593A CN113343521A CN 113343521 A CN113343521 A CN 113343521A CN 202110584593 A CN202110584593 A CN 202110584593A CN 113343521 A CN113343521 A CN 113343521A
- Authority
- CN
- China
- Prior art keywords
- heat source
- focus
- thermal stress
- powder
- gaussian
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22F—WORKING METALLIC POWDER; MANUFACTURE OF ARTICLES FROM METALLIC POWDER; MAKING METALLIC POWDER; APPARATUS OR DEVICES SPECIALLY ADAPTED FOR METALLIC POWDER
- B22F10/00—Additive manufacturing of workpieces or articles from metallic powder
- B22F10/20—Direct sintering or melting
- B22F10/28—Powder bed fusion, e.g. selective laser melting [SLM] or electron beam melting [EBM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22F—WORKING METALLIC POWDER; MANUFACTURE OF ARTICLES FROM METALLIC POWDER; MAKING METALLIC POWDER; APPARATUS OR DEVICES SPECIALLY ADAPTED FOR METALLIC POWDER
- B22F10/00—Additive manufacturing of workpieces or articles from metallic powder
- B22F10/30—Process control
- B22F10/37—Process control of powder bed aspects, e.g. density
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22F—WORKING METALLIC POWDER; MANUFACTURE OF ARTICLES FROM METALLIC POWDER; MAKING METALLIC POWDER; APPARATUS OR DEVICES SPECIALLY ADAPTED FOR METALLIC POWDER
- B22F10/00—Additive manufacturing of workpieces or articles from metallic powder
- B22F10/80—Data acquisition or data processing
- B22F10/85—Data acquisition or data processing for controlling or regulating additive manufacturing processes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y10/00—Processes of additive manufacturing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y40/00—Auxiliary operations or equipment, e.g. for material handling
- B33Y40/10—Pre-treatment
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y50/00—Data acquisition or data processing for additive manufacturing
- B33Y50/02—Data acquisition or data processing for additive manufacturing for controlling or regulating additive manufacturing processes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/25—Process efficiency
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Materials Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Automation & Control Theory (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Plasma & Fusion (AREA)
- Powder Metallurgy (AREA)
Abstract
本申请提供一种基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,所述方法包括如下步骤:S1:基于COMSOL构建三维固体传热和结构力学瞬态模型;S2:确定仿真过程中的参数;S3:确定待熔化的粉末的材料属性;S4:确定移动高斯热源参数;S5:构建粉床的几何模型;S6:实现选区激光熔化的逐层制造;S7:网格划分并计算节点温度;S8:根据步骤S7的结果预测层间热应力分布和残余热应力分布。本方法通过移动的高斯热源模,模拟加工过程中的激光热源效果,并利用均匀的材质粉床来代替粉末床;此外,利用结构力学模块模拟随着热源移动时,该层产生的热应力以及制件变形情况,模拟选区激光熔化技术逐层制造的过程,实现多层间热应力、残余热应力的预测。
Description
技术领域
本发明涉及热应力分布技术领域,尤其涉及一种基于COMSOL预测选区激光熔化技术中层间热应力分布的方法。
背景技术
选区激光熔化技术是基于材料离散-逐渐累加方式制造实体零件的近净成形技术。该技术通常以金属粉末为原料,通过三维模型预分层处理设定激光扫描路径,采用高能量激光束按照设定的扫描路径逐层熔化金属粉末,使其快速凝固、堆积而形成高性能构件。在激光熔化技术过程中,金属材料经历快速加热、凝固和冷却过程,在此过程中会形成较大的热应力以及由固态相变引起的组织应力。这些应力在成形结束后残留在工件内部而成为残余应力。若残余应力超过材料本身的屈服强度,成形件就会发生变形,导致尺寸精度及使用性能的降低。因此,增材制造金属零部件的变形一直是国内外增材制造领域的研究热点之一。现有的预测激光熔化技术制造过程中层间热应力的方法有体积热源方法,即每一层生成一个切片实体,并且带有体积热源,以此来预测层间的热应力的分布,如现有的ANSYS新增的增材制造板块。但是以这种体积热源加载到每一层切片层上,只能预测层与层间的热应力分布,并不能准确预测在每一层上单道扫描累计出现热应力的分布。
因此,亟需一种能准确预测在每一层上单道扫描累计出现热应力的方法。
发明内容
有鉴于此,本发明提供一种基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:所述方法包括如下步骤:
S1:基于COMSOL构建三维固体传热和结构力学瞬态模型;
S2:确定仿真过程中的参数,所述参数包括扫描间隔D_spot、激光扫描速度v_spot、激光功率P_laser、激光半径r_spot、表面辐射率A_Gass、粉末堆积率w_powder和粉末层厚thickness;
S3:确定待熔化的粉末的材料属性,所述材料属性包括热导系数、比热容、材料密度、热膨胀系数、杨氏模量以及泊松比;
S4:确定移动高斯热源参数,所述移动高斯热源参数包括移动高斯热源数量,以及各个移动高斯热源使用与时间相关的差值函数;
S5:构建粉床的几何模型,并确定粉床的几何模型的初始条件、边界热源条件和确定固体力学节点下的边界条件;
S6:实现选区激光熔化的逐层制造;
S7:对粉床的几何模型进行网格划分并确定节点的温度、应力和应变;通过网格划分通过能量守恒和应力平衡确定节点的温度、应力和应变。
S8:根据步骤S7的结果预测层间热应力分布和残余热应力分布。
进一步,步骤S3还包括将待熔化的粉末材料在粉床上堆积的形状进行预处理,即对粉末松装密度进行预处理,所述预处理用于将粉末堆积的粉床近似为一个长方体。
进一步,所述粉末松装密度的取值范围为40%-60%。
进一步,所述粉末松装密度为50%。
进一步,步骤S4中移动高斯热源数量根据切片形状、扫描间距和激光半径决定,使其一次平铺扫描能覆盖目标切片形状。
进一步,步骤S4中的差值函数采用如下方法确定:
x_focus=x_f1(t) (1)
其中,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,x_f1(t)表示使高斯热源移动的x方向的函数;
y_focus=y_f1(t) (2)
其中,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点,y_f1(t)表示使高斯热源移动的y方向的函数,;
r_focus=sqrt((x-x_focus)∧2+(y-y_focus)∧2) (3)
其中,r_focus表示高斯热源的焦点,所述焦点由x向和y向高斯热源的焦点决定,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点;
Flux=((2*A_Gass*P_laser)/(pi*r_spot∧2))*exp(-2*r_focus∧2)/r_spot∧2)(4)
其中,Flux表示高斯热源的热通量,A_Gass表示材料的吸收率,P_laser表示激光热源,pi表示π,r_spot表示激光半径,r_focus表示高斯热源的焦点。
本发明的有益技术效果:本方法通过移动的高斯热源模,考虑材料的相变,材料表面的热辐射与自然对流对于温度场的影响,模拟加工过程中的激光热源效果,并通过对材料属性进行近似处理,利用均匀的材质粉床来代替粉末床;此外,利用结构力学模块模拟随着热源移动时,该层产生的热应力以及制件变形情况,最后利用激活的方法,模拟选区激光熔化技术逐层制造的过程,实现多层间热应力、残余热应力的预测,以及后层的激光作用效果对于前层的影响。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1为本申请的流程图。
图2为本申请的网格剖分图。
图3为本申请的温度场和应力场示意图。
具体实施方式
以下结合说明书附图对本发明做出进一步的说明:
本发明提供一种基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:所述方法包括如下步骤:如图1所示,
S1:基于COMSOL构建三维固体传热和结构力学瞬态模型;
S2:确定仿真过程中的参数,所述参数包括扫描间隔D_spot、激光扫描速度v_spot、激光功率P_laser、激光半径r_spot、表面辐射率A_Gass、粉末堆积率w_powder和粉末层厚thickness;
S3:确定待熔化的粉末的材料属性,所述材料属性包括热导系数、比热容、材料密度、热膨胀系数、杨氏模量以及泊松比;
S4:确定移动高斯热源参数,所述移动高斯热源参数包括移动高斯热源数量,以及各个移动高斯热源使用与时间相关的差值函数;
S5:构建粉床的几何模型,并确定粉床的几何模型的初始条件、边界热源条件和确定固体力学节点下的边界条件;
所述几何模型的初始条件、边界热源条件和确定固体力学节点下的边界条件如下:
固体传热
1.为模拟基板预热效果,所有的初始温度定义为35degC,模型所有表面定义传热系数为80W/(m2·K),对流热通量的边界条件模拟加工过程中风场对于加工的影响。
2.其次,激光作用的材料表面定义为表面对环境辐射的边界条件,模拟在移动的高斯热源作用下,高温材料表面与环境的热交换过程。
3.接着,对于每一层材料定义一个与时间有关的高斯热源,以模拟热源的加载。
固体力学
1.对基板底面施加一个固定约束的边界条件,防止在此方向的变形。
2.对于受激光作用的材料添加热膨胀节点。
3.为了实现逐层制造,使用与时间相关的活化表达式,活化比例因子设为1e-5。
S6:实现选区激光熔化的逐层制造;每一层使用“zigzag”的扫描策略,共三道,来模拟激光的连续加工,其中每一层激光作用的总时长为t_x=0.00116s,并考虑0.2*t_x的冷却时间,而后激活第下一层粉床。
S7:对粉床的几何模型进行网格划分并确定节点的温度、应力和应变;通过网格划分通过能量守恒和应力平衡确定节点的温度、应力和应变。本申请的方法是一种建模,固对于模型进行了简单的网格划分,以实现相应的效果。下图为网格剖分图,以及逐层制造的顺序。如图2所示。
S8:根据步骤S7的结果预测层间热应力分布和残余热应力分布。如图3所示。步骤S8利用COMSOL软件的“后处理”的功能,根据目标物理量来表征结果。步骤S7将利用有限元对粉床进行分割,即化整为零,并确定每个节点的温度、应力和应变,步骤S8将各个节点的计算结果集零为整,从粉床的整体角度通过软件后处理功能获得粉床的温度、应力和应变,从而实现层间热应力和残余热应力的预测。如以应力为例,我们在结果节点下建立三维应力查看节点,输入表达式即可得到想要的结果,同理可查看温度和应变大小,在此不再赘述。
上述技术方案通过移动的高斯热源模,考虑材料的相变,材料表面的热辐射与自然对流对于温度场的影响,模拟加工过程中的激光热源效果,并通过对材料属性进行近似处理,利用均匀的材质粉床来代替粉末床;此外,利用结构力学模块模拟随着热源移动时,该层产生的热应力以及制件变形情况,最后利用激活的方法,模拟选区激光熔化技术逐层制造的过程,实现多层间热应力、残余热应力的预测,以及后层的激光作用效果对于前层的影响。
在本实施例中,步骤S3还包括将待熔化的粉末材料在粉床上堆积的形状进行预处理,即对粉末松装密度进行预处理,所述预处理用于将粉末堆积的粉床近似为一个长方体。因为本申请的方法用于预测层间热应力,不能将熔化过程这一物理现象表征出来,故本申请利用熔化后的较为规则的粉床来预测热应力,将粉床视为规则的长方体是为了近似处理和降低计算难度。
在本实施例中,所述粉末松装密度的取值范围为40%-60%。所述粉末松装密度为50%。在粉床建模过程中,成百上千的粉球堆积在粉床上,粉末球与粉末球之间存在孔隙,松装密度反应了粉床的孔隙率,在本实施例中取松装密度为50%,用以仿真模拟粉末床堆积效果。
在本实施例中,步骤S4中移动高斯热源数量根据切片形状、扫描间距和激光半径决定,使其一次平铺扫描能覆盖目标切片形状。如用一个算例展示高斯热源数量的确定:每一层粉床按照规则模型建模,高斯热源扫描的长度和次数取决于每一层切片形状,本申请的切片为一个规则长方形,在预设的扫描间距、激光半径下,三次扫描即可扫描出本案例中规则的长方体。本领域技术人员可根据实际切片、扫描间距和激光半径来确定高斯热源数量。
在本实施例中,步骤S4中的差值函数采用如下方法确定:
x_focus=x_f1(t) (1)
其中,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,x_f1(t)表示使高斯热源移动的x方向的函数;
y_focus=y_f1(t) (2)
其中,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点,y_f1(t)表示使高斯热源移动的y方向的函数,;
r_focus=sqrt((x-x_focus)∧2+(y-y_focus)∧2) (3)
其中,r_focus表示高斯热源的焦点,所述焦点由x向和y向高斯热源的焦点决定,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点;
Flux=((2*A_Gass*P_laser)/(pi*r_spot∧2))*exp(-2*r_focus∧2)/r_spot∧2)(4)
其中,Flux表示高斯热源的热通量,A_Gass表示材料的吸收率,P_laser表示激光热源,pi表示π,r_spot表示激光半径,r_focus表示高斯热源的焦点。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:所述方法包括如下步骤:
S1:基于COMSOL构建三维固体传热和结构力学瞬态模型;
S2:确定仿真过程中的参数,所述参数包括扫描间隔D_spot、激光扫描速度v_spot、激光功率P_laser、激光半径r_spot、表面辐射率A_Gass、粉末堆积率w_powder和粉末层厚thickness;
S3:确定待熔化的粉末的材料属性,所述材料属性包括热导系数、比热容、材料密度、热膨胀系数、杨氏模量以及泊松比;
S4:确定移动高斯热源参数,所述移动高斯热源参数包括移动高斯热源数量,以及各个移动高斯热源使用与时间相关的差值函数;
S5:构建粉床的几何模型,并确定粉床的几何模型的初始条件、边界热源条件和确定固体力学节点下的边界条件;
S6:实现选区激光熔化的逐层制造;
S7:对粉床的几何模型进行网格划分并确定节点的温度、应力和应变;通过网格划分通过能量守恒和应力平衡确定节点的温度、应力和应变。
S8:根据步骤S7的结果预测层间热应力分布和残余热应力分布。
2.根据权利要求1所述基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:步骤S3还包括将待熔化的粉末材料在粉床上堆积的形状进行预处理,即对粉末松装密度进行预处理,所述预处理用于将粉末堆积的粉床近似为一个长方体。
3.根据权利要求2所述基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:所述粉末松装密度的取值范围为40%-60%。
4.根据权利要求3所述基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:所述粉末松装密度为50%。
5.根据权利要求1所述基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:步骤S4中移动高斯热源数量根据切片形状、扫描间距和激光半径决定,使其一次平铺扫描能覆盖目标切片形状。
6.根据权利要求1所述基于COMSOL预测选区激光熔化过程中层间热应力分布的方法,其特征在于:步骤S4中的差值函数采用如下方法确定:
x_focus=x_f1(t) (1)
其中,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,x_f1(t)表示使高斯热源移动的x方向的函数;
y_focus=y_f1(t) (2)
其中,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点,y_f1(t)表示使高斯热源移动的y方向的函数,;
r_focus=sqrt((x-x_focus)^2+(y-y_focus)^2) (3)
其中,r_focus表示高斯热源的焦点,所述焦点由x向和y向高斯热源的焦点决定,x表示高斯热源移动的x方向,x_focus表示x向高斯热源的焦点,y表示高斯热源移动的y方向,y_focus表示y向高斯热源的焦点;
Flux=((2*A_Gass*P_laser)/(pi*r_spot^2))*exp(-2*r_focus^2)/r_spot^2) (4)
其中,Flux表示高斯热源的热通量,A_Gass表示材料的吸收率,P_laser表示激光热源,pi表示π,r_spot表示激光半径,r_focus表示高斯热源的焦点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110584593.0A CN113343521B (zh) | 2021-05-27 | 2021-05-27 | 基于comsol预测选区激光熔化过程中层间热应力分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110584593.0A CN113343521B (zh) | 2021-05-27 | 2021-05-27 | 基于comsol预测选区激光熔化过程中层间热应力分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113343521A true CN113343521A (zh) | 2021-09-03 |
CN113343521B CN113343521B (zh) | 2023-07-04 |
Family
ID=77471835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110584593.0A Active CN113343521B (zh) | 2021-05-27 | 2021-05-27 | 基于comsol预测选区激光熔化过程中层间热应力分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113343521B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779793A (zh) * | 2021-09-10 | 2021-12-10 | 电子科技大学 | 一种基于射线追踪的激光选区熔化的热源建模方法 |
CN113779906A (zh) * | 2021-11-12 | 2021-12-10 | 北京科技大学 | 镍基高温合金真空感应熔炼过程中脱氮动力学的计算方法 |
CN113976920A (zh) * | 2021-09-27 | 2022-01-28 | 上海交通大学 | 选区激光熔化成形结构残余变形的跨尺度控制方法及*** |
CN116275125A (zh) * | 2023-05-19 | 2023-06-23 | 西北工业大学 | 一种增材制造熔池形成及分布特征预测方法 |
CN116680962A (zh) * | 2023-08-01 | 2023-09-01 | 北京理工大学 | 选区激光熔化技术下成形件表面质量及残余应力预测方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107862163A (zh) * | 2017-12-08 | 2018-03-30 | 南昌航空大学 | 一种激光熔覆纳米YSZ@Ni制备YSZ热障涂层应力场模拟方法 |
CN108038266A (zh) * | 2017-11-17 | 2018-05-15 | 西安铂力特增材技术股份有限公司 | 一种选择性激光修复过程的数值模拟方法 |
CN108062432A (zh) * | 2017-11-17 | 2018-05-22 | 西安铂力特增材技术股份有限公司 | 一种激光选区熔化过程的数值模拟方法 |
CN108717481A (zh) * | 2018-05-08 | 2018-10-30 | 武汉理工大学 | 一种选择性激光熔化过程中温度分布及翘曲变形预测方法 |
CN109317772A (zh) * | 2018-11-20 | 2019-02-12 | 武汉理工大学 | 一种实验表征和数值模拟相结合探究激光钎焊工艺参数的方法 |
CN110899981A (zh) * | 2019-12-04 | 2020-03-24 | 哈尔滨工业大学 | 一种激光改性超精密切削的激光辅助加工硬脆材料方法 |
CN111090937A (zh) * | 2019-12-13 | 2020-05-01 | 北京理工大学 | 基于欧拉网格的增材制造工艺构件尺度的仿真处理方法 |
CN111627503A (zh) * | 2020-05-27 | 2020-09-04 | 燕山大学 | 一种氧化铝陶瓷基复合材料涂层激光熔覆制造过程中应力场的预测方法 |
CN112182921A (zh) * | 2020-08-31 | 2021-01-05 | 北京化工大学 | 一种高性能合金钢选区激光熔化热力耦合行为预测方法 |
CN112651157A (zh) * | 2020-12-24 | 2021-04-13 | 北京星航机电装备有限公司 | 一种航空薄壁件等离子喷涂变形控制方法 |
AU2021100843A4 (en) * | 2021-02-10 | 2021-04-22 | Nanjing Tech University | Simulation And Safety Early-Warning Method For The Structure Failure Of Steel Structure Poles And Towers Induced By Wildfire |
-
2021
- 2021-05-27 CN CN202110584593.0A patent/CN113343521B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038266A (zh) * | 2017-11-17 | 2018-05-15 | 西安铂力特增材技术股份有限公司 | 一种选择性激光修复过程的数值模拟方法 |
CN108062432A (zh) * | 2017-11-17 | 2018-05-22 | 西安铂力特增材技术股份有限公司 | 一种激光选区熔化过程的数值模拟方法 |
CN107862163A (zh) * | 2017-12-08 | 2018-03-30 | 南昌航空大学 | 一种激光熔覆纳米YSZ@Ni制备YSZ热障涂层应力场模拟方法 |
CN108717481A (zh) * | 2018-05-08 | 2018-10-30 | 武汉理工大学 | 一种选择性激光熔化过程中温度分布及翘曲变形预测方法 |
CN109317772A (zh) * | 2018-11-20 | 2019-02-12 | 武汉理工大学 | 一种实验表征和数值模拟相结合探究激光钎焊工艺参数的方法 |
CN110899981A (zh) * | 2019-12-04 | 2020-03-24 | 哈尔滨工业大学 | 一种激光改性超精密切削的激光辅助加工硬脆材料方法 |
CN111090937A (zh) * | 2019-12-13 | 2020-05-01 | 北京理工大学 | 基于欧拉网格的增材制造工艺构件尺度的仿真处理方法 |
CN111627503A (zh) * | 2020-05-27 | 2020-09-04 | 燕山大学 | 一种氧化铝陶瓷基复合材料涂层激光熔覆制造过程中应力场的预测方法 |
CN112182921A (zh) * | 2020-08-31 | 2021-01-05 | 北京化工大学 | 一种高性能合金钢选区激光熔化热力耦合行为预测方法 |
CN112651157A (zh) * | 2020-12-24 | 2021-04-13 | 北京星航机电装备有限公司 | 一种航空薄壁件等离子喷涂变形控制方法 |
AU2021100843A4 (en) * | 2021-02-10 | 2021-04-22 | Nanjing Tech University | Simulation And Safety Early-Warning Method For The Structure Failure Of Steel Structure Poles And Towers Induced By Wildfire |
Non-Patent Citations (8)
Title |
---|
A.BHARATISH ET AL.: "evaluation of thermal residual stresses in laser drilled alumina ceramics using micro-raman spectroscopy and comsol multiphysics", 《OPTICS & LASER TECHNOLOGY》 * |
A.BHARATISH ET AL.: "evaluation of thermal residual stresses in laser drilled alumina ceramics using micro-raman spectroscopy and comsol multiphysics", 《OPTICS & LASER TECHNOLOGY》, no. 70, 15 July 2015 (2015-07-15), pages 76 - 84 * |
姜沐晖: "基于铜膜中间层的聚碳酸酯激光透射焊接残余应力与工艺研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》 * |
姜沐晖: "基于铜膜中间层的聚碳酸酯激光透射焊接残余应力与工艺研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》, no. 2, 15 February 2021 (2021-02-15), pages 022 - 832 * |
孙建华: "钢/铝异种金属激光-电弧复合焊接数值模拟及实验研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》 * |
孙建华: "钢/铝异种金属激光-电弧复合焊接数值模拟及实验研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》, no. 1, 15 January 2021 (2021-01-15), pages 022 - 1082 * |
潘井春: "钢/铝异种金属激光深熔焊过程中熔池行为及机理研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》 * |
潘井春: "钢/铝异种金属激光深熔焊过程中熔池行为及机理研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》, no. 1, 15 January 2019 (2019-01-15), pages 022 - 1212 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113779793A (zh) * | 2021-09-10 | 2021-12-10 | 电子科技大学 | 一种基于射线追踪的激光选区熔化的热源建模方法 |
CN113779793B (zh) * | 2021-09-10 | 2023-09-26 | 电子科技大学 | 一种基于射线追踪的激光选区熔化的热源建模方法 |
CN113976920A (zh) * | 2021-09-27 | 2022-01-28 | 上海交通大学 | 选区激光熔化成形结构残余变形的跨尺度控制方法及*** |
CN113976920B (zh) * | 2021-09-27 | 2022-08-26 | 上海交通大学 | 选区激光熔化成形结构残余变形的跨尺度控制方法及*** |
CN113779906A (zh) * | 2021-11-12 | 2021-12-10 | 北京科技大学 | 镍基高温合金真空感应熔炼过程中脱氮动力学的计算方法 |
CN116275125A (zh) * | 2023-05-19 | 2023-06-23 | 西北工业大学 | 一种增材制造熔池形成及分布特征预测方法 |
CN116275125B (zh) * | 2023-05-19 | 2023-09-01 | 西北工业大学 | 一种增材制造熔池形成及分布特征预测方法 |
CN116680962A (zh) * | 2023-08-01 | 2023-09-01 | 北京理工大学 | 选区激光熔化技术下成形件表面质量及残余应力预测方法 |
CN116680962B (zh) * | 2023-08-01 | 2023-10-17 | 北京理工大学 | 选区激光熔化技术下成形件表面质量及残余应力预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113343521B (zh) | 2023-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113343521A (zh) | 基于comsol预测选区激光熔化过程中层间热应力分布的方法 | |
Moges et al. | A review of model inaccuracy and parameter uncertainty in laser powder bed fusion models and simulations | |
Patil et al. | A generalized feed forward dynamic adaptive mesh refinement and derefinement finite element framework for metal laser sintering—part I: formulation and algorithm development | |
Ganeriwala et al. | A coupled discrete element-finite difference model of selective laser sintering | |
Lavery et al. | A review of Computational Modelling of Additive Layer Manufacturing–multi-scale and multi-physics | |
Patil et al. | A new finite element solver using numerical eigen modes for fast simulation of additive manufacturing processes | |
JP7125764B2 (ja) | 積層造形物の解析方法及び積層造形物の解析装置、並びに積層造形物の製造方法及び積層造形物の製造装置 | |
Zeng et al. | A new dynamic mesh method applied to the simulation of selective laser melting | |
Khan et al. | Probing a novel heat source model and adaptive remeshing technique to simulate laser powder bed fusion with experimental validation | |
CN112182921A (zh) | 一种高性能合金钢选区激光熔化热力耦合行为预测方法 | |
CN111666663B (zh) | 一种slm热应力快速计算方法 | |
CN113976920B (zh) | 选区激光熔化成形结构残余变形的跨尺度控制方法及*** | |
CN115527632A (zh) | 复合材料激光制孔孔型的预测方法和装置 | |
Huo et al. | Influence of different processing parameter on distortion and residual stress of inconel 718 alloys fabricated by selective laser melting (SLM) | |
Samantaray et al. | Finite element simulation of heat transfer in laser additive manufacturing of AlSi10Mg powders | |
Wischeropp et al. | Simulation of the effect of different laser beam intensity profiles on heat distribution in selective laser melting | |
CN115964890A (zh) | 一种3d打印材料属性转换与复杂扫描策略实现方法 | |
WO2023059855A1 (en) | Intelligent scan sequence optimization for powder bed fusion additive manufacturing using linear systems theory | |
Bachy et al. | Simulation and experimental investigation for the 2D and 3D laser direct structuring process | |
Gobal et al. | An adaptive discrete element method for physical modeling of the selective laser sintering process | |
CN111168067B (zh) | 一种基于激光定向能量沉积的孔隙预测与控制方法 | |
Zhang | Nonequilibrium modeling of heat transfer in a gas-saturated powder layer subject to a short-pulsed heat source | |
Duan et al. | Temperature field simulations during selective laser melting process based on fully threaded tree | |
CN115740493A (zh) | 一种电子束选区熔化温度场和应力场的有限元模拟方法 | |
White et al. | Coupling Simulated Annealing and Homogenization to Design Thermally Conductive Hybrid Lattice Support Structures for LPBF |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |