曲嘉部门：航天与建筑工程学院

学科：力学

职务：

职称：副教授

指导
资格：硕士生导师
电话：

传真：**

邮箱：qujia@hrbeu.edu.cn

邮编：150001

地址：黑龙江省哈尔滨市南岗区南通大街145-1号
个人简介
1979年11月出生,男,汉族,中国共产党员,主要从事强度理论、固体实验力学和冲击动力学方面研究，主持完成国家自然科学基金、黑龙江省自然科学基金、哈尔滨市科学基金、教育部博士学科点基金等纵向项目，主持完成横向课题多项。


教育经历
1997.9-2001.7 哈尔滨工程大学 工程力学专业 本科生
2003.9-2005.6 哈尔滨工程大学 固体力学专业 硕士研究生
2004.9-2010.6 哈尔滨工程大学 固体力学专业 博士研究生

工作经历
2001.9-2005.9 哈尔滨工程大学 力学实验中心 助 教
2005.9-2011.9 哈尔滨工程大学 航天与建筑工程学院 讲 师
2011.9至今 哈尔滨工程大学 航天与建筑工程学院 副教授
期间
2015.9-2016.9 美国普渡大学 航天宇航学院 访问学者 联合教授陈为农
研究方向
（1）材料强度理论 ：
主要研究材料应变率相关强度准则，双轴或三轴应力状态对材料强度影响。
（2）实验固体力学 ：
主要研究高速加载技术/高速卸载技术，数字图像相关测量技术，瞬态电测技术等方向。
（3）冲击动力学：
主要研究分离式霍普金森压杆（Split Hopkinson Press Bar，Kolsky Bar），流固耦合冲击加载和测试技术。
（4）CAD与FEM应用
主要应用SolidWorks软件和ANSYS/LS-Dyna软件解决工程问题、分析实验现象。

承担项目
（一）纵向课题
[1] 2011年-2014年，高等学校博士学科点专项科研基金（**016），《钢纤维混凝土动态强度准则的实验研究》，项目负责人。
[2] 2012年-2015年，黑龙江省自然科学基金(No.A201202)，《UHPC动态强度准则的理论与实验研究》，项目负责人。
[3] 2013年-2016年，国家自然科学基金（**），《混凝土应变率相关强度理论研究及实验验证》，项目负责人。
[4] 2013年-2016年，哈尔滨市应用技术研究与开发项目(2013RFQXJ009),《几何光学出平面位移测量装置研究》，项目负责人。
（二）横向课题
[5] 2008年-2009年，渤海造船厂集团有限公司，《主应变测试系统研究》，技术负责人。
[6] 2011年-2012年，哈尔滨核海科技发展有限公司《负荷测量技术研究》，项目负责人。
[7] 2011年-2013年，哈尔滨核海科技发展有限公司《船舶动力设备安装检测技术研究》，项目负责人。
[8] 2012年-2014年，航天科工集团六院41所，《密封圈回弹性能测试系统》，技术负责人。
[9] 2014年-2015年，渤海造船厂集团有限公司，《加载变形对轴系安装影响检测装置》，项目负责人。
[10] 2014年-2015年，航天科工集团六院41所，《二维数字图像相关应变测量系统》，项目负责人。
（三）访学期间研究课题

[11]2015年-2016年，美国普渡大学，《Dynamic Triaxial Compression Experiments on Al Powder》，技术负责人

学术交流
招生信息
本科生授课课程
研究生授课课程
实践性教学
教学研究课题
社会兼职
黑龙江省兵工学会 副秘书长
黑龙江省航空学会 联络员


