| 查看: 489 | 回复: 0 | |||
| 【悬赏金币】回答本帖问题，作者秃头娃将赠送您 15 个金币 | |||
[求助]
正交晶系的三维杨氏模量图怎么作？
|
|||
|
请问一下大家，正交晶系里面的matlab杨氏模量图怎么做出来啊？可以做出立方和四方的(*｀へ?*)，但是正交的怎么都做不出来，附上的我代码，麻烦各位老师同学帮我看一下，谢谢！ %calculate the directional dependence of bulk and elastic moduli (rhombohedral) clear clc S11= 0.001049 ; S22= 0.001038 ; S12= -0.000036; S66= 0.003580 ; S13= -0.000083 ; S55= 0.002496 ; S23= -0.000067; S33= 0.000769; S44= 0.002986; theta=0:pi/36:2*pi; pha=0:pi/36:2*pi; %计算三维图时用 [TH,PH]=meshgrid(theta,pha); ST=sin(TH);CT=cos(TH); SP=sin(PH);CP=cos(PH); L1= ST.*CP; L2= ST.*SP; L3=CP; %the elastic surface of plan (001) dEM=(L1.^4).*S11+(L2.^4).*S22+(L3.^4).*S33+(L2.^2).*(L3.^2).*(2.*S23+S44)+(L1.^2).*(L3.^2).*(2.*S13+S55)+(L1.^2).*(L2.^2).*(2.*S12+S66); EM=1./dEM; XE=EM.*ST.*CP; YE=EM.*ST.*SP; ZE=EM.*CT; figure,surf(XE,YE,ZE),box on,axis equal xlabel('Ex');ylabel('Ey');zlabel('Ez');box on %legend('elastic modul') 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
