comsol与matlab连接使用学习笔记 3 材料随机分布对应力影响
前言
在结构力学中,常常假定材料是均质的,单一域内使用相同材料属性来做相应的结构分析。然而,实际材料很难做到完全均质。在这个案例中,将考察对一个平板简单拉伸的案例,其域内杨氏模量满足正态分布,研究分布的标准差对最大应力的影响。
这个案例中,使用matlab定义一个产生随机分布的外部函数,并且用matlab对标准差进行扫描。
产生随机分布的matlab函数
该函数将定义域内材料杨氏模量随机分布,自变量有点坐标序列 x,杨氏模量分布期望值 E 和标准差 Fe_sigma 。这里用钢的数值,E = 2e11 Pa,泊松比为常数0.3. 函数的返回值是坐标序列对应的杨氏模量序列。
function out = randn_E(x, Fe_E, Fe_sigma)
out = Fe_E + Fe_sigma.*randn(size(x));
保存为 randn_E.m 文件,并与 comsol 模型的.mph文件,放在同一路径下。
comsol 模型设置
打开comsol multiPhysics进行建模
新建:模型向导
选择空间维度:3D
选择物理场:结构力学,稳态
参数定义:
Fe