CN108984887A - 复合材料参数多阶段识别方法 - Google Patents
复合材料参数多阶段识别方法 Download PDFInfo
- Publication number
- CN108984887A CN108984887A CN201810736681.6A CN201810736681A CN108984887A CN 108984887 A CN108984887 A CN 108984887A CN 201810736681 A CN201810736681 A CN 201810736681A CN 108984887 A CN108984887 A CN 108984887A
- Authority
- CN
- China
- Prior art keywords
- parameter
- composite material
- statics
- elastic
- finite element
- 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
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]
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种复合材料确定性参数多阶段识别方法,包括以下步骤:步骤1,在有限元分析软件中建立复合材料的初始有限元模型,并计算其静力学位移响应;步骤2,根据复合材料结构静力学位移响应与复合材料弹性参数之间的力学关系,计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度;步骤3,在步骤1和2的基础上,根据相对灵敏度分析结果使用多阶段方法识别复合材料参数;步骤4,根据试验数据识别复合材料全部弹性参数。本发明使用三点弯位移数据,多阶段模型修正方法识别复合材料全部弹性参数。多阶段识别方法同时解决了识别灵敏度数量级相差较大的参数,导致修正过程中灵敏度矩阵病态的问题。
Description
技术领域
本发明涉及复合材料结构动力学参数识别技术领域,具体涉及一种采用基于静力学位移响应来识别复合材料全部弹性参数的多阶段识别方法。
背景技术
复合材料结构具有高比强度、高比刚度、耐疲劳等优异性能,广泛应用航空航天领域。但由于复合材料结构的不确定性,其材料参数也会随着生产工艺,内部结构不同而发生改变。准确的复合材料的弹性参数是研究复合材料力学性能的基础。复合材料等效参数为正交各向异性,具有9个独立的参数。传统的测量复合材料参数方法通过拉伸,剪切测量应变。再通过应力应变关系计算其刚度。
传统的模型修正方法是通过计算参数对响应的灵敏度,选取灵敏度较大的参数进行识别,放弃了其他的参数。或同时识别所有的参数,对识别结果保留一定的误差。显然这样建立的复合材料模型不会准确,对于结构后续实际工况研究也加入了一定的误差。难以捕捉其实际载荷下的力学性能。
发明内容
本发明的目的在于提供一种复合材料确定性参数多阶段识别方法,以解决现有技术存在的复合材料仿真模拟的精度差的问题。
为实现上述目的,本发明采用的技术方案为:
一种复合材料确定性参数多阶段识别方法,包括以下步骤:
步骤1,在有限元分析软件中建立复合材料的初始有限元模型,并通过有限元分析软件计算其静力学位移响应;
步骤2,根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系,计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度;
步骤3,在步骤1和2的基础上,研究复合材料参数多阶段预测方法,将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集,基于计算静力学位移参数与实验静力学参数的残差建立目标函数,并对待修正参数进行迭代优化,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;并进行仿真验证研究;
步骤4,根据试验数据识别复合材料全部弹性参数。
所述步骤1包括以下步骤:
步骤1.1,在有限元软件中,根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型;其中,9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33,三个泊松比ν12、ν23、ν31,与三个平面内的剪切弹性模量G12、G23、G31;
步骤1.2,根据步骤1.1得到的初始有限元模型,通过有限元分析软件对其进行静力学分析,得到复合材料的静力学位移响应。
所述步骤2包括以下步骤:
步骤2.1,由复合材料力学知识得到材料刚度矩阵C如下式表示:
上式中,E1、E2、E3分别表示主轴方向1,2,3的弹性模量;νij表示单轴应力σ在j 方向作用时引起i方向应变的泊松耦合系数,即泊松比;G12、G23、G31分别表示1-2,2-3,3-1 平面的剪切弹性模量;
刚度矩阵C中各元素与独立弹性参数之间的关系如下所示:
C11=(1-ν23ν32)·E11·B,C12=(ν12+ν13ν32)·E22·B
C22=(1-ν13ν31)·E22·B,C23=(ν23+ν21ν13)·E33·B
C33=(1-ν12ν21)·E33·B,C31=(ν31+ν21ν32)·E11·B
其中
上式中,Ckl分别为6×6刚度矩阵中对应位置上的元素,其中,0≤k≤6,0≤l≤6;
步骤2.2,研究复合材料结构静力学响应与复合材料弹性参数之间的力学关系,并进行相对灵敏度分析,计算公式为:其中,fe为静力试验参数,通过三点弯静力学试验测得;Δf为静力学位移响应与实验参数的残差;p为需要识别的弹性参数; S为相对灵敏度,各个弹性参数对于静力学响应的影响程度,计算公式如下表示:
其中f抽象地代表静力学响应,fi为仿真模型有限元分析出的静力学位移响应,由于刚度矩阵C为对称矩阵故只需要求静力学相响应fi对刚度矩阵中元素C11、C12、C13、 C22、C23、C333、C443、C553、C66的偏导,为刚度矩阵中各元素对待求的独立弹性参数的偏导数。
所述步骤3包括以下步骤:
步骤3.1,根据相对灵敏度分析结果,将相差小于十倍的参数灵敏度划分为相同的集合;则复合材料的9个独立弹性参数划分为n个子集;
步骤3.2,构建以计算响应参数与实验响应参数的加权残差为目标函数;
即在参数p的取值范围p1≤p≤p2内找到一个pa,其中p的取值范围为9个独立弹性参数可能值的上下限;使目标函数J(p)为试验与计算参数的加权残差取极小值;e为特征值的残差,λe,λa(p)分别为试验与计算的参数,加权矩阵W为反映各特征值残差相对权重的对角阵;
步骤3.3,采用分布迭代修正的模式,对步骤3.1中n个集合分别识别,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;从而识别复合材料全部弹性参数。
所述步骤4的根据试验数据识别复合材料全部弹性参数包括以下步骤包括以下步骤:
步骤4.1,开展试验验证研究,采用复合材料板结构三点弯试验位移数据,来识别复合材料的全部的弹性参数。
有益效果:本发明解决了繁琐的传统复合材料参数测量方法,使用三点弯位移数据,多阶段模型修正方法识别复合材料全部弹性参数。多阶段识别方法同时解决了识别灵敏度数量级相差较大的参数,导致修正过程中灵敏度矩阵病态的问题。对于建立准确的复合材料动力学模型具有指导作用,具有工程应用价值。
附图说明
图1是复合材料确定性参数多阶段识别方法具体步骤的流程图;
图2是编织复合材料板仿真模型;
图3是三点弯节点位移对9个独立弹性参数的相对灵敏度绝对值柱状图;
图4是多阶段迭代曲线图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,本发明的一种复合材料确定性参数多阶段识别方法,包括以下步骤:
一种复合材料确定性参数多阶段识别方法,其特征在于:包括以下步骤:
步骤1,在有限元分析软件中建立复合材料的初始有限元模型,并通过有限元分析软件计算其静力学位移响应;具体步骤为:
步骤1.1,在有限元软件中,根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型;其中,9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33,三个泊松比ν12、ν23、ν31,与三个平面内的剪切弹性模量G12、G23、G31;
步骤1.2,根据步骤1.1得到的初始有限元模型,通过有限元分析软件对其进行静力学分析,得到复合材料的静力学位移响应。
步骤2,根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系,计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度;具体步骤为:
步骤2.1,由复合材料力学知识得到材料刚度矩阵C如下式表示:
上式中,E1、E2、E3分别表示主轴方向1,2,3的弹性模量;νij表示单轴应力σ在j 方向作用时引起i方向应变的泊松耦合系数,即泊松比;G12、G23、G31分别表示1-2,2-3,3-1 平面的剪切弹性模量;
刚度矩阵C中各元素与独立弹性参数之间的关系如下所示:
C11=(1-ν23ν32)·E11·B,C12=(ν12+ν13ν32)·E22·B
C22=(1-ν13ν31)·E22·B,C23=(ν23+ν21ν13)·E33·B
C33=(1-ν12ν21)·E33·B,C31=(ν31+ν21ν32)·E11·B
其中
上式中,Ckl分别为6×6刚度矩阵中对应位置上的元素,其中,0≤k≤6,0≤l≤6;
步骤2.2,研究复合材料结构静力学响应与复合材料弹性参数之间的力学关系,并进行相对灵敏度分析,计算公式为:其中,fe为静力试验参数,通过三点弯静力学试验测得;Δf为静力学位移响应与实验参数的残差;p为需要识别的弹性参数; S为相对灵敏度,各个弹性参数对于静力学响应的影响程度,计算公式如下表示:
其中f抽象地代表静力学响应,fi为仿真模型有限元分析出的静力学位移响应,由于刚度矩阵C为对称矩阵故只需要求静力学相响应fi对刚度矩阵中元素C11、C12、C13、 C22、C23、C333、C443、C553、C66的偏导,为刚度矩阵中各元素对待求的独立弹性参数的偏导数。
步骤3,在步骤1和2的基础上,研究复合材料参数多阶段预测方法,将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集,基于计算静力学位移参数与实验静力学参数的残差建立目标函数,并对待修正参数进行迭代优化,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;并进行仿真验证研究;具体步骤为:
步骤3.1,根据相对灵敏度分析结果,将相差小于十倍的参数灵敏度划分为相同的集合;则复合材料的9个独立弹性参数划分为n个子集;
步骤3.2,构建以计算响应参数与实验响应参数的加权残差为目标函数;
即在参数p的取值范围p1≤p≤p2内找到一个pa,其中p的取值范围为9个独立弹性参数可能值的上下限;使目标函数J(p)为试验与计算参数的加权残差取极小值;e为特征值的残差,λe,λa(p)分别为试验与计算的参数,加权矩阵W为反映各特征值残差相对权重的对角阵;
步骤3.3,采用分布迭代修正的模式,对步骤3.1中n个集合分别识别,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;从而识别复合材料全部弹性参数。
步骤4,根据试验数据识别复合材料全部弹性参数,具体步骤为:开展试验验证研究,采用复合材料板结构三点弯试验位移数据,来识别复合材料的全部的弹性参数。
下面结合具体实施例对本发明做进一步说明。
实施例
本实施例为复合材料确定性参数多阶段识别方法,具体包括以下步骤:
步骤1,在有限元分析软件中建立复合材料的初始有限元模型,并计算其静力学位移响应;
研究对象是几何尺寸为100×50×3mm的编织复合材料板,其9个独立弹性参数见表 1,包括三个材料主轴方向的弹性模量E11、E22、E33,三个泊松比ν12、ν23、ν31,与三个平面内的剪切弹性模量G12、G23、G31;在NASTRAN有限元分析软件中根据几何尺寸和受边界条件载荷建立的有限元模型,如图2所示。
表1复合材料板材料参数
步骤2,根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系,计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度;
计算X、Y、Z对这9个系数的相对灵敏度,求得的相对灵敏度结果取绝对值后做柱状图如图3所示。
步骤3,根据相对灵敏度分析结果,将量级相近的参数灵敏度划分为相同的集合;根据上述计算得到的相对灵敏度柱状图,将复合材料9个独立弹性参数分为2组分别为:组1:E11、E22、ν12、G12、G23、G31,组2:E33、ν23、ν31。
构建以计算响应参数与实验响应参数的加权残差为目标函数及分布式优化;
通过matlab编程调用有限元分析软件计算,确定迭代优化目标函数。采用分布迭代修正的模式,对步上步集合分别识别,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;从而识别复合材料全部弹性参数。
步骤4,识别结果仿真验真:通过摄动初值得到实验位移值,摄动的倍数分别为:1.3*E11、1.3*E22、1.3*E33、1.1*ν12、1.1*ν23、1.1*ν31、1.2*G12、1.2*G23、1.2*G31。经过本方法的复合材料确定性参数多阶段识别方法后,上述参数分别为初始值的-1.5× 10-5%、-2.8×10-6%、-0.00267%、-0.00033%、0.225122%、-0.03093%、-0.00015%、-6.2 ×10-5%、-0.00062%。识别后的结果接近参数摄动值。验证了本方法的可行性。图3-图 4为参数迭代结果。
由表2结果可以看出,经过迭代计算得到修正结果,修正后的频率值相对于摄动后的试验频率值,误差绝对值最大仅为0.23%,识别结果较好。
表2参数迭代结果
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种复合材料确定性参数多阶段识别方法,其特征在于:包括以下步骤:
步骤1,在有限元分析软件中建立复合材料的初始有限元模型,并通过有限元分析软件计算其静力学位移响应;
步骤2,根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系,计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度;
步骤3,在步骤1和2的基础上,研究复合材料参数多阶段预测方法,将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集,基于计算静力学位移参数与实验静力学参数的残差建立目标函数,并对待修正参数进行迭代优化,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;并进行仿真验证研究;
步骤4,根据试验数据识别复合材料全部弹性参数。
2.根据权利要求1所述的复合材料确定性参数多阶段识别方法,其特征在于:所述步骤1包括以下步骤:
步骤1.1,在有限元软件中,根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型;其中,9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33,三个泊松比ν12、ν23、ν31,与三个平面内的剪切弹性模量G12、G23、G31;
步骤1.2,根据步骤1.1得到的初始有限元模型,通过有限元分析软件对其进行静力学分析,得到复合材料的静力学位移响应。
3.根据权利要求1所述的复合材料确定性参数多阶段识别方法,其特征在于:所述步骤2包括以下步骤:
步骤2.1,由复合材料力学知识得到材料刚度矩阵C如下式表示:
上式中,E1、E2、E3分别表示主轴方向1,2,3的弹性模量;νij表示单轴应力σ在j方向作用时引起i方向应变的泊松耦合系数,即泊松比;G12、G23、G31分别表示1-2,2-3,3-1平面的剪切弹性模量;
刚度矩阵C中各元素与独立弹性参数之间的关系如下所示:
C11=(1-ν23ν32)·E11·B,C12=(ν12+v13ν32)·E22·B
C22=(1-v13ν31)·E22·B,C23=(ν23+ν21ν13)·E33·B
C33=(1-ν12ν21)·E33·B,C31=(ν31+ν21ν32)·E11·B
其中
上式中,Ckl分别为6×6刚度矩阵中对应位置上的元素,其中,0≤k≤6,0≤l≤6;
步骤2.2,研究复合材料结构静力学响应与复合材料弹性参数之间的力学关系,并进行相对灵敏度分析,计算公式为:其中,fe为静力试验参数,通过三点弯静力学试验测得;Δf为静力学位移响应与实验参数的残差;p为需要识别的弹性参数;S为相对灵敏度,各个弹性参数对于静力学响应的影响程度,计算公式如下表示:
其中f抽象地代表静力学响应,fi为仿真模型有限元分析出的静力学位移响应,由于刚度矩阵C为对称矩阵故只需要求静力学相响应fi对刚度矩阵中元素C11、C12、C13、C22、C23、C333、C443、C553、C66的偏导,为刚度矩阵中各元素对待求的独立弹性参数的偏导数。
4.根据权利要求1所述的复合材料确定性弹性参数多阶段识别方法,其特征在于:所述步骤3包括以下步骤:
步骤3.1,根据相对灵敏度分析结果,将相差小于十倍的参数灵敏度划分为相同的集合;则复合材料的9个独立弹性参数划分为n个子集;
步骤3.2,构建以计算响应参数与实验响应参数的加权残差为目标函数;
即在参数p的取值范围p1≤p≤p2内找到一个pa,其中p的取值范围为9个独立弹性参数可能值的上下限;使目标函数J(p)为试验与计算参数的加权残差取极小值;e为特征值的残差,λe,λa(p)分别为试验与计算的参数,加权矩阵W为反映各特征值残差相对权重的对角阵;
步骤3.3,采用分布迭代修正的模式,对步骤3.1中n个集合分别识别,优先识别灵敏度大的参数集合,直至精度达到要求且所有的参数集合均收敛;从而识别复合材料全部弹性参数。
5.根据权利要求1所述的复合材料确定性弹性参数多阶段识别方法,其特征在于:所述步骤4的根据试验数据识别复合材料全部弹性参数包括以下步骤包括以下步骤:
步骤4.1,开展试验验证研究,采用复合材料板结构三点弯试验位移数据,来识别复合材料的全部的弹性参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810736681.6A CN108984887B (zh) | 2018-07-06 | 2018-07-06 | 复合材料确定性参数多阶段识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810736681.6A CN108984887B (zh) | 2018-07-06 | 2018-07-06 | 复合材料确定性参数多阶段识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108984887A true CN108984887A (zh) | 2018-12-11 |
CN108984887B CN108984887B (zh) | 2023-04-25 |
Family
ID=64536348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810736681.6A Active CN108984887B (zh) | 2018-07-06 | 2018-07-06 | 复合材料确定性参数多阶段识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108984887B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110245455A (zh) * | 2019-06-25 | 2019-09-17 | 南京航空航天大学 | 一种单向碳纤维增强复合材料的设计方法 |
CN114062435A (zh) * | 2021-10-21 | 2022-02-18 | 北京工业大学 | 基于电阻抗成像原理的具有灵敏度修正功能的柔性传感器 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
-
2018
- 2018-07-06 CN CN201810736681.6A patent/CN108984887B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
CN107357992A (zh) * | 2017-07-13 | 2017-11-17 | 东南大学 | 基于聚类分析的复合材料结构有限元模型修正方法 |
Non-Patent Citations (1)
Title |
---|
秦福溶 等: "《复合材料学报》", 12 February 2018 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110245455A (zh) * | 2019-06-25 | 2019-09-17 | 南京航空航天大学 | 一种单向碳纤维增强复合材料的设计方法 |
CN114062435A (zh) * | 2021-10-21 | 2022-02-18 | 北京工业大学 | 基于电阻抗成像原理的具有灵敏度修正功能的柔性传感器 |
Also Published As
Publication number | Publication date |
---|---|
CN108984887B (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021114994A1 (zh) | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 | |
CN108645704B (zh) | 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法 | |
CN103886125B (zh) | 一种钛合金热复合成形数值模拟方法 | |
WO2019011026A1 (zh) | 基于聚类分析的复合材料结构有限元模型修正方法 | |
CN112199823B (zh) | 一种橡胶材料应力松弛及损伤效应的数值预测方法 | |
CN106644711B (zh) | 一种延性材料单轴本构关系测试方法 | |
CN109883825A (zh) | 金属圆棒试样单轴拉伸大应变范围硬化曲线的测量方法 | |
CN109684753A (zh) | 一种管材弯曲回弹角逆向预测和补偿方法 | |
CN105930553A (zh) | 镍基高温合金缺口拉伸强度预测方法 | |
CN106777691B (zh) | 用于结构动力学仿真的橡胶o形圈有限元建模方法 | |
CN103745114A (zh) | 一种钛合金应力松弛数值及回弹计算方法 | |
CN111144016B (zh) | 一种多轴应力状态混凝土变四参数损伤本构 | |
Meng et al. | Study on stress relaxation of membrane structures in the prestress state by considering viscoelastic properties of coated fabrics | |
CN105092374A (zh) | 适用于模型试验中涂层织物类膜材的弹性模量测试方法 | |
CN108984887A (zh) | 复合材料参数多阶段识别方法 | |
Kłosowski et al. | Investigation on rheological properties of technical fabric “Panama” | |
Lina et al. | FEM analysis of spring-backs in age forming of aluminum alloy plates | |
CN110135063B (zh) | 一种串联桁架结构体系非概率失效度计算方法 | |
CN107314938B (zh) | 焊核区域材料塑性参数反演识别的实现方法 | |
CN104392032A (zh) | 基于有限元法的纱线材料参数识别方法 | |
CN108548720A (zh) | I型裂纹弹塑性理论公式获取延性材料j阻力曲线的方法 | |
JP6568782B2 (ja) | ゴムの二軸引張をシミュレーションする方法、装置及びプログラム | |
Naimark et al. | Critical Dynamics of Damage to Composites and Two-Parameter Fracture Criteria | |
Wang et al. | Investigation of the capabilities of yield functions on describing the deformation behavior of 5754O aluminum alloy sheet under combined loading paths | |
Harsono et al. | Material characterization based on simulated spherical-Berkovich indentation tests |
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 |