专利成果
专利名称作者专利
分类申请日期申请号/专利号
一种弯曲承载钢梁邹广平、商磊、张学义、夏兴有、杨丽红、曲嘉发明****.6
一种金属橡胶减震垫的制作方法邹广平、唱忠良、李雨蕾、吴国辉、曲嘉、芦颉发明****.2
一种保护型隔震器五联蠕变试验机曲嘉，高申煣，庞跃钊，孙晓庆，李东昌，夏培秀，薛启超，陈森，黄超发明****.7
一种同时多级加载蠕变试验机曲嘉，孙晓庆，李东昌，黄超，陈森，夏培秀发明****.8
一种三点弯曲疲劳试验机曲嘉，孙晓庆，李东昌，黄超，陈森，夏培秀发明****.1
一种基于拉线式位移传感器的轴系工程构件端面空间位移和角度变化量测量方法曲嘉，李东昌，孙晓庆，黄超，陈森，夏培秀发明****.0
一种基于几何光学的镜面出平面位移测量装置曲嘉，李东昌，孙晓庆，黄超，庞跃钊，高申煣，夏培秀发明****.0
一种基于平行激光的密封圈回弹速率测试装置曲嘉，庞跃钊，陈森，李东昌，孙晓庆，夏培秀发明****.3
一种可调式金属橡胶三向减振装置邹广平；曲嘉；李雨蕾；唱忠良实用
新型****.3
一种弯曲承载钢梁邹广平、商磊、张学义、夏兴有、杨丽红、曲嘉实用
新型****.x
一种保护型隔震器五联蠕变试验机曲嘉，高申煣，庞跃钊，孙晓庆，李东昌，夏培秀，薛启程，陈森，黄超实用
新型****.x
一种同时多级加载蠕变试验机曲嘉，孙晓庆，李东昌，黄超，陈森，夏培秀实用
新型****.3
一种三点弯曲疲劳试验机曲嘉，孙晓庆，李东昌，黄超，陈森，夏培秀实用
新型****.2



出版著作



发表论文
2006年
[1] 邹广平,曲嘉,吴国辉. 大量程扭转试验机标定及扭力扳手设计方法的研究[J]. 计量学报,2006,01:50-52.(EI:**66)
[2] 邹广平,曲嘉. 基于电阻应变计的称重装置[J]. 中国测试技术, 2006,03 :140-142.
[3] 施磊,牛玉梅,曲嘉. 根管预备锥度对根管壁应力影响的三维有限元分析[J]. 口腔医学研究,2006,03:255-257.
2007年
[4] 杨丽红,何蕴增,曲嘉. 实心圆轴扭转问题中的旋率分析[J]. 黑龙江工程学院学报,2007,01:22-25.
[5] Zou Guangping, Qu Jia. The influence of the end of the incidence pole on the KId of the three point bending specimen tested by the SHPB[J]. Key Engineering Materials, v 334-335 I, p 153-156, 2007（EI：**477, WOS:**039）
[6] Zou Guangping, Qu Jia, Qiao Yingjie. A study of the dynamic plane strain fracture toughness on concrete by SHPB[J]. PROCEEDINGS OF THE 3RD WSEAS INTERNATIONAL CONFERENCE ON APPLIED AND THEORETICAL MECHANICS (MECHANICS '07): TOPICS IN ADVANCED THEORETICAL AND APPLIED MECHANICS,p254-257,2007.12 (WOS:**043)
[7] Yang Li-hong, He Yun-zeng, Qu, Jia. Investigation on large deformation constitutive model based on solid cylinder torsion test[J]. PROCEEDINGS OF THE 5TH INTERNATIONAL CONFERENCE ON NONLINEAR MECHANICS,p309-312,2007(WOS:**058)
2008年
[8] Qu Jia, Zou Guang-Ping. A study of the dynamic plane strain fracture toughness of concrete by shpb[J]. Key Engineering Materials, v 385-387, p 21-24, 2008（EI：**510, WOS:**006）
2009年
[9] 曲嘉,邹广平,何蕴增. 扭转冷作硬化对低碳钢高应变率下后屈服强度的影响[A]. 中国力学学会爆炸力学专业委员会冲击动力学专业组.第九届全国冲击动力学学术会议论文集（上册）[C].中国力学学会爆炸力学专业委员会冲击动力学专业组:,2009:5.
2010年
[10] 曲嘉,邹广平,何蕴增. 扭转冷作硬化对低碳钢高低应变率力学性能的影响[J]. 哈尔滨工程大学学报,2010,05:596-600. （EI：**542）
2011年
[11] 曲嘉, 杨鹏飞, 李勇贵, 李雨蕾. 子弹对薄钢板穿透临界速度与入射角度的仿真研究[A]. 中国力学学会爆炸力学专业委员会冲击动力学专业组、北京理工大学爆炸科学与技术国家重点实验室.第十届全国冲击动力学学术会议论文摘要集[C].中国力学学会爆炸力学专业委员会冲击动力学专业组、北京理工大学爆炸科学与技术国家重点实验室,2011:1.
[12] 邹广平, 唱忠良, 曲嘉, 陈思, 王瑞瑞. 平面应变下紧凑拉伸试样的动态断裂韧性的实验研究[J]. 实验力学,2011,03:240-246.
[13] 邹广平, 唱忠良, 曲嘉. 紧凑拉伸试样应力强度因子的动态响应[A]. 中国力学学会、黑龙江省力学学会、辽宁省力学学会.第七届海峡两岸工程力学研讨会论文摘要集[C].中国力学学会、黑龙江省力学学会、辽宁省力学学会:,2011:5.
[14] 杨丽红, 吴林志, 曲嘉. 镍基单晶合金板裂纹尖端塑性区分析[J]. 哈尔滨工程大学学报,2011,11:1434-1438+1500.
[15] Yang Lihong , Wu Linzhi, Qu Jia. Crack tip plastic zone of a nickel-based single crystal plate. Journal of Harbin Engineering University, v 32, n 11, p 1434-1438+1500, November 2011 ,（EI：**040）
[16] Qu Jia, Zou Guang-Ping, Xia Pei-Xiu. Experimental investigation on SFRC split strength under tri-axial stress[J]. Binggong Xuebao/Acta Armamentarii, v 32, n SUPPL. 2, p 346-349, 2011.12（EI：**487）
[17] Yang Lihong , Qu Jia, Zou Guangping. Investigation on finite strain constitutive model with logarithmic strain based on torsion experiments[J]. Key Engineering Materials, v 450, p 522-525, 2011, Advanced Design and Manufacture III, ADM2010（EI：**344, WOS:**129）
[18] Qu Jia, Zou Guangping,Xia Peixiu. Experimentation of Steel Fiber Reinforced Concrete Tri-axial Split Strength under Quasi-Static Loading and High Strain-Rate[J]. ADVANCES IN EXPERIMENTAL MECHANICS VIII/ Applied Mechanics and Materials .v70,p189-194,2011.9 ( WOS:**030)
2012年
[19] Yang Lihong , Qu Jia, He Yunzeng.Stress responses to large simple shear deformation in elasticity based on the logarithmic strain[J]. Key Engineering Materials, v 488-489, p 424-427, 2012, Advances in Fracture and Damage Mechanics X（EI：**083,WOS:**105）
[20] Yang Lihong , Zou Guangping, Qu Jia. Yield criterion and crack tip plastic zone of nickel-based single crystal[J]. Key Engineering Materials, v 525-526, p 341-344, 2012, Advances in Fracture and Damage Mechanics XI（EI：**122,WOS:**086）
2013年
[21] 曲嘉,陈庚,杨育伟. 橡胶密封圈回弹性能的仿真分析[A]. 中国力学学会、西安交通大学.中国力学大会——2013论文摘要集[C].中国力学学会、西安交通大学, 2013:1.
[22] Qu Jia , Chen Geng, Yang Yu Wei. Finite element analysis of rubber sealing ring resilience behavior[J]. Advanced Materials Research, v 705, p 410-414, 2013, MEMS and Mechanics（EI：**555）.
[23] Chen Geng, Yang Yu-Wei, Qu Jia. The influence of stress relaxation on resilience of seal ring[J]. Dongbei Daxue Xuebao/Journal of Northeastern University, v 34, n S2, p 130-133, 2013.12 ,（EI：**785）.
2014年
[24] Chen, Geng, Qu Jia, Tian ChunYu. Numerical analysis of ring beam effect on the modal of LNG storage tank[J]. Applied Mechanics and Materials, v 638-640, p 172-177, 2014（EI：**875）
2015年
[25] 黄超, 李东昌, 曲嘉. 数字图像相关技术在处理橡胶类材料埋藏缺陷失效现象的研究[A]. 中国力学学会、上海交通大学.中国力学大会-2015论文摘要集[C].中国力学学会、上海交通大学,2015:1.
[26] 孙晓庆, 黄超, 曲嘉. 周期结构对SHPB波形整形效果的仿真研究[A]. 中国力学学会、上海交通大学.中国力学大会-2015论文摘要集[C].中国力学学会、上海交通大学,2015:1.
[27] 曲嘉, 邹广平, 夏培秀. 钢纤维混凝土准静态和冲击三轴劈拉强度实验研究[J].中国科技论文在线精品论文,v 7,No.15,p1501-1507,2014.8.
2016年
[28] Li Dongchang, Chao Huang, Jia Qu. Unable to detect burial defective rubber sheet in Digital Image Correlation method[C]. Key Engineering Materials, v 665, p 209-212, 2016(EI:**703)
[29] 曲嘉, 李东昌, 黄超. 动态断裂韧性实验中DIC技术应用[J],中国测试,v42,p44-48,2016.10
[30] 曲嘉, 庞跃钊, 孙晓庆. 316L三点弯曲试样动静态断裂韧性对比实验研究[J], 已录用待发表.
[31] 曲嘉, 高申煣, 陈森. 高应变率下SFRC三轴劈拉强度实验研究[J],已录用待发表.
[32] 曲嘉, 高申煣, 孙晓庆. 混凝土强度动态增强因子的线性计算模型研究[J],已投稿.
[33] 曲嘉, 庞跃钊, 李东昌. 316L不锈钢不同应变率下的冲击拉伸力学性能研究[J],已投稿.
荣誉

奖励



率相关强度准则
应变率相关强度准则计算程序MATLAB源代码(面向科研人员和学生免费开放)
function h = ezimplot3(varargin)
if ishandle(varargin{1})
cax = varargin{1};
axes(cax);
args{:} = varargin{2:end};
else
args = varargin;
end
[fun domain n color] = argcheck(args{:});
xm = linspace(domain(1),domain(2),n);
ym = linspace(domain(3),domain(4),n);
zm = linspace(domain(5),domain(6),n);
[x,y,z] = meshgrid(xm,ym,zm);
[f_handle f_text] = fix_fun(fun);
try
fvalues = f_handle(x,y,z);
catch ME
error('Ezimplot3:Functions', 'FUN must have no more than 3 arguments');
end
h = patch(isosurface(x,y,z,fvalues,0));
isonormals(x,y,z,fvalues,h)
set(h,'FaceColor',color,'EdgeColor','none');
xlabel('x');ylabel('y');zlabel('z');
alpha(0.7)
grid on; view([1,1,1]); axis equal; camlight; lighting gouraud
title([f_text,' = 0']);
function [f dom n color] = argcheck(varargin)
switch nargin
case 0
error('Ezimplot3:Arguments',...
'At least "fun" argument must be given');
case 1
f = varargin{1};
dom = [-2*pi, 2*pi];
n = 60;
color = 'red';
case 2
f = varargin{1};
if isa(varargin{2},'double') && length(varargin{2})>1
dom = varargin{2};
n = 60;
color = 'red';
elseif isa(varargin{2},'double') && length(varargin{2})==1
n = varargin{2};
dom = [-2*pi, 2*pi];
color = 'red';
elseif isa(varargin{2},'char')
dom = [-2*pi, 2*pi];
n = 60;
color = varargin{2};
end
case 3
f = varargin{1};
dom = varargin{2};
n = varargin{3};
color = 'red';
case 4
f = varargin{1};
dom = varargin{2};
n = varargin{3};
color = varargin{4};
otherwise
warning('Ezimplot3:Arguments', ...
'Attempt will be made only with the 4 first arguments');
f = varargin{1};
dom = varargin{2};
n = varargin{3};
color = varargin{4};
end
if length(dom) == 2
dom = repmat(dom,1,3);
elseif length(dom) ~= 6
error('Ezimplot3:Arguments',...
'Input argument "domain" must be a row vector of size 2 or size 6');
end
function [f_hand f_text] = fix_fun(fun)
types = {'char','sym','function_handle','inline'};
type = '';
for i=1:size(types,2)
if isa(fun,types{i})
type = types{i};
break;
end
end
switch type



计算程序(续1)
case 'char'
if exist([fun,'.m'],'file')
syms x y z;
if nargin(str2func(fun)) == 3
f_sym = eval([fun,'(x,y,z)']); else
error('Ezimplot3:Arguments',...
'%s must be a function of 3 arguments or unknown function',fun);
end
f_text = strrep(char(f_sym),' ','');
f_hand = eval(['@(x,y,z)',vectorize(f_text),';']);
else
f_hand = eval(['@(x,y,z)',vectorize(fun),';']);
f_text = strrep(fun,'.',''); f_text = strrep(f_text,' ',''); end
case 'sym'
f_hand = eval(['@(x,y,z)',vectorize(fun),';']);
f_text = strrep(char(fun),' ','');
case {'function_handle', 'inline'}
syms x y z;
if nargin(fun) == 3 && numel(symvar(char(fun))) == 3
f_sym = fun(x,y,z);
else
error('Ezimplot3:Arguments',...
'%s must be function of 3 arguments or unknown function',char(fun));
end
f_text = strrep(char(f_sym),' ','');
f_hand = eval(['@(x,y,z)',vectorize(f_text),';']);
otherwise
error('First argument "fun" must be of type character, simbolic, function handle or inline');
end
k=1; %下面分别计算含参数的两个动态模型
syms x y z;
fk=(1/22877.17).*((x-y+k.*(3.*(x-y)/11438.59)./((2*x-y-z+8.77).^2./11438.59+(2*y-x-z+8.77).^2./11438.59+(2*z-y-x+8.77).^2./11438.59).^0.5)^2+...
(x-z+k.*(3.*(x-z)./11438.59)./((2*x-y-z+8.77).^2./11438.59+(2*y-x-z+8.77)^2./11438.59+(2*z-y-x+8.77).^2./11438.59).^0.5).^2+...
(y-z+k.*(3.*(y-z)./11438.59)./((2*x-y-z+8.77).^2./11438.59+(2*y-x-z+8.77).^2./11438.59+(2*z-y-x+8.77).^2./11438.59).^0.5).^2)+...
(1/8.77).*(x+y+z+26.31*k./((2*x-y-z+8.77).^2./11438.59+(2*y-x-z+8.77).^2./11438.59+(2*z-y-x+8.77).^2./11438.59).^0.5)-1;
f2=@(x,y,z)fk;
ezimplot3(f2,[-20000,20000]) %带参数的扩张模型
clear;
f1=@(x,y,z)((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-1/105.36;
ezimplot3(f1,[-20000,20000]); %带参数的随动模型
clear %求解随动模型中单轴拉压值
k=1e-6; %10-6应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('f=0','x')
clear
k=1e-3; %10-3应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('f=0','x')
clear




(续2)
k=1; %100应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f；
solve('f=0','x')
clear
k=10; %101应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('f=0','x')
clear
k=100; %102应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')
clear
k=1000; %103应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')
clear
k=10000; %104应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')
fsolve('f=0','x')
clear
k=100; %102应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')
clear
k=1000; %103应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')
clear
k=10000; %104应变率下的单轴拉压值
syms x;
y=0;z=0;
f=((x-y).^2+(x-z).^2+(y-z).^2)/(6*20086.16)+(x+y+z)/92.4-(1+0.14*k)/105.36;
f
solve('=0','x')